Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 279,334 C→A 23.3% I141I (ATC→ATA)  yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273279,3340CA23.3% 759.9 / 163.0 290I141I (ATC→ATA) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base C (115/107);  new base A (29/39);  total (144/146)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 296 total aligned reads displayed.

ATCGGCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCT  >  CP009273/279185‑279480
                                                                                                                                                     |                                                                                                                                                  
 tCGGCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCg                                                                                                                                                   <  1:302108/150‑1 (MQ=255)
  cGGCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCgt                                                                                                                                                  >  1:1941073/1‑150 (MQ=255)
     cATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGtt                                                                                                                                               >  2:3111329/1‑150 (MQ=255)
     cATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGtt                                                                                                                                               >  2:3733488/1‑150 (MQ=255)
      atgcagatgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGttt                                                                                                                                              <  1:3283113/150‑1 (MQ=255)
         cagatgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTcgc                                                                                                                                           >  2:2762137/1‑150 (MQ=255)
          agatgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTGTcgcc                                                                                                                                          >  2:3846470/1‑150 (MQ=255)
             tgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCgccgcc                                                                                                                                       <  2:2528573/150‑1 (MQ=255)
               cagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTg                                                                                                                                     >  2:3343164/1‑150 (MQ=255)
               cagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTg                                                                                                                                     >  1:2033546/1‑150 (MQ=255)
                 ggcggcggcGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGAt                                                                                                                                   >  1:702544/1‑150 (MQ=255)
                  gcggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATc                                                                                                                                  >  2:1804600/1‑150 (MQ=255)
                  gcggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATc                                                                                                                                  <  2:3908323/150‑1 (MQ=255)
                  gcggcggcgAAGGAGATCACCCGCAAGGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATc                                                                                                                                  <  2:3908406/150‑1 (MQ=255)
                   cggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATcc                                                                                                                                 <  2:3166407/150‑1 (MQ=255)
                    ggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCg                                                                                                                                <  2:2893030/150‑1 (MQ=255)
                    ggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCg                                                                                                                                <  1:3897123/150‑1 (MQ=255)
                      cggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCt                                                                                                                              >  2:2330647/1‑150 (MQ=255)
                         cgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTccc                                                                                                                           <  1:1650448/150‑1 (MQ=255)
                         cgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATACCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTccc                                                                                                                           >  2:4502463/1‑150 (MQ=255)
                          gAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCt                                                                                                                          >  1:1299999/1‑150 (MQ=255)
                          gAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCt                                                                                                                          >  1:1663169/1‑150 (MQ=255)
                           aaGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTg                                                                                                                         >  1:3685135/1‑150 (MQ=255)
                                gaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGAc                                                                                                                    >  2:4551167/1‑150 (MQ=255)
                                gaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGAc                                                                                                                    <  2:1995415/150‑1 (MQ=255)
                                gaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGAc                                                                                                                    <  1:3568725/150‑1 (MQ=255)
                                   cACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGAcgcg                                                                                                                 >  1:49794/1‑150 (MQ=255)
                                     cccGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcggc                                                                                                               >  2:2540590/1‑150 (MQ=255)
                                      ccGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcggcg                                                                                                              <  2:3478588/150‑1 (MQ=255)
                                      ccGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcggcg                                                                                                              <  2:4054539/150‑1 (MQ=255)
                                         cAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCggcgggc                                                                                                           <  2:207634/150‑1 (MQ=255)
                                         cAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTAGCCGCCTGATCCGCTCCCTGCCGACGCggcgggc                                                                                                           <  1:569654/150‑1 (MQ=255)
                                          aaTGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCggcgggcg                                                                                                          >  2:2519063/1‑150 (MQ=255)
                                               cGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGAt                                                                                                     >  2:2183614/1‑150 (MQ=255)
                                                gggATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATc                                                                                                    <  2:4360732/150‑1 (MQ=255)
                                                    tCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCg                                                                                                >  2:1914986/1‑150 (MQ=255)
                                                      ccGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTa                                                                                              <  2:849140/150‑1 (MQ=255)
                                                        gTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGc                                                                                            >  1:1523803/1‑150 (MQ=255)
                                                        gTTCGCGGCCTTCGTCAGCGATCCGTACGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGc                                                                                            >  2:3610914/1‑150 (MQ=255)
                                                         ttCGCGGCCTTCGTCAGCTATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCTGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAgag                                                                                           <  2:179801/150‑3 (MQ=255)
                                                         ttCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCg                                                                                           <  2:2492284/150‑1 (MQ=255)
                                                         ttCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCg                                                                                           <  2:4252320/150‑1 (MQ=255)
                                                              ggCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTg                                                                                      >  1:4494069/1‑150 (MQ=255)
                                                                 cTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGGATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGAGATTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGa                                                                                   >  2:3949107/1‑150 (MQ=255)
                                                                   tCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATa                                                                                 <  2:4076869/150‑1 (MQ=255)
                                                                    cGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATaa                                                                                >  2:258064/1‑150 (MQ=255)
                                                                    cGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATaa                                                                                <  1:3441544/150‑1 (MQ=255)
                                                                       cAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAgg                                                                             <  1:3400027/150‑1 (MQ=255)
                                                                        aGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAggg                                                                            >  2:3076633/1‑150 (MQ=255)
                                                                        aGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAggg                                                                            >  2:3596688/1‑150 (MQ=255)
                                                                         gCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGc                                                                           <  1:4304511/150‑1 (MQ=255)
                                                                          cGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCt                                                                          >  2:1647152/1‑150 (MQ=255)
                                                                          cGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGGTTCGCCGCCGGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCAGCGATAAAGGGCt                                                                          >  2:878210/1‑150 (MQ=255)
                                                                           gATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGAGAAAGGGCTg                                                                         >  2:1221053/1‑150 (MQ=255)
                                                                            aTCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGc                                                                        <  2:4342013/150‑1 (MQ=255)
                                                                               cGTGCGACGGGCGCTCGCAGGGCACGCACGGTCTGTTCGATTCCCTGCCGTACCGCAACGACGCGGAGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCg                                                                     <  2:191685/150‑1 (MQ=255)
                                                                                     aCGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGccacca                                                               >  2:577113/1‑150 (MQ=255)
                                                                                      cGGGCGCTCTCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCAt                                                              <  2:1011418/150‑1 (MQ=255)
                                                                                       gggCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATg                                                             >  2:3723610/1‑150 (MQ=255)
                                                                                       gggCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATg                                                             >  1:1009525/1‑150 (MQ=255)
                                                                                           gcTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTg                                                         >  1:3224731/1‑150 (MQ=255)
                                                                                           gcTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCAGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTg                                                         >  1:3224726/1‑150 (MQ=255)
                                                                                             tCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTgcg                                                       <  1:1395207/150‑1 (MQ=255)
                                                                                               gCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCt                                                     >  1:1721746/1‑150 (MQ=255)
                                                                                                cAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTg                                                    >  1:3988558/1‑150 (MQ=255)
                                                                                                  gggCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGc                                                  <  2:698962/150‑1 (MQ=255)
                                                                                                    gcacgcacGGTATGTTCGATTCCCTGCCTTACAGCAAAGACGCGGCGATAGTGATTCGCCGCCTGATCCGCTCCCTGCAGAAGCGGCGGGCGGTGATCGGCGTAGCGACCTGAGATAAATGGCTGCCCGCCAACATTATTGCGCTGGcct                                                >  1:825444/1‑149 (MQ=255)
                                                                                                    gcacgcacGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcg                                                >  1:1787295/1‑150 (MQ=255)
                                                                                                     cacgcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcgc                                               >  1:124562/1‑150 (MQ=255)
                                                                                                     cacgcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcgc                                               <  1:254463/150‑1 (MQ=255)
                                                                                                       cgcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGa                                             <  1:2829851/150‑1 (MQ=255)
                                                                                                        gcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGAt                                            <  2:4161175/150‑1 (MQ=255)
                                                                                                        gcacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGAt                                            <  2:645441/150‑1 (MQ=255)
                                                                                                         cacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATg                                           >  1:2084872/1‑150 (MQ=255)
                                                                                                              tATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGa                                      >  1:2698724/1‑150 (MQ=255)
                                                                                                               aTGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGAc                                     >  1:2481010/1‑150 (MQ=255)
                                                                                                                 gTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCt                                   <  1:1070330/150‑1 (MQ=255)
                                                                                                                   tCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGc                                 >  2:1821232/1‑150 (MQ=255)
                                                                                                                     gATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCg                               >  2:3946010/1‑150 (MQ=255)
                                                                                                                          ccTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTAt                          >  1:2276395/1‑150 (MQ=255)
                                                                                                                                gTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGt                    <  2:4040020/150‑1 (MQ=255)
                                                                                                                                gTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACATGCGATAAAGGGCTGCCCGCCACCATGATTGCGCGGGACGCGATGCACGACCTGCCGACTATTCTGGt                    <  1:4286439/150‑1 (MQ=255)
                                                                                                                                gTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCTTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGt                    <  2:1046663/150‑1 (MQ=255)
                                                                                                                                 tACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTg                   <  1:96866/150‑1 (MQ=255)
                                                                                                                                 tACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTg                   >  2:2467776/1‑150 (MQ=255)
                                                                                                                                  aCCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGc                  <  2:134459/150‑1 (MQ=255)
                                                                                                                                  aCCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGc                  >  2:796387/1‑150 (MQ=255)
                                                                                                                                  aCCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGc                  >  1:2729634/1‑150 (MQ=255)
                                                                                                                                     gCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgg               <  2:3013073/150‑1 (MQ=255)
                                                                                                                                      cAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcggg              >  1:437434/1‑150 (MQ=255)
                                                                                                                                       aaCGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATCAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggc             >  2:40950/1‑150 (MQ=255)
                                                                                                                                       aaCGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggc             >  1:1323006/1‑150 (MQ=255)
                                                                                                                                       aaCGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggc             >  2:3320281/1‑150 (MQ=255)
                                                                                                                                       aaCGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCAGGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggc             >  1:2317085/1‑150 (MQ=255)
                                                                                                                                             gcggcgATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCgggcggggcg       <  1:1112010/150‑1 (MQ=255)
                                                                                                                                              cggcgATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGa      >  1:25307/1‑150 (MQ=255)
                                                                                                                                              cggcgATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGa      >  2:1500801/1‑150 (MQ=255)
                                                                                                                                               ggcgATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGAc     >  2:2289734/1‑150 (MQ=255)
                                                                                                                                               ggcgATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGAc     <  1:2552116/150‑1 (MQ=255)
                                                                                                                                                  gATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCt  <  2:3359563/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATCGGCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCT  >  CP009273/279185‑279480

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: ‑