Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,362,578 G→T 50.4% noncoding (49/76 nt) pheU ← tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,362,5780GT50.4% 11.7 / 86.0 84noncoding (49/76 nt)pheUtRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base G (25/17);  new base T (18/24);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTG  >  NC_000913/4362430‑4362724
                                                                                                                                                    |                                                                                                                                                  
ttGCGGGTGATGTGGGGAGAATCTGGTTGAGGTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGGCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTAAATGAATTTGGTGCCCGGACTCGGAATCGAACCAAg                                                                                                                                                    <  2:380059/149‑1 (MQ=255)
  gCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGTCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCCAGGa                                                                                                                                                  >  1:619179/1‑149 (MQ=255)
    tGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGaca                                                                                                                                                >  2:325800/1‑149 (MQ=255)
               ggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTc                                                                                                                                     <  1:1151694/149‑1 (MQ=255)
               ggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTc                                                                                                                                     <  1:36931/149‑1 (MQ=255)
                   aaTCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAATCGAAAAACGGACCTCCGTGGAGGTCGGTTTATATGAATTTGGTGCCCGGACTCGGAATCGACCCAAGGACACGCGGATTTTCAAta                                                                                                                                 >  1:676686/1‑148 (MQ=255)
                   aaTCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATc                                                                                                                                 >  1:118150/1‑149 (MQ=255)
                        ggTTGAGTTCGGTAGAATTGATTTGGAGTCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGAGGAGGTCCGTTTATATGAATTTGGTGCCCGTACTCTGAATAGAACCAATGACACTGGGATTTTCAATCCCCTg                                                                                                                            >  2:778221/1‑149 (MQ=255)
                         gTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGGGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGc                                                                                                                           >  2:803004/1‑149 (MQ=255)
                              gTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTAc                                                                                                                      <  2:1014194/149‑1 (MQ=255)
                                 cGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACACGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGa                                                                                                                   >  2:546526/1‑149 (MQ=255)
                                  ggTAGACTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGAc                                                                                                                  <  1:1021647/149‑1 (MQ=255)
                                  ggTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGGCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGAc                                                                                                                  <  1:1026353/149‑1 (MQ=255)
                                    tAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTg                                                                                                                >  1:56946/1‑149 (MQ=255)
                                     aGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa                                                                                                               <  1:853283/149‑1 (MQ=255)
                                      gAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAg                                                                                                              <  1:138836/149‑1 (MQ=255)
                                      cttcagATCTGGAGGCAGAACGCTTAAATCGTGGAGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGTGATTTTCAATCCCCTGCTCTACCGACTGAg                                                                                                              <  1:97054/144‑1 (MQ=255)
                                         ttGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATAAGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCGTCCGACTGAGCTa                                                                                                           >  1:1177294/1‑149 (MQ=255)
                                             tttGGAGGCATAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCg                                                                                                       >  2:1130446/1‑149 (MQ=255)
                                              ttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCgg                                                                                                      <  2:90908/149‑1 (MQ=255)
                                                    gCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACg                                                                                                <  1:457616/149‑1 (MQ=255)
                                                     cAGAACGCTTAAATCGTGGCGTCCTGAAACTAAAAACGGACCTCCGTGGAGGTCCGTTTATATGTATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACgg                                                                                               <  2:85801/149‑1 (MQ=255)
                                                        aaagCTTAACTCGTGGCGTCCTGAAACGAAAAACGGACCTACGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGgc                                                                                            <  2:191300/146‑1 (MQ=255)
                                                        aaCGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGgc                                                                                            >  1:97449/1‑149 (MQ=255)
                                                                aaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTaaa                                                                                    <  2:633546/149‑1 (MQ=255)
                                                                aaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCGCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTaaa                                                                                    <  1:660579/149‑1 (MQ=255)
                                                                 aTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAGTAAAc                                                                                   <  1:614799/149‑1 (MQ=255)
                                                                   cGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCg                                                                                 <  1:235628/149‑1 (MQ=255)
                                                                         gTCCTGAAACGAAAAACCGACCGCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGTACCAAGGCCACGGGGATTTTAAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCt                                                                           <  1:677389/149‑1 (MQ=21)
                                                                           ccTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGc                                                                         <  1:430443/149‑1 (MQ=255)
                                                                            cTGAAACGAAAAACGGACCTCCGTGGAGCTCCGTTTATAGGAATTAGGTGCCCCGACTCGGAATCGACCGACGTACACGTGGCTTTTCAATCCCCTGCTATACCGACTGACCTATCCGGGCATCGCGGCGCATCAAACCGTAATCAGca                                                                        >  2:81109/1‑149 (MQ=255)
                                                                             tGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACAGAGCTATCCGGGCAACGGGGCGCATTAAGCCGTAATCTGcac                                                                       >  2:223542/1‑149 (MQ=255)
                                                                               aaaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGCCACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACAt                                                                     >  1:597111/1‑149 (MQ=255)
                                                                               aaaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACCCGGGGATGTTCAATCCCCTGCTCTCCCGACTGCGCTATCCGGGCAACGGGGCTCATTAAACCGTAATCTGCACAt                                                                     >  1:670349/1‑149 (MQ=21)
                                                                                aaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGCGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAAGCGTAATCTGCACAtg                                                                    >  2:436956/1‑148 (MQ=255)
                                                                                 aCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCCCAtct                                                                   >  2:129483/1‑149 (MQ=255)
                                                                                 aCGAAAAACGGACCTCCGTGGAGGTCCGGTTATATGAATATGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGAAATCTGCACAtct                                                                   >  1:374450/1‑149 (MQ=255)
                                                                                   gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCg                                                                 <  1:789791/149‑1 (MQ=255)
                                                                                   gAAAAACGGACCTCCGAGGAGGTCAGTTTATATGAATTTGGTGACCGGACTCGGAATCGAACCAATGACACGGGGAATTTCAATCCCCTGCTCTACCGACTGAGCTATCCCGCCAACGGTGCGCATTAAACCGTAATCTGAACATCTCg                                                                 >  1:1031752/1‑149 (MQ=21)
                                                                                    aaaaaCGGACCTCCGTGGAGGTCCGTTTATATGATTTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCAGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAGCCGGAATCTGCAGATGTCGt                                                                >  2:584945/1‑149 (MQ=18)
                                                                                     taaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTc                                                               >  1:404721/2‑149 (MQ=255)
                                                                                     aaaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTc                                                               <  1:803004/149‑1 (MQ=255)
                                                                                      aaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCa                                                              >  2:1005430/1‑149 (MQ=255)
                                                                                        aCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCacc                                                            >  1:199611/1‑147 (MQ=255)
                                                                                          ggACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCt                                                          >  1:969188/1‑149 (MQ=255)
                                                                                             ccTCCGTGGAGGTCCGGTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCACCTGCTCTACCGACTGAGCTATCCGGGCAACGGGTCGCATTAAACCGTAATCTGCACATCTCGGCAACCTAAt                                                       >  2:666937/1‑149 (MQ=21)
                                                                                               tCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAttt                                                     <  1:546526/149‑1 (MQ=255)
                                                                                                 cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGCCTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa                                                   <  2:48334/149‑1 (MQ=255)
                                                                                                 cGTGGAGGTCCGTTTATAGGAATTTGGTGCCCTGACTCGGAATCGAACCAATGACACGGTGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCa                                                   <  1:785646/149‑1 (MQ=255)
                                                                                                     gAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTaa                                               >  2:95868/1‑149 (MQ=255)
                                                                                                     gAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAAACCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCTCATAAAACCGTAATCTGCACATCTCGTCAACCTAATTTCATTaa                                               >  2:1069779/1‑149 (MQ=21)
                                                                                                       ggTCGGTTTATATGAAATTGGTGCCCGGTCTGGGAATCTAACCAATGACACGGGGTTGTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTaaaa                                             <  2:1137817/149‑1 (MQ=21)
                                                                                                          ccGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCg                                          >  2:16424/1‑149 (MQ=255)
                                                                                                           cGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGa                                         >  2:488928/1‑149 (MQ=255)
                                                                                                           cGTTTATATGAATTTGGTGCACGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGa                                         >  1:580456/1‑149 (MQ=255)
                                                                                                            gTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGAt                                        >  1:325452/1‑149 (MQ=255)
                                                                                                             tttATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGAtt                                       <  1:481209/149‑1 (MQ=255)
                                                                                                             tttATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCAATAAACCGTAAACAGCACATCACCTCAACCTAATTTAAGAAAAAGCGagt                                       >  1:766876/1‑147 (MQ=18)
                                                                                                              ttATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCACCCTAATTTCAGTAAAAGCGATTc                                      >  2:404981/1‑149 (MQ=255)
                                                                                                                 tatGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAAc                                   >  2:973821/1‑149 (MQ=255)
                                                                                                                     aaTTTGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCGATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCt                               <  2:943806/149‑1 (MQ=255)
                                                                                                                      aTTTGGTGCCCTGAGTCTGAATCTAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGAATGAGCTATCCGGGCAACGGGGCGCTTTAAACCGTAATCTGCTCATCTCTTCAAGCTAATTTCAGTAAAAGCGATTCAACTGCtt                              <  2:44051/149‑1 (MQ=18)
                                                                                                                       tttGGTGCCCGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTa                             <  1:325800/149‑1 (MQ=255)
                                                                                                                        ttGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTaa                            <  1:45599/149‑1 (MQ=255)
                                                                                                                          ggTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGa                          <  2:179471/149‑1 (MQ=255)
                                                                                                                            tGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGAtt                        >  1:570472/1‑149 (MQ=255)
                                                                                                                             gCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTg                       >  1:442646/1‑149 (MQ=255)
                                                                                                                              cccGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGc                      <  2:969188/149‑1 (MQ=255)
                                                                                                                               ccGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGTCTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCg                     <  1:95868/149‑1 (MQ=255)
                                                                                                                               ccGGACTCGGAATCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCg                     <  1:1069779/149‑1 (MQ=255)
                                                                                                                               ccGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCGGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCAGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCg                     <  2:274988/149‑1 (MQ=255)
                                                                                                                                     tCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAAt               >  1:793493/1‑149 (MQ=255)
                                                                                                                                       ggAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGCGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCACCCTAATTTCAGTAAACGCGATTCAACTGCTTAACATTGCGGCAAATCg             >  2:994770/1‑149 (MQ=255)
                                                                                                                                       ggAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAGCCGTAATCTGCACATCTCGTCGACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCTCATCg             >  1:839781/1‑149 (MQ=255)
                                                                                                                                        gAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCACCCTAATTTCTGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGc            <  1:404981/149‑1 (MQ=255)
                                                                                                                                         aaTCGAACCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGct           <  2:544668/149‑1 (MQ=255)
                                                                                                                                              aaCCAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCttt      >  2:592208/1‑149 (MQ=255)
                                                                                                                                               aCCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCtttt     <  2:619179/149‑1 (MQ=255)
                                                                                                                                                ccAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTGTTTTc    >  2:591775/1‑149 (MQ=255)
                                                                                                                                                 cAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAACGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCt   <  2:1068220/149‑1 (MQ=255)
                                                                                                                                                 cAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCt   <  2:379770/149‑1 (MQ=255)
                                                                                                                                                 cAATGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAAATTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCt   <  2:1072738/149‑1 (MQ=255)
                                                                                                                                                  aaGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCACCCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCtt  >  2:275801/1‑148 (MQ=255)
                                                                                                                                                  aaGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCtg  <  1:485463/149‑1 (MQ=255)
                                                                                                                                                  aaGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCtg  >  1:76976/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTGCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAACGGGGCGCATTAAACCGTAATCTGCACATCTCGTCAACCTAATTTCAGTAAAAGCGATTCAACTGCTTAAGATTGCGGCAAATCGCTTCTTTTCTG  >  NC_000913/4362430‑4362724

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: