Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 107,553 (T)9→8 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.95.3% 441.9 / 14.5 106intergenic (+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 . (55/46);  total (59/47)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 129 total aligned reads displayed.

TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107446‑107640
                                                                                                   |                                                                                               
tCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGt                                                                                                 >  2:3267011/1‑100 (MQ=255)
    aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                             <  1:1751762/100‑1 (MQ=255)
    aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                             >  2:758689/1‑100 (MQ=255)
     gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                            >  1:2521158/1‑100 (MQ=255)
     gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                            <  1:362688/100‑1 (MQ=255)
     gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                            <  1:1126941/100‑1 (MQ=255)
      cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  1:304480/100‑1 (MQ=255)
      cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  1:2431574/100‑1 (MQ=255)
      cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  1:1710470/100‑1 (MQ=255)
      cgcCTTCAGCTGTACTGGCATAACCACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  2:2850381/100‑1 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                         <  2:2076370/100‑2 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                         <  1:2679203/100‑2 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                         <  2:1013885/100‑2 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                         <  2:1013973/100‑2 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                         >  2:3048016/1‑99 (MQ=255)
         cTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttga                                                                                        >  1:1347040/1‑98 (MQ=255)
          ttCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttgaa                                                                                       >  2:1136444/1‑97 (MQ=255)
          ttCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttattttgaa                                                                                       >  2:1806780/1‑97 (MQ=255)
           tCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttgaat                                                                                      <  1:831453/100‑5 (MQ=255)
            cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAtt                                                                                    >  1:3116815/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAtttt                                                                                  >  1:2442738/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAtttt                                                                                  >  1:1419849/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAtttt                                                                                  >  2:3314044/1‑100 (MQ=255)
                 gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtata                                                                               >  2:13781/1‑100 (MQ=255)
                 gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtata                                                                               <  2:1579512/100‑1 (MQ=255)
                 gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtata                                                                               <  1:2241871/100‑1 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtatacg                                                                             >  1:1737830/1‑98 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAag                                                                             <  1:121437/100‑1 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAag                                                                             <  2:730987/100‑1 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAag                                                                             <  1:1702404/100‑1 (MQ=255)
                    cTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAaga                                                                            >  1:173050/1‑100 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                          <  2:2251155/100‑1 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                          >  1:770878/1‑100 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                          >  1:965192/1‑100 (MQ=255)
                       gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCt                                                                         >  2:3143013/1‑100 (MQ=255)
                        cATAACGACATTTATACTGTCGTATAAAATTCGACTGGCATATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCta                                                                        <  2:1661735/100‑1 (MQ=255)
                          tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCtata                                                                      <  2:1492754/100‑1 (MQ=255)
                           aaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaa                                                                     <  2:106625/100‑1 (MQ=255)
                           aaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaa                                                                     <  1:1572306/100‑1 (MQ=255)
                                cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATGTTATAAGAGCTATAAAAAAc                                                                <  2:2048295/100‑1 (MQ=255)
                                  tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACgg                                                              <  1:1445945/100‑1 (MQ=255)
                                    tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            >  2:1001118/1‑100 (MQ=255)
                                    tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            <  2:890640/100‑1 (MQ=255)
                                     ataCTGTCGGATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            <  1:2294422/99‑1 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:1837134/1‑100 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:3069840/1‑100 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                >  2:49392/1‑100 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                <  1:2334551/100‑1 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                >  1:418070/1‑100 (MQ=255)
                                                  aaaTTCGACTGGCAAATCTGGCCCTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                              <  1:1001120/100‑1 (MQ=255)
                                                  aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                              >  2:261367/1‑100 (MQ=255)
                                                  aaaTGCGACTGGCAAATCTGGCATTCTCTCCGGCCATGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtttt                                               <  2:432603/99‑1 (MQ=255)
                                                    aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                            >  2:807556/1‑100 (MQ=255)
                                                      tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                          <  1:2472670/100‑1 (MQ=255)
                                                       cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         >  2:1841921/1‑100 (MQ=255)
                                                       cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  2:1783699/100‑1 (MQ=255)
                                                       cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  2:3200173/100‑1 (MQ=255)
                                                        gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                        >  2:56152/1‑100 (MQ=255)
                                                          cTGTCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      <  1:2949815/100‑1 (MQ=255)
                                                           tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGTATTTTTTTCGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                     >  2:1223556/1‑100 (MQ=255)
                                                           tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAATAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                     <  2:336486/100‑1 (MQ=255)
                                                           tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                     >  2:2717530/1‑100 (MQ=255)
                                                           tGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAAAAAAAAAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACt                                     <  2:557655/100‑1 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:2032091/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  1:1500722/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  1:536903/1‑100 (MQ=255)
                                                             gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                   >  1:1298692/1‑100 (MQ=255)
                                                             gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                   >  2:1549448/1‑100 (MQ=255)
                                                                aaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                <  2:1500718/100‑1 (MQ=255)
                                                                   cTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                              <  1:108921/99‑1 (MQ=255)
                                                                    tGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                            >  2:509124/1‑100 (MQ=255)
                                                                    tGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                            >  1:1972029/1‑100 (MQ=255)
                                                                    tGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                            >  1:2239637/1‑100 (MQ=255)
                                                                    tGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                            >  1:2694439/1‑100 (MQ=255)
                                                                      gCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaa                          <  1:3143030/100‑1 (MQ=255)
                                                                       cACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGTACTTTTCCGCACAAc                         >  2:2563713/1‑100 (MQ=255)
                                                                         ctctctCCGGCCAGGTGAACCAGTCG‑GTTTTTTTGAATTTTATAAGAGCTATACCAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACt                        <  2:2800585/99‑1 (MQ=255)
                                                                          gagagagtgCCAGATGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                      >  1:1027164/9‑100 (MQ=255)
                                                                           ctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  1:1169602/1‑100 (MQ=255)
                                                                                cGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCa                >  2:3020799/1‑100 (MQ=255)
                                                                                cGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCa                >  1:2864593/1‑100 (MQ=255)
                                                                                cGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCa                <  2:524565/100‑1 (MQ=255)
                                                                                  gCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAtt              <  2:921078/100‑1 (MQ=255)
                                                                                    cAGGTGATCCAGTCGTTTTTTTTT‑AATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg            <  2:41035/100‑1 (MQ=255)
                                                                                    cAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg            >  1:2142695/1‑100 (MQ=255)
                                                                                    cAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg            <  1:1806785/100‑1 (MQ=255)
                                                                                     aGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGt           <  1:82755/100‑1 (MQ=255)
                                                                                      ggTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg          >  2:2600314/1‑100 (MQ=255)
                                                                                      ggTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg          >  1:1283047/1‑100 (MQ=255)
                                                                                      ggTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg          >  1:1711797/1‑100 (MQ=255)
                                                                                       gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         <  1:509125/100‑1 (MQ=255)
                                                                                       gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         <  1:13781/100‑1 (MQ=255)
                                                                                        tGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCt        >  2:850408/1‑100 (MQ=255)
                                                                                         gAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtg       >  1:2644140/1‑100 (MQ=255)
                                                                                           aCCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     <  2:1626582/100‑1 (MQ=255)
                                                                                           aCCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     >  1:2925740/1‑100 (MQ=255)
                                                                                           aCCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     <  1:2032099/100‑1 (MQ=255)
                                                                                            ccAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg    >  1:121828/1‑100 (MQ=255)
                                                                                              aGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAc  <  1:2472174/100‑1 (MQ=255)
                                                                                               gTCGTTTTTTTTTTTTTTTTTTAAAAGCTAGAAAAAACGGAGCTAACGCTGTTTTCTTAAGCTCTTTTCGGCACAACTTCTCTTCATTCTTgtttttttt  >  2:2996244/1‑92 (MQ=255)
                                                                                                   |                                                                                               
TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107446‑107640

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