Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,412,205 C→A 56.4% D151Y (GAT→TAT)  yfcD ← putative NUDIX hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,412,2050CA56.4% 21.5 / inf 258D151Y (GAT→TAT) yfcDputative NUDIX hydrolase
Reads supporting (aligned to +/- strand):  ref base C (55/57);  new base A (86/59);  total (141/117)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

CGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGCCCGTGCTCG  >  CP009273/2412056‑2412353
                                                                                                                                                     |                                                                                                                                                    
cgcTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATc                                                                                                                                                      >  2:1764619/1‑150 (MQ=255)
 gcTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCg                                                                                                                                                     >  2:3740928/1‑150 (MQ=255)
   tCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCa                                                                                                                                                   <  2:1363977/150‑1 (MQ=255)
   tCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCa                                                                                                                                                   <  2:2963754/150‑1 (MQ=255)
   tCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCa                                                                                                                                                   <  2:2668773/150‑1 (MQ=255)
   tCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCa                                                                                                                                                   >  2:2870824/1‑150 (MQ=255)
   tCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCa                                                                                                                                                   >  2:223529/1‑150 (MQ=255)
     aGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGc                                                                                                                                                 >  1:652605/1‑150 (MQ=255)
      gTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCg                                                                                                                                                >  2:142850/1‑150 (MQ=255)
           tGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCgg                                                                                                                                           >  1:1991543/1‑150 (MQ=255)
             ggCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTg                                                                                                                                         >  1:2691048/1‑150 (MQ=255)
              gCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGa                                                                                                                                        <  1:3237188/150‑1 (MQ=255)
               cTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGAt                                                                                                                                       <  1:1306163/150‑1 (MQ=255)
                tACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGAtt                                                                                                                                      >  1:2515228/1‑150 (MQ=255)
                   cACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATttct                                                                                                                                   <  1:2275818/150‑1 (MQ=255)
                    aCTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATttctt                                                                                                                                  <  2:2566605/150‑1 (MQ=255)
                     cTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAAATCATAGCAGCGTGCGGTGATttcttc                                                                                                                                 >  1:780495/1‑150 (MQ=255)
                        aaaaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCgg                                                                                                                              >  2:559602/1‑150 (MQ=255)
                        aaaaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCgg                                                                                                                              >  2:1962407/1‑150 (MQ=255)
                          aaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCg                                                                                                                            <  1:712454/150‑1 (MQ=255)
                          aaTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCg                                                                                                                            >  1:546973/1‑150 (MQ=255)
                           aTGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGt                                                                                                                           >  1:1658038/1‑150 (MQ=255)
                             gTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTca                                                                                                                         >  1:4119203/1‑150 (MQ=255)
                              ttAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTcag                                                                                                                        >  1:1854649/1‑150 (MQ=255)
                              ttAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTcag                                                                                                                        >  1:1854659/1‑150 (MQ=255)
                                 agagCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTcagcca                                                                                                                     <  1:1222428/150‑1 (MQ=255)
                                 agagCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTcagcca                                                                                                                     >  2:3706692/1‑150 (MQ=255)
                                   agCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTcagccagc                                                                                                                   <  2:270178/150‑1 (MQ=255)
                                   agCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTcagccagc                                                                                                                   >  1:2454615/1‑150 (MQ=255)
                                      cGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCaaa                                                                                                                >  1:1953955/1‑150 (MQ=255)
                                        cAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACt                                                                                                              >  2:3373730/1‑150 (MQ=255)
                                         taTTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACtt                                                                                                             <  2:3363033/149‑1 (MQ=255)
                                         aagtGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACtt                                                                                                             <  2:1207786/147‑1 (MQ=255)
                                         aaTTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACtt                                                                                                             >  2:3660333/1‑150 (MQ=255)
                                         aaTTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACtt                                                                                                             >  1:4011238/1‑150 (MQ=255)
                                               ggCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGa                                                                                                       >  2:1664412/1‑150 (MQ=255)
                                                gCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGAc                                                                                                      >  2:3982229/1‑150 (MQ=255)
                                                gCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGAc                                                                                                      >  2:4102666/1‑150 (MQ=255)
                                                gCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGAc                                                                                                      >  2:3965688/1‑150 (MQ=255)
                                                gCTCTTTTTTAATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGATAGCGCTTTCAGCGAGACTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGAGTCAGCCAGCAAACTTCACTGAc                                                                                                      >  2:2196960/1‑150 (MQ=255)
                                                 ctctTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACt                                                                                                     >  1:1677366/1‑150 (MQ=255)
                                                   ctTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTc                                                                                                   >  1:2376698/1‑150 (MQ=255)
                                                   ctTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTc                                                                                                   <  2:274085/150‑1 (MQ=255)
                                                     tttttCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCAt                                                                                                 >  1:1553650/1‑150 (MQ=255)
                                                     tttttCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCAt                                                                                                 <  2:1134079/150‑1 (MQ=255)
                                                      ttttCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATc                                                                                                <  1:3826155/150‑1 (MQ=255)
                                                        ttCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCtt                                                                                              <  2:3610872/150‑1 (MQ=255)
                                                             tCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGt                                                                                         >  2:4129692/1‑150 (MQ=255)
                                                              cTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTa                                                                                        >  1:130244/1‑150 (MQ=255)
                                                               tGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAg                                                                                       >  1:2869619/1‑150 (MQ=255)
                                                               tGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAg                                                                                       <  1:3512470/150‑1 (MQ=255)
                                                               tGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAg                                                                                       >  1:3006912/1‑150 (MQ=255)
                                                                gCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGt                                                                                      >  2:84289/1‑150 (MQ=255)
                                                                 ccGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTg                                                                                     <  1:3471662/150‑1 (MQ=255)
                                                                   gTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCg                                                                                   >  2:3778196/1‑150 (MQ=255)
                                                                   gTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCg                                                                                   >  1:2952811/1‑150 (MQ=255)
                                                                   gTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCg                                                                                   >  1:2742167/1‑150 (MQ=255)
                                                                      tCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGAGATTTCTTCCGGCGTCAGCCAGCACACTACACTGACTTCATCATCCTGTAGTGCGAAg                                                                                >  2:2185055/1‑150 (MQ=255)
                                                                        aGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAggg                                                                              >  1:2198414/1‑150 (MQ=255)
                                                                         gTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGa                                                                             <  2:1048167/150‑1 (MQ=255)
                                                                          tctcTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGAc                                                                            <  2:4081960/150‑1 (MQ=255)
                                                                          tctcTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGAc                                                                            <  2:2415001/150‑1 (MQ=255)
                                                                          tctcTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGAc                                                                            <  1:3190976/150‑1 (MQ=255)
                                                                             ctACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGt                                                                         <  1:2578212/150‑1 (MQ=255)
                                                                             ctACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGt                                                                         >  1:3621028/1‑150 (MQ=255)
                                                                               aCGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTga                                                                       >  1:2038803/1‑150 (MQ=255)
                                                                               aCGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTga                                                                       >  1:2038928/1‑150 (MQ=255)
                                                                                 ggCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTgaga                                                                     >  1:3395544/1‑150 (MQ=255)
                                                                                   cTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTgagaga                                                                   >  2:3355651/1‑150 (MQ=255)
                                                                                   cTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTgagaga                                                                   <  1:3901711/150‑1 (MQ=255)
                                                                                    ttCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGAc                                                                  <  2:401744/150‑1 (MQ=255)
                                                                                       aTTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCa                                                               >  1:3249289/1‑150 (MQ=255)
                                                                                        tttttGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAg                                                              >  1:3719038/1‑150 (MQ=255)
                                                                                          tttGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCt                                                            <  1:1231503/150‑1 (MQ=255)
                                                                                               cATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACa                                                       >  2:4161025/1‑150 (MQ=255)
                                                                                                  tGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATg                                                    <  2:3791434/150‑1 (MQ=255)
                                                                                                   gcgcTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATTGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATgc                                                   <  1:3311306/150‑1 (MQ=255)
                                                                                                    cgcTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATgcg                                                  >  1:830224/1‑150 (MQ=255)
                                                                                                      cTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGcc                                                >  1:689501/1‑150 (MQ=255)
                                                                                                            ccACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGac                                          <  2:3759659/150‑1 (MQ=255)
                                                                                                              acaACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCTTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGacac                                        <  2:161346/150‑1 (MQ=255)
                                                                                                                aaCGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACgg                                      >  2:263269/1‑150 (MQ=255)
                                                                                                                  cGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCa                                    >  2:349431/1‑150 (MQ=255)
                                                                                                                   gCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCaa                                   <  1:3949889/150‑1 (MQ=255)
                                                                                                                      aGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAttt                                <  2:3625642/150‑1 (MQ=255)
                                                                                                                       gcgcTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAtttt                               <  2:2585762/150‑1 (MQ=255)
                                                                                                                       gcgcTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAtttt                               <  2:3088233/150‑1 (MQ=255)
                                                                                                                       gcgcTTTCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAAtttt                               >  2:2043912/1‑150 (MQ=255)
                                                                                                                         gcTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTa                             <  2:1256651/150‑1 (MQ=255)
                                                                                                                           tttCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATc                           <  2:3856510/150‑1 (MQ=255)
                                                                                                                             tCAGCGAGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCtt                         <  2:716879/150‑1 (MQ=255)
                                                                                                                                   aGTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAAt                   <  1:3866599/150‑1 (MQ=255)
                                                                                                                                    gTCTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATa                  >  2:3332248/1‑150 (MQ=255)
                                                                                                                                      cTGGAGTGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGa                >  2:2317425/1‑150 (MQ=255)
                                                                                                                                            tGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGcc          <  2:2786170/150‑1 (MQ=255)
                                                                                                                                            tGAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGcc          >  1:2255828/1‑150 (MQ=255)
                                                                                                                                             gAACTCATAGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGccc         >  2:2135812/1‑150 (MQ=255)
                                                                                                                                                cTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGCCCGTg      >  2:661772/1‑150 (MQ=255)
                                                                                                                                                  cATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGCCCGTGCt    >  1:3012789/1‑150 (MQ=255)
                                                                                                                                                    tCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGCCCGTGCTCg  <  1:3070394/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCTCAGTTGATGGGCTACCACTAAAAATGTTAAGAGCCGCAATTGCGGCTCTTTTTTCATTCTGCCGTTTCAGTCTCTACGGCTTCATTTTTGGCATTGCGCTTCATCCACAACGCTAGCGCTTTCAGCGAGTCTGGAGTGAACTCATCGCAGCGTGCGGTGATTTCTTCCGGCGTCAGCCAGCAAACTTCACTGACTTCATCTTCCTGTAGTGCGAAGGGACCGTGAGAGACGCAGCTGAACAATGCGCCCCAGACACGGCAATTTTTATCTTCGAAATAGAACTGCCCGTGCTCG  >  CP009273/2412056‑2412353

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: