Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,867,151 A→G 27.2% F134S (TTT→TCT)  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,1510AG27.2% 252.4 / 117.3 169F134S (TTT→TCT) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base A (64/59);  new base G (24/22);  total (88/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 174 total aligned reads displayed.

AAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCA  >  NC_000913/2867052‑2867250
                                                                                                   |                                                                                                   
     aaTAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTtc                                                                                                >  1:1190893/1‑100 (MQ=255)
     aaTAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTtc                                                                                                >  1:2535663/1‑100 (MQ=255)
       tAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACGtctc                                                                                              <  2:131696/100‑1 (MQ=255)
        aGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTctct                                                                                             >  2:1206757/1‑100 (MQ=255)
         gTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTa                                                                                            <  1:1084811/100‑1 (MQ=255)
         gTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTa                                                                                            >  1:2330700/1‑100 (MQ=255)
            cGGGTTTGGTTCATAATCGCCCGTTCGATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcg                                                                                         >  2:1696860/1‑100 (MQ=255)
             gggTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgt                                                                                        >  1:837798/1‑99 (MQ=255)
              ggTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACcgcg                                                                                       <  2:3195484/100‑1 (MQ=255)
                  tGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGa                                                                                   <  2:284201/100‑1 (MQ=255)
                  tGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCgg                                                                                    >  1:2503593/1‑99 (MQ=255)
                   ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGAt                                                                                  >  2:1716963/1‑100 (MQ=255)
                   ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGAt                                                                                  >  1:2782866/1‑100 (MQ=255)
                   ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGAt                                                                                  >  1:2852988/1‑100 (MQ=255)
                   ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGAt                                                                                  >  1:88606/1‑100 (MQ=255)
                   ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGAt                                                                                  >  2:1270890/1‑100 (MQ=255)
                     ttCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                 <  1:2108127/99‑1 (MQ=255)
                      tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAg                                                                               >  2:1679134/1‑100 (MQ=255)
                       cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                              <  1:137584/100‑1 (MQ=255)
                        aTAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGcc                                                                             >  2:725081/1‑100 (MQ=255)
                        aTAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGcc                                                                             >  1:218947/1‑100 (MQ=255)
                               cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAgagatcggaa                                                                      >  2:3314399/1‑91 (MQ=255)
                                ccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTg                                                                     >  2:202942/1‑100 (MQ=255)
                                ccGTTCAAGCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCGCTACCGTGCGGATCAGCCCCAGGTTg                                                                     <  1:1995100/100‑1 (MQ=255)
                                 cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGc                                                                    <  1:633758/100‑1 (MQ=255)
                                 cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGc                                                                    <  1:3324903/100‑1 (MQ=255)
                                 cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTgn                                                                    <  1:127582/100‑2 (MQ=255)
                                  gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCCCGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGcc                                                                   >  1:1504298/1‑100 (MQ=255)
                                  gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGcc                                                                   >  2:3005239/1‑100 (MQ=255)
                                   ttCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGccc                                                                  <  1:2358858/100‑1 (MQ=255)
                                      aaTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGccn                                                                  <  1:528169/97‑2 (MQ=255)
                                        tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTc                                                             <  1:592342/100‑1 (MQ=255)
                                        tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTc                                                             <  2:3024958/100‑1 (MQ=255)
                                        tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTc                                                             <  1:808582/100‑1 (MQ=255)
                                        tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTc                                                             >  2:2180704/1‑100 (MQ=255)
                                          gTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                           >  2:413832/1‑100 (MQ=255)
                                            cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATa                                                         >  2:2939032/1‑100 (MQ=255)
                                            cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATa                                                         <  1:3231164/100‑1 (MQ=255)
                                              ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                       >  2:300649/1‑100 (MQ=255)
                                               gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                      <  2:2330717/100‑1 (MQ=255)
                                                cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                     >  1:1462020/1‑100 (MQ=255)
                                                cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                     >  2:855406/1‑100 (MQ=255)
                                                cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                     >  2:3059106/1‑100 (MQ=255)
                                                  aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTcc                                                   <  1:444562/100‑1 (MQ=255)
                                                   aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCa                                                  >  2:616310/1‑100 (MQ=255)
                                                   aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCa                                                  >  2:2102328/1‑100 (MQ=255)
                                                    tCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCa                                                  >  2:1339202/1‑99 (MQ=255)
                                                    tCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAg                                                 <  2:3327739/100‑1 (MQ=255)
                                                    tCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAg                                                 <  1:67644/100‑1 (MQ=255)
                                                     ctaccaGGTTGCGCATGTTGAGAAGCGGAAAACACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                >  2:2283037/3‑100 (MQ=255)
                                                     ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                >  1:2603406/1‑100 (MQ=255)
                                                     ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                >  2:1874578/1‑100 (MQ=255)
                                                      caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                               >  2:3088900/1‑100 (MQ=255)
                                                          aatttGCGTATGATGAGAGGCGGAGGCCGCGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCACGTGGCACTCTTCGATAAGGTcc                                                   <  1:2860631/89‑1 (MQ=255)
                                                          aGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGc                                           >  1:2044566/1‑100 (MQ=255)
                                                          aGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGGCAATCTTCTCTACCGTGCGGATCAGCCCCTGGTTGCTCTCTTCGATAAGGTCCAGCAACg                                            <  1:621424/99‑1 (MQ=255)
                                                           ggTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGcc                                          <  2:1548037/100‑1 (MQ=255)
                                                               gCGTATGTTGGCAAGCGGACACCGCCTGCCGAGTTAAACCGGTTAACTGGGGGGATCAGCCCCAGGTTACCCTCGTCCATAAGGTCCAGCTTCCCCAGAc                                      <  1:1756345/100‑1 (MQ=255)
                                                                cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                     <  2:442945/100‑1 (MQ=255)
                                                                cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                     >  2:2565448/1‑100 (MQ=255)
                                                                cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                     <  1:2411796/100‑1 (MQ=255)
                                                                 gTATGTTGAGAAGCGGAAACCAAGCTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGACAGACCa                                    >  1:1707653/1‑100 (MQ=255)
                                                                  tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCAc                                   <  1:1716952/100‑1 (MQ=255)
                                                                  tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCAc                                   >  2:16004/1‑100 (MQ=255)
                                                                    tGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                 >  1:2703017/1‑100 (MQ=255)
                                                                    tGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                 >  1:1675073/1‑100 (MQ=255)
                                                                       tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                              >  1:1601986/1‑100 (MQ=255)
                                                                       tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                              <  2:1177939/100‑1 (MQ=255)
                                                                       tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                              >  2:1206439/1‑100 (MQ=255)
                                                                        gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                             >  2:1827706/1‑100 (MQ=255)
                                                                         ngaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGcc                            >  1:559503/2‑100 (MQ=255)
                                                                          gaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTACCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                           >  1:1692510/1‑100 (MQ=255)
                                                                          gaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                           <  1:1472071/100‑1 (MQ=255)
                                                                            aGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATa                         >  1:1895773/1‑100 (MQ=255)
                                                                             gCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATCAGGTCCAGCAACGCCAGACCACGATTGCCATaa                        >  1:2401653/1‑100 (MQ=255)
                                                                               ggAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                      >  1:1830117/1‑100 (MQ=255)
                                                                                gAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcgg                     >  2:193047/1‑100 (MQ=255)
                                                                                 aaaCCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                    >  2:745716/1‑100 (MQ=255)
                                                                                 aaaCCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                    <  1:927756/100‑1 (MQ=255)
                                                                                 aaaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                    <  2:1982165/100‑1 (MQ=255)
                                                                                  aaCCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                   <  2:2367549/100‑1 (MQ=255)
                                                                                  aaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                   <  2:1074734/100‑1 (MQ=255)
                                                                                  aaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                    >  1:102097/1‑99 (MQ=255)
                                                                                   aCCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATTAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcgg                  <  2:3224683/100‑1 (MQ=255)
                                                                                    ccACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                 >  2:936948/1‑100 (MQ=255)
                                                                                     cACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                <  2:2528480/100‑1 (MQ=255)
                                                                                      aCGTTCCGGGTCAAACTTCTATACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACAACGATTGCCATAACGGCGGGc                <  1:918491/99‑1 (MQ=255)
                                                                                        gTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAt             >  2:1355255/1‑100 (MQ=255)
                                                                                          tCCGGGTTAAACTTCTCTACCGCGCGGTTCAGCCACAGGTTGCCCCCTTCGATAAGGTCCAGCAACGCCAGACCACGCTTGCCATCACGGCGGGCAAtt            <  1:2196813/99‑1 (MQ=255)
                                                                                           ccGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAtttt          >  1:1582745/1‑100 (MQ=255)
                                                                                           ccGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAt             <  1:1178380/97‑1 (MQ=255)
                                                                                            cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt         <  2:1902664/100‑1 (MQ=255)
                                                                                             gggTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa        <  1:1134375/100‑1 (MQ=255)
                                                                                             gggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa        <  2:2603426/100‑1 (MQ=255)
                                                                                             gggTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa        >  1:330927/1‑100 (MQ=255)
                                                                                               gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc      <  1:847833/100‑1 (MQ=255)
                                                                                               gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc      <  1:3088874/100‑1 (MQ=255)
                                                                                                 cAAACTTCTCTACCGTGCTGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCAGAACGGCGGGCAATTTTTaccac    <  1:1530338/100‑1 (MQ=255)
                                                                                                  aaaCTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc   >  1:292966/1‑100 (MQ=255)
                                                                                                   aaCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAccacca  >  1:3347634/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCA  >  NC_000913/2867052‑2867250

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