Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,192,525 C→T 11.5% 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,5250CT11.5% 986.0 / 89.1 322S155N (AGC→AAC) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (169/113);  new base T (13/24);  total (183/139)
Fisher's exact test for biased strand distribution p-value = 4.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 322 total aligned reads displayed.

GACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTC  >  CP009273/1192378‑1192673
                                                                                                                                                   |                                                                                                                                                    
  cTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTtttt                                                                                                                                                  <  1:4062351/150‑1 (MQ=255)
  attttATATTTCATGGCCTTGTGAAGAAGATATTTCTTGCCCTGCGTGTTCATCGAGTATTCCATAGTTAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCtttt                                                                                                                                                  <  1:4592609/149‑1 (MQ=255)
    ttAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTtat                                                                                                                                                <  1:707665/150‑1 (MQ=255)
     tAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTtata                                                                                                                                               >  1:1905604/1‑150 (MQ=255)
      aaTATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTatat                                                                                                                                              <  2:39016/150‑1 (MQ=255)
      aaTATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCATGTTTCGTATTTTTTGAATTATGAGCTTTTatat                                                                                                                                              >  2:2555921/1‑150 (MQ=255)
       atatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATAtt                                                                                                                                             >  1:2718968/1‑150 (MQ=255)
        tatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTg                                                                                                                                            <  1:5200220/150‑1 (MQ=255)
           ttCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGAtt                                                                                                                                         <  1:2040376/150‑1 (MQ=255)
               aGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACa                                                                                                                                     <  2:1376910/150‑1 (MQ=255)
                 gCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAAt                                                                                                                                   >  2:4896447/1‑150 (MQ=255)
                   cATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATaa                                                                                                                                 >  1:5645194/1‑150 (MQ=255)
                   cATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATaa                                                                                                                                 >  1:1461762/1‑150 (MQ=255)
                       tgAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTATTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGt                                                                                                                             <  2:2344287/150‑1 (MQ=255)
                          agaaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGtctc                                                                                                                          <  1:4943086/150‑1 (MQ=255)
                          agaaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGtctc                                                                                                                          <  1:5073801/150‑1 (MQ=255)
                          agaaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGtctc                                                                                                                          >  1:4842559/1‑150 (MQ=255)
                          agaaCATATTTCTTGCACTGCGAGTTAATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGtctc                                                                                                                          <  2:3440511/150‑1 (MQ=255)
                           gaaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCa                                                                                                                         >  2:1752196/1‑150 (MQ=255)
                            aaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAc                                                                                                                        >  2:756073/1‑150 (MQ=255)
                                 atTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAcctcct                                                                                                                   >  2:312211/1‑150 (MQ=255)
                                  tttCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTg                                                                                                                  >  1:2647794/1‑150 (MQ=255)
                                     cTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAtt                                                                                                               <  1:1156255/150‑1 (MQ=255)
                                     cTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAtt                                                                                                               >  1:5593964/1‑150 (MQ=255)
                                      ttGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAttt                                                                                                              >  2:1577886/1‑150 (MQ=255)
                                      ttGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATATGAGTCTCACCTCCTGAttt                                                                                                              >  2:4263017/1‑150 (MQ=255)
                                        gCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTa                                                                                                            <  1:2078767/150‑1 (MQ=255)
                                         cACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAg                                                                                                           >  2:4625479/1‑150 (MQ=255)
                                         cACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAg                                                                                                           >  1:4710619/1‑150 (MQ=255)
                                          aCTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGt                                                                                                          >  2:4740524/1‑150 (MQ=255)
                                          aCTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGt                                                                                                          <  2:2215752/150‑1 (MQ=255)
                                           cTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTa                                                                                                         <  1:4435880/150‑1 (MQ=255)
                                             gCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAt                                                                                                       <  1:3801432/150‑1 (MQ=255)
                                               gAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAtct                                                                                                     >  1:268736/1‑150 (MQ=255)
                                                    cATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGa                                                                                                >  2:920161/1‑150 (MQ=255)
                                                      tCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATg                                                                                              >  1:1100442/1‑150 (MQ=255)
                                                           tAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAg                                                                                         <  2:3270529/150‑1 (MQ=255)
                                                               cgaTAGTAAGAAACACTTTTTGAATATTTGTATATGTATCCTTAAAGATGTTATCTTACTGGGTCGTATATATTGTAGTATGAGCTTTTATATTGATTTTCAATAAGAGACTCACCTCCTGATTTTAGTACTCTCGTGTTGCTAAGAACt                                                                                     >  2:2978588/3‑150 (MQ=255)
                                                               ccATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACt                                                                                     <  2:5211446/150‑1 (MQ=255)
                                                               ccATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACt                                                                                     >  1:1733310/1‑150 (MQ=255)
                                                                 aTAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                                                   >  1:686104/1‑150 (MQ=255)
                                                                    gTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAtttt                                                                                <  1:2349362/150‑1 (MQ=255)
                                                                    gTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAtttt                                                                                >  2:3287333/1‑150 (MQ=255)
                                                                    gTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAtttt                                                                                >  2:1083320/1‑150 (MQ=255)
                                                                     tAAGAAACACTTTTTTAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTTTGAGCTTTTATATTGATTTACTATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTg                                                                               <  1:328844/150‑1 (MQ=255)
                                                                     tAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTg                                                                               >  1:3635170/1‑150 (MQ=255)
                                                                      aaGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGt                                                                              <  2:189735/150‑1 (MQ=255)
                                                                        gAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCg                                                                            <  2:4331751/150‑1 (MQ=255)
                                                                        gAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCg                                                                            >  2:2425217/1‑150 (MQ=255)
                                                                        gAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCg                                                                            >  1:2228623/1‑150 (MQ=255)
                                                                         aaaCACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGa                                                                           >  2:3586089/1‑150 (MQ=255)
                                                                          aaCACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATGTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGAt                                                                          >  1:3969784/1‑150 (MQ=255)
                                                                           acacTTTTTGAATTTTTTTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGAtt                                                                         <  2:550883/150‑1 (MQ=255)
                                                                            cacTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTACTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTg                                                                        >  1:825132/1‑150 (MQ=255)
                                                                             acTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGt                                                                       >  2:3144946/1‑150 (MQ=255)
                                                                               tttttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTAt                                                                     >  1:2852799/1‑150 (MQ=255)
                                                                               tttttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTAt                                                                     <  2:2028523/150‑1 (MQ=255)
                                                                               tttttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTAt                                                                     >  1:4577325/1‑150 (MQ=255)
                                                                                      aTTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACa                                                              >  1:422031/1‑150 (MQ=255)
                                                                                      aTTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACa                                                              <  1:3297434/150‑1 (MQ=255)
                                                                                           tGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGaa                                                         <  2:3682143/150‑1 (MQ=255)
                                                                                            gTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAAc                                                        >  2:1645223/1‑150 (MQ=255)
                                                                                             tAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACg                                                       >  1:3144928/1‑150 (MQ=255)
                                                                                             tAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACg                                                       >  1:4090804/1‑150 (MQ=255)
                                                                                              aaaTGTATCCGTAAAGATGTTTTCTTCATGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAATTTTTCGATTGTATCCCGACAAGTAACgg                                                      >  2:4465604/1‑150 (MQ=255)
                                                                                               aaTGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGc                                                     >  1:5031542/1‑150 (MQ=255)
                                                                                                 tGTATCCGTTAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCag                                                   <  2:2899431/150‑1 (MQ=255)
                                                                                                 tGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCAGACAAGGAACGGCag                                                   <  1:4963690/150‑1 (MQ=255)
                                                                                                    aTCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCagaga                                                <  2:4224395/150‑1 (MQ=255)
                                                                                                    aTCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCagaga                                                >  2:652900/1‑150 (MQ=255)
                                                                                                          aaaGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACAt                                          >  1:3002214/1‑150 (MQ=255)
                                                                                                          aaaGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACAt                                          <  1:1935235/150‑1 (MQ=255)
                                                                                                            aGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTa                                        >  2:3202431/1‑150 (MQ=255)
                                                                                                            aGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTa                                        >  1:967924/1‑150 (MQ=255)
                                                                                                             gATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTag                                       >  1:5326377/1‑150 (MQ=255)
                                                                                                                gTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGAATTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGc                                    >  2:4574310/1‑150 (MQ=255)
                                                                                                                 ttttCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAAAATTAGAGCa                                   >  2:4408452/1‑150 (MQ=255)
                                                                                                                   ttcttcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCaaa                                 >  1:5487922/1‑150 (MQ=255)
                                                                                                                     cttcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAgg                               >  2:2186925/1‑150 (MQ=255)
                                                                                                                     cttcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAgg                               >  2:5013350/1‑150 (MQ=255)
                                                                                                                      ttcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGa                              >  1:3542752/1‑150 (MQ=255)
                                                                                                                        ccTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATg                            <  1:1566332/150‑1 (MQ=255)
                                                                                                                         cTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGa                           >  1:1095202/1‑150 (MQ=255)
                                                                                                                          tGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGaa                          <  1:4639133/150‑1 (MQ=255)
                                                                                                                            tttCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAAt                        <  2:5362244/150‑1 (MQ=255)
                                                                                                                            tttCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAAt                        >  1:930275/1‑150 (MQ=255)
                                                                                                                             ttCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATc                       <  1:1283754/150‑1 (MQ=255)
                                                                                                                             ttCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATc                       >  1:215224/1‑150 (MQ=255)
                                                                                                                              tCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCg                      <  1:2453465/150‑1 (MQ=255)
                                                                                                                               cGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGt                     >  1:3115834/1‑150 (MQ=255)
                                                                                                                                gTATTTTTTGAAGTATGAGCTTTTATATTGATTGACAATAAGAGTCTCTCCTCCTGATTTTAGTACTCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTa                    <  1:1028574/150‑1 (MQ=255)
                                                                                                                                  aTTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAc                  <  2:3768130/150‑1 (MQ=255)
                                                                                                                                   ttttttGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAcc                 <  2:5645452/150‑1 (MQ=255)
                                                                                                                                     ttttGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAcccc               <  2:5308778/150‑1 (MQ=255)
                                                                                                                                          aaGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAatt          >  2:1794862/1‑150 (MQ=255)
                                                                                                                                              aTGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTAttt      <  2:2718968/150‑1 (MQ=255)
                                                                                                                                               tGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTAtttt     >  2:2597898/1‑150 (MQ=255)
                                                                                                                                               tGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTAattt     >  2:2639202/1‑150 (MQ=255)
                                                                                                                                                gAGCTTTTATATTTATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTg    <  2:4072818/150‑1 (MQ=255)
                                                                                                                                                 aGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGt   >  2:126018/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTC  >  CP009273/1192378‑1192673

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