Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,428,697 A→G 5.7% 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 5.7% 280.8 / 6.0 88intergenic (‑71/+24)ubiX/purF3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (41/41);  new base G (3/2);  total (44/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 = 3.18e-01

TCGCGCAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAG  >  NC_000913/2428550‑2428844
                                                                                                                                                   |                                                                                                                                                   
tCGCGCAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAg                                                                                                                                                    >  3:668850/1‑149 (MQ=255)
    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                                                                                            <  2:710970/143‑1 (MQ=255)
                                                              aaTGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACt                                                                                            >  1:710970/1‑143 (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                                                                                                                                       <  2:908203/64‑1 (MQ=255)
                                                                                                  tGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCa                                                                                                                                       >  1:908203/1‑64 (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                                  >  2:73172/1‑148 (MQ=255)
                                                                                                                   ggATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGa                                  <  1:21593/148‑1 (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                                 <  3:71613/142‑1 (MQ=255)
                                                                                                                          ggcggGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAg                                 >  4:71613/1‑142 (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)
                                                                                                                                                   |                                                                                                                                                   
TCGCGCAGAACCTGTAATAAGCGCACGCCATAAATCGCGCCGCTGGCACCGCTGATGCCTACAATGAGTCGTTTCATAAATAGTTGCCCTGTCAGACTTGCAGGCAGACTTTGCAGGATTTCGGCGGGAGTTGCAAGTCAGGGTGCCAGACCGGCACCCTCAGCGAAGGCATCATCCTTCGTTATGCATTTCGAGATTTTCCACTTCGTTCTGACGTTGCACTGCTTTGGCGTCATCATTACGTAACGTATCGAGGAAATCGAGGTAGCCCTGATCAACATCTTTGGTGACGTAG  >  NC_000913/2428550‑2428844

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