Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,428,697 A→G 8.6% intergenic (‑71/+24) 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,6970AG 8.6% 214.8 / 9.8 70intergenic (‑71/+24)ubiX/purF3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/36);  new base G (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 = 5.62e-01

ATCGCGCAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGA  >  NC_000913/2428549‑2428845
                                                                                                                                                    |                                                                                                                                                    
aTCGCGCAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCa                                                                                                                                                      <  4:158089/149‑1 (MQ=255)
 tCGCGCAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAg                                                                                                                                                     >  1:169290/1‑149 (MQ=255)
      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                                                                                                                         >  2:268145/1‑149 (MQ=255)
                             cATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATcc                                                                                                                         <  2:378618/149‑1 (MQ=255)
                             cATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATcc                                                                                                                         >  3:374644/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:378153/149‑1 (MQ=255)
                                              gCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCga                                                                                                        <  4:317637/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                                            <  1:161115/149‑1 (MQ=255)
                                                                                                          agACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGa                                            >  3:206642/1‑149 (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                    <  1:270615/149‑1 (MQ=255)
                                                                                                                                  gTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCaa                    <  4:461281/149‑1 (MQ=255)
                                                                                                                                     gCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACAt                 >  3:469792/1‑149 (MQ=255)
                                                                                                                                          tgatggTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                  >  2:122993/5‑127 (MQ=255)
                                                                                                                                          tgaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                  >  4:378632/3‑127 (MQ=255)
                                                                                                                                          tgaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                  <  1:122993/125‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
ATCGCGCAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGA  >  NC_000913/2428549‑2428845

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: