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,1640CT94.8% 286.4 / 3.9 96A130T (GCG→ACG) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (47/44);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAC  >  NC_000913/2867067‑2867253
                                                                                                 |                                                                                         
gTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgc                                                                                         >  2:836058/1‑100 (MQ=255)
 tttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgcg                                                                                        <  2:2013397/100‑1 (MQ=255)
 tttGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACcgtgcg                                                                                        <  1:2508030/100‑1 (MQ=255)
    ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGAt                                                                                     >  1:1118316/1‑100 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                    >  2:1076205/1‑100 (MQ=255)
     gTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATc                                                                                    >  2:66775/1‑100 (MQ=255)
       tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAg                                                                                  <  1:272847/100‑1 (MQ=255)
       tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAg                                                                                  <  1:430047/100‑1 (MQ=255)
       tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAg                                                                                  >  1:1321335/1‑100 (MQ=255)
       tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAg                                                                                  >  1:485278/1‑100 (MQ=255)
       tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                                 <  2:575425/100‑1 (MQ=255)
        cATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGc                                                                                 <  1:2377563/100‑1 (MQ=255)
          tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGccc                                                                               >  2:306205/1‑100 (MQ=255)
          tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGccc                                                                               <  1:2126303/100‑1 (MQ=255)
           aaTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGcccc                                                                              <  2:2189816/100‑1 (MQ=255)
               gaccGTTCAATCGTCTGGCGAATCCAACAGATTGCGTATGTTGAGACGCGGAAACCACGTTCCGAGTCAAACTTCTCTACCGTGCGGATCAGCCCCAg                                                                            >  2:1490636/3‑98 (MQ=255)
                cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGtt                                                                         >  1:670846/1‑100 (MQ=255)
                  atctCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGc                                                                       <  1:848716/97‑1 (MQ=255)
                    ttCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGccc                                                                     <  2:2239923/100‑1 (MQ=255)
                    ttCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGccc                                                                     <  1:1258116/100‑1 (MQ=255)
                      cAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGcagat                                                                   >  2:848709/1‑96 (MQ=255)
                      cAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCctc                                                                   >  2:1069570/1‑100 (MQ=255)
                       aaTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCctct                                                                  >  1:1191413/1‑100 (MQ=255)
                        aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCtt                                                                 <  2:1699111/100‑1 (MQ=255)
                         tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCctctcc                                                                <  1:1734269/100‑3 (MQ=255)
                          cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                               <  1:2492070/100‑1 (MQ=255)
                          cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                               >  2:1667980/1‑100 (MQ=255)
                          cGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCg                                                               <  1:1383017/100‑1 (MQ=255)
                           gTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAATTTCTCTACCGTGCGGATCAGCCCCAGGTTGCACTCTTCGa                                                              >  2:1232757/1‑100 (MQ=255)
                             cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGAGGAAACCACGTTCCGGGTCAAACTGCACTCACGTTCGGAACAGCCCCAGTTTGCCCTCTTCGAt                                                             >  2:52549/1‑99 (MQ=255)
                              tGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATaa                                                           <  2:1947900/100‑1 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                          >  1:563792/1‑100 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                          <  1:2560324/100‑1 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                          >  1:174045/1‑100 (MQ=255)
                               ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAg                                                          >  2:1649053/1‑100 (MQ=255)
                                gCGACTCCACCAGGTTGCGTCTGTTGAGAAGCGGAGACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                         <  2:1592585/100‑1 (MQ=255)
                                gCGAATCCATCAGTTTGCGTATGTTGAGAAGCGGGAGCCGCGTTCCGGGTCAAACTTCTCGGCGGTGCGGATCAGCCGCAGGTTGCCCTCTTCGATAAgg                                                         <  2:1641614/100‑1 (MQ=255)
                                gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                         <  1:836065/100‑1 (MQ=255)
                                gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                         <  1:1198038/100‑1 (MQ=255)
                                gCGAAGCAACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                         <  2:1288584/100‑1 (MQ=255)
                                 cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                        >  1:43455/1‑100 (MQ=255)
                                   aaTCCACCAGGTTGCATATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTcc                                                      <  2:1631119/100‑1 (MQ=255)
                                    aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCa                                                     <  1:818862/100‑1 (MQ=255)
                                     tCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAg                                                    >  2:1442464/1‑100 (MQ=255)
                                     tCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAg                                                    <  2:2282718/100‑1 (MQ=255)
                                      ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                   >  1:597930/1‑100 (MQ=255)
                                      ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                   >  1:781168/1‑100 (MQ=255)
                                       caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                  >  2:724866/1‑100 (MQ=255)
                                       caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                  >  2:998883/1‑100 (MQ=255)
                                        accaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCaa                                                 >  1:1393254/1‑100 (MQ=255)
                                         ccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAAc                                                >  2:2541191/1‑100 (MQ=255)
                                             gTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCa                                            <  2:70145/100‑1 (MQ=255)
                                              ttGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAg                                           <  2:1789173/100‑1 (MQ=255)
                                              ttGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAg                                           >  2:1139558/1‑100 (MQ=255)
                                              ttGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAg                                           >  2:1438244/1‑100 (MQ=255)
                                                 cGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGAcc                                        <  1:46759/100‑1 (MQ=255)
                                                  gTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCa                                       >  1:559961/1‑100 (MQ=255)
                                                   tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCAc                                      >  2:570614/1‑100 (MQ=255)
                                                    aTGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACg                                     >  1:2251796/1‑100 (MQ=255)
                                                     tGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                    >  2:727257/1‑100 (MQ=255)
                                                     tGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                    >  1:1652706/1‑100 (MQ=255)
                                                       ttGAGAAGCGGAAACCACGTTCCGGGCCAAACTTGTCTACCGTGCGGAGCAGTCCCAGGGTGCCCTCTTCGATAAGGTCGCGCAACGCAAGAGCACGAtt                                  >  2:743265/1‑100 (MQ=255)
                                                        tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                                 <  2:2067901/100‑1 (MQ=255)
                                                        tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTg                                 >  2:2418227/1‑100 (MQ=255)
                                                         gagaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                >  1:1653879/1‑100 (MQ=255)
                                                          tctaGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGcc                               <  1:1759578/97‑1 (MQ=255)
                                                           naaGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                              >  2:2229944/2‑100 (MQ=255)
                                                           gaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                              <  2:2334852/100‑1 (MQ=255)
                                                           gaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                              <  1:1649065/100‑1 (MQ=255)
                                                            aaGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCAt                             <  2:1720556/100‑1 (MQ=255)
                                                             aGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCaga                            >  2:1759565/1‑98 (MQ=255)
                                                                ggAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                         >  1:1686804/1‑100 (MQ=255)
                                                                 gAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGcccgccgtt                        >  1:1869887/1‑93 (MQ=255)
                                                                 gAAACCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcgg                        >  1:973613/1‑100 (MQ=255)
                                                                  aaaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                       <  1:747457/100‑1 (MQ=255)
                                                                   aaCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                      <  1:1657441/100‑1 (MQ=255)
                                                                    aCCACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcgg                     <  1:451815/100‑1 (MQ=255)
                                                                     ccACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCggg                    <  1:1708426/100‑1 (MQ=255)
                                                                      cACTTTCCGTGTCAAACGTCTCTACCGTGCGGATCAGGCCCGAGTTGCCCTCTTCGAGAAGGTCCAGCAACGCCAGACCACGATTTCCATAACGGCggg                    <  2:2548111/99‑1 (MQ=255)
                                                                      cACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                   <  2:284117/100‑1 (MQ=255)
                                                                      cACGTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                   >  1:727308/1‑100 (MQ=255)
                                                                         gTTCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAt                >  1:234808/1‑100 (MQ=255)
                                                                          ttCCGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAtn               <  1:2229958/100‑2 (MQ=255)
                                                                            ccGGGTCAAACTTCTCTACCGTGCGGATCAGCCGCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGTGGGCAAtttt             >  2:1684715/1‑100 (MQ=255)
                                                                             cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt            <  1:1069578/100‑1 (MQ=255)
                                                                             cGGGTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt            >  1:2055071/1‑100 (MQ=255)
                                                                             cGGGTCAAACTTCTCTACCGTGCGGAACAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt            >  2:788474/1‑100 (MQ=255)
                                                                              gggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa           >  2:204612/1‑100 (MQ=255)
                                                                              gggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa           >  2:1259808/1‑100 (MQ=255)
                                                                              gggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa           <  2:1088945/100‑1 (MQ=255)
                                                                               ggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac          <  2:519849/100‑1 (MQ=255)
                                                                               ggTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac          <  1:2294773/100‑1 (MQ=255)
                                                                                gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCGAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc         <  2:2050004/100‑1 (MQ=255)
                                                                                gTCAAACTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc         <  1:2248055/100‑1 (MQ=255)
                                                                                   aaaCTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc      >  1:1726610/1‑100 (MQ=255)
                                                                                    aaCTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAccacca     >  1:1574067/1‑100 (MQ=255)
                                                                                     aCTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAg    <  1:2418241/100‑1 (MQ=255)
                                                                                     aCTTCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAg    <  1:502151/100‑1 (MQ=255)
                                                                                       ttCTCTACCGTGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc  <  1:868523/100‑1 (MQ=255)
                                                                                                 |                                                                                         
GTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAC  >  NC_000913/2867067‑2867253

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