Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,935,100 G→T 100% intergenic (‑413/‑68) yieP ← / → rrsC putative transcriptional regulator/16S ribosomal RNA of rrnC operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,935,1000GT98.8% 402.3 / ‑1.4 169intergenic (‑413/‑68)yieP/rrsCputative transcriptional regulator/16S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (61/106);  total (61/108)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 169 total aligned reads displayed.

AAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAA  >  CP009273/3934953‑3935162
                                                                                                                                                   |                                                              
aaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTaa                                                              <  1:5073241/150‑1 (MQ=21)
aaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTaa                                                              >  2:2208176/1‑150 (MQ=21)
aaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTaa                                                              <  1:1995487/150‑1 (MQ=21)
aaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTaa                                                              <  2:5015182/150‑1 (MQ=21)
 aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAAt                                                             <  1:1157862/150‑1 (MQ=21)
 aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAAt                                                             >  1:3227584/1‑150 (MQ=21)
  atGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATa                                                            >  1:852555/1‑150 (MQ=21)
   tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATAc                                                           <  2:1176155/150‑1 (MQ=21)
   tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATAc                                                           <  2:1156269/150‑1 (MQ=21)
    gCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATAcc                                                          <  1:3661909/150‑1 (MQ=21)
     cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCa                                                         >  2:5141378/1‑150 (MQ=21)
     cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCa                                                         >  1:4488722/1‑150 (MQ=21)
     cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCa                                                         >  1:4999956/1‑150 (MQ=21)
       tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAACAATTAATACCAAg                                                       >  1:3652413/1‑150 (MQ=21)
       tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAg                                                       <  2:2114716/150‑1 (MQ=21)
       tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAg                                                       <  2:3659114/150‑1 (MQ=21)
       tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAg                                                       <  1:1007420/150‑1 (MQ=21)
        gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGt                                                      <  1:3868797/150‑1 (MQ=21)
        gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGt                                                      <  1:2650497/150‑1 (MQ=21)
         aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGtc                                                     >  2:3352422/1‑150 (MQ=21)
         aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGtc                                                     >  2:1102496/1‑150 (MQ=21)
         aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGtc                                                     >  1:1918781/1‑150 (MQ=21)
             tGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCaa                                                 >  1:548487/1‑150 (MQ=21)
              gTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAag                                                >  2:110704/1‑150 (MQ=21)
               tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAaga                                               >  1:865570/1‑150 (MQ=21)
               tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAaga                                               >  2:801797/1‑150 (MQ=21)
               tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAaga                                               >  2:4755532/1‑150 (MQ=21)
                aGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAagag                                              <  2:4300846/150‑1 (MQ=21)
                aGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAagag                                              >  1:4517719/1‑150 (MQ=21)
                 gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGt                                             <  1:1842290/150‑1 (MQ=21)
                 gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGt                                             <  1:2590541/150‑1 (MQ=21)
                 gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGt                                             <  1:4801374/150‑1 (MQ=21)
                 gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGt                                             <  2:566987/150‑1 (MQ=21)
                   gggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                           >  2:4888065/1‑150 (MQ=21)
                   gggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                           <  1:2965019/150‑1 (MQ=21)
                   gggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                           >  1:5024820/1‑150 (MQ=21)
                    ggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCTAGAGTGaa                                          >  2:4841452/1‑150 (MQ=21)
                    ggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGaa                                          <  2:1041781/150‑1 (MQ=21)
                     gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAac                                         <  2:1112524/150‑1 (MQ=21)
                     gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAac                                         >  1:4949497/1‑150 (MQ=21)
                     gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAac                                         <  2:4785428/150‑1 (MQ=21)
                       aaGCGTATTATGCACACCCCGCGCCGCTGANAAAAANCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTNTCAAGAGTGAacac                                       <  2:2751328/150‑1 (MQ=18)
                       aaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAacac                                       >  1:492801/1‑150 (MQ=18)
                       aaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAacac                                       >  1:3735581/1‑150 (MQ=18)
                       aaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAacac                                       >  2:2976251/1‑150 (MQ=18)
                       aaGCGTATTATGCACACCCCGCGCCGATGAGAAAAAGAAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAacac                                       <  2:1526236/150‑1 (MQ=16)
                        aGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACg                                      <  2:848181/150‑1 (MQ=17)
                         gCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGt                                     <  2:1618065/150‑1 (MQ=17)
                          cGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATATGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTa                                    <  2:2653942/150‑1 (MQ=17)
                            tattatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAt                                  <  1:91290/150‑1 (MQ=17)
                            tattatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAt                                  <  2:4344972/150‑1 (MQ=17)
                            tattatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAt                                  <  2:3855903/150‑1 (MQ=17)
                             attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAtt                                 <  2:752755/150‑1 (MQ=17)
                             attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAtt                                 <  2:752698/150‑1 (MQ=17)
                             attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAtt                                 <  2:2868375/150‑1 (MQ=17)
                             attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAtt                                 <  2:4233532/150‑1 (MQ=17)
                              ttatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTc                                <  2:248212/150‑1 (MQ=17)
                              ttatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTc                                >  2:2111379/1‑150 (MQ=17)
                              ttatGCACACCCCGCGCCGATGATAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTc                                <  1:188381/150‑1 (MQ=12)
                                 tGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                             <  1:1758519/150‑1 (MQ=17)
                                 tGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                             <  2:4237262/150‑1 (MQ=17)
                                  gCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa                            >  1:863280/1‑150 (MQ=17)
                                  gCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa                            <  2:4898568/150‑1 (MQ=17)
                                   cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTAc                           <  2:4299001/150‑1 (MQ=17)
                                   cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTAc                           <  1:671042/150‑1 (MQ=17)
                                     cacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGa                         >  2:207961/1‑150 (MQ=17)
                                      acCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                        <  1:3107870/150‑1 (MQ=17)
                                      acCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                        <  1:118360/150‑1 (MQ=17)
                                       ccccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAg                       <  1:3152603/150‑1 (MQ=17)
                                        cccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt                      <  1:4008442/150‑1 (MQ=17)
                                         ccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGtt                     <  2:5307367/150‑1 (MQ=17)
                                         ccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGtt                     <  1:657519/150‑1 (MQ=32)
                                          cgcgCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt                    <  2:1040469/150‑1 (MQ=17)
                                          cgcgCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt                    <  2:3286227/150‑1 (MQ=17)
                                          cgcgCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt                    <  2:1329191/150‑1 (MQ=17)
                                           gcgcCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTa                   <  1:3433302/150‑1 (MQ=17)
                                             gccgcTTAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAAt                 <  2:2534867/150‑1 (MQ=17)
                                             gccgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAAt                 >  1:4306602/1‑150 (MQ=17)
                                              ccgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAAtt                <  1:1410684/150‑1 (MQ=17)
                                               cgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTc               <  2:2333320/150‑1 (MQ=17)
                                                gcTGAGAAAAACCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCt              <  1:61797/150‑1 (MQ=17)
                                                 cTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAg                       <  1:219950/140‑1 (MQ=17)
                                                  tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt            <  1:3110972/150‑1 (MQ=17)
                                                  tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt            >  1:4573279/1‑150 (MQ=17)
                                                  tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt            <  2:4617889/150‑1 (MQ=17)
                                                   gagaAAAAGCAACGCGGCACTGCTCTTTAACAAGTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg           <  2:2156295/150‑1 (MQ=12)
                                                   gagaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg           <  1:3466412/150‑1 (MQ=17)
                                                   gagaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg           >  1:2390972/1‑150 (MQ=17)
                                                     gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTc               >  1:4299505/1‑144 (MQ=17)
                                                     gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAg         <  1:5312975/150‑1 (MQ=17)
                                                       aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCa       <  1:2825610/150‑1 (MQ=17)
                                                       aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCa       >  1:4648444/1‑150 (MQ=17)
                                                         aagcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGACCATc     <  1:1620813/150‑1 (MQ=17)
                                                          agcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCa    >  1:4180051/1‑150 (MQ=17)
                                                          agcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCa    >  2:2430898/1‑150 (MQ=17)
                                                          agcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGACCATCa    >  2:2551017/1‑150 (MQ=17)
                                                            caaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCaaa  >  1:2950145/1‑150 (MQ=14)
                                                            caaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCaaa  >  1:5225730/1‑150 (MQ=14)
                                                            caaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCaaa  <  2:4365612/150‑1 (MQ=14)
                                                            caaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCaaa  <  2:1853544/150‑1 (MQ=14)
                                                                                                                                                   |                                                              
AAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAA  >  CP009273/3934953‑3935162

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: