Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,401,343 A→T 55.8% S75S (TCA→TCT yjfK → DUF2491 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,401,3430AT55.8% 86.2 / inf 236S75S (TCA→TCTyjfKDUF2491 family protein
Reads supporting (aligned to +/- strand):  ref base A (45/59);  new base T (63/69);  total (108/128)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 241 total aligned reads displayed.

CTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGCGATGGGAGCAAT  >  CP009273/4401194‑4401492
                                                                                                                                                     |                                                                                                                                                     
ctcAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCa                                                                                                                                                       <  2:4381772/150‑1 (MQ=255)
ctcAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCa                                                                                                                                                       <  2:1375962/150‑1 (MQ=255)
atcAATAGTGGCTTTACGCTCGACACGGTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGATGAATTTACGGTAGCCGCCGTCAGCCACATCGATATGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCa                                                                                                                                                       >  2:3060882/2‑150 (MQ=255)
    aTAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCg                                                                                                                                                   <  2:1191861/150‑1 (MQ=255)
    aTAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCg                                                                                                                                                   >  2:2916203/1‑150 (MQ=255)
    aTAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCg                                                                                                                                                   >  2:1400820/1‑150 (MQ=255)
      aGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGAt                                                                                                                                                 >  1:466223/1‑150 (MQ=255)
         ggCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAg                                                                                                                                              >  1:726694/1‑150 (MQ=255)
          gCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGt                                                                                                                                             >  2:733800/1‑150 (MQ=255)
            tttACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGttt                                                                                                                                           <  2:3732617/150‑1 (MQ=255)
               aCGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTg                                                                                                                                        <  1:3856569/150‑1 (MQ=255)
               aCGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTg                                                                                                                                        <  1:571919/150‑1 (MQ=255)
               aCGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTg                                                                                                                                        <  1:3626633/150‑1 (MQ=255)
                  cTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCaa                                                                                                                                     <  2:7021/150‑1 (MQ=255)
                    cGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAAt                                                                                                                                   <  2:2242294/150‑1 (MQ=255)
                    cGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAAt                                                                                                                                   >  2:4037347/1‑150 (MQ=255)
                     gACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATc                                                                                                                                  <  1:4240117/150‑1 (MQ=255)
                      acacGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACGTCTGGCGATGAGTTTCTGCAAATCa                                                                                                                                 >  1:4559202/1‑150 (MQ=255)
                        acGGTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGGAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAAt                                                                                                                               >  1:214638/1‑150 (MQ=255)
                           ttAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATacc                                                                                                                            <  2:4431682/150‑1 (MQ=255)
                            tAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATacca                                                                                                                           >  2:848100/1‑150 (MQ=255)
                            tAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATacca                                                                                                                           <  2:4231599/150‑1 (MQ=255)
                             aGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATaccac                                                                                                                          >  1:4039476/1‑150 (MQ=255)
                             aGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATaccac                                                                                                                          >  1:1007629/1‑150 (MQ=255)
                             aGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATaccac                                                                                                                          >  1:1802878/1‑150 (MQ=255)
                             aGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATaccac                                                                                                                          <  2:3618882/150‑1 (MQ=255)
                                gtttcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACcgg                                                                                                                       <  1:4063884/150‑1 (MQ=255)
                                gtttcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACcgg                                                                                                                       >  2:3300516/1‑150 (MQ=255)
                                gtttcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACcgg                                                                                                                       <  1:993918/150‑1 (MQ=255)
                                      tttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGa                                                                                                                 >  1:132936/1‑150 (MQ=255)
                                      tttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGa                                                                                                                 >  2:2015596/1‑150 (MQ=255)
                                       ttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGANTTTCTGCAAATCAATACCACCGGCGGCGaa                                                                                                                <  1:4286712/150‑1 (MQ=255)
                                        tGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAg                                                                                                               <  1:1757807/150‑1 (MQ=255)
                                         gCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGa                                                                                                              >  2:2598050/1‑150 (MQ=255)
                                          cTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGat                                                                                                             >  1:800668/1‑150 (MQ=255)
                                           tGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGata                                                                                                            <  2:1125244/150‑1 (MQ=255)
                                             gAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATAtt                                                                                                          <  2:1643363/150‑1 (MQ=255)
                                               aGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATtga                                                                                                        >  1:377608/1‑150 (MQ=255)
                                                gATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATtgat                                                                                                       <  2:676135/150‑1 (MQ=255)
                                                 aTGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATtgatg                                                                                                      <  2:2431516/150‑1 (MQ=255)
                                                 aTGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATtgatg                                                                                                      >  1:828586/1‑150 (MQ=255)
                                                  tGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATtgatga                                                                                                     >  1:4552693/1‑150 (MQ=255)
                                                     aTTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACAt                                                                                                  >  2:2973649/1‑150 (MQ=255)
                                                         cTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATtgat                                                                                              <  2:3610712/150‑1 (MQ=255)
                                                           gATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATtgatga                                                                                            >  2:3773725/1‑150 (MQ=255)
                                                            aTTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTgatgat                                                                                           >  1:3029998/1‑150 (MQ=255)
                                                            aTTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTgatgat                                                                                           >  2:130600/1‑150 (MQ=255)
                                                            aTTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTgatgat                                                                                           >  1:2056344/1‑150 (MQ=255)
                                                             ttGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGata                                                                                          <  1:160280/150‑1 (MQ=255)
                                                             ttGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGata                                                                                          <  1:1670583/150‑1 (MQ=255)
                                                             ttGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGata                                                                                          <  1:2366761/150‑1 (MQ=255)
                                                               gcgcTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATc                                                                                        <  2:873126/150‑1 (MQ=255)
                                                                cgcTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCa                                                                                       >  2:4095467/1‑150 (MQ=255)
                                                                  cTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAg                                                                                     >  1:2406085/1‑150 (MQ=255)
                                                                    gCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGct                                                                                   <  1:3749079/150‑1 (MQ=255)
                                                                      cGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGctct                                                                                 <  2:1552346/150‑1 (MQ=255)
                                                                        ggTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCttt                                                                               <  1:2068122/150‑1 (MQ=255)
                                                                        ggTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCttt                                                                               >  1:1518929/1‑150 (MQ=255)
                                                                           gAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTc                                                                            >  1:1270419/1‑150 (MQ=255)
                                                                            aGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCt                                                                           <  1:391477/150‑1 (MQ=255)
                                                                               aaTTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATg                                                                        <  2:3965459/150‑1 (MQ=255)
                                                                               aaTTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATg                                                                        >  1:4114061/1‑150 (MQ=255)
                                                                               aaTTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATg                                                                        >  1:783742/1‑150 (MQ=255)
                                                                                   tACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAaga                                                                    <  2:192237/150‑1 (MQ=255)
                                                                                     cGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAagaga                                                                  <  2:3339963/150‑1 (MQ=255)
                                                                                      ggTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAagagag                                                                 <  1:511643/150‑1 (MQ=255)
                                                                                      ggTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAagagag                                                                 <  2:4427632/150‑1 (MQ=255)
                                                                                       gTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGc                                                                <  2:1725340/150‑1 (MQ=255)
                                                                                          gccgccGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTAc                                                             >  1:4191959/1‑150 (MQ=255)
                                                                                           ccgccgTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACg                                                            <  1:3860504/150‑1 (MQ=255)
                                                                                           ccgccgTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACg                                                            <  1:1988216/150‑1 (MQ=255)
                                                                                             gccgTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAATCTCTTTGTCTATGAAGAGAGCTACGGt                                                          >  1:2627320/1‑150 (MQ=255)
                                                                                               cgTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTAt                                                        >  2:1432325/1‑150 (MQ=255)
                                                                                                gTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATc                                                       <  1:330135/150‑1 (MQ=255)
                                                                                                 tCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATGGATGAGATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCa                                                      <  2:3569973/150‑1 (MQ=255)
                                                                                                   aGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGt                                                    <  2:182417/150‑1 (MQ=255)
                                                                                                            gATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGt                                           >  2:798284/1‑150 (MQ=255)
                                                                                                             aTCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTc                                          >  1:2167362/1‑150 (MQ=255)
                                                                                                               cTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTCTGAAGAGAGCTACGGTATCAGTAAAGAGAGTCAc                                        >  2:632690/1‑150 (MQ=255)
                                                                                                                tGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACt                                       >  1:4356871/1‑150 (MQ=255)
                                                                                                                    cggcggTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGgcg                                   >  1:2506314/1‑150 (MQ=255)
                                                                                                                     ggcggTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGgcgc                                  >  1:2863706/1‑150 (MQ=255)
                                                                                                                      gcggTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGcgcg                                 <  2:3140702/150‑1 (MQ=255)
                                                                                                                      gcggTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGcgcg                                 <  2:1333093/150‑1 (MQ=255)
                                                                                                                      gcggTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGcgcg                                 <  2:4054592/150‑1 (MQ=255)
                                                                                                                       cggTAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGa                                <  2:466223/150‑1 (MQ=255)
                                                                                                                          tAGTCAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGc                             <  1:3919227/150‑1 (MQ=255)
                                                                                                                              cAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATc                         <  1:1847846/150‑1 (MQ=255)
                                                                                                                              cAGATTTTCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATc                         >  1:1092121/1‑150 (MQ=255)
                                                                                                                                   tttCCGCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGc                    >  1:2540723/1‑150 (MQ=255)
                                                                                                                                        gCTACTACACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCTGCGAAGATATTGATGACATTTATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAg               <  1:3219586/150‑1 (MQ=255)
                                                                                                                                         ctactaCACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGc              <  2:709322/150‑1 (MQ=255)
                                                                                                                                         ctactaCACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGc              >  2:2734326/1‑150 (MQ=255)
                                                                                                                                         ctactaCACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGc              >  2:2260358/1‑150 (MQ=255)
                                                                                                                                            ctacACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGCGAt           <  2:3422663/150‑1 (MQ=255)
                                                                                                                                            ctacACTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGCGAt           >  1:786470/1‑150 (MQ=255)
                                                                                                                                              acacTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGCGATgg         >  2:1365000/1‑150 (MQ=255)
                                                                                                                                              acacTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGCGATgg         <  2:3893570/150‑1 (MQ=255)
                                                                                                                                                 cTTCTGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGCGATGGGAg      >  2:393502/1‑150 (MQ=255)
                                                                                                                                                     aGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGCGATGGGAGCAAt  <  1:2719834/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGCGGCGAAGATATTGATGACATTGATGATATCAAGCTCTTTGTCTATGAAGAGAGCTACGGTATCAGTAAAGAGAGTCACTGGCGCGAGGCGATCAACGCCAAAGCGATGGGAGCAAT  >  CP009273/4401194‑4401492

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: