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.93.8% 314.9 / 14.4 81intergenic (+71/‑160)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base . (40/36);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGA  >  NC_000913/107446‑107639
                                                                                                   |                                                                                              
tCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGt                                                                                                >  2:825683/1‑100 (MQ=255)
tCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGt                                                                                                <  2:1553289/100‑1 (MQ=255)
tCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGt                                                                                                <  2:943254/100‑1 (MQ=255)
 cAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCCTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtt                                                                                               <  2:269627/100‑1 (MQ=255)
  aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttt                                                                                              >  2:353152/1‑100 (MQ=255)
    aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                            >  2:1960910/1‑100 (MQ=255)
     gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                           >  2:747960/1‑100 (MQ=255)
     gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                           <  2:1177082/100‑1 (MQ=255)
      cgcCTTCAGCTCTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                            <  2:110726/98‑1 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCTCTCTCTCCGGCCAGGTGAACCAGTCGttttttttt                                                                                        <  1:2239333/100‑1 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                        <  2:951575/100‑2 (MQ=255)
         cTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttga                                                                                       >  2:2234110/1‑98 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCCGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                                 >  2:775190/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                                 >  2:2304693/1‑100 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAag                                                                            <  1:2507474/100‑1 (MQ=255)
                    cTGGCATAACGACATTTAGACTGGCGTATAAAATTGGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAag                                                                            <  1:458163/99‑1 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGATTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                         <  2:1818882/100‑1 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                         >  2:1850301/1‑100 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                         <  2:751124/100‑1 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                         <  2:415083/100‑1 (MQ=255)
                       gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCt                                                                        <  2:241870/100‑1 (MQ=255)
                        cATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCta                                                                       >  2:1998982/1‑100 (MQ=255)
                            aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGACCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGAtata                                                                     >  1:1900394/1‑98 (MQ=255)
                             cGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaaa                                                                  <  2:1335591/100‑1 (MQ=255)
                                cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTATTTTTTAAGAGCTATAAAAAAc                                                               >  1:68909/1‑100 (MQ=255)
                                cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGTCTTTTATAATAGCTATAAAAAAc                                                               >  1:1269166/1‑100 (MQ=255)
                                cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                               >  2:142256/1‑100 (MQ=255)
                                cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGGCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                               >  1:2425089/1‑100 (MQ=255)
                                 aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                               <  2:874893/99‑1 (MQ=255)
                                 aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACg                                                              <  2:76198/100‑1 (MQ=255)
                                   ttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGt                                                            <  2:2011980/100‑1 (MQ=255)
                                     ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTTTTAGAGCTATAAAAAACGGTGc                                                          >  1:1468485/1‑100 (MQ=255)
                                     aaaatgaagTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTAAAAAAAAAAAAAAAAACgg                                                             <  1:2139081/89‑1 (MQ=255)
                                      taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCg                                                         >  2:343527/1‑100 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                       >  2:1927406/1‑100 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                       >  1:2065777/1‑100 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                       >  2:1641476/1‑100 (MQ=255)
                                         tGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAc                                                      >  2:145963/1‑100 (MQ=255)
                                         tGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAAc                                                      <  2:415896/100‑1 (MQ=255)
                                          gTCGTATGAAATTCGTCTTGTAATTCTTTTACTCTCTCCGCGCACGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGTTGCGaa                                                       <  1:1258529/98‑1 (MQ=255)
                                             gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGc                                                    <  2:2065790/98‑1 (MQ=255)
                                             gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTg                                                  >  2:2322797/1‑100 (MQ=255)
                                               ataAAATTCGACTGGCAAATCTGACACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                                >  1:874863/1‑100 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAATAGCTATAAAAAACGGTGCGAACGCTGttt                                               >  1:144396/1‑100 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                               <  2:1836510/100‑1 (MQ=255)
                                                  aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                             >  2:1985184/1‑100 (MQ=255)
                                                  aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                             >  2:776873/1‑100 (MQ=255)
                                                  aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                             >  2:1652665/1‑100 (MQ=255)
                                                    aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                           <  2:1513435/100‑1 (MQ=255)
                                                       cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                        <  1:521818/100‑1 (MQ=255)
                                                       cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                        >  2:2485959/1‑100 (MQ=255)
                                                        gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                       >  2:1836192/1‑100 (MQ=255)
                                                          cTGGCATATTTTGCATTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGATCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                     <  1:820197/100‑1 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                     >  2:67086/1‑100 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGGGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                     <  1:316667/100‑1 (MQ=255)
                                                           tGGCAAATCTGGCATTCTTTCCGGTCAGGTGGATTAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACTGTGCGAACGCTGTTTTCTTAAGc                                       <  1:1367442/97‑1 (MQ=255)
                                                           tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                    <  2:944945/100‑1 (MQ=255)
                                                             gCAAATCTGGCACTCCCTCCGGCCAGGTGAACCAGTCG‑‑TTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtttt                                 >  2:1059763/1‑100 (MQ=255)
                                                              cAAAACTGTCACTCTCCCCGGCCAGGTGAACCAGTCGTTTTTTTTGATTTTTATAAAAAAAAAAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                  <  1:172393/99‑1 (MQ=255)
                                                               agaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAATAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                >  2:1752696/3‑100 (MQ=255)
                                                               aaaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                <  2:143574/100‑1 (MQ=255)
                                                                 aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                              >  2:1189931/1‑100 (MQ=255)
                                                                 aTCTGGCACTCTCTCCCGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                              <  1:1516514/100‑1 (MQ=255)
                                                                  tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                             <  1:879872/100‑1 (MQ=255)
                                                                   cTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGca                            >  2:2280931/1‑100 (MQ=255)
                                                                   cTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGca                            >  2:2280962/1‑100 (MQ=255)
                                                                    tGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                           <  1:1788189/100‑1 (MQ=255)
                                                                      gcccccgctcCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcaca                          <  2:819221/93‑1 (MQ=255)
                                                                        aCTCTCTCCGGCCAGGTGAACCAGTCGGTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                           <  1:338480/96‑1 (MQ=255)
                                                                          tctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGATTTTTATAAGAACTATCAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                        >  1:572351/1‑97 (MQ=255)
                                                                          tctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                     >  2:2590339/1‑100 (MQ=255)
                                                                           ctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                    >  2:2390941/1‑100 (MQ=255)
                                                                            tctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                   <  1:1856911/100‑1 (MQ=255)
                                                                            tccccGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                   <  2:2345593/97‑1 (MQ=255)
                                                                             ctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                  >  1:2376272/1‑100 (MQ=255)
                                                                              tcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCtt                 <  2:2210622/100‑1 (MQ=255)
                                                                              tcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCtt                 >  1:50807/1‑100 (MQ=255)
                                                                              tcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAACCGCTTTTCCGCACAACTTATCtt                 >  1:785533/1‑100 (MQ=255)
                                                                               ccGGCCAGGTGAACCTGTCGTTTTTTTTTTTTTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                >  1:659911/1‑100 (MQ=255)
                                                                               ccGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                >  2:162167/1‑100 (MQ=255)
                                                                               ccGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                >  1:780439/1‑100 (MQ=255)
                                                                                cGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                >  2:2297216/1‑99 (MQ=255)
                                                                                    cAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg           >  2:2040253/1‑100 (MQ=255)
                                                                                    cAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg           <  1:1483931/100‑1 (MQ=255)
                                                                                    cAGGTGAACCAGTAGTTTTTTTTGAATTTTAAAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg           >  2:1744070/1‑100 (MQ=255)
                                                                                      ggTGAACCAGTCGTTTTTTTTTATTTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTGCGCACAACTTATCTTCATTCGTg         >  1:1162612/1‑100 (MQ=255)
                                                                                      ggTGAACCAGTCG‑‑TTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc        <  1:1059759/100‑1 (MQ=255)
                                                                                       gTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc        <  1:470546/100‑1 (MQ=255)
                                                                                         gAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTGCGCACAACTTATCTTCATTCGTGCtg      <  2:1578927/100‑1 (MQ=255)
                                                                                          aaCCAATCGTTTTTTTTGAATTTTATAAGATCTATAAAAAACCGTCCGAACTCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg         <  1:2154360/96‑1 (MQ=255)
                                                                                            ccAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg   <  2:2301343/100‑1 (MQ=255)
                                                                                             cAGTCGTTTTTTTTGAATTTTATGAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGa  >  2:1600538/1‑100 (MQ=255)
                                                                                             cAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAGACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGATCACCTTATCTTCATTCGTGTtgtg    >  1:975774/1‑98 (MQ=255)
                                                                                                   |                                                                                              
TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGA  >  NC_000913/107446‑107639

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