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.90.2% 301.6 / 25.6 82intergenic (+71/‑160)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base . (34/40);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

AGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107450‑107640
                                                                                               |                                                                                               
aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                             >  2:366238/1‑100 (MQ=255)
 gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                            >  2:579417/1‑100 (MQ=255)
 gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttt                                                                                              <  1:2148741/98‑1 (MQ=255)
   gcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           <  2:1984882/99‑1 (MQ=255)
    ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                         <  1:218119/100‑2 (MQ=255)
    ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                           >  2:1442699/1‑98 (MQ=255)
       tCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttgaat                                                                                      <  2:1068772/100‑5 (MQ=255)
        cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAtt                                                                                    <  1:929823/100‑1 (MQ=255)
          gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTTTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAtttt                                                                                  >  2:2120269/1‑100 (MQ=255)
            tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtat                                                                                <  1:395569/100‑1 (MQ=255)
             gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtcta                                                                               >  2:795300/1‑100 (MQ=255)
                cTGGCATAACGACCTTTATACTGTCGTTTAAAATTCGACTGGCAAATCTGGCACTCTCTCCGACCATTTGAACCAGTCG‑TTTTTTTTGAATTTTATAaga                                                                            <  1:477291/100‑1 (MQ=255)
                cTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAaga                                                                            <  1:477293/100‑1 (MQ=255)
                  ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                          <  1:1920519/100‑1 (MQ=255)
                  ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                          <  2:7983/100‑1 (MQ=255)
                  ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                          >  1:1130106/1‑100 (MQ=255)
                   gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCt                                                                         >  1:782687/1‑100 (MQ=255)
                    cATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCta                                                                        <  2:2093970/100‑1 (MQ=255)
                     aTAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCtat                                                                       >  2:2250159/1‑100 (MQ=255)
                      tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCtata                                                                      <  1:1975376/100‑1 (MQ=255)
                      tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCtat                                                                       <  1:1402432/99‑1 (MQ=255)
                       aaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaa                                                                     >  1:1831225/1‑100 (MQ=255)
                          gacgcttATACTATCGTATAAAATTCGACTGTAAAATTTTCCACGCTTTCCGCCCAGATGAACCAGT‑GGTTTTTTTTGAATTTTATAAGAGCTATaaaaa                                                                  <  2:2106470/95‑1 (MQ=255)
                           aCATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaaaaa                                                                 >  2:324281/1‑100 (MQ=255)
                            cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                >  2:483253/1‑100 (MQ=255)
                             aTTTATACTGTCGTGTAAAAGTCTACTGGCAAATTTGGCATTCTTTCCGGCCAGGTGAACCAGGCGTTGTTTTTTGAATTTTATAAAAGCTATaaaaa                                                                  <  2:332435/98‑1 (MQ=255)
                              tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACgg                                                              >  1:2218258/1‑100 (MQ=255)
                              tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACgg                                                              >  1:392813/1‑100 (MQ=255)
                                tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            <  1:857379/100‑1 (MQ=255)
                                 ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                           <  2:617287/100‑1 (MQ=255)
                                  taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCg                                                          <  1:1759532/100‑1 (MQ=255)
                                   aCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                           >  2:2160721/1‑98 (MQ=255)
                                   aCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCAATAAAAAACGGTGc                                                           <  2:105982/98‑1 (MQ=255)
                                    cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATaaaaagatcggaaga                                                        >  1:1359923/1‑90 (MQ=255)
                                    cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:2044940/1‑100 (MQ=255)
                                    cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:1825855/1‑100 (MQ=255)
                                    cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:2283921/1‑100 (MQ=255)
                                      gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                      >  2:1435081/1‑100 (MQ=255)
                                            taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                >  1:1908611/1‑100 (MQ=255)
                                            taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGCGCGAACGCTGttt                                                <  1:1169902/100‑1 (MQ=255)
                                            ctcttccgatcTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                <  1:522899/90‑1 (MQ=255)
                                              aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                              <  2:1757481/100‑1 (MQ=255)
                                              aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                              <  2:29934/100‑1 (MQ=255)
                                               aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  1:2253502/100‑1 (MQ=255)
                                               aaTTCGACTGGCAAATCTGGCACTCGTTCCGGTCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  2:404440/100‑1 (MQ=255)
                                                aTTCGACTGGCAAATCTGGCACTCACTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                            <  2:2034878/100‑1 (MQ=255)
                                                  tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                          <  2:1206279/100‑1 (MQ=255)
                                                   caaCTGGCAAATCGGGCACTCTCTACGGCCAGGTGAACCAGTCGTTTTTTTTTTAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  1:899921/97‑1 (MQ=255)
                                                   cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  2:1106534/100‑1 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCTCCGGGAAGGTGAACCGGTCG‑GTTTTTTTGGAGTTTATAAGAGCTAGAAAAAACgg                                                              >  2:2355372/1‑78 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGCGCTATAACAAACGGTGCGAACGCTGTTTTCTTAAg                                         >  2:737532/1‑99 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  1:1147204/99‑1 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                        >  1:1562782/1‑100 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCACCGGCCAAGTGAACCAGTCG‑TTTTTTTTGAATTTTAAAAAAAAAAAAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  2:2076699/99‑1 (MQ=255)
                                                     aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCa                                       <  2:1085270/100‑1 (MQ=255)
                                                      cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  1:720063/1‑100 (MQ=255)
                                                        ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAGAAAGCGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:914102/1‑100 (MQ=255)
                                                        ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  1:1043076/1‑100 (MQ=255)
                                                        ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  1:1215069/1‑100 (MQ=255)
                                                          cAAATCTGGTACTCTGTTCGGCCAGGTGAAGCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGTACGCTGTTTTCTTAAGCAc                                      <  2:2255228/96‑1 (MQ=255)
                                                           aaaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                 <  2:1600694/100‑1 (MQ=255)
                                                             aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTTAGCACTTTTCCg                               >  2:168405/1‑100 (MQ=255)
                                                             aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               <  1:1784134/100‑1 (MQ=255)
                                                             aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               <  1:962782/100‑1 (MQ=255)
                                                                 ggCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcaca                           <  2:607646/100‑1 (MQ=255)
                                                                  gCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaa                          >  1:2119546/1‑100 (MQ=255)
                                                                   cACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         >  1:2319938/1‑100 (MQ=255)
                                                                     ctctctCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         <  1:90196/98‑1 (MQ=255)
                                                                      tctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                      <  1:1557278/100‑1 (MQ=255)
                                                                      tccctcCGGCCAGGTGAACAAGT‑GGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                      <  2:1599865/97‑1 (MQ=255)
                                                                       ctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  2:710669/1‑100 (MQ=255)
                                                                       ctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     <  1:638949/100‑1 (MQ=255)
                                                                        tctcCGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTATAACAGCTATAAAAAACGGTGCGTACGCTGTTTTCTTAAGCACTTTTCCTCACAACTTAt                     >  2:1963981/1‑99 (MQ=255)
                                                                        tctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                    >  1:413431/1‑100 (MQ=255)
                                                                        tctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                    >  2:612914/1‑100 (MQ=255)
                                                                        tctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGATGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  2:460340/1‑99 (MQ=255)
                                                                         ctcCGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTGAGAACGCTTTTTTCTTAAGCACTTTTACGCACAACTtctc                    >  2:354646/1‑99 (MQ=255)
                                                                         ctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                   <  1:751898/100‑1 (MQ=255)
                                                                          tcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCtt                  <  2:1805258/100‑1 (MQ=255)
                                                                                cAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAtt              <  2:344621/98‑1 (MQ=255)
                                                                                 aGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTc             >  2:135354/1‑98 (MQ=255)
                                                                                  ggTGAACGAGTAGTTTTTTTTTTTTTTTTTTTA‑AGCTATAAAAAACGGTTTGATTTTTGTTTTTTATAGTACTTTTCTGCATTACTTTTCTTCTTTAGt           >  2:1252843/1‑99 (MQ=255)
                                                                                   gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg          <  2:2003981/99‑1 (MQ=255)
                                                                                   gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         >  1:292057/1‑100 (MQ=255)
                                                                                     gAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtg       <  1:442597/100‑1 (MQ=255)
                                                                                       aCCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     >  2:456600/1‑100 (MQ=255)
                                                                                       aCCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     <  2:203924/100‑1 (MQ=255)
                                                                                        ccAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg    >  2:2362084/1‑100 (MQ=255)
                                                                                          aGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAc  <  1:483255/100‑1 (MQ=255)
                                                                                               |                                                                                               
AGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107450‑107640

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: ‑