Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,536,360 C→A 11.1% 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,3600CA11.1% 1025.9 / 65.7 286intergenic (+37/‑166)feoC/yhgAputative DNA‑binding transcriptional regulator/transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base C (136/118);  new base A (20/12);  total (156/130)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

TATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTCAAACAA  >  CP009273/3536212‑3536509
                                                                                                                                                    |                                                                                                                                                     
 aTGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGt                                                                                                                                                     <  1:112182/150‑1 (MQ=255)
   gggCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCt                                                                                                                                                   <  1:3993318/150‑1 (MQ=255)
    ggCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTg                                                                                                                                                  >  1:2674795/1‑150 (MQ=255)
      cAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCa                                                                                                                                                >  2:3348603/1‑150 (MQ=255)
       aaaGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCAt                                                                                                                                               >  2:2169233/1‑150 (MQ=255)
       aaaGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCAt                                                                                                                                               <  1:2532386/150‑1 (MQ=255)
       aaaGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCAt                                                                                                                                               >  1:342700/1‑150 (MQ=255)
        aaGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATc                                                                                                                                              >  1:168367/1‑150 (MQ=255)
        aaGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTAACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATc                                                                                                                                              >  1:4560002/1‑150 (MQ=255)
           ccGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCgg                                                                                                                                           <  2:2589776/150‑1 (MQ=255)
             gTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGc                                                                                                                                         >  1:3873740/1‑150 (MQ=255)
              tGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCt                                                                                                                                        <  1:844044/150‑1 (MQ=255)
                cGGAGTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCttt                                                                                                                                      <  1:4400472/150‑1 (MQ=255)
                 ggATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtttt                                                                                                                                     <  2:373431/150‑1 (MQ=255)
                 ggATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtttt                                                                                                                                     >  2:1926618/1‑150 (MQ=255)
                 ggATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtttt                                                                                                                                     >  1:3479044/1‑150 (MQ=255)
                 ggATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtttt                                                                                                                                     >  2:4008985/1‑150 (MQ=255)
                   aTTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtttttt                                                                                                                                   >  2:409867/1‑150 (MQ=255)
                   aTTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtttttt                                                                                                                                   >  1:1633656/1‑150 (MQ=255)
                   aTTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCtttttt                                                                                                                                   >  1:4174604/1‑150 (MQ=255)
                    ttCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTtg                                                                                                                                  <  2:1990338/150‑1 (MQ=255)
                    ttCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTtg                                                                                                                                  >  1:2494563/1‑150 (MQ=255)
                      cAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGGTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTtgcg                                                                                                                                >  1:987539/1‑150 (MQ=255)
                       aGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTtgcgt                                                                                                                               >  2:3874359/1‑150 (MQ=255)
                          aagaaCCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTtgcgtgcg                                                                                                                            >  1:971189/1‑150 (MQ=255)
                             aaCCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGACGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCt                                                                                                                         >  2:2345439/1‑150 (MQ=255)
                             aaCCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCt                                                                                                                         <  1:658344/150‑1 (MQ=255)
                              aCCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCtt                                                                                                                        <  2:3511341/150‑1 (MQ=255)
                               ccTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCttt                                                                                                                       <  1:3035506/150‑1 (MQ=255)
                               ccTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCttt                                                                                                                       <  1:1772601/150‑1 (MQ=255)
                                  gACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCa                                                                                                                    >  1:525708/1‑150 (MQ=255)
                                  gACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCa                                                                                                                    >  1:3231781/1‑150 (MQ=255)
                                    cGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCGTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTATTTTGCGTGCGGCTTTCCATa                                                                                                                  <  2:3623140/150‑1 (MQ=255)
                                     ggCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATaa                                                                                                                 >  2:3669035/1‑150 (MQ=255)
                                      gCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATaaa                                                                                                                >  2:1822970/1‑150 (MQ=255)
                                        tGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATaaaaa                                                                                                              >  2:2728880/1‑150 (MQ=255)
                                        tGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATaaaaa                                                                                                              >  2:2728798/1‑150 (MQ=255)
                                             ctctGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCaa                                                                                                         >  2:3824809/1‑150 (MQ=255)
                                              tctGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAAc                                                                                                        >  1:4637292/1‑150 (MQ=255)
                                                  gCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCtt                                                                                                    <  1:1905628/150‑1 (MQ=255)
                                                   cAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTg                                                                                                   <  2:764516/150‑1 (MQ=255)
                                                   cAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTg                                                                                                   >  1:3302682/1‑150 (MQ=255)
                                                    aGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgc                                                                                                  <  1:680165/150‑1 (MQ=255)
                                                    aGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgc                                                                                                  >  2:1133564/1‑150 (MQ=255)
                                                     gTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgca                                                                                                 <  1:1932820/150‑1 (MQ=255)
                                                      ttGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgcag                                                                                                >  2:4386064/1‑150 (MQ=255)
                                                         tAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCAc                                                                                             <  2:4574778/150‑1 (MQ=255)
                                                          aaaaGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACg                                                                                            >  1:2088268/1‑150 (MQ=255)
                                                           aaaGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACgg                                                                                           >  1:861472/1‑150 (MQ=255)
                                                             aGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCg                                                                                         <  1:402294/150‑1 (MQ=255)
                                                                      gaaggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCt                                                                                <  1:3729851/150‑1 (MQ=255)
                                                                      gaaggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCt                                                                                <  1:1324411/150‑1 (MQ=255)
                                                                       aaggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCGGCATCCGGGCTTGTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGGTCCtt                                                                               >  2:3771675/1‑150 (MQ=255)
                                                                        aggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCttt                                                                              >  2:3080909/1‑150 (MQ=255)
                                                                        aggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCATAAGTTCCttt                                                                              <  2:4153555/150‑1 (MQ=255)
                                                                         ggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTg                                                                             >  2:481286/1‑150 (MQ=255)
                                                                         ggaaAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCATTg                                                                             <  1:2006749/150‑1 (MQ=255)
                                                                           aaaaGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGaa                                                                           <  2:3165161/150‑1 (MQ=255)
                                                                             aaGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAg                                                                         <  1:250673/150‑1 (MQ=255)
                                                                             aaGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAg                                                                         >  2:3593896/1‑150 (MQ=255)
                                                                             aaGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAg                                                                         <  1:3112628/150‑1 (MQ=255)
                                                                               gCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCa                                                                       >  2:4459425/1‑150 (MQ=255)
                                                                                ccTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCAt                                                                      <  1:2741168/150‑1 (MQ=255)
                                                                                ccTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCAt                                                                      <  2:2190692/150‑1 (MQ=255)
                                                                                    tatGCGCGTGTGGTGGGCGGTGCGTTAAACGTACTCCATCGCCTGTTTTGGTAAGGCCGGTATGCGTCTGCTTCCGGGCTTTTTTTCTTGCGGCTTTCAATAAAAATGCATCTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCg                                                                  <  1:3088861/148‑1 (MQ=255)
                                                                                     cTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGc                                                                 >  1:1051025/1‑150 (MQ=255)
                                                                                       gcgcGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAg                                                               >  2:2409024/1‑150 (MQ=255)
                                                                                           gAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGAt                                                           >  1:2398910/1‑150 (MQ=255)
                                                                                            aGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATg                                                          >  1:1961202/1‑150 (MQ=255)
                                                                                              tggtggGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGaa                                                        <  2:1686665/150‑1 (MQ=255)
                                                                                               ggtggGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGaaa                                                       <  2:2835629/150‑1 (MQ=255)
                                                                                               ggtggGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGaaa                                                       <  2:3889779/150‑1 (MQ=255)
                                                                                                gtggGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGaaaa                                                      <  2:4463124/150‑1 (MQ=255)
                                                                                                   ggCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTc                                                   >  1:2602764/1‑150 (MQ=255)
                                                                                                      gcTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCt                                                >  2:3585157/1‑150 (MQ=255)
                                                                                                          cGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATa                                            >  2:3910118/1‑150 (MQ=255)
                                                                                                            ttAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATaca                                          >  1:1190913/1‑150 (MQ=255)
                                                                                                            ttAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATaca                                          >  1:2801590/1‑150 (MQ=255)
                                                                                                             tAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGAGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATacac                                         >  2:1571345/1‑150 (MQ=255)
                                                                                                              aaCCTTACTCCATCGCCTGTTTTGGAAAGCCCGNTATGCGTCTNCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATAcaca                                        <  1:257000/150‑1 (MQ=255)
                                                                                                              aaCCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTACCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATAcaca                                        >  2:1048021/1‑150 (MQ=255)
                                                                                                               aCCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAg                                       >  2:3873966/1‑150 (MQ=255)
                                                                                                                ccTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAgg                                      <  2:1457006/150‑1 (MQ=255)
                                                                                                                 cTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGgt                                     <  1:100199/150‑1 (MQ=255)
                                                                                                                     cTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTgg                                 >  1:1786978/1‑150 (MQ=255)
                                                                                                                         aTCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTg                             <  2:3395906/150‑1 (MQ=255)
                                                                                                                         aTCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTg                             <  2:3202788/150‑1 (MQ=255)
                                                                                                                           cGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGgc                           <  1:3808709/150‑1 (MQ=255)
                                                                                                                            gCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGgcg                          <  2:857997/150‑1 (MQ=255)
                                                                                                                             ccTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGgcgc                         >  1:610964/1‑150 (MQ=255)
                                                                                                                              cTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGcgcg                        <  2:4315464/150‑1 (MQ=255)
                                                                                                                               tGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGt                       >  2:2023798/1‑150 (MQ=255)
                                                                                                                                gTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTa                      >  2:2108626/1‑150 (MQ=255)
                                                                                                                                gTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTa                      <  1:3779330/150‑1 (MQ=255)
                                                                                                                                    tGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGt                  <  2:1219552/150‑1 (MQ=255)
                                                                                                                                       aaaGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTcgc               >  1:4296336/1‑150 (MQ=255)
                                                                                                                                         aGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCgg             <  1:2849997/150‑1 (MQ=255)
                                                                                                                                         aGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCgg             <  2:927044/150‑1 (MQ=255)
                                                                                                                                           cccGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCGGCa           >  1:2329567/1‑150 (MQ=255)
                                                                                                                                              ggTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTc        >  1:2793684/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAGTGGCGCGTAGAGTCGCGGCATTCAAACAA  >  CP009273/3536212‑3536509

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: