Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,867,164 C→T 100% A130T (GCG→ACG)  rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,1640CT97.4% 363.6 / ‑0.3 117A130T (GCG→ACG) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (65/49);  total (66/51)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

Only 100 of 120 total aligned reads displayed.

GTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTT  >  NC_000913/2867067‑2867260
                                                                                                 |                                                                                                
gTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgc                                                                                                >  1:1662531/1‑100 (MQ=255)
 tttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgcg                                                                                               >  2:2499680/1‑100 (MQ=255)
    ggTTCATAATCGCCCAGTCAATCGTCTGGCGAATCCACTAGGTTACGTATGTTGAGAAGCGGAAACCAAATTCCGGGTAAAACTTCTATAAcgggcg                                                                                               >  2:287359/1‑97 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                           <  1:1052106/100‑1 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                           >  2:907530/1‑100 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                           <  2:2213192/100‑1 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                           <  2:1120395/100‑1 (MQ=255)
      ttCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCa                                                                                          <  1:1873197/100‑1 (MQ=255)
       tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGAt                                                                                            >  2:2436487/1‑97 (MQ=255)
       tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAg                                                                                         >  2:910941/1‑100 (MQ=255)
        cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                                        <  1:2542976/100‑1 (MQ=255)
        cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                                        >  2:2463229/1‑100 (MQ=255)
         aTAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGcc                                                                                       <  1:713480/100‑1 (MQ=255)
          tAATCTCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGccc                                                                                      >  2:54334/1‑100 (MQ=255)
          tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAAAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGccc                                                                                      >  1:812356/1‑100 (MQ=255)
          tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATCTTGAGAAGCGGACACCACGTTCCGGGTCAAACTTCTCTGCCGTGCGGATCAGccc                                                                                      >  2:1769172/1‑100 (MQ=255)
            aTCGCCCGTTCAATCGTCTGGTGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCa                                                                                    <  1:614934/100‑1 (MQ=255)
             tCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAg                                                                                   >  1:2904896/1‑100 (MQ=255)
             tCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAg                                                                                   >  1:3044469/1‑100 (MQ=255)
              cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAgg                                                                                  >  1:2861377/1‑100 (MQ=255)
              cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAgg                                                                                  >  1:2012147/1‑100 (MQ=255)
              cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAgg                                                                                  <  2:1109032/100‑1 (MQ=255)
               gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGt                                                                                 >  2:2493622/1‑100 (MQ=255)
               gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGt                                                                                 >  1:2928860/1‑100 (MQ=255)
               gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGGATGTTGAGAAGCGGAAACCACGTTCGGGGGCAAACTTCTCTACCGTGCGGATCAGCCCCAGGt                                                                                 >  2:2271967/1‑100 (MQ=255)
                  cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGc                                                                              <  2:1077153/100‑1 (MQ=255)
                    ttCAATCGTCTGGCGAATCCTCCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGccc                                                                            <  2:531519/100‑1 (MQ=255)
                     tCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCct                                                                           <  1:1165536/100‑1 (MQ=255)
                      cAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCctc                                                                          >  1:2247430/1‑100 (MQ=255)
                        aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCtt                                                                        >  2:772470/1‑100 (MQ=255)
                          cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                                      <  1:1769330/100‑1 (MQ=255)
                             cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCACAGGTTGCCCTCTTCGa                                                                     >  2:1118611/1‑98 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                                 >  2:3024554/1‑100 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                                 >  2:11253/1‑100 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                                 >  1:2741500/1‑100 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                                 <  1:2709004/100‑1 (MQ=255)
                                gCGAATCCATCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                                <  1:1901856/100‑1 (MQ=255)
                                gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                                <  2:885614/100‑1 (MQ=255)
                                 cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCCCTTCGATAAGGt                                                               >  2:1226783/1‑100 (MQ=255)
                                   aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTcc                                                             <  1:1610456/100‑1 (MQ=255)
                                   aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTcc                                                             <  2:1849303/100‑1 (MQ=255)
                                   aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTcc                                                             <  2:2012124/100‑1 (MQ=255)
                                       caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                          <  1:1128136/99‑1 (MQ=255)
                                       caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                         >  1:1022520/1‑100 (MQ=255)
                                       caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                         >  2:123131/1‑100 (MQ=255)
                                       caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                         >  1:1373187/1‑100 (MQ=255)
                                        nccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCaa                                                        >  2:1506308/2‑100 (MQ=255)
                                         ccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAAc                                                       <  1:2660010/100‑1 (MQ=255)
                                           aGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTACGGATCAGCCCCAGGTCGCCCTCTTCGATAAGGTCCAGCAACGc                                                     >  1:1910097/1‑100 (MQ=255)
                                            ggTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGcc                                                    <  1:1578792/100‑1 (MQ=255)
                                            ggTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGcc                                                    <  1:496008/100‑1 (MQ=255)
                                            ggTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATCAGGTCCAGCAACGcc                                                    <  2:2156086/100‑1 (MQ=255)
                                            ggTTGCGTATGTTGAGAAGCAGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGAGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGc                                                     >  2:1097909/1‑99 (MQ=255)
                                                gCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAc                                                <  1:2385978/100‑1 (MQ=255)
                                                 cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                               >  2:2031081/1‑100 (MQ=255)
                                                   tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCAc                                             <  1:1460626/100‑1 (MQ=255)
                                                     tGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTTCTGATCAGCCCAAAGTTACCCTCTTCGATAAAGTCCAGCAACGCCAGACCACGa                                           >  1:2624794/1‑100 (MQ=255)
                                                     tGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                           >  2:2724226/1‑100 (MQ=255)
                                                      gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTGCAGCAACGCCAGACCACGAt                                          >  2:40602/1‑100 (MQ=255)
                                                       ttGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAtt                                         >  2:1552979/1‑100 (MQ=255)
                                                       ttGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAtt                                         >  2:1640044/1‑100 (MQ=255)
                                                         gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                       <  1:2001045/100‑1 (MQ=255)
                                                         gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                       <  1:195991/100‑1 (MQ=255)
                                                         gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                       <  1:2192655/100‑1 (MQ=255)
                                                              gCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCTCGATTGCCATaa                                  >  2:2735864/1‑100 (MQ=255)
                                                              gCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATaa                                  >  1:2795052/1‑100 (MQ=255)
                                                              gCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATaa                                  <  1:2526125/100‑1 (MQ=255)
                                                               cGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAc                                 <  1:123132/100‑1 (MQ=255)
                                                                 gAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                                >  1:1803290/1‑99 (MQ=255)
                                                                 gAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGATCCAGCAACGCCAGACCACGATTGCCATAAcgg                               >  1:3619/1‑100 (MQ=255)
                                                                    aCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcgg                            >  2:2691110/1‑100 (MQ=255)
                                                                     ccACGTTCCGGGTCAAACTTCTCTACCGTGCGTATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                           <  2:2325570/100‑1 (MQ=255)
                                                                     ccACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                           >  2:505503/1‑100 (MQ=255)
                                                                     ccACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                           >  1:2352133/1‑100 (MQ=255)
                                                                     ccACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                           <  2:2763386/100‑1 (MQ=255)
                                                                        cGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCaa                        >  1:1877274/1‑100 (MQ=255)
                                                                             cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt                   >  2:2496160/1‑100 (MQ=255)
                                                                             cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt                   >  2:1707619/1‑100 (MQ=255)
                                                                              gggTCAAACTTCTCTACCGTGCTGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                  <  1:413151/100‑1 (MQ=255)
                                                                              gggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                  >  1:457584/1‑100 (MQ=255)
                                                                                gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACAACGACTGCCATAACGGCGGGCAATTTTTacc                >  2:416270/1‑100 (MQ=255)
                                                                                 tCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacca               >  2:795232/1‑100 (MQ=255)
                                                                                  cAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCTGGCAATTTTTaccac              >  1:2996015/1‑100 (MQ=255)
                                                                                  cAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccac              >  2:2908155/1‑100 (MQ=255)
                                                                                   aaaCTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc             <  1:1500477/100‑1 (MQ=255)
                                                                                       ttCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc         >  2:1781826/1‑100 (MQ=255)
                                                                                       ttCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc         >  2:2630059/1‑100 (MQ=255)
                                                                                       ttCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc         <  1:2553104/100‑1 (MQ=255)
                                                                                        tctcTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACg        <  2:2150399/100‑1 (MQ=255)
                                                                                        tctcTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACg        <  2:1143252/100‑1 (MQ=255)
                                                                                         ntctACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGc       >  1:893078/2‑100 (MQ=255)
                                                                                           ctACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCaa     >  1:2260941/1‑100 (MQ=255)
                                                                                           ctACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCaa     >  1:608732/1‑100 (MQ=255)
                                                                                            tACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg    <  1:1763192/100‑1 (MQ=255)
                                                                                             aCCGTGCGGATCAGCGCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg    <  2:1337591/99‑1 (MQ=255)
                                                                                             aCCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt   <  2:608728/100‑1 (MQ=255)
                                                                                             aCCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt   <  1:505506/100‑1 (MQ=255)
                                                                                             aCCGTGCGGATAAGCACCAAGTTAACATCTTGGATAAGGTCCAGCAACGCCAGACCACGCTTGCTATCACGGCCGGTAATTTTTACTACCATGCGCaa     >  2:2634279/1‑98 (MQ=255)
                                                                                              ccGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGtt  >  2:1527415/1‑100 (MQ=255)
                                                                                              ccGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGtt  <  1:110188/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTT  >  NC_000913/2867067‑2867260

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