Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,867,164 C→T 62.9% 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,1640CT62.9% 120.8 / 173.0 165A130T (GCG→ACG) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base C (37/24);  new base T (52/51);  total (89/76)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 168 total aligned reads displayed.

GGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTA  >  NC_000913/2867066‑2867261
                                                                                                  |                                                                                                 
 gTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgc                                                                                                 <  2:2231889/100‑1 (MQ=255)
   ttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCgg                                                                                               <  1:1315314/100‑1 (MQ=255)
    tGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCgg                                                                                               >  1:2503593/1‑99 (MQ=255)
    tGGTTAATAATCGCCCGTCTAATCGACTGGCGAATCCACCAGGTTGCGTATGTTGAGAGGCGGAAGCCACGTTCCGTGTCAAACTTCTCTACTGTGCGGa                                                                                              <  1:2014727/100‑1 (MQ=255)
     ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGAt                                                                                             >  2:1716963/1‑100 (MQ=255)
     ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGAt                                                                                             >  1:88606/1‑100 (MQ=255)
      gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                            <  1:2788977/100‑1 (MQ=255)
      gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                            >  1:1692616/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)
         cATAATCGCCCGTTCAATCGTCTGGCAAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                                         <  1:3347139/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)
                 cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGtt                                                                                 >  1:1752389/1‑100 (MQ=255)
                 cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGtt                                                                                 >  1:176015/1‑100 (MQ=255)
                  ccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTg                                                                                >  2:202942/1‑100 (MQ=255)
                   cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGc                                                                               <  1:633758/100‑1 (MQ=255)
                   cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGc                                                                               <  1:3324903/100‑1 (MQ=255)
                    gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCCCGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGcc                                                                              >  1:1504298/1‑100 (MQ=255)
                    gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGcc                                                                              <  1:234228/100‑1 (MQ=255)
                    gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGcc                                                                              >  2:3021593/1‑100 (MQ=255)
                    gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGcc                                                                              >  2:2546849/1‑100 (MQ=255)
                    gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGcc                                                                              >  2:3005239/1‑100 (MQ=255)
                       cAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCCAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGccc                                                                             >  1:2729814/1‑98 (MQ=255)
                        aaTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGccn                                                                             <  1:528169/97‑2 (MQ=255)
                          tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTc                                                                        <  2:3024958/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)
                                ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCCTGTTCCGGGGCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTTCCCTCTTCCATAAg                                                                  >  1:3017844/1‑100 (MQ=255)
                                ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                                  >  1:700798/1‑100 (MQ=255)
                                 gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                                 <  2:2330717/100‑1 (MQ=255)
                                  cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                                >  2:855406/1‑100 (MQ=255)
                                  cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                                >  1:1462020/1‑100 (MQ=255)
                                  cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                                <  2:1315017/100‑1 (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)
                                       ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                           >  1:2603406/1‑100 (MQ=255)
                                        caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                          >  2:2354742/1‑100 (MQ=255)
                                        caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                          >  2:3088900/1‑100 (MQ=255)
                                        caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                          >  1:897527/1‑100 (MQ=255)
                                           caGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACAGTGCGCATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCCAc                                                        >  1:1483499/1‑99 (MQ=255)
                                            aatttGCGTATGATGAGAGGCGGAGGCCGCGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCACGTGGCACTCTTCGATAAGGTcc                                                              <  1:2860631/89‑1 (MQ=255)
                                            aGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGc                                                      <  1:2906161/100‑1 (MQ=255)
                                            aGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGc                                                      >  1:2044566/1‑100 (MQ=255)
                                               ttGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAg                                                   <  1:616307/100‑1 (MQ=255)
                                                tGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGa                                                  >  2:1925021/1‑100 (MQ=255)
                                                  cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                                <  1:2411796/100‑1 (MQ=255)
                                                  cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                                >  2:2565448/1‑100 (MQ=255)
                                                  cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                                <  2:2875735/100‑1 (MQ=255)
                                                   gTATGTTGAGAAGCGGAAACCAAGCTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGACAGACCa                                               >  1:1707653/1‑100 (MQ=255)
                                                    tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCAc                                              <  2:1026408/100‑1 (MQ=255)
                                                     aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACAGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACg                                             >  1:109835/1‑100 (MQ=255)
                                                      tGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                            >  1:1675073/1‑100 (MQ=255)
                                                      tGTTGAGAAGCGGAAACCACGTTCCGGGTCAGACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                            >  1:2703017/1‑100 (MQ=255)
                                                         tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                                         >  1:1601986/1‑100 (MQ=255)
                                                         tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                                         >  2:1206439/1‑100 (MQ=255)
                                                         tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                                         <  2:1177939/100‑1 (MQ=255)
                                                          gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                        >  2:1827706/1‑100 (MQ=255)
                                                          gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                        <  2:2640757/100‑1 (MQ=255)
                                                           ngaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGcc                                       >  1:559503/2‑100 (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)
                                                                    aaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                              <  2:1567161/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)
                                                                      ccACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                            <  2:1863083/100‑1 (MQ=255)
                                                                       cACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                           <  2:2436282/100‑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)
                                                                              cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCaa                         >  1:3119378/1‑95 (MQ=255)
                                                                              cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt                    <  2:1902664/100‑1 (MQ=255)
                                                                               gggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                   <  2:1604489/100‑1 (MQ=255)
                                                                               gggTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                   >  1:330927/1‑100 (MQ=255)
                                                                               gggTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                   <  2:417673/100‑1 (MQ=255)
                                                                                 gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc                 <  1:847833/100‑1 (MQ=255)
                                                                                 gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc                 >  2:1241671/1‑100 (MQ=255)
                                                                                   cAAACTTCTCTACCGTGCTGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCAGAACGGCGGGCAATTTTTaccac               <  1:1530338/100‑1 (MQ=255)
                                                                                     gaCTTCTCTCCCGCTCGGATCAGCCCCAGGTTGCCATCTCCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTaccacc              <  1:814844/98‑1 (MQ=255)
                                                                                        ttCTCTACCGTGCGGATGAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc          <  1:3380223/100‑1 (MQ=255)
                                                                                           tctACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCa       >  2:3317867/1‑100 (MQ=255)
                                                                                            ctACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCaa      >  1:735878/1‑100 (MQ=255)
                                                                                              aCCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt    >  2:319215/1‑100 (MQ=255)
                                                                                              aCCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt    >  1:600552/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
GGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTA  >  NC_000913/2867066‑2867261

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: ‑