Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,428,699 A→T 8.6% intergenic (‑73/+22) ubiX ← / ← purF 3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase/amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,428,6990AT 8.6% 232.7 / 11.2 70intergenic (‑73/+22)ubiX/purF3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/36);  new base T (4/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACG  >  NC_000913/2428555‑2428847
                                                                                                                                                |                                                                                                                                                    
cAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCACGGTGCCAGACCgg                                                                                                                                                  >  4:527985/1‑149 (MQ=255)
  gAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCa                                                                                                                                                >  3:489607/1‑149 (MQ=255)
          aataaGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGc                                                                                                                                        <  3:259612/149‑1 (MQ=255)
             aaGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGaa                                                                                                                                     >  1:187757/1‑149 (MQ=255)
              aGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGaat                                                                                                                                    <  2:528146/149‑2 (MQ=255)
              aGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAg                                                                                                                                    <  4:448133/149‑1 (MQ=255)
              aGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAg                                                                                                                                    <  1:422436/149‑1 (MQ=255)
                       cATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATcc                                                                                                                           >  1:109283/1‑149 (MQ=255)
                       cATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATcc                                                                                                                           >  3:374644/1‑149 (MQ=255)
                       cATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATcc                                                                                                                           <  2:378618/149‑1 (MQ=255)
                       cATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATcc                                                                                                                           >  2:268145/1‑149 (MQ=255)
                        aTAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCt                                                                                                                          >  1:419643/1‑149 (MQ=255)
                        aTAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCt                                                                                                                          <  2:143026/149‑1 (MQ=255)
                            aTCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGt                                                                                                                      >  2:516054/1‑149 (MQ=255)
                             tCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGtt                                                                                                                     <  3:576770/149‑1 (MQ=255)
                               gcgcCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTTCAGGCAGACTTTGCAGGATGTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTAt                                                                                                                   <  3:381614/149‑1 (MQ=255)
                                   cgcTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCAt                                                                                                               >  4:352774/1‑149 (MQ=255)
                                       ggCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCg                                                                                                           >  4:213737/1‑149 (MQ=255)
                                        gCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCga                                                                                                          <  4:317637/149‑1 (MQ=255)
                                        gCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCga                                                                                                          <  4:378153/149‑1 (MQ=255)
                                        gCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCga                                                                                                          >  3:461281/1‑149 (MQ=255)
                                          aCCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCgaga                                                                                                        >  3:205030/1‑149 (MQ=255)
                                             gCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGAttt                                                                                                     >  3:289203/1‑149 (MQ=255)
                                                     cTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTc                                                                                             >  2:593042/1‑149 (MQ=255)
                                                                      aTAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTg                                                                            >  3:397231/1‑149 (MQ=255)
                                                                          aTAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCttt                                                                        <  2:158988/149‑1 (MQ=255)
                                                                           tAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTg                                                                       <  4:397231/149‑1 (MQ=255)
                                                                                gCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTCTGCATTTCCAGATTTTCCACTTCGTTCTGCCGTCGCACTGCTTTGCCGtc                                                                  <  4:133496/149‑1 (MQ=255)
                                                                                gCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGtc                                                                  <  4:374644/149‑1 (MQ=255)
                                                                                 cccTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCt                                                                                        <  2:105738/126‑1 (MQ=255)
                                                                                 cccTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCt                                                                                        >  1:105738/1‑126 (MQ=255)
                                                                                 cccTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGtca                                                                 >  1:544281/1‑149 (MQ=255)
                                                                                    tgacAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATATCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGtcatca                                                              >  4:415358/4‑149 (MQ=255)
                                                                                    tGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGtcatca                                                              >  1:448844/1‑149 (MQ=255)
                                                                                       cAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTa                                                           <  2:419643/149‑1 (MQ=255)
                                                                                        aGACTTGCAGGCAGACTTTGCAGGATTTCGGGGGGAGTTGCAAGTCAGGGTGCCAGACCGGCAGACTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTac                                                          <  2:187757/149‑1 (MQ=255)
                                                                                        aGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGcatcat                                                                                                                             <  2:177550/82‑1 (MQ=255)
                                                                                        aGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGcatcat                                                                                                                             >  1:177550/1‑82 (MQ=255)
                                                                                        aGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTac                                                          >  1:50870/1‑149 (MQ=255)
                                                                                        aGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTac                                                          <  4:289203/149‑1 (MQ=255)
                                                                                            ttGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTacgtaa                                                      <  1:516054/149‑1 (MQ=255)
                                                                                             tGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTacgtaac                                                     <  3:213737/149‑1 (MQ=255)
                                                                                                    agACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGa                                              >  3:206642/1‑149 (MQ=255)
                                                                                                    agACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGa                                              <  1:161115/149‑1 (MQ=255)
                                                                                                      aCTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAgg                                            >  2:66132/1‑149 (MQ=255)
                                                                                                        tttGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGa                                              >  3:435645/1‑145 (MQ=255)
                                                                                                        tttGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGa                                              <  4:435645/145‑1 (MQ=255)
                                                                                                        tttGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGaa                                          >  4:119575/1‑149 (MQ=255)
                                                                                                         ttGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCCTTGGCGTCATCATTACGTAACGTATCGAGGaaa                                         <  3:527985/149‑1 (MQ=255)
                                                                                                               gATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAgg                                   <  4:489607/149‑1 (MQ=255)
                                                                                                                        gggAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGa                          <  2:308725/149‑1 (MQ=255)
                                                                                                                        gggAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGa                          <  4:514203/149‑1 (MQ=255)
                                                                                                                         ggAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGATGTAGCCCTGAt                         <  3:352774/149‑1 (MQ=255)
                                                                                                                         ggAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGAt                         <  1:207296/149‑1 (MQ=255)
                                                                                                                           aGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCa                       <  2:169290/149‑1 (MQ=255)
                                                                                                                            gTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCaa                      <  4:461281/149‑1 (MQ=255)
                                                                                                                            gTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCaa                      <  1:270615/149‑1 (MQ=255)
                                                                                                                               gCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACAt                   >  3:469792/1‑149 (MQ=255)
                                                                                                                                    tgatggTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                    >  2:122993/5‑127 (MQ=255)
                                                                                                                                    tgaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                    <  1:122993/125‑1 (MQ=255)
                                                                                                                                    tgaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                    >  4:378632/3‑127 (MQ=255)
                                                                                                                                    tgaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                    <  3:378632/125‑1 (MQ=255)
                                                                                                                                    tgaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTg              >  4:51738/3‑149 (MQ=255)
                                                                                                                                    tCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTg              <  4:469792/149‑1 (MQ=255)
                                                                                                                                     gaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTgg             >  2:221195/2‑149 (MQ=255)
                                                                                                                                     cAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTgg             <  4:573082/149‑1 (MQ=255)
                                                                                                                                     cAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCATTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCGTTgg             >  4:305132/1‑149 (MQ=255)
                                                                                                                                      aGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt            <  2:109283/149‑1 (MQ=255)
                                                                                                                                        ggTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGa          <  1:593042/149‑1 (MQ=255)
                                                                                                                                              aGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGa    >  3:170570/1‑149 (MQ=255)
                                                                                                                                                aCCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACg  >  1:543144/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGACG  >  NC_000913/2428555‑2428847

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: