Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,681,932 T→C P33P (CCA→CCG)  ydgC ← GlpM family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,681,9320TC99.3% 500.9 / ‑5.0 144P33P (CCA→CCG) ydgCGlpM family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (64/78);  total (66/78)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 144 total aligned reads displayed.

TCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCC  >  NC_000913/1681835‑1682027
                                                                                                 |                                                                                               
tCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCgg                                                                                               <  2:2195045/100‑1 (MQ=255)
  gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                             >  1:2634260/1‑100 (MQ=255)
  gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                             <  2:65977/100‑1 (MQ=255)
  gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                             >  2:599610/1‑100 (MQ=255)
  gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                             <  2:3136561/100‑1 (MQ=255)
    ccACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATc                                                                                           <  2:1263148/100‑1 (MQ=255)
    ccACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAACGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATc                                                                                           >  2:1819006/1‑100 (MQ=255)
      acaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAg                                                                                         <  2:483899/100‑1 (MQ=255)
       caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAg                                                                                         >  1:2547239/1‑99 (MQ=255)
       caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                        >  2:290259/1‑100 (MQ=255)
        aTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGcc                                                                                       <  2:3396641/100‑1 (MQ=255)
         tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                      >  2:1957326/1‑100 (MQ=255)
         tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGCCGGGAAAAGCGGAATCAGctc                                                                                      >  1:2108996/1‑98 (MQ=255)
          aCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCg                                                                                     >  1:1429761/1‑100 (MQ=255)
          aCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCg                                                                                     >  1:2638503/1‑100 (MQ=255)
           cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                    >  1:49236/1‑100 (MQ=255)
            tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                   >  2:2268759/1‑100 (MQ=255)
            tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                   <  1:2746524/100‑1 (MQ=255)
             aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg                                                                                  >  2:112093/1‑100 (MQ=255)
             aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg                                                                                  >  1:2989626/1‑100 (MQ=255)
              aaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGa                                                                                 <  2:2762929/100‑1 (MQ=255)
               aaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGat                                                                                <  1:2814117/100‑1 (MQ=255)
                aaTGAAGGATGCGCGTAAGGCTTCAATGCCGCGATCGCTGGCAACAATATAATGCGCGATAAGAAAAAAAAAAAAAAAAAGCGGAATCAGCCCGGCGata                                                                               <  1:1163769/100‑1 (MQ=255)
                       gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                        >  1:557684/1‑100 (MQ=255)
                       gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                        <  1:1191017/100‑1 (MQ=255)
                           ngcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                    >  2:255627/2‑100 (MQ=255)
                           cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                    >  2:2223201/1‑100 (MQ=255)
                            gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                   <  2:1994283/100‑1 (MQ=255)
                            gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                   <  1:2147828/100‑1 (MQ=255)
                             cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTcn                                                                  <  2:933825/100‑2 (MQ=255)
                             cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                  <  2:678894/100‑1 (MQ=255)
                              gTAAGGCTTCAATTCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGt                                                                 >  2:2354485/1‑100 (MQ=255)
                               tAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtt                                                                <  2:1459824/100‑1 (MQ=255)
                                aaGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttt                                                               >  1:1161825/1‑100 (MQ=255)
                                 aGGCCTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtttt                                                              >  1:30672/1‑100 (MQ=255)
                                  ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                             >  2:1077806/1‑100 (MQ=255)
                                  ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                             >  1:2159731/1‑100 (MQ=255)
                                  ggCTTCAATGCCCCGTTCGCGGGCAACAATATAATGCGCGATAGGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                             <  1:1025087/100‑1 (MQ=255)
                                   gcgtCAATGCCGGGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGGCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                            <  1:1883768/97‑1 (MQ=255)
                                   gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                            <  2:3306910/100‑1 (MQ=255)
                                   gCTTCAATGCCGCGTTAGCAGGCAACAATAAAATCAGCAATAAGCGCAAAGGACAGGAAAAGCGAAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                            <  1:1267485/100‑1 (MQ=255)
                                    cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGTTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           <  1:330501/100‑1 (MQ=255)
                                    cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttttc                                                           <  1:330305/100‑3 (MQ=255)
                                    cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           <  2:722977/100‑1 (MQ=255)
                                    cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           <  1:1245900/100‑1 (MQ=255)
                                    cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           <  2:1803038/100‑1 (MQ=255)
                                     ttCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCt                                                          <  2:2875203/100‑1 (MQ=255)
                                       cAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGTGATATAATAATTTTTCGTTTTTGCTaa                                                        >  2:2738793/1‑100 (MQ=255)
                                        aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                       >  2:1369807/1‑100 (MQ=255)
                                        aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                       >  2:483568/1‑100 (MQ=255)
                                        aaTGCCGCGTTCACTGGCAACAATATAATGGGCGATAAGCGCAAAGGTCGTGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                       <  1:3131372/100‑1 (MQ=255)
                                          tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaaan                                                     <  1:933965/100‑2 (MQ=255)
                                           gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAaca                                                    >  2:1412677/1‑100 (MQ=255)
                                            ccGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAacac                                                   >  1:195689/1‑100 (MQ=255)
                                            ccGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAacac                                                   <  2:2109027/100‑1 (MQ=255)
                                             cccgTTCGCTGGCAACAATATAATGCGCGATAAGCGAAAAGTTCGGTAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                  <  1:1633645/98‑1 (MQ=255)
                                                 ttCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                              <  2:2839189/100‑1 (MQ=255)
                                                 ttCGCTGGCAACAATATAATGCGCCATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                              <  1:1204299/100‑1 (MQ=255)
                                                  tCGCTGGCAACAATATAATGCGCGATAAGCGGAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATACTTTTTCGTTTTTGCTAAAACACCAATc                                              <  1:1454727/99‑1 (MQ=255)
                                                  tCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCa                                             <  1:1292990/100‑1 (MQ=255)
                                                    gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                           >  1:2089675/1‑100 (MQ=255)
                                                    gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                           <  2:3027201/100‑1 (MQ=255)
                                                    gCTGGCAACAATAGAATGCGCGATAAGCGCAAAGGTCGGGAAACGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                           <  1:385680/100‑1 (MQ=255)
                                                    gCTGGCAAAAATATAATGCGAGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATGTTTCGTTTTTGCTAAAACACCAATCAAc                                           <  1:1550770/100‑1 (MQ=255)
                                                     cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                          >  1:2579419/1‑100 (MQ=255)
                                                       ggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc                                        >  1:2937757/1‑100 (MQ=255)
                                                       ggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAAGCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc                                        <  1:1803191/100‑1 (MQ=255)
                                                        gCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCa                                       <  2:3161287/100‑1 (MQ=255)
                                                         caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                      >  2:820909/1‑100 (MQ=255)
                                                             aaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAcc                                  <  2:653010/100‑1 (MQ=255)
                                                                ataATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgc                               <  1:2618774/100‑1 (MQ=255)
                                                                ataATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgc                               <  2:1161840/100‑1 (MQ=255)
                                                                 taATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgcg                              >  2:831757/1‑100 (MQ=255)
                                                                  aaTGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgcgc                             >  1:2004597/1‑100 (MQ=255)
                                                                   aTGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCa                            >  2:2492981/1‑100 (MQ=255)
                                                                   aTGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCa                            >  2:900101/1‑100 (MQ=255)
                                                                    tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCa                            <  2:2817229/99‑1 (MQ=255)
                                                                    tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                           <  1:1772867/100‑1 (MQ=255)
                                                                    tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                           <  2:2792860/100‑1 (MQ=255)
                                                                    tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                           >  2:1390020/1‑100 (MQ=255)
                                                                     gcgcGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAcc                          <  2:1180524/100‑1 (MQ=255)
                                                                        cgATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAg                       >  1:2606406/1‑100 (MQ=255)
                                                                         gATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAgc                      <  2:2489288/100‑1 (MQ=255)
                                                                         gATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAgc                      >  2:302051/1‑100 (MQ=255)
                                                                            aaGCTCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACCGCACTACCAGCGCACCGAGcgcg                   <  2:148093/100‑1 (MQ=255)
                                                                            aaGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCAATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGcgcg                   <  2:511445/100‑1 (MQ=255)
                                                                             aGCGCAACGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAACCACCAATCCACAGCACTACCAGCGCACCGCGCGCgg                  <  1:1378008/100‑1 (MQ=255)
                                                                              gcgaAAAGGTCGGGAAAAGAGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGc                 >  1:3250150/1‑100 (MQ=255)
                                                                               cgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCt                <  1:78075/100‑1 (MQ=255)
                                                                                  aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg             >  2:3109122/1‑100 (MQ=255)
                                                                                    aGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt           <  1:2268727/100‑1 (MQ=255)
                                                                                     gttctggAAAAGCAGAATCACCCCCGCGATATAATAATTTTTCGTCTTTGCTAAAACACGAATCAACAGCACTACCAGCGCACCGAGCGCGGCTATGAtt          <  1:72125/95‑1 (MQ=255)
                                                                                     ggTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAtt          <  2:2634298/100‑1 (MQ=255)
                                                                                         gggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATCACCa      >  2:1912936/1‑100 (MQ=255)
                                                                                          ggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAg     >  2:504314/1‑100 (MQ=255)
                                                                                          ggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAg     >  1:964441/1‑100 (MQ=255)
                                                                                          ggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAg     <  2:1773411/100‑1 (MQ=255)
                                                                                           gAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAg     <  1:1197408/99‑1 (MQ=255)
                                                                                           gAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGc    >  2:1257674/1‑100 (MQ=255)
                                                                                             aaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGccc  <  2:2370185/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCC  >  NC_000913/1681835‑1682027

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