Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,673,034 T→C 13.6% F330L (TTT→CTT)  hcaD → phenylpropionate dioxygenase, ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,673,0340TC13.6% 1615.9 / 200.4 618F330L (TTT→CTT) hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base T (245/286);  new base C (39/45);  total (287/331)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

Only 100 of 644 total aligned reads displayed.

TCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGAT  >  NC_000913/2672886‑2673182
                                                                                                                                                    |                                                                                                                                                    
tCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGt                                                                                                                                                      >  1:1405077/1‑149 (MQ=255)
  ttGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGttt                                                                                                                                                    <  4:1243565/149‑1 (MQ=255)
  ttGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGttt                                                                                                                                                    >  4:1849697/1‑149 (MQ=255)
  ttGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGttt                                                                                                                                                    >  1:784636/1‑149 (MQ=255)
   tGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTa                                                                                                                                                   <  4:1710202/149‑1 (MQ=255)
     aTAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTAtt                                                                                                                                                 <  2:6991756/149‑1 (MQ=255)
      tAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTg                                                                                                                                                <  3:7004698/149‑1 (MQ=255)
           gTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGat                                                                                                                                           <  2:4846144/149‑1 (MQ=255)
               aCTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGc                                                                                                                                       <  4:2301690/149‑1 (MQ=255)
                 tACACCGCTGCGAAAGCTGGGGAAACGCCAATAACCAGGCGCAAATTGCCGCTGACGCAATGTTAGGGCTAACGCTACAGCTACTGCCGTCGCAGTGGTTCTGGAGCGCTCAGTACAGTGCTAACTTCCAGTTTATTGCCGATATGCGt                                                                                                                                     <  2:5974815/149‑1 (MQ=255)
                 tACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTAtt                                                                                                                                                 >  2:7156927/1‑137 (MQ=255)
                   cacCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTgg                                                                                                                                   <  3:3946843/149‑1 (MQ=255)
                    acCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTgg                                                                                                                                   <  2:2786521/148‑1 (MQ=255)
                    acCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGc                                                                                                                                  >  4:3469783/1‑149 (MQ=255)
                     ccGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCg                                                                                                                                 <  2:1547415/149‑1 (MQ=255)
                            gAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTg                                                                                                                          <  1:497698/149‑1 (MQ=255)
                                gCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCtt                                                                                                                      <  3:5419789/149‑1 (MQ=255)
                                 cTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCttt                                                                                                                     >  1:6229465/1‑149 (MQ=255)
                                 cTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCttt                                                                                                                     <  3:1475370/149‑1 (MQ=255)
                                      aaaaCGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGt                                                                                                                <  2:2219754/149‑1 (MQ=255)
                                       aaaCGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTg                                                                                                               <  1:6989291/149‑1 (MQ=255)
                                        aaCGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTgg                                                                                                              >  3:2383671/1‑149 (MQ=255)
                                         aCGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGc                                                                                                             >  4:6642122/1‑149 (MQ=255)
                                           gCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCaa                                                                                                           >  3:6491543/1‑149 (MQ=255)
                                             cAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAAcc                                                                                                         >  1:2094093/1‑149 (MQ=255)
                                             cAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAAcc                                                                                                         >  1:5854851/1‑149 (MQ=255)
                                              aaTAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAAccc                                                                                                        >  4:742592/1‑149 (MQ=255)
                                              aaTAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAAccc                                                                                                        >  2:5501915/1‑149 (MQ=255)
                                              aaTAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAAccc                                                                                                        >  2:2634262/1‑149 (MQ=255)
                                                          cAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGaa                                                                                            <  1:5407000/149‑1 (MQ=255)
                                                          cAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGaa                                                                                            <  4:2450683/149‑1 (MQ=255)
                                                          cAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGaa                                                                                            <  2:3634036/149‑1 (MQ=255)
                                                          cAAATTGCCACTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGaa                                                                                            <  4:4572374/149‑1 (MQ=255)
                                                             aTTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGc                                                                                         <  4:58219/149‑1 (MQ=255)
                                                               tgccgctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCCGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGa                                                                                       <  4:4039902/149‑1 (MQ=255)
                                                                gccgctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGAt                                                                                      >  2:5070271/1‑149 (MQ=255)
                                                                 ccgctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGTCATGGTTCTGGAGCGATCAGTACAGTGATAACTTGCAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGAtt                                                                                     >  4:4385202/1‑149 (MQ=255)
                                                                 ccgctgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAAccc                                                                                                        >  1:1211796/1‑130 (MQ=255)
                                                                     tgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                                                                 <  3:5712495/149‑1 (MQ=255)
                                                                     tgccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                                                                 <  4:6491543/149‑1 (MQ=255)
                                                                      gccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGtt                                                                                <  4:4220487/149‑1 (MQ=255)
                                                                      ccgatcctGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGtt                                                                                <  1:3091387/142‑1 (MQ=255)
                                                                       ccgcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGttt                                                                               >  1:6736790/1‑149 (MQ=255)
                                                                         gcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTaa                                                                             >  3:5015496/1‑149 (MQ=255)
                                                                         gcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTaa                                                                             >  1:1581059/1‑149 (MQ=255)
                                                                           aaTGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATc                                                                           <  1:3875945/149‑1 (MQ=255)
                                                                              gTTAGGGCTACCGCTACCCCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTTGGGATTTGCGTGGCGATGACTGGGTTTGTCGTGGGCACCCGGAAACTCAGAAGGCGGTTTGGTTTagattgg                                                                        >  2:3091387/1‑143 (MQ=255)
                                                                                  gggCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTa                                                                              >  4:2974885/1‑139 (MQ=255)
                                                                                  gggCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCaaaa                                                                    <  4:3180517/149‑1 (MQ=255)
                                                                                   ggCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTCGTCGTGGCAACCCGGAAACGCAGAAGGCGATTTGGTTTAATCTGCAAAAc                                                                   <  2:3321769/149‑1 (MQ=255)
                                                                                    gctaccgctaccGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACg                                                                  <  3:1727651/149‑1 (MQ=255)
                                                                                      taccgctaccgcTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGc                                                                <  3:1929085/149‑1 (MQ=255)
                                                                                          gctaccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAAt                                                                            <  1:4162048/133‑1 (MQ=255)
                                                                                          gctaccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAAt                                                                            >  2:4162048/1‑133 (MQ=255)
                                                                                          gctaccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGt                                                              <  4:7034442/147‑1 (MQ=255)
                                                                                          gctaccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGc                                                            >  1:6427991/1‑149 (MQ=255)
                                                                                             accgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTa                                                         <  3:3198204/149‑1 (MQ=255)
                                                                                              ccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGTCGTGCTTAt                                                        <  1:3496572/149‑1 (MQ=255)
                                                                                               cgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATc                                                       >  2:6677174/1‑149 (MQ=255)
                                                                                                gctacTGCCGCCGCCGTGGTTCTGGAGCTATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTTGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCTGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCg                                                      <  4:1402273/149‑1 (MQ=255)
                                                                                                 ctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCgg                                                     >  4:5032588/1‑149 (MQ=255)
                                                                                                    cTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGc                                                  <  3:3469783/149‑1 (MQ=255)
                                                                                                    cTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGc                                                  <  3:951233/149‑1 (MQ=255)
                                                                                                          ccgccgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCAc                                            >  1:604779/1‑149 (MQ=255)
                                                                                                            gccgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                                                                 >  3:1833659/1‑110 (MQ=255)
                                                                                                            gccgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGc                                          >  2:3953496/1‑149 (MQ=255)
                                                                                                              cgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTg                                                                                                               <  1:2431628/78‑1 (MQ=255)
                                                                                                              cgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTg                                        >  3:6144141/1‑149 (MQ=255)
                                                                                                               gTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTTCTTATCGGTGCGGTCACGCTGa                                       >  3:4378101/1‑149 (MQ=255)
                                                                                                                tGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGaa                                      >  2:1419491/1‑149 (MQ=255)
                                                                                                                  gTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGc                                                                >  3:5786910/1‑121 (MQ=255)
                                                                                                                     cTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAgg                                 >  1:2915281/1‑149 (MQ=255)
                                                                                                                     cTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAgg                                 <  2:1957481/149‑1 (MQ=255)
                                                                                                                     cTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAgg                                 <  3:797897/149‑1 (MQ=255)
                                                                                                                           cGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTCTGTCGTGGCAACCCTGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTAAATCAGGGGCGta                           <  4:6167184/149‑2 (MQ=255)
                                                                                                                            gATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTga                          <  3:1284455/149‑1 (MQ=255)
                                                                                                                             aTCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCttt                                                                                                                     >  2:6438897/1‑57 (MQ=255)
                                                                                                                             aTCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGt                                                              >  4:2050967/1‑112 (MQ=255)
                                                                                                                             aTCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTgag                         >  3:4500918/1‑149 (MQ=255)
                                                                                                                             aTCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTgag                         >  3:2988908/1‑149 (MQ=255)
                                                                                                                                 gTACAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTc                     <  1:4448585/149‑1 (MQ=255)
                                                                                                                                   aCAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGc                   >  2:1031594/1‑149 (MQ=255)
                                                                                                                                    cAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTAt                                                        >  4:2719113/1‑111 (MQ=255)
                                                                                                                                    cAGTGATAACTTACAGCTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGcc                  <  3:2858392/149‑1 (MQ=255)
                                                                                                                                     aGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGccc                 >  4:533559/1‑149 (MQ=255)
                                                                                                                                       tGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGAttt                                                                                    <  1:4626093/80‑1 (MQ=255)
                                                                                                                                       tGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGAttt                                                                                    >  2:4626093/1‑80 (MQ=255)
                                                                                                                                         aTAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAAtt             >  3:2167423/1‑149 (MQ=255)
                                                                                                                                          tAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTc            <  4:5647686/149‑1 (MQ=255)
                                                                                                                                          tAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCAGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTc            >  1:2742835/1‑149 (MQ=255)
                                                                                                                                            aCTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTTGCCCAATTCGc          <  3:4166971/149‑1 (MQ=255)
                                                                                                                                              ttACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCaa        >  1:3326527/1‑149 (MQ=255)
                                                                                                                                               tACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCATAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCaaa       <  4:4882351/149‑1 (MQ=255)
                                                                                                                                               tACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCaaa       <  2:6093313/149‑1 (MQ=255)
                                                                                                                                                 cAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATg     >  1:1008101/1‑149 (MQ=255)
                                                                                                                                                 cAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATg     >  2:1699256/1‑149 (MQ=255)
                                                                                                                                                   gTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGa   >  4:3916101/1‑149 (MQ=255)
                                                                                                                                                    tttATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAggg                                <  2:847983/119‑1 (MQ=255)
                                                                                                                                                    tttATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGAt  >  4:3676549/1‑149 (MQ=255)
                                                                                                                                                    tttATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGAt  <  3:5296302/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCTTGATAATGGTGCACTACACCGCTGCGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGAT  >  NC_000913/2672886‑2673182

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: