Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,990,937 A→G 9.4% intergenic (+45/+17) 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,9370AG 9.4% 238.4 / 11.7 87intergenic (+45/+17)tyrP/yecAtyrosine:H(+) symporter/UPF0149 family protein YecA
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base G (5/3);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATG  >  NC_000913/1990788‑1991085
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                                     <  1:615377/150‑1 (MQ=255)
                 ggCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCctct                                                                                                                                     <  2:282252/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                                                                           >  1:1441555/1‑150 (MQ=255)
                                                                           gCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTAc                                                                           >  2:518019/1‑150 (MQ=255)
                                                                            cGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACg                                                                          >  2:130252/1‑150 (MQ=255)
                                                                                  ggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCgg                                                                    <  1:897680/150‑1 (MQ=255)
                                                                                  ggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCgg                                                                    >  1:343393/1‑150 (MQ=255)
                                                                                      gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                                >  2:1421075/1‑150 (MQ=255)
                                                                                      gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                                <  2:657602/150‑1 (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)
                                                                                                                                                    gACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATg  >  2:383060/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATG  >  NC_000913/1990788‑1991085

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: