Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,990,933 G→T 8.9% intergenic (+41/+21) tyrP → / ← yecA tyrosine:H(+) symporter/UPF0149 family protein YecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,990,9330GT 8.9% 289.3 / 13.8 90intergenic (+41/+21)tyrP/yecAtyrosine:H(+) symporter/UPF0149 family protein YecA
Reads supporting (aligned to +/- strand):  ref base G (41/41);  new base T (5/3);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

GGGTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGC  >  NC_000913/1990786‑1991082
                                                                                                                                                   |                                                                                                                                                     
gggTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCt                                                                                                                                                     >  1:489690/1‑150 (MQ=255)
 ggTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTg                                                                                                                                                    <  2:599246/150‑1 (MQ=255)
 ggTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTg                                                                                                                                                    >  1:1547601/1‑150 (MQ=255)
  gTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGa                                                                                                                                                   >  1:685369/1‑150 (MQ=255)
   tCAAAGGTGGTCGTCCGGCGCTGGTCGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGac                                                                                                                                                  >  2:527968/1‑150 (MQ=255)
    cAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGaca                                                                                                                                                 >  2:1167305/1‑150 (MQ=255)
          tggtCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAt                                                                                                                                           <  1:722107/150‑1 (MQ=255)
          tggtCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAt                                                                                                                                           <  1:508486/150‑1 (MQ=255)
          tggtCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAAt                                                                                                                                           <  1:816978/150‑1 (MQ=255)
           ggtCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAAtt                                                                                                                                          >  1:1635301/1‑150 (MQ=255)
                gccGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGcc                                                                                                                                     <  1:846963/149‑1 (MQ=255)
                 ccGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAtaaggaat                                                                                                                                    <  2:1614506/150‑8 (MQ=38)
                   ggCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCctct                                                                                                                                  <  2:282252/150‑1 (MQ=255)
                   ggCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCctct                                                                                                                                  <  1:615377/150‑1 (MQ=255)
                      gcTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTaa                                                                                                                               >  1:657602/1‑150 (MQ=255)
                        gggtggtggtGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATg                                                                                                                             <  1:1581220/149‑1 (MQ=255)
                         ggtggtggtGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgc                                                                                                                            <  1:1167305/150‑1 (MQ=255)
                          gtggtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgca                                                                                                                           <  2:685369/150‑1 (MQ=255)
                           tggtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGTAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgcag                                                                                                                          >  2:1457049/1‑150 (MQ=255)
                            ggtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgcagg                                                                                                                         >  1:1108865/1‑150 (MQ=255)
                             gtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgcaggc                                                                                                                        <  2:446764/150‑1 (MQ=255)
                             ctcgagattCTCTTTGGTATAGCTGTGATTGGCGTGCAATTTTTGATTGCCGCAGGGTTGTTACCAAAAGTCGTGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgcaggc                                                                                                                        >  2:233621/8‑150 (MQ=255)
                               ggtgTTTCTCTGTGTTATTACTGTCATTGACCTGCAAATTTTGATTACGGCAGGATTGTTACAAGAACTCGGGTGATCAGATAGCCTCAAATTCCTTAATGGGTGACAGAATTAACGCTGACACACAATTTGGCCTCTTAATgcaggcaa                                                                                                                      >  2:716090/1‑149 (MQ=255)
                               ggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGTAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgcaggcag                                                                                                                      <  1:1457049/150‑1 (MQ=255)
                               ggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgcaggcag                                                                                                                      <  2:489690/150‑1 (MQ=255)
                                  gtTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCAc                                                                                                                   >  2:897680/1‑150 (MQ=255)
                                      ctctGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCt                                                                                                               <  2:622207/150‑1 (MQ=255)
                                         tgtgGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTaa                                                                                                            >  1:1239003/1‑150 (MQ=255)
                                          gtgGTATTGCTGTGATTGTCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTTCCATAATTAACGCTGACACCCAATTTGGCCTATTAATTCAGGCAGCACTTCTTaaa                                                                                                           <  2:810430/150‑1 (MQ=255)
                                            ggTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAAtt                                                                                                         <  1:1492795/150‑1 (MQ=255)
                                            ggTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAAtt                                                                                                         >  1:1421736/1‑150 (MQ=255)
                                               aTTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCt                                                                                                      >  1:1470520/1‑150 (MQ=255)
                                                 tGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCtt                                                                                                     >  1:212415/1‑149 (MQ=255)
                                                 tGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCtt                                                                                                     <  2:212415/149‑1 (MQ=255)
                                                 tGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTa                                                                                                    >  2:981772/1‑150 (MQ=255)
                                                   cTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTAcc                                                                                                  <  2:1404343/150‑1 (MQ=255)
                                                    tgtgATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCa                                                                                                 <  1:702901/150‑1 (MQ=255)
                                                      tgATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACt                                                                                               >  1:71497/1‑150 (MQ=255)
                                                       taTTGGCTTGCAATTTTTGATTGCGGCAGGGTTGTTACGAGAAGTGGGGTGATCAGATAGCGTCAAATTCCTTATTGGGTGCCAGAATTAACGGTGACAGCCAATTTGGACTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTa                                                                                              <  2:495993/149‑1 (MQ=255)
                                                       gATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTa                                                                                              >  2:1061331/1‑150 (MQ=255)
                                                            gCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCa                                                                                         >  2:1270356/1‑150 (MQ=255)
                                                              gTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACg                                                                                       >  1:1186406/1‑150 (MQ=255)
                                                                  aaTTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCa                                                                                   >  2:611500/1‑150 (MQ=255)
                                                                    tttttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAg                                                                                 <  2:79347/150‑1 (MQ=255)
                                                                    tttttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAg                                                                                 <  1:187509/150‑1 (MQ=255)
                                                                       ttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGAt                                                                              <  1:735270/150‑1 (MQ=255)
                                                                           ttGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACTACTACCGCACGGGCAAGGATCGtt                                                                          <  2:512992/150‑1 (MQ=255)
                                                                            tGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTa                                                                         >  2:1353649/1‑150 (MQ=255)
                                                                             gCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTAc                                                                        >  2:518019/1‑150 (MQ=255)
                                                                             gCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTAc                                                                        >  1:1441555/1‑150 (MQ=255)
                                                                              cGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACg                                                                       >  2:130252/1‑150 (MQ=255)
                                                                                    ggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCgg                                                                 >  1:343393/1‑150 (MQ=255)
                                                                                    ggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCgg                                                                 <  1:897680/150‑1 (MQ=255)
                                                                                        gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                             <  2:657602/150‑1 (MQ=255)
                                                                                        gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                             >  2:1421075/1‑150 (MQ=255)
                                                                                        gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                             >  1:1555544/1‑150 (MQ=255)
                                                                                          tACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGGCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtctc                                                           >  1:1115847/1‑150 (MQ=255)
                                                                                          tACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtctc                                                           >  1:286038/1‑150 (MQ=255)
                                                                                                 aGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGATTTAACGCTGACACCCAATTTGGCCTCTTAATGCATGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCt                                                    <  2:1342692/150‑1 (MQ=255)
                                                                                                     gggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGa                                                <  2:1635301/150‑1 (MQ=255)
                                                                                                            cAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTg                                         <  2:343393/150‑1 (MQ=255)
                                                                                                            cAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTg                                         >  2:742673/1‑150 (MQ=255)
                                                                                                            cAGATAGCATCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACAAAATTTGGACTCTTAATGCAGGCAGCACACCTTAAATTTCTTACCACTACAGCACGAGCAAGGAGCGTTACGCCCCCGTTTCTCTTATTCTTTGATCGGttt                                         >  2:1066690/1‑149 (MQ=255)
                                                                                                             aGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGc                                        <  1:233621/150‑1 (MQ=255)
                                                                                                               aTAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg                                      <  2:1441555/150‑1 (MQ=255)
                                                                                                                tAGCCTCAAATTCCTTATTGTGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGa                                     <  1:1298833/150‑1 (MQ=255)
                                                                                                                tAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGa                                     <  1:122856/150‑1 (MQ=255)
                                                                                                                 aGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGaa                                    >  2:1298131/1‑150 (MQ=255)
                                                                                                                  gCCTCAAATTCCTTATTGGGTGCCAGAATTACCGCTGACCCCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTCCCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCGGCTTTGATCGGTTGCTGacc                                   >  1:761400/1‑148 (MQ=255)
                                                                                                                  gCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTATGCTTTGATCGGTTGCTGAAc                                   >  2:1632084/1‑150 (MQ=255)
                                                                                                                    cTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg                                 >  2:1142454/1‑150 (MQ=255)
                                                                                                                      cAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCt                               >  2:296496/1‑150 (MQ=255)
                                                                                                                          ttCCTTATTGGGGGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCttttt                           <  2:195378/150‑1 (MQ=255)
                                                                                                                            ccaaaTTGGGTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc                         >  1:1081564/5‑150 (MQ=39)
                                                                                                                              ttATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTg                       <  2:71497/150‑1 (MQ=255)
                                                                                                                                aTTGGGTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCg                     >  2:1301313/1‑150 (MQ=255)
                                                                                                                                  tGGGTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTAAACAGCTTTTTCCTGCGGa                   >  1:142333/1‑150 (MQ=255)
                                                                                                                                   gggTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGAt                  >  2:1149305/1‑150 (MQ=255)
                                                                                                                                    ggTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATAGGTTGCTGAACAGAGGTTTCATGCGGATg                 >  1:733627/1‑150 (MQ=255)
                                                                                                                                    ggTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATg                 >  2:702391/1‑150 (MQ=255)
                                                                                                                                      tGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgcg               <  2:831438/150‑1 (MQ=255)
                                                                                                                                        ccAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGcc             <  2:286038/150‑1 (MQ=255)
                                                                                                                                        ccAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGcc             <  2:631103/150‑1 (MQ=255)
                                                                                                                                          aGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCAt           >  1:128555/1‑150 (MQ=255)
                                                                                                                                           gAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATc          <  1:1353649/150‑1 (MQ=255)
                                                                                                                                           gAATTAACGATGACAACCAATTTGGCCTCTTAATGCAGCCAGCACTGCTTAAATTTCTTACCAATAAAGCACGGCCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGCTTGCTGAACACCTTTTTCCTGCGTATCCGCCATc          >  2:465202/1‑150 (MQ=255)
                                                                                                                                            aaTTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATcc         <  1:1061331/150‑1 (MQ=255)
                                                                                                                                               tAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTTTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGt      <  2:1186406/150‑1 (MQ=255)
                                                                                                                                                aaCGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTa     <  1:518019/150‑1 (MQ=255)
                                                                                                                                                 aCGCTGTCACCGAATTTGGCCTCTTAATGCTGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAATGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAAGAGCTTTTTCCTGCGGATGCGCCATCCAGTAc    <  1:716090/150‑1 (MQ=255)
                                                                                                                                                   gCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGc  <  2:1163296/150‑1 (MQ=255)
                                                                                                                                                   gCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTATTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTACTGCGGATGCGCCATCCAGTACGc  >  1:18205/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGGTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGC  >  NC_000913/1990786‑1991082

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