Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 263,172 A→C 100% pseudogene (41/63 nt) ykfN ← pseudogene, YdiA family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913263,1720AC96.7% 631.3 / 7.4 185pseudogene (41/63 nt)ykfNpseudogene, YdiA family
Reads supporting (aligned to +/- strand):  ref base A (1/5);  new base C (94/82);  total (97/88)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

Only 100 of 185 total aligned reads displayed.

GGCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAC  >  NC_000913/263075‑263267
                                                                                                 |                                                                                               
ggCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGc                                                                                               >  1:1005453/1‑100 (MQ=255)
ggCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGc                                                                                               <  1:2190789/100‑1 (MQ=255)
 gCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCt                                                                                              >  2:3849920/1‑100 (MQ=255)
   ttCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCa                                                                                            <  2:3319891/100‑1 (MQ=255)
      gCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTa                                                                                         >  1:3597064/1‑100 (MQ=255)
      gCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTa                                                                                         <  2:3464966/100‑1 (MQ=255)
       cTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTAt                                                                                        <  2:2734119/100‑1 (MQ=255)
        tCACTATACATCCTTGGCTGCAGGTTTAGTTGTCCACCACTCCTAAATTTAATGTGTTGGCACTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATg                                                                                       >  1:305582/1‑100 (MQ=255)
           cgataCATCCTTGGCTGCAGCTTTAGTGGTACACCACTCCGAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATc                                                                                    <  2:1677293/98‑1 (MQ=255)
           cTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGAt                                                                                     >  2:3199777/1‑99 (MQ=255)
           cTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATc                                                                                    <  1:1040059/100‑1 (MQ=255)
             cttccgatctGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCgg                                                                                  <  1:685653/91‑1 (MQ=255)
             ataCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCgg                                                                                  <  1:2518222/100‑1 (MQ=255)
              taCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGt                                                                                 <  2:858374/100‑1 (MQ=255)
              taCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGt                                                                                 <  1:3334405/100‑1 (MQ=255)
               aCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGtt                                                                                <  1:1716682/100‑1 (MQ=255)
               aCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGtt                                                                                <  1:3191746/100‑1 (MQ=255)
                naTCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTa                                                                               >  2:2113295/2‑100 (MQ=255)
                cATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTa                                                                               >  2:3141214/1‑100 (MQ=255)
                 aTCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTaa                                                                              >  2:104314/1‑100 (MQ=255)
                  tCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATATGTTTGCAATGTATTCAATGAAGCAGGAACAAATTTGCGCATTATGAGCggagaat                                                                             >  1:3386863/1‑95 (MQ=255)
                    cTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATa                                                                            >  1:264281/1‑99 (MQ=255)
                    cTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATAc                                                                           <  1:1613966/100‑1 (MQ=255)
                    cTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATAc                                                                           <  1:3346438/100‑1 (MQ=255)
                    cTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATAc                                                                           <  2:644300/100‑1 (MQ=255)
                     ttGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACt                                                                          <  2:2935485/100‑1 (MQ=255)
                     cttgCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACt                                                                          <  2:1192621/97‑1 (MQ=255)
                       ggCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCggagatcggaag                                                                        >  2:685658/1‑90 (MQ=255)
                       ggCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTc                                                                        <  2:527525/100‑1 (MQ=255)
                        gCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGTCAATGTGTCCAATAAAGCTCGAACAAACTCGCTCATTATGATAGGTTAATAc                                                                           >  1:2695867/1‑96 (MQ=255)
                        gCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCa                                                                       <  2:3676818/100‑1 (MQ=255)
                         cTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCaa                                                                      >  2:1937567/1‑100 (MQ=255)
                         cTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCaa                                                                      >  2:245392/1‑100 (MQ=255)
                          tGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTTGCTCATTATGATCGGTTAATACTTCAAc                                                                     >  1:1247122/1‑100 (MQ=255)
                          tGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAAc                                                                     >  2:2679273/1‑100 (MQ=255)
                               gTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTg                                                                <  2:2172619/100‑1 (MQ=255)
                               gTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTg                                                                >  2:3145862/1‑100 (MQ=255)
                               gTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTg                                                                >  2:924242/1‑100 (MQ=255)
                                tttAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTgg                                                               <  2:2219960/100‑1 (MQ=255)
                                 ttAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGt                                                              <  1:638524/100‑1 (MQ=255)
                                  tAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCCGGtt                                                             >  1:3497264/1‑100 (MQ=255)
                                   aGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTg                                                            <  1:3103293/100‑1 (MQ=255)
                                   aGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTg                                                            <  1:2364394/100‑1 (MQ=255)
                                   aGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTg                                                            >  1:3411896/1‑100 (MQ=255)
                                   aGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAACCTCAAACAAATTCGCTCATTATGATCGGTTAATACTTCAACATCTGGTTg                                                            >  1:722415/1‑100 (MQ=255)
                                     ttGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCa                                                          <  2:1308900/100‑1 (MQ=255)
                                      tGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCAt                                                         <  2:1343306/100‑1 (MQ=255)
                                        tACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGa                                                       >  1:2456192/1‑100 (MQ=255)
                                        tACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGa                                                       <  2:1296731/100‑1 (MQ=255)
                                         acacCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAt                                                      >  2:2388092/1‑100 (MQ=255)
                                          caccacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAtt                                                     >  2:3393135/1‑100 (MQ=255)
                                          caccacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAtt                                                     >  1:312029/1‑100 (MQ=255)
                                            ccacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgt                                                   >  2:327860/1‑100 (MQ=255)
                                             cacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtt                                                  <  1:1433381/100‑1 (MQ=255)
                                             cacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtt                                                  <  1:2171113/100‑1 (MQ=255)
                                               cTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtttg                                                <  1:968819/100‑1 (MQ=255)
                                                tCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtttgt                                               <  1:1346628/100‑1 (MQ=255)
                                                 ccTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTc                                              >  2:623757/1‑100 (MQ=255)
                                                 ccTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTc                                              >  1:1037610/1‑100 (MQ=255)
                                                  cTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTcc                                             <  1:2408411/100‑1 (MQ=255)
                                                   tAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCg                                            >  2:1821449/1‑100 (MQ=255)
                                                   tAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCg                                            <  1:3393108/100‑1 (MQ=255)
                                                          aaTGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAg                                     >  2:2535056/1‑100 (MQ=255)
                                                                ttGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAcg                               <  1:2801604/100‑1 (MQ=255)
                                                                   gCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACgcgc                            >  2:3509941/1‑100 (MQ=255)
                                                                    cAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGcc                           <  1:3769898/100‑1 (MQ=255)
                                                                     aaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCt                          >  2:3632260/1‑100 (MQ=255)
                                                                       tgtgTTCAACAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGc                        >  2:1367770/1‑100 (MQ=255)
                                                                         tgtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCg                      >  2:1744486/1‑100 (MQ=255)
                                                                          gtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCgg                     >  1:3840341/1‑100 (MQ=255)
                                                                          gtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCgg                     <  2:3553148/100‑1 (MQ=255)
                                                                          gtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCgg                     >  2:2588364/1‑100 (MQ=255)
                                                                          gtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCgg                     <  2:2268083/100‑1 (MQ=255)
                                                                          gtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCgg                     >  1:2143411/1‑100 (MQ=255)
                                                                          gtTCAATAAAGCTCGAACAAATTCGCCCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCgg                     >  1:3284539/1‑100 (MQ=255)
                                                                           ttCAATAAAGCTCGAACAAATTCGCTCATTATGATGGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCggg                    <  1:969176/100‑1 (MQ=255)
                                                                           ttCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCggg                    <  2:1648921/100‑1 (MQ=255)
                                                                           ttCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCggg                    <  2:890330/100‑1 (MQ=255)
                                                                            tCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgt                   <  1:3692829/100‑1 (MQ=255)
                                                                            tCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgt                   <  1:1139366/100‑1 (MQ=255)
                                                                             cAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgta                  >  1:141458/1‑100 (MQ=255)
                                                                             cAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgta                  >  2:199701/1‑100 (MQ=255)
                                                                               aTAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgtagt                <  1:3297568/100‑1 (MQ=255)
                                                                                tAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGctactacccg               <  1:1775068/100‑10 (MQ=255)
                                                                                 aaaGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCAGGGtagtag              >  1:2885916/1‑100 (MQ=255)
                                                                                   aGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCa            <  2:1464056/100‑1 (MQ=255)
                                                                                     cTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAgg          >  2:3866188/1‑100 (MQ=255)
                                                                                     cTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAgg          >  2:480480/1‑100 (MQ=255)
                                                                                     cTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAgg          >  1:1856318/1‑100 (MQ=255)
                                                                                      tCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGgc         <  2:2438265/100‑1 (MQ=255)
                                                                                      tCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGgc         >  2:1087211/1‑100 (MQ=255)
                                                                                      tCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGgc         <  2:130596/100‑1 (MQ=255)
                                                                                       cGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCa            >  1:491633/1‑96 (MQ=255)
                                                                                       cGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGgcg        >  2:3709885/1‑100 (MQ=255)
                                                                                       cGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGgcg        >  2:1173582/1‑100 (MQ=255)
                                                                                        gAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCAAGATTGTGTGTCCGTAAAAAGATAACGCGCCCGCCGGGTAATAGCAgg          >  1:977400/1‑97 (MQ=255)
                                                                                         aaCAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTACCATGCGCa      >  1:2584258/1‑100 (MQ=255)
                                                                                          aCAAATTCGCTCATTATGATCGGTTAATACTTCAACTTTTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCAt     >  1:2386896/1‑100 (MQ=255)
                                                                                          aCAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCAt     >  1:260725/1‑100 (MQ=255)
                                                                                           cAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCAtt    >  2:837125/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GGCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAC  >  NC_000913/263075‑263267

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