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)
        cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtt                                                                                    <  1:929823/100‑1 (MQ=255)
          gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTTTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                                  >  2:2120269/1‑100 (MQ=255)
            tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtat                                                                                <  1:395569/100‑1 (MQ=255)
             gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtcta                                                                               >  2:795300/1‑100 (MQ=255)
                cTGGCATAACGACCTTTATACTGTCGTTTAAAATTCGACTGGCAAATCTGGCACTCTCTCCGACCATTTGAACCAGTCGTTTTTTTTGAATTTTATAaga                                                                            <  1:477291/100‑1 (MQ=255)
                cTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAaga                                                                            <  1:477293/100‑1 (MQ=255)
                  ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                          <  1:1920519/100‑1 (MQ=255)
                  ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                          <  2:7983/100‑1 (MQ=255)
                  ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                          >  1:1130106/1‑100 (MQ=255)
                   gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCt                                                                         >  1:782687/1‑100 (MQ=255)
                    cATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCta                                                                        <  2:2093970/100‑1 (MQ=255)
                     aTAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtat                                                                       >  2:2250159/1‑100 (MQ=255)
                      tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtata                                                                      <  1:1975376/100‑1 (MQ=255)
                      tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtat                                                                       <  1:1402432/99‑1 (MQ=255)
                       aaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaa                                                                     >  1:1831225/1‑100 (MQ=255)
                          gacgcttATACTATCGTATAAAATTCGACTGTAAAATTTTCCACGCTTTCCGCCCAGATGAACCAGTGGTTTTTTTTGAATTTTATAAGAGCTATaaaaa                                                                  <  2:2106470/95‑1 (MQ=255)
                           aCATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaaaaa                                                                 >  2:324281/1‑100 (MQ=255)
                            cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                >  2:483253/1‑100 (MQ=255)
                             aTTTATACTGTCGTGTAAAAGTCTACTGGCAAATTTGGCATTCTTTCCGGCCAGGTGAACCAGGCGTTGTTTTTTGAATTTTATAAAAGCTATaaaaa                                                                  <  2:332435/98‑1 (MQ=255)
                              tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACgg                                                              >  1:2218258/1‑100 (MQ=255)
                              tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACgg                                                              >  1:392813/1‑100 (MQ=255)
                                tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                            <  1:857379/100‑1 (MQ=255)
                                 ataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                           <  2:617287/100‑1 (MQ=255)
                                  taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCg                                                          <  1:1759532/100‑1 (MQ=255)
                                   aCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGc                                                           >  2:2160721/1‑98 (MQ=255)
                                   aCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCAATAAAAAACGGTGc                                                           <  2:105982/98‑1 (MQ=255)
                                    cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaaaagatcggaaga                                                        >  1:1359923/1‑90 (MQ=255)
                                    cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:2044940/1‑100 (MQ=255)
                                    cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:1825855/1‑100 (MQ=255)
                                    cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:2283921/1‑100 (MQ=255)
                                      gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                      >  2:1435081/1‑100 (MQ=255)
                                            taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                >  1:1908611/1‑100 (MQ=255)
                                            taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGCGCGAACGCTGttt                                                <  1:1169902/100‑1 (MQ=255)
                                            ctcttccgatcTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                <  1:522899/90‑1 (MQ=255)
                                              aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                              <  2:1757481/100‑1 (MQ=255)
                                              aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTc                                              <  2:29934/100‑1 (MQ=255)
                                               aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  1:2253502/100‑1 (MQ=255)
                                               aaTTCGACTGGCAAATCTGGCACTCGTTCCGGTCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  2:404440/100‑1 (MQ=255)
                                                aTTCGACTGGCAAATCTGGCACTCACTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                            <  2:2034878/100‑1 (MQ=255)
                                                  tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                          <  2:1206279/100‑1 (MQ=255)
                                                   caaCTGGCAAATCGGGCACTCTCTACGGCCAGGTGAACCAGTCGTTTTTTTTTTAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  1:899921/97‑1 (MQ=255)
                                                   cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  2:1106534/100‑1 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCTCCGGGAAGGTGAACCGGTCGGTTTTTTTGGAGTTTATAAGAGCTAGAAAAAACgg                                                              >  2:2355372/1‑78 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGCGCTATAACAAACGGTGCGAACGCTGTTTTCTTAAg                                         >  2:737532/1‑99 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  1:1147204/99‑1 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                        >  1:1562782/1‑100 (MQ=255)
                                                    gACTGGCAAATCTGGCACTCTCACCGGCCAAGTGAACCAGTCGTTTTTTTTGAATTTTAAAAAAAAAAAAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         <  2:2076699/99‑1 (MQ=255)
                                                     aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCa                                       <  2:1085270/100‑1 (MQ=255)
                                                      cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  1:720063/1‑100 (MQ=255)
                                                        ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAGAAAGCGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:914102/1‑100 (MQ=255)
                                                        ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  1:1043076/1‑100 (MQ=255)
                                                        ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  1:1215069/1‑100 (MQ=255)
                                                          cAAATCTGGTACTCTGTTCGGCCAGGTGAAGCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGTACGCTGTTTTCTTAAGCAc                                      <  2:2255228/96‑1 (MQ=255)
                                                           aaaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTc                                 <  2:1600694/100‑1 (MQ=255)
                                                             aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTTAGCACTTTTCCg                               >  2:168405/1‑100 (MQ=255)
                                                             aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               <  1:1784134/100‑1 (MQ=255)
                                                             aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               <  1:962782/100‑1 (MQ=255)
                                                                 ggCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcaca                           <  2:607646/100‑1 (MQ=255)
                                                                  gCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaa                          >  1:2119546/1‑100 (MQ=255)
                                                                   cACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         >  1:2319938/1‑100 (MQ=255)
                                                                     ctctctCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         <  1:90196/98‑1 (MQ=255)
                                                                      tctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                      <  1:1557278/100‑1 (MQ=255)
                                                                      tccctcCGGCCAGGTGAACAAGTGGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                      <  2:1599865/97‑1 (MQ=255)
                                                                       ctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  2:710669/1‑100 (MQ=255)
                                                                       ctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     <  1:638949/100‑1 (MQ=255)
                                                                        tctcCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTATAACAGCTATAAAAAACGGTGCGTACGCTGTTTTCTTAAGCACTTTTCCTCACAACTTAt                     >  2:1963981/1‑99 (MQ=255)
                                                                        tctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                    >  1:413431/1‑100 (MQ=255)
                                                                        tctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                    >  2:612914/1‑100 (MQ=255)
                                                                        tctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGATGTTTTCTTAAGCACTTTTCCGCACAACTTAt                     >  2:460340/1‑99 (MQ=255)
                                                                         ctcCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTTTAAGAGCTATAAAAAACGGTGAGAACGCTTTTTTCTTAAGCACTTTTACGCACAACTtctc                    >  2:354646/1‑99 (MQ=255)
                                                                         ctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                   <  1:751898/100‑1 (MQ=255)
                                                                          tcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCtt                  <  2:1805258/100‑1 (MQ=255)
                                                                                cAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAtt              <  2:344621/98‑1 (MQ=255)
                                                                                 aGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTc             >  2:135354/1‑98 (MQ=255)
                                                                                  ggTGAACGAGTAGTTTTTTTTTTTTTTTTTTTAAGCTATAAAAAACGGTTTGATTTTTGTTTTTTATAGTACTTTTCTGCATTACTTTTCTTCTTTAGt           >  2:1252843/1‑99 (MQ=255)
                                                                                   gTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg          <  2:2003981/99‑1 (MQ=255)
                                                                                   gTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         >  1:292057/1‑100 (MQ=255)
                                                                                     gAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtg       <  1:442597/100‑1 (MQ=255)
                                                                                       aCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     >  2:456600/1‑100 (MQ=255)
                                                                                       aCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     <  2:203924/100‑1 (MQ=255)
                                                                                        ccAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg    >  2:2362084/1‑100 (MQ=255)
                                                                                          aGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAc  <  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: