Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 461,499 G→T 100% R424S (CGC→AGC)  ybaE ← putative ABC superfamily transporter periplasmic‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273461,4990GT98.7% 994.9 / 2.0 312R424S (CGC→AGC) ybaEputative ABC superfamily transporter periplasmic‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (156/152);  total (159/153)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 314 total aligned reads displayed.

CCAGAGTATATTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGAC  >  CP009273/461352‑461643
                                                                                                                                                   |                                                                                                                                                
 cAGAGTATATTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTc                                                                                                                                               >  1:358118/1‑150 (MQ=255)
 cAGAGTATATTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTc                                                                                                                                               <  2:2570342/150‑1 (MQ=255)
      tataTTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCa                                                                                                                                          >  2:5152435/1‑150 (MQ=255)
      tataTTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCa                                                                                                                                          >  2:3680849/1‑150 (MQ=255)
       atatTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCAt                                                                                                                                         <  2:3037742/150‑1 (MQ=255)
       agagtCCGGGGCTTCGCCAAATACTCTGTCGCCCATCTTGAGGGCTGCGTGTGCCTGTAGGGTCGTGTCGCCCCAGTTTATTGCGTTCTGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCAt                                                                                                                                         <  2:2970657/146‑1 (MQ=255)
          ttCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGt                                                                                                                                      <  1:1019298/150‑1 (MQ=255)
           tCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTa                                                                                                                                     >  2:3265695/1‑150 (MQ=255)
           tCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTa                                                                                                                                     >  2:3439479/1‑150 (MQ=255)
            ccGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTAt                                                                                                                                    >  1:1377620/1‑150 (MQ=255)
             cGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATg                                                                                                                                   <  1:4534/150‑1 (MQ=255)
             cGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGGTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATg                                                                                                                                   <  1:1762914/150‑1 (MQ=255)
              ggTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGa                                                                                                                                  >  1:2917549/1‑150 (MQ=255)
                tGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAg                                                                                                                                >  2:2159618/1‑150 (MQ=255)
                  cTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGtt                                                                                                                              <  1:3143161/150‑1 (MQ=255)
                  cTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGtt                                                                                                                              <  1:2517464/150‑1 (MQ=255)
                   ttCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTc                                                                                                                             >  2:1984420/1‑150 (MQ=255)
                     cGCCAATTACTCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTa                                                                                                                           <  1:2622786/150‑1 (MQ=255)
                     cGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTa                                                                                                                           <  2:3537869/150‑1 (MQ=255)
                        cAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCg                                                                                                                        <  1:270050/150‑1 (MQ=255)
                         aaTTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCgg                                                                                                                       >  2:1862139/1‑150 (MQ=255)
                          aTTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGt                                                                                                                      >  2:1932852/1‑150 (MQ=255)
                           ttAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTa                                                                                                                     <  1:2871245/150‑1 (MQ=255)
                            tAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAg                                                                                                                    >  1:608871/1‑150 (MQ=255)
                            tAATCTGTCGCCCATCATCAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTCGGCTTTCTGCCATGGAATGAAGTTCTATCGGTAg                                                                                                                    >  2:3726860/1‑150 (MQ=255)
                                   tCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATaa                                                                                                             >  2:155758/1‑150 (MQ=255)
                                        catcatGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAg                                                                                                        >  2:1267316/1‑150 (MQ=255)
                                          tcatGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCg                                                                                                      <  1:3660716/150‑1 (MQ=255)
                                           catGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGt                                                                                                     >  1:130138/1‑150 (MQ=255)
                                             tGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCa                                                                                                   <  1:489870/150‑1 (MQ=255)
                                             tGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCATCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCa                                                                                                   >  2:3493056/1‑150 (MQ=255)
                                              gAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATACTTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCAGGGTCTGAAGTTCTATCGGTCGGTGATAAACCAGCGTCaa                                                                                                  >  2:2216161/1‑150 (MQ=255)
                                              gAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCaa                                                                                                  <  1:1076151/150‑1 (MQ=255)
                                              gAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCaa                                                                                                  >  2:437854/1‑150 (MQ=255)
                                              gAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCaa                                                                                                  <  2:32278/150‑1 (MQ=255)
                                              gAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTCTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCaa                                                                                                  >  2:3211462/1‑150 (MQ=255)
                                                 gTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGt                                                                                               >  1:3674604/1‑150 (MQ=255)
                                                    tGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGtttt                                                                                            >  2:762983/1‑150 (MQ=255)
                                                     gCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTc                                                                                           <  1:99607/150‑1 (MQ=255)
                                                      cGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCg                                                                                          <  2:4242320/150‑1 (MQ=255)
                                                         gtgCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTa                                                                                       <  2:3404383/150‑1 (MQ=255)
                                                          tgCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAg                                                                                      <  2:4475840/150‑1 (MQ=255)
                                                           gCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGt                                                                                     >  1:92897/1‑150 (MQ=255)
                                                             cTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGttt                                                                                   >  2:232928/1‑150 (MQ=255)
                                                              tGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTg                                                                                  <  2:1665569/150‑1 (MQ=255)
                                                               gTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGa                                                                                 >  2:60324/1‑150 (MQ=255)
                                                                  gggTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACtt                                                                              <  1:4337858/150‑1 (MQ=255)
                                                                          tcgtcCCAGTTTTTTTCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGt                                                                      <  1:4823774/150‑1 (MQ=255)
                                                                          tcgtcCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGt                                                                      >  2:2547720/1‑150 (MQ=255)
                                                                           cgtcCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTa                                                                     <  2:641213/150‑1 (MQ=255)
                                                                             tcCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTAcc                                                                   <  1:3439479/150‑1 (MQ=255)
                                                                              cccAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCt                                                                  <  1:2316541/150‑1 (MQ=255)
                                                                               ccAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTg                                                                 <  2:1477618/150‑1 (MQ=255)
                                                                               ccAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTg                                                                 >  2:3358514/1‑150 (MQ=255)
                                                                                 aGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGcc                                                               >  1:4572393/1‑150 (MQ=255)
                                                                                 aGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGcc                                                               <  2:3960390/150‑1 (MQ=255)
                                                                                  gTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCa                                                              <  1:2965576/150‑1 (MQ=255)
                                                                                  gTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCa                                                              >  2:4554437/1‑150 (MQ=255)
                                                                                     ttttGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATg                                                           <  2:4683388/150‑1 (MQ=255)
                                                                                       ttGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCg                                                         <  1:3994433/150‑1 (MQ=255)
                                                                                          cGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGaa                                                      >  2:4157805/1‑150 (MQ=255)
                                                                                               tGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTc                                                 >  2:4936518/1‑150 (MQ=255)
                                                                                                gAAAAATAATTGTGAGTTCACAGCCGTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTcc                                                <  2:1175201/150‑1 (MQ=255)
                                                                                                  aaaaTAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAg                                              >  2:4982570/1‑150 (MQ=255)
                                                                                                     ataatTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCt                                           <  1:3807484/150‑1 (MQ=255)
                                                                                                         ttGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCa                                       >  1:3900468/1‑150 (MQ=255)
                                                                                                              aGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAAt                                  >  2:1611616/1‑150 (MQ=255)
                                                                                                               gTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATg                                 >  1:1523475/1‑150 (MQ=255)
                                                                                                                ttCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGc                                <  2:840205/150‑1 (MQ=255)
                                                                                                                ttCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGc                                <  1:1411479/150‑1 (MQ=255)
                                                                                                                 tCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCa                               <  1:2369673/150‑1 (MQ=255)
                                                                                                                   acaGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATg                             >  1:3205357/1‑150 (MQ=255)
                                                                                                                     aGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGAc                           <  2:4888470/150‑1 (MQ=255)
                                                                                                                     aGCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGAc                           >  1:4453049/1‑150 (MQ=255)
                                                                                                                      gCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACt                          <  1:1153521/150‑1 (MQ=255)
                                                                                                                      gCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACt                          <  2:3025237/150‑1 (MQ=255)
                                                                                                                      gCCTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACt                          <  1:5052888/150‑1 (MQ=255)
                                                                                                                        cTTCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTgg                        >  1:1261847/1‑150 (MQ=255)
                                                                                                                         ttCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGc                       <  1:1984420/150‑1 (MQ=255)
                                                                                                                         ttCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGc                       >  2:753585/1‑150 (MQ=255)
                                                                                                                         ttCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGc                       >  2:1687914/1‑150 (MQ=255)
                                                                                                                         ttCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGc                       <  2:2057507/150‑1 (MQ=255)
                                                                                                                         ttCCGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGGAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGc                       <  2:4614191/150‑1 (MQ=255)
                                                                                                                           ccGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCgg                     >  1:4073534/1‑150 (MQ=255)
                                                                                                                            cGCTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGt                    <  1:2943008/150‑1 (MQ=255)
                                                                                                                              cTGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGa                  >  1:226958/1‑150 (MQ=255)
                                                                                                                               tGCCAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGAt                 <  1:4624951/150‑1 (MQ=255)
                                                                                                                                 ccAGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCa               <  2:3860101/150‑1 (MQ=255)
                                                                                                                                   aGTGTCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTGGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAgg             <  2:717058/150‑1 (MQ=255)
                                                                                                                                     tgtCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGtt           <  1:2046859/150‑1 (MQ=255)
                                                                                                                                       tCGCCTGTAGGCTTTCTTCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGtttt         <  1:1926887/150‑1 (MQ=255)
                                                                                                                                        cGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTtc        <  2:2789865/150‑1 (MQ=255)
                                                                                                                                        cGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTtc        <  1:762983/150‑1 (MQ=255)
                                                                                                                                        cGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTtc        <  1:243291/150‑1 (MQ=255)
                                                                                                                                         gCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTtct       >  1:4708898/1‑150 (MQ=255)
                                                                                                                                          ccTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTtctc      >  1:4592589/1‑150 (MQ=255)
                                                                                                                                           cTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTcc     <  1:5090618/150‑1 (MQ=255)
                                                                                                                                            tGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCg    >  1:3128916/1‑150 (MQ=255)
                                                                                                                                             gTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGa   >  1:4208556/1‑150 (MQ=255)
                                                                                                                                              tAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACAAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCTTGACTGGCGGTGATCAGGTTTTCTCCGAc  <  1:1831111/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCAGAGTATATTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGAC  >  CP009273/461352‑461643

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: