Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,197,734 C→T 5.1% R219K (AGA→AAA)  cohE ← e14 prophage; repressor protein phage e14

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,197,7340CT 5.1% 641.6 / 14.8 178R219K (AGA→AAA) cohEe14 prophage; repressor protein phage e14
Reads supporting (aligned to +/- strand):  ref base C (78/90);  new base T (6/3);  total (85/93)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 187 total aligned reads displayed.

TCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTGTA  >  CP009273/1197586‑1197883
                                                                                                                                                    |                                                                                                                                                     
tCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtct                                                                                                                                                      <  1:2016877/150‑1 (MQ=255)
  gtacgtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtcttc                                                                                                                                                    >  2:3454090/1‑150 (MQ=255)
   tacgtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtcttct                                                                                                                                                   <  1:1858837/150‑1 (MQ=255)
    acgtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATcttctt                                                                                                                                                  <  1:80305/150‑1 (MQ=255)
         ccGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGttt                                                                                                                                             <  2:556796/150‑1 (MQ=255)
          cgggcgggCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGtttt                                                                                                                                            <  1:2736649/150‑1 (MQ=255)
           gggcgggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTc                                                                                                                                           >  1:1019767/1‑150 (MQ=255)
           gggcgggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTc                                                                                                                                           >  1:3405858/1‑150 (MQ=255)
                 gcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTgaga                                                                                                                                     <  1:1089476/150‑1 (MQ=255)
                 gcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTgaga                                                                                                                                     <  2:4638817/150‑1 (MQ=255)
                      ttCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATc                                                                                                                                >  2:4671217/1‑150 (MQ=255)
                      ttCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATc                                                                                                                                >  2:2206752/1‑150 (MQ=255)
                       tCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATTTTCTTGGTTTTCCTGAGAAAATca                                                                                                                               >  2:874451/1‑150 (MQ=255)
                        cGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATTTTCTTGGTTTTCCTGAGAAAATcac                                                                                                                              >  1:2507050/1‑150 (MQ=255)
                         gTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATcaca                                                                                                                             >  2:3397685/1‑150 (MQ=255)
                         gTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATcaca                                                                                                                             >  2:1014661/1‑150 (MQ=255)
                         gTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATcaca                                                                                                                             <  1:454602/150‑1 (MQ=255)
                          tGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAg                                                                                                                            <  1:3780800/150‑1 (MQ=255)
                          tGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAg                                                                                                                            >  2:2028992/1‑150 (MQ=255)
                          tGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAg                                                                                                                            >  2:2029013/1‑150 (MQ=255)
                          tGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCAAAg                                                                                                                            >  2:1944917/1‑150 (MQ=255)
                                tCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCaa                                                                                                                      <  2:3927756/150‑1 (MQ=255)
                                    gcctgTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATtat                                                                                                                  >  1:2900092/1‑150 (MQ=255)
                                      ctgTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATag                                                                                                                <  1:3932229/150‑1 (MQ=255)
                                        gtttgttGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATagag                                                                                                              >  1:2267288/1‑150 (MQ=255)
                                          ttgttgTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCa                                                                                                            <  2:178798/150‑1 (MQ=255)
                                           tgttgtTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCaa                                                                                                           <  2:4328147/150‑1 (MQ=255)
                                           tgttgtTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCaa                                                                                                           <  1:4195118/150‑1 (MQ=255)
                                            gttgttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAAt                                                                                                          <  1:3755241/150‑1 (MQ=255)
                                              tgttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTa                                                                                                        <  1:3689905/150‑1 (MQ=255)
                                              tgttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTa                                                                                                        <  1:3784051/150‑1 (MQ=255)
                                               gttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTAc                                                                                                       >  2:4700998/1‑150 (MQ=255)
                                                tttCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTAcc                                                                                                      <  2:3783879/150‑1 (MQ=255)
                                                 ttCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCg                                                                                                     >  1:17416/1‑150 (MQ=255)
                                                  tctcTTGGGTACATTATTTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATATTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGt                                                                                                    >  1:2295615/1‑150 (MQ=255)
                                                    tctTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTg                                                                                                  >  1:2314612/1‑150 (MQ=255)
                                                     ctTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGa                                                                                                 <  1:1269968/150‑1 (MQ=255)
                                                        gggTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCt                                                                                              <  2:3050928/150‑1 (MQ=255)
                                                         ggTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCtt                                                                                             <  2:302968/150‑1 (MQ=255)
                                                          gTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTa                                                                                            <  1:2218145/150‑1 (MQ=255)
                                                           tACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTaa                                                                                           <  2:1749707/150‑1 (MQ=255)
                                                            aCATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAAt                                                                                          <  2:495661/150‑1 (MQ=255)
                                                            aCATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAAt                                                                                          <  2:3542610/150‑1 (MQ=255)
                                                                tATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTaa                                                                                      >  2:1904244/1‑150 (MQ=255)
                                                                 aTGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAg                                                                                     <  2:308994/150‑1 (MQ=255)
                                                                 aTGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAg                                                                                     <  1:179678/150‑1 (MQ=255)
                                                                    tATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATTTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCt                                                                                  >  1:3695489/1‑150 (MQ=255)
                                                                      tctcAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCa                                                                                <  1:3291700/150‑1 (MQ=255)
                                                                        tcAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAgg                                                                              <  1:2519956/150‑1 (MQ=255)
                                                                          aaaGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGcc                                                                            <  2:3037293/150‑1 (MQ=255)
                                                                          aaaGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGcc                                                                            >  2:1552629/1‑150 (MQ=255)
                                                                            aGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAg                                                                          >  2:3136225/1‑150 (MQ=255)
                                                                               tACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTATGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGttt                                                                       <  1:3848033/150‑1 (MQ=255)
                                                                                     gTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGttt                                                                 >  1:150671/1‑150 (MQ=255)
                                                                                        aaGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAAc                                                              <  1:3589864/150‑1 (MQ=255)
                                                                                           tataAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCt                                                           >  2:4189960/1‑150 (MQ=255)
                                                                                             taAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCttt                                                         >  1:1898125/1‑150 (MQ=255)
                                                                                             taAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCttt                                                         <  1:3458132/150‑1 (MQ=255)
                                                                                              aaaaaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTg                                                        >  2:4345164/1‑150 (MQ=255)
                                                                                              aaaaaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTg                                                        <  1:3323216/150‑1 (MQ=255)
                                                                                               aaaaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTga                                                       <  2:3098769/150‑1 (MQ=255)
                                                                                                 aaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTGgaga                                                     >  1:4593739/1‑150 (MQ=255)
                                                                                                   aaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTGTGGGTTTAACGCTTTGAGATa                                                   >  2:1398590/1‑150 (MQ=255)
                                                                                                   aaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATATTGTATCTTCTTGGGTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCCTTGATCATAATGTAAGTCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATa                                                   >  2:2681138/1‑150 (MQ=255)
                                                                                                    aCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATaa                                                  >  2:1905939/1‑150 (MQ=255)
                                                                                                     ccTGCCGAAGCAGGTTCATACACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAAc                                                 <  2:4483394/150‑1 (MQ=255)
                                                                                                         ccGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCtg                                             >  2:3377265/1‑150 (MQ=255)
                                                                                                         ccGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCtg                                             <  1:3348726/150‑1 (MQ=255)
                                                                                                          cGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCtgt                                            >  1:2863911/1‑150 (MQ=255)
                                                                                                          cGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTGGAGATAACGCtgt                                            >  1:1008146/1‑150 (MQ=255)
                                                                                                            aaGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTgtgt                                          <  1:1618287/150‑1 (MQ=255)
                                                                                                            aaGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTgtgt                                          >  1:3556885/1‑150 (MQ=255)
                                                                                                             aGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTc                                         >  2:3611185/1‑150 (MQ=255)
                                                                                                               cAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATTTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTccc                                       >  2:2456784/1‑150 (MQ=255)
                                                                                                                aGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCa                                      >  1:1802792/1‑150 (MQ=255)
                                                                                                                 ggTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAt                                     >  1:3901456/1‑150 (MQ=255)
                                                                                                                      aTAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAtcttct                                <  2:3263504/150‑1 (MQ=255)
                                                                                                                      aTAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAtcttct                                <  1:2686744/150‑1 (MQ=255)
                                                                                                                      aTAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAtcttct                                <  1:2685441/150‑1 (MQ=255)
                                                                                                                         aaCATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGATTTGAGATAACGCTGTGTCCCATCTTCTATc                             >  1:4730684/1‑150 (MQ=255)
                                                                                                                           cATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCACTTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCaa                           >  2:2983300/1‑150 (MQ=255)
                                                                                                                           cATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTAACAATTATAGAGAAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGATTGGGTTTAACGATTTGAGATAACGCTGTGTCCCATCTTCTATCaa                           >  2:708005/1‑150 (MQ=255)
                                                                                                                                  tAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTg                    >  2:4628533/1‑150 (MQ=255)
                                                                                                                                  tAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTg                    >  2:345173/1‑150 (MQ=255)
                                                                                                                                   aGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGa                   >  1:279983/1‑150 (MQ=255)
                                                                                                                                   aGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGa                   >  1:1770016/1‑150 (MQ=255)
                                                                                                                                     gCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAg                 <  2:331611/150‑1 (MQ=255)
                                                                                                                                     gCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAg                 <  2:331892/150‑1 (MQ=255)
                                                                                                                                     gCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTGTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAg                 >  2:1137118/1‑150 (MQ=255)
                                                                                                                                      cTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAgg                >  1:1854540/1‑150 (MQ=255)
                                                                                                                                       tttGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAggt               <  1:2374551/150‑1 (MQ=255)
                                                                                                                                        ttGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAggtg              <  2:2835379/150‑1 (MQ=255)
                                                                                                                                         tGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAggtgg             >  2:2766945/1‑150 (MQ=255)
                                                                                                                                              ttGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCa        <  2:4436945/150‑1 (MQ=255)
                                                                                                                                               tGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATGATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCAc       <  1:2461113/150‑1 (MQ=255)
                                                                                                                                                 tATCGTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCt     <  2:1536090/150‑1 (MQ=255)
                                                                                                                                                  aTCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTg    <  2:403998/150‑1 (MQ=255)
                                                                                                                                                   tcttctTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTGt   >  1:2472644/1‑150 (MQ=255)
                                                                                                                                                   tcttctTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTGt   >  1:3170236/1‑150 (MQ=255)
                                                                                                                                                    cttcttGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTGTa  <  2:1039144/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTGTA  >  CP009273/1197586‑1197883

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