Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,373,985 C→A 100% intergenic (+37/‑117) yhcB → / → degQ DUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,373,9850CA99.7% 1135.5 / ‑4.5 367intergenic (+37/‑117)yhcB/degQDUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (147/217);  total (147/220)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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 382 total aligned reads displayed.

AATCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTGCA  >  CP009273/3373838‑3374131
                                                                                                                                                   |                                                                                                                                                  
 aTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCCTTGCGCCCTCCTCTtatct                                                                                                                                                 >  1:1705866/1‑150 (MQ=255)
 aTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtatct                                                                                                                                                 <  2:687690/150‑1 (MQ=255)
  tCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtatctc                                                                                                                                                <  2:3665545/150‑1 (MQ=255)
   cGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctcc                                                                                                                                               <  1:5053242/150‑1 (MQ=255)
   cGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctcc                                                                                                                                               >  1:2052714/1‑150 (MQ=255)
    gTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccc                                                                                                                                              >  1:2108382/1‑150 (MQ=255)
      cTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccctc                                                                                                                                            >  2:31522/1‑150 (MQ=255)
          cAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCg                                                                                                                                        <  2:1304975/150‑1 (MQ=255)
          cAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAAGTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCg                                                                                                                                        >  2:2670629/1‑150 (MQ=255)
                cTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATc                                                                                                                                  >  1:586192/1‑150 (MQ=255)
                  gAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAt                                                                                                                                >  2:416925/1‑150 (MQ=255)
                        ataaaataTCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCGGCAGCCTACAGGCGCGAAGCGCGACTAATATATTATACGGGCGCAGCCATAGCGCCCACCACTAATCTCCCTCCCCGACAAGCAAATAATc                                                                                                                          >  2:1326764/8‑150 (MQ=255)
                         gCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCt                                                                                                                         <  1:3056853/150‑1 (MQ=255)
                                  aGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAAGTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCa                                                                                                                >  2:409397/1‑150 (MQ=255)
                                   ggCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCAt                                                                                                               <  1:4078767/150‑1 (MQ=255)
                                   ggCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCAt                                                                                                               >  1:3355584/1‑150 (MQ=255)
                                    gCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCAtt                                                                                                              <  2:4447020/150‑1 (MQ=255)
                                       ccGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGtt                                                                                                           <  1:1089644/150‑1 (MQ=255)
                                       ccGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGtt                                                                                                           <  1:3128609/150‑1 (MQ=255)
                                       ccGGTACAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGtt                                                                                                           <  2:1608626/150‑1 (MQ=255)
                                        cGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCATATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTa                                                                                                          >  2:249829/1‑150 (MQ=255)
                                         ggTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAg                                                                                                         <  2:4940255/150‑1 (MQ=255)
                                         ggTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAg                                                                                                         <  1:2431798/150‑1 (MQ=255)
                                         ggTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAg                                                                                                         <  2:3321875/150‑1 (MQ=255)
                                           tGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGcc                                                                                                       <  2:2835149/150‑1 (MQ=255)
                                             cAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGt                                                                                                     >  2:650939/1‑150 (MQ=255)
                                               gATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCt                                                                                                   <  1:5195002/150‑1 (MQ=255)
                                                aTGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTg                                                                                                  <  2:2291409/150‑1 (MQ=255)
                                                aTGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTg                                                                                                  <  1:183870/150‑1 (MQ=255)
                                                 tGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGa                                                                                                 >  1:680291/1‑150 (MQ=255)
                                                    cTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGaaaa                                                                                              >  2:3268317/1‑150 (MQ=255)
                                                       gcgACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTc                                                                                           <  1:2080259/150‑1 (MQ=255)
                                                       gcgACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTCTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTc                                                                                           >  1:4727168/1‑150 (MQ=255)
                                                         gACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCaa                                                                                         >  2:1526235/1‑150 (MQ=255)
                                                         gACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCaa                                                                                         >  1:4886871/1‑150 (MQ=255)
                                                            tATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATaa                                                                                      <  1:964742/150‑1 (MQ=255)
                                                            tATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATaa                                                                                      <  2:1500870/150‑1 (MQ=255)
                                                            tATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATaa                                                                                      >  2:3306709/1‑150 (MQ=255)
                                                               tCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACAt                                                                                   <  1:4555713/150‑1 (MQ=255)
                                                                cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATc                                                                                  <  2:2814211/150‑1 (MQ=255)
                                                                 tGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCa                                                                                 <  2:2164263/150‑1 (MQ=255)
                                                                 tGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCa                                                                                 <  2:586192/150‑1 (MQ=255)
                                                                  gAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaa                                                                                >  1:4668009/1‑150 (MQ=255)
                                                                  gAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaa                                                                                >  2:2893433/1‑150 (MQ=255)
                                                                     ggCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACt                                                                             <  2:401861/150‑1 (MQ=255)
                                                                     ggCCCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACt                                                                             <  1:729209/150‑1 (MQ=255)
                                                                        gcATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGtt                                                                          >  2:2312512/1‑150 (MQ=255)
                                                                          aTCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGtttt                                                                        <  1:279062/150‑1 (MQ=255)
                                                                          aTCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTACCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGtttt                                                                        <  2:2331503/150‑1 (MQ=255)
                                                                           tCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTg                                                                       >  1:1228530/1‑150 (MQ=255)
                                                                            ccGGCCCGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGa                                                                      >  2:3674529/1‑150 (MQ=255)
                                                                              ggCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                    <  1:1873756/150‑1 (MQ=255)
                                                                              ggCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                    >  2:911314/1‑150 (MQ=255)
                                                                              ggCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                    >  2:4712926/1‑150 (MQ=255)
                                                                              ggCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                    <  1:4470493/150‑1 (MQ=255)
                                                                              ggCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                    >  2:2059319/1‑150 (MQ=255)
                                                                              ggCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAt                                                                    >  2:1167805/1‑150 (MQ=255)
                                                                                ccTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtct                                                                  >  2:30503/1‑150 (MQ=255)
                                                                                 ctgctgCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtctc                                                                 >  2:81426/1‑150 (MQ=255)
                                                                                 ctgctgCGTACTGGCGCGAAGAGCGACTAATTTAATTTTCGGGCGCAGCCATAGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtctc                                                                 <  2:3538461/150‑1 (MQ=255)
                                                                                  tgctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATctct                                                                >  2:1912043/1‑150 (MQ=255)
                                                                                  tgctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATctct                                                                <  2:3790912/150‑1 (MQ=255)
                                                                                     tgcGTACTGGCGCGAAGCGCGACTAATTTTTTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCGCCCTCCCCGACTATCATTTAATCTGGTGTCTCAGTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCtttt                                                             <  1:1035375/150‑1 (MQ=255)
                                                                                     tgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCtttt                                                             >  1:2897934/1‑150 (MQ=255)
                                                                                      gcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTc                                                            <  2:3926564/150‑1 (MQ=255)
                                                                                      gcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTc                                                            <  1:1491529/150‑1 (MQ=255)
                                                                                       cGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCt                                                           <  2:259457/150‑1 (MQ=255)
                                                                                        gTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCtt                                                          <  1:2686796/150‑1 (MQ=255)
                                                                                        gTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCtt                                                          >  1:3551743/1‑150 (MQ=255)
                                                                                            tGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCa                                                      <  2:4607795/150‑1 (MQ=255)
                                                                                              gcgcGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCAtt                                                    <  1:1644485/150‑1 (MQ=255)
                                                                                               cgcgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTc                                                   >  2:5090022/1‑150 (MQ=255)
                                                                                                 cgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAg                                                 >  2:478727/1‑150 (MQ=255)
                                                                                                   aaGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGt                                               <  1:4501442/150‑1 (MQ=255)
                                                                                                   aaGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGt                                               >  2:5150633/1‑150 (MQ=255)
                                                                                                    aGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTa                                              >  1:413536/1‑150 (MQ=255)
                                                                                                    aGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTa                                              <  1:1776906/150‑1 (MQ=255)
                                                                                                      cgcgACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACg                                            <  2:467528/150‑1 (MQ=255)
                                                                                                       gcgACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACga                                           >  2:1812104/1‑150 (MQ=255)
                                                                                                       gcgACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCACTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACga                                           <  1:262241/150‑1 (MQ=255)
                                                                                                              atttatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAgg                                    >  1:1350370/1‑150 (MQ=255)
                                                                                                                 tatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaat                                 <  1:2270426/150‑1 (MQ=255)
                                                                                                                 tatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaat                                 <  2:427923/150‑1 (MQ=255)
                                                                                                                 tatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaat                                 <  1:2358429/150‑1 (MQ=255)
                                                                                                                    ttttCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaataat                              <  2:1680497/150‑1 (MQ=255)
                                                                                                                     tttCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATg                             <  2:1441627/150‑1 (MQ=255)
                                                                                                                       tCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaa                           <  1:837364/150‑1 (MQ=255)
                                                                                                                            cgcAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaac                      >  1:4394019/1‑150 (MQ=255)
                                                                                                                               aGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaa                   >  2:3703669/1‑150 (MQ=255)
                                                                                                                                 ccATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAcc                 <  1:4946191/150‑1 (MQ=255)
                                                                                                                                  cATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAccc                >  1:528724/1‑150 (MQ=255)
                                                                                                                                   aTTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCa               >  1:115531/1‑150 (MQ=255)
                                                                                                                                    ttGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAATCACAg              >  1:4234664/1‑150 (MQ=255)
                                                                                                                                    ttGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTATTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAg              >  1:249433/1‑150 (MQ=255)
                                                                                                                                        gcCCTCCTCTTATCTCCCTCCCCGACTATCAGTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGt          >  2:452839/1‑150 (MQ=255)
                                                                                                                                          cctcctCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTg        >  2:2453374/1‑150 (MQ=255)
                                                                                                                                           ctcctcTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGa       >  1:3626122/1‑150 (MQ=255)
                                                                                                                                             cctcTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGt     >  2:611012/1‑150 (MQ=255)
                                                                                                                                                cttatCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTGCa  >  1:3474207/1‑150 (MQ=255)
                                                                                                                                                cttatCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTGCa  >  1:3733548/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTGCA  >  CP009273/3373838‑3374131

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: