Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,835,133 T→A 6.6% intergenic (+88/+40) norW → / ← hypF NADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,835,1330TA 6.6% 188.8 / 3.6 61intergenic (+88/+40)norW/hypFNADH:flavorubredoxin oxidoreductase/carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base A (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GTGGTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAA  >  NC_000913/2834986‑2835277
                                                                                                                                                   |                                                                                                                                                
gtgGTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGtt                                                                                                                                                  <  6:128071/148‑1 (MQ=255)
   gTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                              >  5:267005/1‑149 (MQ=255)
            gATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                     <  6:321032/149‑1 (MQ=255)
            gATCGGATGAAAGAGGCATGTGGATTGTTGAAAACATTGCCGATGTAGGTGCGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                     <  5:8878/149‑1 (MQ=255)
                ggATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                   <  2:277071/147‑1 (MQ=255)
                  aTGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATg                                                                                                                               <  5:183317/149‑1 (MQ=255)
                         agGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTcc                                                                                                                         >  5:141186/1‑148 (MQ=255)
                          ggCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCaa                                                                                                                       <  8:104880/149‑1 (MQ=255)
                             aTTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                    >  8:4889/1‑133 (MQ=255)
                             aTTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                    <  7:4889/133‑1 (MQ=255)
                             aTTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATg                                                                                                                    >  8:12442/1‑149 (MQ=255)
                             aTTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATg                                                                                                                    <  8:313649/149‑1 (MQ=255)
                             aTTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATg                                                                                                                    >  4:11285/1‑149 (MQ=255)
                                     tgttgAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGc                                                                                                                   <  6:233510/142‑1 (MQ=255)
                                     tgttgAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGc                                                                                                                   >  5:233510/1‑142 (MQ=255)
                                       ttgAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAt                                                                                                                     >  7:244135/1‑138 (MQ=255)
                                       ttgAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAt                                                                                                                     <  8:244135/138‑1 (MQ=255)
                                         gAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCaaa                                                                                                                      <  2:313758/135‑1 (MQ=255)
                                         gAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCaaa                                                                                                                      >  1:313758/1‑135 (MQ=255)
                                            aaCATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCAGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATcc                                                                                                     >  4:271350/1‑149 (MQ=255)
                                                   ccGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAt                                                                                                                     <  4:176915/126‑1 (MQ=255)
                                                   ccGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAAt                                                                                                                     >  3:176915/1‑126 (MQ=255)
                                                   ccGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTg                                                                                               <  3:11285/148‑1 (MQ=255)
                                                    cGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGa                                                                                             <  3:1545/149‑1 (MQ=255)
                                                           ggTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACCCGTTTTATCCGTTCTGGACTTCAc                                                                                       <  1:35229/148‑1 (MQ=255)
                                                               ggCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCt                                                                                  <  8:26318/149‑1 (MQ=255)
                                                                  tACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAAc                                                                               >  4:102914/1‑149 (MQ=255)
                                                                    cTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCa                                                                             <  7:97404/149‑1 (MQ=255)
                                                                          cTAAAATGTCGGATGCGACGCTGGCGCGTATTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCATATTCAACACGTTTTATCCTTTCTGGACTTCACCCGCTAACCAACgcgc                                                                       >  4:75106/1‑149 (MQ=255)
                                                                           tAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGcgcc                                                                      >  5:110266/1‑149 (MQ=255)
                                                                              aaTGGCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAACTGCAACACGTTTTATCCGTTCTGGACTTCACCCGCGAACCACCGCGCCgca                                                                   <  6:42271/149‑1 (MQ=25)
                                                                                  tCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTATACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCaa                                                               <  8:31577/149‑1 (MQ=255)
                                                                                   cGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCaa                                                               >  6:267491/1‑148 (MQ=255)
                                                                                    ggATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCCACACGTTTTATCCGTTCTGGACTTCACCCACTAACCAAcgcg                                                                        >  1:63197/1‑138 (MQ=255)
                                                                                    ggATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCACTAACCAAcgcg                                                                        <  2:63197/138‑1 (MQ=255)
                                                                                       tGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAAcc                                                          <  6:141186/149‑1 (MQ=255)
                                                                                         cGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAAcccc                                                        <  2:267279/149‑1 (MQ=255)
                                                                                          caCGCTCGCGCTTCTAATCCGACCTACTTGGACTCATGTGTATTCCGGATAAGGCGTTTACGCAGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAAccccc                                                       >  8:275966/2‑149 (MQ=38)
                                                                                             gCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGc                                                    <  3:102914/149‑1 (MQ=255)
                                                                                                ggCGCGTCTTATCCGCCCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGccc                                                  >  2:349590/1‑148 (MQ=255)
                                                                                                  cgcgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGccc                                                  >  4:99199/1‑146 (MQ=255)
                                                                                                     gTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAgaga                                            <  5:267687/149‑1 (MQ=255)
                                                                                                      tCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAgaga                                            >  3:160532/1‑148 (MQ=255)
                                                                                                                 cTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATc                                <  6:267005/149‑1 (MQ=255)
                                                                                                                  tACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCa                               >  2:12536/1‑149 (MQ=255)
                                                                                                                  tACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCa                               >  7:293695/1‑149 (MQ=255)
                                                                                                                           gCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCgg                       >  1:215959/1‑148 (MQ=255)
                                                                                                                           gCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCgg                       >  4:250691/1‑148 (MQ=255)
                                                                                                                           gCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCgg                       >  4:24628/1‑148 (MQ=255)
                                                                                                                             aTGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTa                     >  4:116260/1‑148 (MQ=255)
                                                                                                                                   ttgCCGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCtg              >  7:159471/3‑149 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCtgt             >  7:184856/1‑149 (MQ=255)
                                                                                                                                     gCCGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGTCCAAAAGCAACACGTATTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCCGCAATACCCCCCTGCCCCAGAGACAAACCGCCATCCCCCGCCGGAAAACTCtgtg            >  2:239685/1‑149 (MQ=255)
                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGGCTATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGa          <  7:131242/149‑1 (MQ=255)
                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACTCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTgg           >  8:137737/1‑148 (MQ=255)
                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTCCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTgg           >  7:296610/1‑148 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAAGCGCCATCACCCGCCGGTAAACTCTGTGGaaa        <  7:137781/149‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGaaa        <  7:156192/149‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGaaa        <  2:215959/149‑1 (MQ=255)
                                                                                                                                           tgcggCGTAAACGCCTTATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAaga      <  8:159471/146‑1 (MQ=255)
                                                                                                                                              ggCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAACCTCTGTGGAAAGAGc    >  2:56831/1‑148 (MQ=255)
                                                                                                                                               gCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCaa  >  5:29250/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTGGTCAGTGAGGATCGGATGAAAGAGGCATTTGGATTGTTGAAAACATTGCCGATGTAGGTGGGCTACTGTGCCTAAAATGTCGGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAA  >  NC_000913/2834986‑2835277

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑