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)
            cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGCTGAACAAGTCG‑TTTTTTTTGAAtt                                                                                     >  2:82186/1‑100 (MQ=255)
             aGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAttt                                                                                    >  1:1920246/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGCCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCTGTCGttttttttgcattt                                                                                     >  2:1450588/1‑93 (MQ=255)
               cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTta                                                                                  <  1:943572/100‑1 (MQ=255)
                tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtat                                                                                 <  1:925601/100‑1 (MQ=255)
                tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtat                                                                                 >  1:1972480/1‑100 (MQ=255)
                 gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTtat                                                                                 <  1:999817/100‑1 (MQ=255)
                 gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtata                                                                                <  1:1090161/100‑1 (MQ=255)
                  tACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATaa                                                                               >  2:1142501/1‑100 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAag                                                                              >  1:1449055/1‑100 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAag                                                                              >  1:215929/1‑100 (MQ=255)
                     tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAagag                                                                            <  2:1403044/100‑1 (MQ=255)
                       gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCt                                                                          <  1:245266/100‑1 (MQ=255)
                          tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCtata                                                                       <  1:178998/100‑1 (MQ=255)
                          tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCAGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAATAGCtat                                                                        >  2:1093754/1‑99 (MQ=255)
                            aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaa                                                                     >  2:1464538/1‑100 (MQ=255)
                            aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaa                                                                     >  1:564758/1‑100 (MQ=255)
                            aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaa                                                                     >  2:636103/1‑100 (MQ=255)
                                cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                 >  1:1707337/1‑100 (MQ=255)
                                 aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACg                                                                >  1:1352/1‑100 (MQ=255)
                                    tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                             <  2:1930772/100‑1 (MQ=255)
                                     ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                            <  1:1398170/100‑1 (MQ=255)
                                     ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                            <  1:1424469/100‑1 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                         >  1:1873957/1‑100 (MQ=255)
                                         cgatctataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAc                                                        <  1:968560/95‑1 (MQ=255)
                                          gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                       <  2:2282362/100‑1 (MQ=255)
                                          gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                       <  1:1285486/100‑1 (MQ=255)
                                            cGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCt                                                     >  1:574943/1‑100 (MQ=255)
                                              tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGt                                                   >  2:1869171/1‑100 (MQ=255)
                                               ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAcagatcg                                                  >  2:968552/1‑94 (MQ=255)
                                               ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                  >  1:2227552/1‑100 (MQ=255)
                                               ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                  <  1:1279513/100‑1 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                 <  1:1401724/100‑1 (MQ=255)
                                                ctcttccgatcTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                 <  1:2392852/90‑1 (MQ=255)
                                                    agtCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑GTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                             <  2:279747/98‑1 (MQ=255)
                                                    aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                             <  2:2505811/100‑1 (MQ=255)
                                                    aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                             <  1:1329606/100‑1 (MQ=255)
                                                    aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAAACAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                             <  2:2505810/100‑1 (MQ=255)
                                                     ttCGACTGGCAAATCTGGCACTCTCTCCAGCCAGGTGAACCAGGCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                            <  2:329474/100‑1 (MQ=255)
                                                      tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATGAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                           <  2:1371766/100‑1 (MQ=255)
                                                        gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                         <  1:876320/100‑1 (MQ=255)
                                                        gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                         <  1:712933/100‑1 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAACACGGTGCGAAAGCTGTTGTCTTAAGCAc                                       >  1:2091580/1‑100 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                       <  1:2485477/100‑1 (MQ=255)
                                                           tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                      <  1:320212/100‑1 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                     >  1:2026154/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                     >  1:1661813/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                     >  1:1035127/1‑100 (MQ=255)
                                                             gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                    >  1:723673/1‑100 (MQ=255)
                                                             gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                    >  1:837864/1‑100 (MQ=255)
                                                              cAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtttt                                   >  2:2286521/1‑100 (MQ=255)
                                                              cAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑AATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtttt                                   >  2:996514/1‑100 (MQ=255)
                                                               aaaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                  >  1:1182469/1‑100 (MQ=255)
                                                                aaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                 <  1:2226104/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtttt                                   >  1:322129/1‑96 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                               <  2:1076425/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                               <  2:442148/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                               <  1:437374/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGACAGGTTAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                                <  1:2081051/99‑1 (MQ=255)
                                                                   cTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGca                              >  1:1005767/1‑100 (MQ=255)
                                                                       cACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                          >  1:592456/1‑100 (MQ=255)
                                                                       cACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                          <  1:2455754/100‑1 (MQ=255)
                                                                       cACTCTCTCCGGCCAGGTGAACCAGTCG‑TTGTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                          <  2:2032251/100‑1 (MQ=255)
                                                                         ctctctCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACtt                        >  1:120751/1‑100 (MQ=255)
                                                                         ctctctCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACtt                        >  1:24896/1‑100 (MQ=255)
                                                                            tctcCGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                     >  2:958367/1‑100 (MQ=255)
                                                                             ctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGTATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                    >  2:998506/1‑100 (MQ=255)
                                                                             ctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                    >  1:1477273/1‑100 (MQ=255)
                                                                             ctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAACAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTGTCt                    >  2:344545/1‑100 (MQ=255)
                                                                               ccGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATTAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTATTCCGCTCAAATTCTCTTc                  >  2:2227536/1‑100 (MQ=255)
                                                                               ccGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                  >  1:973019/1‑100 (MQ=255)
                                                                                 ggCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAACTATAAAAAACAGTGTAAACACTTTTTTCTTAAGCTCTTTTCCGCACAACTTctct                    >  2:2318351/1‑96 (MQ=255)
                                                                                    cAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg             <  2:215927/100‑1 (MQ=255)
                                                                                      ggTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg           >  1:1401605/1‑100 (MQ=255)
                                                                                       gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc          >  2:653424/1‑100 (MQ=255)
                                                                                       gTGAACCAGTCG‑TGTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc          <  2:2021632/100‑1 (MQ=255)
                                                                                        tGAACCAGTCG‑TTTTTTTTGATTTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGTTCTTTTCCGTACAACTTATCTTCATTCTTgtt         >  2:676554/1‑98 (MQ=255)
                                                                                           aCCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg      <  1:31553/100‑1 (MQ=255)
                                                                                           aCCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg      <  2:1666945/100‑1 (MQ=255)
                                                                                           aCCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg      <  2:2178101/100‑1 (MQ=255)
                                                                                            ccAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg     <  1:653431/100‑1 (MQ=255)
                                                                                               gTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTTCGATCGCTGTTTTCTTAAGCACTTTTCTGCACCACTTATCTTCATTCGTGCTGTGTACt  >  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: ‑