Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 107,553 (T)9→8 100% intergenic (+79/‑152) lpxC → / → secM UDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913107,5450T.92.6% 315.4 / 19.0 81intergenic (+71/‑160)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base T (5/1);  new base . (37/38);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACT  >  NC_000913/107446‑107641
                                                                                                   |                                                                                                
tCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGt                                                                                                  >  2:2179918/1‑100 (MQ=255)
 cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtt                                                                                                 >  2:437370/1‑100 (MQ=255)
  aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttt                                                                                                >  1:536820/1‑100 (MQ=255)
    aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                              >  2:2216066/1‑100 (MQ=255)
      cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                            <  1:445469/100‑1 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                          >  1:689461/1‑99 (MQ=255)
          ttCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttgaa                                                                                        >  2:2195096/1‑97 (MQ=255)
           tCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttgaat                                                                                       >  1:1379430/1‑96 (MQ=255)
            cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGCTGAACAAGTCGTTTTTTTTGAAtt                                                                                     >  2:82186/1‑100 (MQ=255)
             aGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAttt                                                                                    >  1:1920246/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGCCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCTGTCGttttttttgcattt                                                                                     >  2:1450588/1‑93 (MQ=255)
               cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTta                                                                                  <  1:943572/100‑1 (MQ=255)
                tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtat                                                                                 <  1:925601/100‑1 (MQ=255)
                tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtat                                                                                 >  1:1972480/1‑100 (MQ=255)
                 gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTtat                                                                                 <  1:999817/100‑1 (MQ=255)
                 gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtata                                                                                <  1:1090161/100‑1 (MQ=255)
                  tACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATaa                                                                               >  2:1142501/1‑100 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAag                                                                              >  1:1449055/1‑100 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAag                                                                              >  1:215929/1‑100 (MQ=255)
                     tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAagag                                                                            <  2:1403044/100‑1 (MQ=255)
                       gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCt                                                                          <  1:245266/100‑1 (MQ=255)
                          tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtata                                                                       <  1:178998/100‑1 (MQ=255)
                          tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCAGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAATAGCtat                                                                        >  2:1093754/1‑99 (MQ=255)
                            aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaa                                                                     >  2:1464538/1‑100 (MQ=255)
                            aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaa                                                                     >  1:564758/1‑100 (MQ=255)
                            aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaa                                                                     >  2:636103/1‑100 (MQ=255)
                                cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                 >  1:1707337/1‑100 (MQ=255)
                                 aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACg                                                                >  1:1352/1‑100 (MQ=255)
                                    tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                             <  2:1930772/100‑1 (MQ=255)
                                     ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                            <  1:1398170/100‑1 (MQ=255)
                                     ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                            <  1:1424469/100‑1 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                         >  1:1873957/1‑100 (MQ=255)
                                         cgatctataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAc                                                        <  1:968560/95‑1 (MQ=255)
                                          gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                       <  2:2282362/100‑1 (MQ=255)
                                          gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                       <  1:1285486/100‑1 (MQ=255)
                                            cGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCt                                                     >  1:574943/1‑100 (MQ=255)
                                              tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGt                                                   >  2:1869171/1‑100 (MQ=255)
                                               ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAcagatcg                                                  >  2:968552/1‑94 (MQ=255)
                                               ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                  >  1:2227552/1‑100 (MQ=255)
                                               ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                  <  1:1279513/100‑1 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                 <  1:1401724/100‑1 (MQ=255)
                                                ctcttccgatcTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                 <  1:2392852/90‑1 (MQ=255)
                                                    agtCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGGTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                             <  2:279747/98‑1 (MQ=255)
                                                    aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                             <  2:2505811/100‑1 (MQ=255)
                                                    aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                             <  1:1329606/100‑1 (MQ=255)
                                                    aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAAACAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                             <  2:2505810/100‑1 (MQ=255)
                                                     ttCGACTGGCAAATCTGGCACTCTCTCCAGCCAGGTGAACCAGGCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                            <  2:329474/100‑1 (MQ=255)
                                                      tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATGAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                           <  2:1371766/100‑1 (MQ=255)
                                                        gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                         <  1:876320/100‑1 (MQ=255)
                                                        gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                         <  1:712933/100‑1 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAACACGGTGCGAAAGCTGTTGTCTTAAGCAc                                       >  1:2091580/1‑100 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                       <  1:2485477/100‑1 (MQ=255)
                                                           tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                      <  1:320212/100‑1 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                     >  1:2026154/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                     >  1:1661813/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                     >  1:1035127/1‑100 (MQ=255)
                                                             gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                    >  1:723673/1‑100 (MQ=255)
                                                             gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                    >  1:837864/1‑100 (MQ=255)
                                                              cAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtttt                                   >  2:2286521/1‑100 (MQ=255)
                                                              cAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtttt                                   >  2:996514/1‑100 (MQ=255)
                                                               aaaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                  >  1:1182469/1‑100 (MQ=255)
                                                                aaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                 <  1:2226104/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtttt                                   >  1:322129/1‑96 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                               <  2:1076425/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                               <  2:442148/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                               <  1:437374/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGACAGGTTAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                                <  1:2081051/99‑1 (MQ=255)
                                                                   cTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGca                              >  1:1005767/1‑100 (MQ=255)
                                                                       cACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                          >  1:592456/1‑100 (MQ=255)
                                                                       cACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                          <  1:2455754/100‑1 (MQ=255)
                                                                       cACTCTCTCCGGCCAGGTGAACCAGTCGTTGTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                          <  2:2032251/100‑1 (MQ=255)
                                                                         ctctctCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACtt                        >  1:120751/1‑100 (MQ=255)
                                                                         ctctctCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACtt                        >  1:24896/1‑100 (MQ=255)
                                                                            tctcCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                     >  2:958367/1‑100 (MQ=255)
                                                                             ctcCGGCCAGGTGAACCAGTCGTTTTTTTTGTATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                    >  2:998506/1‑100 (MQ=255)
                                                                             ctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                    >  1:1477273/1‑100 (MQ=255)
                                                                             ctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAACAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTGTCt                    >  2:344545/1‑100 (MQ=255)
                                                                               ccGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATTAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTATTCCGCTCAAATTCTCTTc                  >  2:2227536/1‑100 (MQ=255)
                                                                               ccGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                  >  1:973019/1‑100 (MQ=255)
                                                                                 ggCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTTTAAGAACTATAAAAAACAGTGTAAACACTTTTTTCTTAAGCTCTTTTCCGCACAACTTctct                    >  2:2318351/1‑96 (MQ=255)
                                                                                    cAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg             <  2:215927/100‑1 (MQ=255)
                                                                                      ggTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg           >  1:1401605/1‑100 (MQ=255)
                                                                                       gTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc          >  2:653424/1‑100 (MQ=255)
                                                                                       gTGAACCAGTCGTGTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc          <  2:2021632/100‑1 (MQ=255)
                                                                                        tGAACCAGTCGTTTTTTTTGATTTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGTTCTTTTCCGTACAACTTATCTTCATTCTTgtt         >  2:676554/1‑98 (MQ=255)
                                                                                           aCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg      <  1:31553/100‑1 (MQ=255)
                                                                                           aCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg      <  2:1666945/100‑1 (MQ=255)
                                                                                           aCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg      <  2:2178101/100‑1 (MQ=255)
                                                                                            ccAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg     <  1:653431/100‑1 (MQ=255)
                                                                                               gTCGTTTTTTTTTTTTTTTTTAAGAGCTATAAAAAACGGTTCGATCGCTGTTTTCTTAAGCACTTTTCTGCACCACTTATCTTCATTCGTGCTGTGTACt  >  2:1244878/1‑100 (MQ=255)
                                                                                                   |                                                                                                
TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGACT  >  NC_000913/107446‑107641

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: