Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,192,525 C→T 9.9% S155N (AGC→AAC)  ymfD ← e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,192,5250CT 9.9% 848.4 / 59.7 263S155N (AGC→AAC) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (129/108);  new base T (13/13);  total (142/121)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 268 total aligned reads displayed.

CGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGT  >  CP009273/1192377‑1192672
                                                                                                                                                    |                                                                                                                                                   
cGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGc                                                                                                                                                     <  1:4385414/149‑1 (MQ=255)
cGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGc                                                                                                                                                     <  2:2479593/149‑1 (MQ=255)
cGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGc                                                                                                                                                     <  2:3871416/149‑1 (MQ=255)
  aCTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCtt                                                                                                                                                   <  1:4676656/149‑1 (MQ=255)
    tttAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCtttt                                                                                                                                                 <  1:300701/149‑1 (MQ=255)
     ttAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTta                                                                                                                                                <  1:3836771/149‑1 (MQ=255)
       aaTATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTtata                                                                                                                                              >  2:3493538/1‑149 (MQ=255)
        atatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTatat                                                                                                                                             >  2:3128263/1‑149 (MQ=255)
        atatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTatat                                                                                                                                             >  2:2661857/1‑149 (MQ=255)
         tatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATAtt                                                                                                                                            >  1:17502/1‑149 (MQ=255)
           tttCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATAttta                                                                                                                                          >  2:2148725/1‑147 (MQ=255)
            ttCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGAt                                                                                                                                         >  2:3298683/1‑149 (MQ=255)
              cAAGGCGATGTGTAGAACATATTTCTTGCTCTGCGAGTTCATCGAGTTATCCATATTTAGTTAGACTTTTTGAATTTTTGTATATGTTTCCGTAAAGATGTTTTCTTCCTGTTTCGTTTTTTTTTAAGTATGAGCTTTTATATTGtttt                                                                                                                                       <  1:4349746/149‑1 (MQ=255)
                aGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTAc                                                                                                                                     >  1:2920362/1‑149 (MQ=255)
                 ggCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACa                                                                                                                                    >  2:4300843/1‑149 (MQ=255)
                 ggCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACa                                                                                                                                    >  2:749577/1‑149 (MQ=255)
                   ccATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAAt                                                                                                                                  <  1:1306293/149‑1 (MQ=255)
                     aTGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATaa                                                                                                                                <  1:3634152/149‑1 (MQ=255)
                      tgtgAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAag                                                                                                                               <  2:4091804/149‑1 (MQ=255)
                       gtgAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAaga                                                                                                                              >  1:4359429/1‑149 (MQ=255)
                        tgAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAagag                                                                                                                             <  1:2407947/149‑1 (MQ=255)
                                  atTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAcctcc                                                                                                                   <  2:3008663/149‑1 (MQ=255)
                                     tCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGa                                                                                                                <  2:3480031/149‑1 (MQ=255)
                                         gCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAtttt                                                                                                            <  1:1899335/149‑1 (MQ=255)
                                            cTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGt                                                                                                         >  1:3710783/1‑149 (MQ=255)
                                              gCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTaa                                                                                                       >  1:3081151/1‑149 (MQ=255)
                                                 aGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAtct                                                                                                    <  2:2415522/149‑1 (MQ=255)
                                                       tCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGAt                                                                                              >  2:3401113/1‑149 (MQ=255)
                                                       tCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGAt                                                                                              <  1:4540962/149‑1 (MQ=255)
                                                         gAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGc                                                                                            >  1:3661641/1‑149 (MQ=255)
                                                             aaTCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAg                                                                                        >  1:2126298/1‑149 (MQ=255)
                                                              aTCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGa                                                                                       >  1:2236605/1‑149 (MQ=255)
                                                               tCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGaa                                                                                      >  2:572197/1‑149 (MQ=255)
                                                                 cATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACt                                                                                    <  1:841244/149‑1 (MQ=255)
                                                                  aTAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTa                                                                                   <  2:2574344/149‑1 (MQ=255)
                                                                  aTAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTa                                                                                   >  2:2238325/1‑149 (MQ=255)
                                                                  aTAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTa                                                                                   <  2:2277055/149‑1 (MQ=255)
                                                                    aGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAtt                                                                                 >  2:4402843/1‑149 (MQ=255)
                                                                    aGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAtt                                                                                 >  1:2572003/1‑149 (MQ=255)
                                                                    aGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAtt                                                                                 >  2:3341087/1‑149 (MQ=255)
                                                                        aGAAACAGTTTTTGTATTTTTTTTAATGTATCCGTAAAGATGTTTTCTTCCTGTTTAGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTATGAACTATTTTGt                                                                             <  2:219182/149‑1 (MQ=255)
                                                                        aGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGt                                                                             <  2:148131/149‑1 (MQ=255)
                                                                         gAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTc                                                                            >  2:1788180/1‑149 (MQ=255)
                                                                         gAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTc                                                                            >  2:1788218/1‑149 (MQ=255)
                                                                         gAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTc                                                                            >  2:2544448/1‑149 (MQ=255)
                                                                          aaaCACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCg                                                                           >  1:1535467/1‑149 (MQ=255)
                                                                            acacTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGAt                                                                         >  1:3378443/1‑149 (MQ=255)
                                                                             cacTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGAtt                                                                        <  2:3768497/149‑1 (MQ=255)
                                                                              acTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTg                                                                       <  1:2682478/149‑1 (MQ=255)
                                                                               cTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGt                                                                      <  1:4219433/149‑1 (MQ=255)
                                                                                tttttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTa                                                                     <  2:69767/149‑1 (MQ=255)
                                                                                tttttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTa                                                                     <  2:1116885/149‑1 (MQ=255)
                                                                                tttttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTa                                                                     >  2:1251407/1‑149 (MQ=255)
                                                                                tttttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTa                                                                     >  2:420812/1‑149 (MQ=255)
                                                                                  tttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATc                                                                   <  1:4300843/149‑1 (MQ=255)
                                                                                   ttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATcc                                                                  >  2:130010/1‑149 (MQ=255)
                                                                                       aTTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGAc                                                              >  1:3167852/1‑149 (MQ=255)
                                                                                         ttttGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACaa                                                            <  1:849528/149‑1 (MQ=255)
                                                                                         ttttGTAAATGTATCCGTAAAGATGATTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACaa                                                            <  1:4485871/149‑1 (MQ=255)
                                                                                            tGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATGGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGa                                                         <  2:3977650/149‑1 (MQ=255)
                                                                                                aaTGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACgg                                                     >  2:2207813/1‑149 (MQ=255)
                                                                                                aaTGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACgg                                                     >  1:630941/1‑149 (MQ=255)
                                                                                                 aTGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACACGGAACGGc                                                    >  1:3602647/1‑149 (MQ=255)
                                                                                                    tATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCaga                                                 <  1:4182193/149‑1 (MQ=255)
                                                                                                    tATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCaga                                                 >  1:1424594/1‑149 (MQ=255)
                                                                                                      tCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAATATTCTCACCTCCTGATTTTAGTAATCTCTTTATTCTAAGAAATATTTTGTCGATTGTATCCCGACAAGGAAAGGCacaga                                               >  1:2866461/1‑149 (MQ=255)
                                                                                                        cGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAgagaga                                             <  1:983011/149‑1 (MQ=255)
                                                                                                        cGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAgagaga                                             <  1:705665/149‑1 (MQ=255)
                                                                                                          tAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAAc                                           <  1:2371750/149‑1 (MQ=255)
                                                                                                           aaaGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACa                                          <  2:3033614/149‑1 (MQ=255)
                                                                                                            aaGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACAt                                         <  2:2791729/149‑1 (MQ=255)
                                                                                                             aGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACAtt                                        >  1:3505727/1‑149 (MQ=255)
                                                                                                             aGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTAATGATTTTATTAATATCTTGATGCTAAGAACTATTTTTTCGATTGTATCACGACAATGAAAGGCAGAGAGAACAtt                                        >  1:1737522/1‑149 (MQ=255)
                                                                                                              gATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTa                                       >  1:1098148/1‑149 (MQ=255)
                                                                                                              gATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTa                                       >  2:531047/1‑149 (MQ=255)
                                                                                                               aTGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTag                                      >  2:1888822/1‑149 (MQ=255)
                                                                                                                tGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTaga                                     >  2:3848645/1‑149 (MQ=255)
                                                                                                                     tcttcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCaaa                                >  1:4447870/1‑149 (MQ=255)
                                                                                                                      cttcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAg                               >  1:546774/1‑149 (MQ=255)
                                                                                                                       ttcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAgg                              >  1:1593229/1‑149 (MQ=255)
                                                                                                                        tcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGa                             <  2:4569775/149‑1 (MQ=255)
                                                                                                                          cTGTTTCGTATTTTTTGAAGTTTGAGTTTTTATATTGATGTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATg                           <  2:2355058/149‑1 (MQ=255)
                                                                                                                            gTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGaa                         <  1:501083/149‑1 (MQ=255)
                                                                                                                              ttCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAAt                       <  2:850635/149‑1 (MQ=255)
                                                                                                                                cGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCg                     <  1:1159956/149‑1 (MQ=255)
                                                                                                                                    ttttttGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAc                 >  2:1553370/1‑149 (MQ=255)
                                                                                                                                     tttttGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAcc                <  2:2700870/149‑1 (MQ=255)
                                                                                                                                        ttGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAccccc             <  1:4660229/149‑1 (MQ=255)
                                                                                                                                          gAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCaa           >  2:2915254/1‑149 (MQ=255)
                                                                                                                                            aGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAatt         <  2:4377607/149‑1 (MQ=255)
                                                                                                                                            aGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAatt         <  1:2200276/149‑1 (MQ=255)
                                                                                                                                             gTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAatta        >  2:161041/1‑149 (MQ=255)
                                                                                                                                             gTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAatta        >  2:4638271/1‑149 (MQ=255)
                                                                                                                                              tATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAattat       >  2:2500769/1‑149 (MQ=255)
                                                                                                                                              tATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAattat       >  2:886817/1‑149 (MQ=255)
                                                                                                                                                  aGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTg   <  1:3128263/149‑1 (MQ=255)
                                                                                                                                                  aGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTg   <  1:2661857/149‑1 (MQ=255)
                                                                                                                                                   gCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGt  >  2:2088187/1‑149 (MQ=255)
                                                                                                                                                   gCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGt  >  1:3975013/1‑149 (MQ=255)
                                                                                                                                                   gCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTATTGATGCTAAGAACTATTTTGTCTATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGt  <  2:1737522/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGT  >  CP009273/1192377‑1192672

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: