Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,536,360 C→A 15.0% intergenic (+37/‑166) feoC → / → yhgA putative DNA‑binding transcriptional regulator/transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,536,3600CA15.0% 816.0 / 84.7 247intergenic (+37/‑166)feoC/yhgAputative DNA‑binding transcriptional regulator/transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base C (98/112);  new base A (17/20);  total (115/132)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 254 total aligned reads displayed.

GTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTCAAACA  >  CP009273/3536211‑3536508
                                                                                                                                                     |                                                                                                                                                    
gTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGc                                                                                                                                                      <  1:2663712/150‑1 (MQ=255)
gTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGc                                                                                                                                                      <  1:2285048/150‑1 (MQ=255)
gTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGc                                                                                                                                                      <  1:1008025/150‑1 (MQ=255)
   tGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTc                                                                                                                                                   >  2:901281/1‑150 (MQ=255)
    gggCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCt                                                                                                                                                  <  1:423753/150‑1 (MQ=255)
    gggCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGAATGAGTCt                                                                                                                                                  <  2:4090694/150‑1 (MQ=255)
     ggCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTg                                                                                                                                                 <  1:2683544/150‑1 (MQ=255)
     ggCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTg                                                                                                                                                 <  2:2162503/150‑1 (MQ=255)
       cAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCa                                                                                                                                               >  1:2414141/1‑150 (MQ=255)
        aaaGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCAt                                                                                                                                              >  2:2256918/1‑150 (MQ=255)
        aaaGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGAAAAAGCTGCCCGGAAGGAAAAGCCTGTCAGCGCGAGTGGTGGGCGCTGCGAGAACCTTACTCCATCGCCAGGTTTGGAAAGCCCGGGATGCGTCTGCAt                                                                                                                                              >  1:4258880/1‑150 (MQ=255)
             cGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTGGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCggg                                                                                                                                         <  2:2046326/150‑1 (MQ=255)
              gTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGc                                                                                                                                        >  1:4183608/1‑150 (MQ=255)
                gCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtt                                                                                                                                      <  1:2194449/150‑1 (MQ=255)
                 cGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCttt                                                                                                                                     <  1:4569999/150‑1 (MQ=255)
                  ggATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtttt                                                                                                                                    <  1:3958613/150‑1 (MQ=255)
                  ggATACAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCCTCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCtttt                                                                                                                                    >  2:4545177/1‑150 (MQ=255)
                   gATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCttttt                                                                                                                                   <  2:1741592/150‑1 (MQ=255)
                   gATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTGGGAAAGCCCGGTATGCGTCTGCATCCGGGCttttt                                                                                                                                   >  2:1041380/1‑150 (MQ=255)
                    aTTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtttttt                                                                                                                                  <  1:4533034/150‑1 (MQ=255)
                       cAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTtgcg                                                                                                                               >  1:1736476/1‑150 (MQ=255)
                        aGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTtgcgt                                                                                                                              >  2:65465/1‑150 (MQ=255)
                         ggAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCCTCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTtgcgtg                                                                                                                             <  2:1480401/150‑1 (MQ=255)
                             gaaCCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGc                                                                                                                         <  2:2534444/150‑1 (MQ=255)
                              aaCCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCt                                                                                                                        <  1:3190093/150‑1 (MQ=255)
                                   gACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCa                                                                                                                   <  1:2936838/150‑1 (MQ=255)
                                    aCGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCAt                                                                                                                  >  1:2743508/1‑150 (MQ=255)
                                    aCGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCAt                                                                                                                  <  1:3375237/150‑1 (MQ=255)
                                     cGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATa                                                                                                                 <  2:324499/150‑1 (MQ=255)
                                           ccTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATg                                                                                                           <  2:4292754/150‑1 (MQ=255)
                                               tctGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAAc                                                                                                       <  1:2348334/150‑1 (MQ=255)
                                                ctGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAAct                                                                                                      >  2:1351318/1‑150 (MQ=255)
                                                   gCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCtt                                                                                                   <  1:3849511/150‑1 (MQ=255)
                                                    cAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTg                                                                                                  >  2:4455968/1‑150 (MQ=255)
                                                     aGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgc                                                                                                 >  1:4454652/1‑150 (MQ=255)
                                                      gTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgca                                                                                                <  2:299421/150‑1 (MQ=255)
                                                      gTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGCTGGGCGCTGCGTTAACATTACTCCATCGCCTGGTTTGGAAAGCCCGGAATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTGgca                                                                                                >  2:3828293/1‑150 (MQ=255)
                                                       ttGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgcag                                                                                               <  1:2002827/150‑1 (MQ=255)
                                                         gTAACAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgcagca                                                                                             <  2:3584298/150‑1 (MQ=255)
                                                           aaaaGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACg                                                                                           >  1:2103084/1‑150 (MQ=255)
                                                            aaaGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACgg                                                                                          >  2:3204299/1‑150 (MQ=255)
                                                            aaaGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACgg                                                                                          <  2:2852700/150‑1 (MQ=255)
                                                             aaGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGc                                                                                         >  1:1267777/1‑150 (MQ=255)
                                                              aGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCg                                                                                        <  1:2617222/150‑1 (MQ=255)
                                                                 tGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTaa                                                                                     <  1:3126317/150‑1 (MQ=255)
                                                                 tGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTaa                                                                                     <  1:1775958/150‑1 (MQ=255)
                                                                  gCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAg                                                                                    <  2:1814886/150‑1 (MQ=255)
                                                                   cccGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAACAATGCAACTCTTGCAGCACGGCGTAAGt                                                                                   <  2:1713180/150‑1 (MQ=255)
                                                                    ccGGAAGGAAAATCCTGTCTGCGCGAGTGGTGGGCGCTGAGTTAAACTTACTCCATCGCCTGTTTTGTAAAGCCCGGTATGCGTCTGCATACGGGCATTTTTGCGTGCGGCTTTACATAAAAATGCAACTCTTGCAGCACGGCGTAAGtt                                                                                  <  2:1180001/150‑1 (MQ=255)
                                                                       gaaggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCt                                                                               >  1:1649978/1‑150 (MQ=255)
                                                                       gaaggaaAAGCCTGTCTGCGCGAGTGGAGGGCGCATCGTTAACCTTACTCCATCGCCTCTTTAGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGAGGCATTACATAAAAATGAAACACTTGAAGAACGGCGTAAGTAcca                                                                               >  1:1650195/1‑149 (MQ=255)
                                                                        aaggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCtt                                                                              <  2:3657358/150‑1 (MQ=255)
                                                                        aaggaaAAGACTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTATGCATCCGGGCTTTTTTGAGTGCGTCTTTCCATAATATTGCAACTCTTGCAGCACGGCGTAAGTTCCtt                                                                              <  1:3855569/150‑1 (MQ=255)
                                                                         aggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCttt                                                                             >  1:4540318/1‑150 (MQ=255)
                                                                            aaaaGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGTTAAGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTTaa                                                                          >  1:371146/1‑150 (MQ=255)
                                                                              aaGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAg                                                                        >  1:4351840/1‑150 (MQ=255)
                                                                              aaGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAg                                                                        <  2:1520796/150‑1 (MQ=255)
                                                                                   tGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCAtct                                                                   >  2:1587448/1‑150 (MQ=255)
                                                                                    gTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCAGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTTAAAGCAtctc                                                                  >  1:152247/1‑150 (MQ=255)
                                                                                     tCTGCGCGTGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCg                                                                 <  1:160757/150‑1 (MQ=255)
                                                                                     tCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCATTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCg                                                                 <  2:4559012/150‑1 (MQ=255)
                                                                                      cTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGc                                                                >  1:498174/1‑150 (MQ=255)
                                                                                       tGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCa                                                               >  1:1909441/1‑150 (MQ=255)
                                                                                            gAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGAt                                                          >  1:1865656/1‑150 (MQ=255)
                                                                                             aGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATg                                                         <  1:193796/150‑1 (MQ=255)
                                                                                             aGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATg                                                         <  2:4127577/150‑1 (MQ=255)
                                                                                              gtggtgGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGa                                                        <  2:1140704/150‑1 (MQ=255)
                                                                                               tggtggGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGaa                                                       >  1:2735796/1‑150 (MQ=255)
                                                                                               tggtggGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGaa                                                       >  1:3570951/1‑150 (MQ=255)
                                                                                                  tggGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAAc                                                    >  2:3542442/1‑150 (MQ=255)
                                                                                                   gggCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACt                                                   >  1:810666/1‑150 (MQ=255)
                                                                                                    ggCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTc                                                  <  1:2538213/150‑1 (MQ=255)
                                                                                                    ggCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTc                                                  <  1:733207/150‑1 (MQ=255)
                                                                                                       gcTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCt                                               >  2:396197/1‑150 (MQ=255)
                                                                                                       gcTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCt                                               <  2:2367667/150‑1 (MQ=255)
                                                                                                        cTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTa                                              <  2:1611119/150‑1 (MQ=255)
                                                                                                         tGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTaa                                             <  2:4343179/150‑1 (MQ=255)
                                                                                                          gCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAAt                                            <  2:1761790/150‑1 (MQ=255)
                                                                                                           cGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATa                                           >  2:3290912/1‑150 (MQ=255)
                                                                                                            gTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATac                                          <  1:2855014/150‑1 (MQ=255)
                                                                                                               aaCCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATAcaca                                       >  1:1684517/1‑150 (MQ=255)
                                                                                                                   ttACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTATTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGgtg                                   <  1:1167025/150‑1 (MQ=255)
                                                                                                                      cTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTgg                                >  2:2937949/1‑150 (MQ=255)
                                                                                                                         cATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGt                             >  2:3352318/1‑150 (MQ=255)
                                                                                                                         cATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGt                             >  2:766812/1‑150 (MQ=255)
                                                                                                                          aTCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTg                            <  1:3700274/150‑1 (MQ=255)
                                                                                                                            cGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGgc                          <  2:3804446/150‑1 (MQ=255)
                                                                                                                               cTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGAGGcgcg                       >  2:2672321/1‑150 (MQ=255)
                                                                                                                                tGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGt                      <  1:3581114/150‑1 (MQ=255)
                                                                                                                                 gTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTa                     <  2:2712708/150‑1 (MQ=255)
                                                                                                                                      ggAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTc                <  1:3468857/150‑1 (MQ=255)
                                                                                                                                        aaaGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATCCACAGGTGTGGAGTGGCGCGTAGAGTcgc              >  2:3725622/1‑150 (MQ=255)
                                                                                                                                        aaaGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTcgc              <  1:3983309/150‑1 (MQ=255)
                                                                                                                                         aaGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTcgcg             >  2:4366019/1‑150 (MQ=255)
                                                                                                                                          aGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCgg            <  1:853199/150‑1 (MQ=255)
                                                                                                                                               ggTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACCCAGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTc       <  1:2256283/150‑1 (MQ=255)
                                                                                                                                                gTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTCa      <  1:932831/150‑1 (MQ=255)
                                                                                                                                                gTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCGCGCAGGGATGAAACCTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTcc      >  2:1370395/1‑149 (MQ=255)
                                                                                                                                                    gCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACCGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTCAaaca  <  2:3182542/150‑1 (MQ=255)
                                                                                                                                                    gCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTCAaaca  <  2:3182509/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTCAAACA  >  CP009273/3536211‑3536508

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: