Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,164 C→T 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.9% 315.1 / ‑2.4 93A130T (GCG→ACG) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (52/39);  total (54/39)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAG  >  NC_000913/2867067‑2867258
                                                                                                 |                                                                                              
gTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgc                                                                                              <  1:672651/100‑1 (MQ=255)
gTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgc                                                                                              >  1:502971/1‑100 (MQ=255)
 tttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgcg                                                                                             <  1:1266829/100‑1 (MQ=255)
  ttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCgg                                                                                            <  2:1221615/100‑1 (MQ=255)
  ttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCgg                                                                                            >  2:1186838/1‑100 (MQ=255)
  ttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCgg                                                                                            >  2:428375/1‑100 (MQ=255)
  ttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCgg                                                                                            <  2:186353/100‑1 (MQ=255)
   tGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGa                                                                                           >  1:1192484/1‑100 (MQ=255)
   tGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGa                                                                                           <  1:794061/100‑1 (MQ=255)
    ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGAt                                                                                          <  1:229180/100‑1 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGCATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATc                                                                                         >  2:1544008/1‑100 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                         >  2:1057673/1‑100 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                         >  1:939868/1‑100 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                         >  1:8484/1‑100 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                         >  2:619541/1‑100 (MQ=255)
       tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAg                                                                                       >  2:732550/1‑100 (MQ=255)
        cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                                      >  1:1095132/1‑100 (MQ=255)
        cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                                      >  1:693159/1‑100 (MQ=255)
        cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                                      >  1:61610/1‑100 (MQ=255)
         aTAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGcc                                                                                     <  1:1046569/100‑1 (MQ=255)
           aaTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGcccc                                                                                   >  2:275117/1‑100 (MQ=255)
           aaTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGcccc                                                                                   >  1:577394/1‑100 (MQ=255)
            tctgCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCa                                                                                  <  1:603162/97‑1 (MQ=255)
              cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAgg                                                                                >  2:300536/1‑100 (MQ=255)
              cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAgg                                                                                >  1:615198/1‑100 (MQ=255)
               gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCaaga                                                                               >  2:603160/1‑97 (MQ=255)
               gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGt                                                                               >  2:513348/1‑100 (MQ=255)
               gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGt                                                                               >  1:1330017/1‑100 (MQ=255)
                 ccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTg                                                                             <  2:1495993/100‑1 (MQ=255)
                   gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGcc                                                                           <  1:917254/100‑1 (MQ=255)
                    ttCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGccc                                                                          >  1:1325641/1‑100 (MQ=255)
                        aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCtt                                                                      >  1:1086689/1‑100 (MQ=255)
                         tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTc                                                                     >  2:261763/1‑100 (MQ=255)
                         cgatCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTc                                                                     <  1:620575/97‑1 (MQ=255)
                           gTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                                   >  1:500208/1‑100 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTcagatcg                                                               >  2:620573/1‑94 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                               >  2:370846/1‑100 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                               <  2:1440616/100‑1 (MQ=255)
                                gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                              <  2:426323/100‑1 (MQ=255)
                                  gAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTc                                                            >  2:162174/1‑100 (MQ=255)
                                   aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTcc                                                           >  1:1310914/1‑100 (MQ=255)
                                    aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCa                                                          <  1:693772/100‑1 (MQ=255)
                                    aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCa                                                          <  1:276979/100‑1 (MQ=255)
                                      ncaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGTCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                        >  2:729346/2‑100 (MQ=255)
                                       caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                       <  2:510250/100‑1 (MQ=255)
                                       caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                       >  1:113970/1‑100 (MQ=255)
                                       caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                       >  2:456171/1‑100 (MQ=255)
                                        accaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCaa                                                      >  1:33527/1‑100 (MQ=255)
                                        accaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGCTGCCCTCTTCGATAAGGTCCAGCaa                                                      >  2:1392385/1‑100 (MQ=255)
                                            ggTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGcc                                                  >  1:625257/1‑100 (MQ=255)
                                              ttGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAg                                                >  1:589886/1‑100 (MQ=255)
                                               tGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGa                                               <  1:1142828/100‑1 (MQ=255)
                                                gCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAc                                              <  1:75538/100‑1 (MQ=255)
                                                gCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAc                                              <  2:763157/100‑1 (MQ=255)
                                                 cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                             >  1:877545/1‑100 (MQ=255)
                                                  gTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCa                                            >  2:1092604/1‑100 (MQ=255)
                                                   tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCAc                                           <  1:24866/100‑1 (MQ=255)
                                                   tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCAc                                           >  1:1089069/1‑100 (MQ=255)
                                                    aTGTTGAGAAGCGGCAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACg                                          >  2:1242889/1‑100 (MQ=255)
                                                     tGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                         >  2:1314573/1‑100 (MQ=255)
                                                     tGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                         >  1:1424591/1‑100 (MQ=255)
                                                        tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                                      >  2:1543917/1‑100 (MQ=255)
                                                         gagaAGCGGAAACCACGTTCCGGGTCAAAGTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                     <  2:787537/100‑1 (MQ=255)
                                                         gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                     >  2:856998/1‑100 (MQ=255)
                                                         gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                     >  2:1244502/1‑100 (MQ=255)
                                                            aaGCGGAAACCACGTTCCGGGTAAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCATCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGACAt                                  <  2:1319328/100‑1 (MQ=255)
                                                             aGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATa                                 <  1:732552/100‑1 (MQ=255)
                                                                ngAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                              >  1:1509295/2‑100 (MQ=255)
                                                                   aaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                           <  2:1095127/100‑1 (MQ=255)
                                                                   aaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                           >  2:1035166/1‑100 (MQ=255)
                                                                   aaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                           <  1:332846/100‑1 (MQ=255)
                                                                    cgtggtttCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcgg                          >  2:1470009/7‑100 (MQ=255)
                                                                      cACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                        <  1:194115/100‑1 (MQ=255)
                                                                            ccGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAtttt                  <  1:1314580/100‑1 (MQ=255)
                                                                             cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt                 <  2:813713/100‑1 (MQ=255)
                                                                             cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt                 <  2:939863/100‑1 (MQ=255)
                                                                             cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt                 <  2:577392/100‑1 (MQ=255)
                                                                             cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCAATAACGGCGGACAAttttt                 >  2:50924/1‑100 (MQ=255)
                                                                               ggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac               <  2:1086684/100‑1 (MQ=255)
                                                                               ggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac               <  1:508251/100‑1 (MQ=255)
                                                                               ggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac               <  2:765838/100‑1 (MQ=255)
                                                                                 tCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacca             <  1:306599/100‑1 (MQ=255)
                                                                                   aaaCTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc           >  1:1134505/1‑100 (MQ=255)
                                                                                    aaCTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAccacca          >  2:1300185/1‑100 (MQ=255)
                                                                                     aCTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAg         <  2:8484/100‑1 (MQ=255)
                                                                                      cTTCTCTACCGTGCGGCTCAGCCCCAGGTTGCCCTCTTCGCTAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTATTACCACCAg         >  2:454946/1‑99 (MQ=255)
                                                                                      cTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGa        >  1:1299073/1‑100 (MQ=255)
                                                                                      cTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGa        <  2:189386/100‑1 (MQ=255)
                                                                                      cTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGa        <  2:697200/100‑1 (MQ=255)
                                                                                       ttCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc       <  2:695656/100‑1 (MQ=255)
                                                                                       ttCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc       <  2:50755/100‑1 (MQ=255)
                                                                                        tctcTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAcn      <  1:729348/100‑2 (MQ=255)
                                                                                        tctcTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACg      >  2:432846/1‑100 (MQ=255)
                                                                                         ctctACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGc     >  2:228381/1‑100 (MQ=255)
                                                                                          tctACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGc     <  1:1558745/99‑1 (MQ=255)
                                                                                            tACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCTACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg  >  2:876777/1‑100 (MQ=255)
                                                                                                 |                                                                                              
GTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAG  >  NC_000913/2867067‑2867258

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