Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,192,525 C→T 19.8% 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,5250CT19.8% 198.2 / 42.9 86S155N (AGC→AAC) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/44);  new base T (12/5);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 1.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTC  >  CP009273/1192377‑1192673
                                                                                                                                                    |                                                                                                                                                    
cGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCt                                                                                                                                                     <  2:2948698/150‑1 (MQ=255)
  aCTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCttt                                                                                                                                                   >  2:3405316/1‑150 (MQ=255)
     ttAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTtat                                                                                                                                                <  2:3793082/150‑1 (MQ=255)
       aaTATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTatat                                                                                                                                              >  2:1799510/1‑150 (MQ=255)
       aaTATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTatat                                                                                                                                              <  1:1990390/150‑1 (MQ=255)
        atatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATAtt                                                                                                                                             >  1:2790994/1‑150 (MQ=255)
         tatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTg                                                                                                                                            <  1:3696879/150‑1 (MQ=255)
         tatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTg                                                                                                                                            <  1:3129304/150‑1 (MQ=255)
         tatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTg                                                                                                                                            <  1:3129083/150‑1 (MQ=255)
         tatTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGCTTCGTATTTTTTGAAGTATGAGCTTTTATATTg                                                                                                                                            <  1:503058/150‑1 (MQ=255)
                aGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTaaa                                                                                                                                     >  2:1791508/1‑148 (MQ=255)
                aGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACa                                                                                                                                     >  2:720632/1‑150 (MQ=255)
                 ggacATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACaa                                                                                                                                    <  2:3516871/147‑1 (MQ=255)
                  gCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAAt                                                                                                                                   <  2:3613753/150‑1 (MQ=255)
                  gCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAAt                                                                                                                                   <  2:1152697/150‑1 (MQ=255)
                     aTGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTACGAGCTTTTATATTGATTTACAATAag                                                                                                                                <  1:2121666/150‑1 (MQ=255)
                      tgtgAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAaga                                                                                                                               >  2:4227461/1‑150 (MQ=255)
                        tgAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGt                                                                                                                             >  1:4685789/1‑150 (MQ=255)
                             aaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAc                                                                                                                        <  1:135111/150‑1 (MQ=255)
                             aaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAc                                                                                                                        <  1:4411574/150‑1 (MQ=255)
                             aaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAc                                                                                                                        <  1:135600/150‑1 (MQ=255)
                             aaCATATTTCTGGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATCTGGATTTACAATAAGAGTCTCAc                                                                                                                        <  2:3392524/150‑1 (MQ=255)
                                atatTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAcctc                                                                                                                     <  2:2372884/150‑1 (MQ=255)
                                atatTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAcctc                                                                                                                     <  2:4822727/150‑1 (MQ=255)
                                 tatTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAcctcc                                                                                                                    <  2:1305450/150‑1 (MQ=255)
                                  atTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCAcctcct                                                                                                                   >  1:1708276/1‑150 (MQ=255)
                                  atTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAcctcct                                                                                                                   >  1:735827/1‑150 (MQ=255)
                                    ttCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGTTGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGa                                                                                                                 <  1:4371195/150‑1 (MQ=255)
                                    ttCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGa                                                                                                                 <  1:1919193/150‑1 (MQ=255)
                                     tCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAt                                                                                                                <  1:3427000/150‑1 (MQ=255)
                                      cTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAtt                                                                                                               <  1:3577741/150‑1 (MQ=255)
                                      cTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAtt                                                                                                               >  1:144185/1‑150 (MQ=255)
                                          cACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAg                                                                                                           >  2:1174706/1‑150 (MQ=255)
                                           aCTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGt                                                                                                          <  1:1162653/150‑1 (MQ=255)
                                                gAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAtct                                                                                                     <  2:469396/150‑1 (MQ=255)
                                                gAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAtct                                                                                                     <  2:3297269/150‑1 (MQ=255)
                                                gAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAtct                                                                                                     <  2:468803/150‑1 (MQ=255)
                                                   ttCATCGATTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCtt                                                                                                  <  2:3615773/150‑1 (MQ=255)
                                                   ttCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCtt                                                                                                  <  2:3615669/150‑1 (MQ=255)
                                                   ttCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCtt                                                                                                  <  2:3615762/150‑1 (MQ=255)
                                                   ttCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCtt                                                                                                  <  2:187430/150‑1 (MQ=255)
                                                    tCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTg                                                                                                 >  2:1361497/1‑150 (MQ=255)
                                                     cATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGa                                                                                                >  1:2995447/1‑150 (MQ=255)
                                                      aTCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGAt                                                                                               <  1:3306340/150‑1 (MQ=255)
                                                       tCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATg                                                                                              <  1:2717433/150‑1 (MQ=255)
                                                           gTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTaa                                                                                          <  1:3044125/150‑1 (MQ=255)
                                                             aaTCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGa                                                                                        >  2:863498/1‑150 (MQ=255)
                                                                ccATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACt                                                                                     >  2:3195326/1‑150 (MQ=255)
                                                                 cATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAGTAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTa                                                                                    <  2:2419914/150‑1 (MQ=255)
                                                                 cATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTa                                                                                    <  2:2419894/150‑1 (MQ=255)
                                                                  aTAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                                                   <  1:3762567/150‑1 (MQ=255)
                                                                      tAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTg                                                                               >  1:1822495/1‑150 (MQ=255)
                                                                        aGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTc                                                                             >  1:1094190/1‑150 (MQ=255)
                                                                        aGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTc                                                                             >  1:3221090/1‑150 (MQ=255)
                                                                         gAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCg                                                                            >  2:1585732/1‑150 (MQ=255)
                                                                         gAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCg                                                                            >  2:1586244/1‑150 (MQ=255)
                                                                          aaaCACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGa                                                                           >  1:1526425/1‑150 (MQ=255)
                                                                           aaCACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGAt                                                                          >  2:2417522/1‑150 (MQ=255)
                                                                           aaCACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGAt                                                                          >  2:1892221/1‑150 (MQ=255)
                                                                                    tGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCg                                                                 >  1:1762296/1‑150 (MQ=255)
                                                                                    tGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCg                                                                 >  1:1763332/1‑150 (MQ=255)
                                                                                     gAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGa                                                                <  1:2785169/150‑1 (MQ=255)
                                                                                        tttttGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATATCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACaa                                                             >  2:3341337/1‑150 (MQ=255)
                                                                                         ttttGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACCAg                                                            >  2:3293280/1‑150 (MQ=255)
                                                                                          tttGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAgg                                                           >  2:587875/1‑150 (MQ=255)
                                                                                              tAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACg                                                       >  2:757111/1‑150 (MQ=255)
                                                                                              tAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACg                                                       <  2:4309483/150‑1 (MQ=255)
                                                                                                      tCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCagagag                                               <  1:1656175/150‑1 (MQ=255)
                                                                                                      tCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCagagag                                               >  1:651976/1‑150 (MQ=255)
                                                                                                        cGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGaa                                             >  2:153083/1‑150 (MQ=255)
                                                                                                           aaaGATGTTTTCTTCCTGTTTCGTATTTTTTTAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACAt                                          <  1:916269/150‑1 (MQ=255)
                                                                                                            aaGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTTTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACAtt                                         >  2:3467911/1‑150 (MQ=255)
                                                                                                            aaGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACAtt                                         >  2:3467851/1‑150 (MQ=255)
                                                                                                                   tttCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCaa                                  >  2:1494526/1‑150 (MQ=255)
                                                                                                                     tcttcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAg                                <  2:3870985/150‑1 (MQ=255)
                                                                                                                        tcCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGAt                             <  1:1361497/150‑1 (MQ=255)
                                                                                                                           tGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGaa                          >  1:129849/1‑150 (MQ=255)
                                                                                                                                cGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGt                     <  2:3425022/150‑1 (MQ=255)
                                                                                                                                  tATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTaa                   <  2:4803848/150‑1 (MQ=255)
                                                                                                                                   aTTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAc                  <  2:2759989/150‑1 (MQ=255)
                                                                                                                                   aTTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAc                  <  2:3095422/150‑1 (MQ=255)
                                                                                                                                    ttttttGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAcc                 <  1:304801/150‑1 (MQ=255)
                                                                                                                                     tttttGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAccc                >  2:4480811/1‑150 (MQ=255)
                                                                                                                                          gAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAat           <  2:4107282/150‑1 (MQ=255)
                                                                                                                                                tGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTAtttt     <  1:3291745/150‑1 (MQ=255)
                                                                                                                                                   gCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTc  >  1:580233/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTC  >  CP009273/1192377‑1192673

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: