Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,428,699 A→T 5.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 5.6% 307.0 / 7.2 89intergenic (‑73/+22)ubiX/purF3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (42/42);  new base T (3/2);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

GCAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGAC  >  NC_000913/2428554‑2428846
                                                                                                                                                 |                                                                                                                                                   
gcAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCg                                                                                                                                                  <  4:102058/149‑1 (MQ=255)
   gAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCa                                                                                                                                               >  4:780341/1‑149 (MQ=255)
      ccTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCGGTCTGGCAcc                                                                                                                                             <  1:889539/148‑1 (MQ=255)
       cTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCt                                                                                                                                           >  4:20178/1‑149 (MQ=255)
        tGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCt                                                                                                                                           >  2:320748/1‑148 (MQ=255)
        tGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCt                                                                                                                                           >  1:740219/1‑148 (MQ=255)
        tGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCt                                                                                                                                           >  1:657633/1‑148 (MQ=255)
         gTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTc                                                                                                                                          >  1:88769/1‑148 (MQ=255)
         gTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCa                                                                                                                                         <  3:622378/149‑1 (MQ=255)
           aataaGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAg                                                                                                                                        >  2:141178/1‑148 (MQ=255)
              aaGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGa                                                                                                                                     <  2:319948/148‑1 (MQ=255)
               aGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGaa                                                                                                                                    >  1:417502/1‑148 (MQ=255)
               aGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAg                                                                                                                                   <  3:358310/149‑1 (MQ=255)
                       ccATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGcatcat                                                                                                                            <  1:50694/148‑1 (MQ=255)
                         aTAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATcc                                                                                                                          >  2:161495/1‑148 (MQ=255)
                           aaaTCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCtt                                                                                                                        <  1:703745/148‑1 (MQ=255)
                               cgcgCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGTCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTa                                                                                                                   <  3:335973/149‑1 (MQ=255)
                               cgcgCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTa                                                                                                                   <  3:297260/149‑1 (MQ=255)
                                       tGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCAttt                                                                                                            <  1:320748/148‑1 (MQ=255)
                                         gCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCga                                                                                                         >  4:761969/1‑149 (MQ=255)
                                          cACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCgag                                                                                                        >  4:830111/1‑149 (MQ=255)
                                           aCCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCgag                                                                                                        <  2:96774/148‑1 (MQ=255)
                                            ccGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGAt                                                                                                      <  3:20178/149‑1 (MQ=255)
                                                    gCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACt                                                                                              >  3:679438/1‑149 (MQ=255)
                                                     ccTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACtt                                                                                             <  4:724832/149‑1 (MQ=255)
                                                        aCAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCg                                                                                           >  1:669010/1‑148 (MQ=255)
                                                         cAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGtt                                                                                         <  3:547476/149‑1 (MQ=255)
                                                          aaTGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACt                                                                                              >  1:710970/1‑143 (MQ=255)
                                                          aaTGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACt                                                                                              <  2:710970/143‑1 (MQ=255)
                                                               gtcatgtCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGCCTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACg                                                                                   >  4:607943/7‑149 (MQ=255)
                                                                    ttCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCa                                                                              <  3:122102/149‑1 (MQ=255)
                                                                        tAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGTGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTc                                                                                        <  2:675233/135‑1 (MQ=255)
                                                                        tAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTc                                                                                        >  1:675233/1‑135 (MQ=255)
                                                                        tAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTg                                                                           >  2:21593/1‑148 (MQ=255)
                                                                         aaaTAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGc                                                                          >  1:115288/1‑148 (MQ=255)
                                                                           aTAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCttt                                                                       <  4:800456/149‑1 (MQ=255)
                                                                            tAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCttt                                                                       >  2:933511/1‑148 (MQ=255)
                                                                                tGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGt                                                                  >  3:51021/1‑149 (MQ=255)
                                                                                 gCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGtc                                                                 <  3:780341/149‑1 (MQ=255)
                                                                                 gCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGt                                                                  <  2:417502/148‑1 (MQ=255)
                                                                                  cccTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGtc                                                                 <  1:161495/148‑1 (MQ=255)
                                                                                   ccTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGtca                                                                >  2:962514/1‑148 (MQ=255)
                                                                                       tCAGACTTGCAGGCAGACTTTGCAGGATTACGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCAtt                                                           >  3:569235/1‑149 (MQ=255)
                                                                                             ttGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTacgta                                                      <  1:560944/148‑1 (MQ=255)
                                                                                              tGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCa                                                                                                                                         >  1:908203/1‑64 (MQ=255)
                                                                                              tGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCa                                                                                                                                         <  2:908203/64‑1 (MQ=255)
                                                                                              tGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGa                                                                                                                                     >  2:635552/1‑68 (MQ=255)
                                                                                              tGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGa                                                                                                                                     <  1:635552/68‑1 (MQ=255)
                                                                                              tGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGtc                                                                 <  2:53106/136‑1 (MQ=255)
                                                                                              tGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGtc                                                                 >  1:53106/1‑136 (MQ=255)
                                                                                              tGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTacgtaac                                                    <  4:51021/149‑1 (MQ=255)
                                                                                                   gcagACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTAt                                                <  1:440941/148‑1 (MQ=255)
                                                                                                    cagACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATc                                               >  2:148037/1‑148 (MQ=255)
                                                                                                    cagACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATc                                               >  2:948744/1‑148 (MQ=255)
                                                                                                     agACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGa                                             >  4:847268/1‑149 (MQ=255)
                                                                                                         tttGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGa                                          >  1:921753/1‑148 (MQ=255)
                                                                                                          ttGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGaa                                         >  2:225705/1‑148 (MQ=255)
                                                                                                          ttGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGaa                                         >  1:117161/1‑148 (MQ=255)
                                                                                                              aGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                     <  2:657633/148‑1 (MQ=255)
                                                                                                              aGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCg                                     <  2:740219/148‑1 (MQ=255)
                                                                                                               ggATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGa                                    <  1:21593/148‑1 (MQ=255)
                                                                                                               ggATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGa                                    >  2:73172/1‑148 (MQ=255)
                                                                                                               ggATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGa                                    <  1:868462/148‑1 (MQ=255)
                                                                                                               ggATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGa                                    <  1:73172/148‑1 (MQ=255)
                                                                                                               ggATTGCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                   <  4:302341/149‑1 (MQ=255)
                                                                                                                   ttCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTa                                >  1:648306/1‑148 (MQ=255)
                                                                                                                      ggcggGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                   >  4:71613/1‑142 (MQ=255)
                                                                                                                      ggcggGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                   <  3:71613/142‑1 (MQ=255)
                                                                                                                      ggcggGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGGCATCATTACGTAACGTATCGAGGAAATCGAGGTAGcc                             <  2:117161/148‑1 (MQ=255)
                                                                                                                       gcggGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGccc                            <  1:619767/148‑1 (MQ=255)
                                                                                                                       gcggGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCt                           <  4:742312/149‑1 (MQ=255)
                                                                                                                         gggAGTTGCAAGTCAGGGTGCCTGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGGGGTAGCCCTGa                         <  3:761969/149‑1 (MQ=255)
                                                                                                                          ggAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGAt                        >  4:863308/1‑149 (MQ=255)
                                                                                                                              ttGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCaa                     <  1:728228/148‑1 (MQ=255)
                                                                                                                              ttGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAAc                    >  4:843972/1‑149 (MQ=255)
                                                                                                                                  aaGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCt                >  4:407493/1‑149 (MQ=255)
                                                                                                                                   gctgaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCtt               >  4:533078/5‑149 (MQ=255)
                                                                                                                                   aGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCtt               <  4:720954/149‑1 (MQ=255)
                                                                                                                                   aGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCt                <  1:49044/148‑1 (MQ=255)
                                                                                                                                      gaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGaa                                         >  1:731361/2‑120 (MQ=255)
                                                                                                                                      gaGGGTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGaa                                         <  2:731361/119‑1 (MQ=255)
                                                                                                                                        gggTGCCGGTCTGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGt           >  2:545712/1‑148 (MQ=255)
                                                                                                                                        gggTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTg          >  3:153712/1‑149 (MQ=255)
                                                                                                                                          gTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGa         <  1:962514/148‑1 (MQ=255)
                                                                                                                                           tGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGAc        <  1:29131/148‑1 (MQ=255)
                                                                                                                                            gCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACg       >  1:83475/1‑148 (MQ=255)
                                                                                                                                              cAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAg    >  3:856028/1‑149 (MQ=255)
                                                                                                                                               aGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAg    >  1:630952/1‑148 (MQ=255)
                                                                                                                                                gACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGa   >  2:817473/1‑148 (MQ=255)
                                                                                                                                                 aCCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGAc  >  2:462330/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAGAC  >  NC_000913/2428554‑2428846

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