Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,909,342 T→A 100% E150D (GAA→GAT proQ ← RNA chaperone, putative ProP translation regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,909,3420TA95.5% 1064.7 / 26.1 309E150D (GAA→GATproQRNA chaperone, putative ProP translation regulator
Reads supporting (aligned to +/- strand):  ref base T (5/9);  new base A (170/125);  total (175/134)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 312 total aligned reads displayed.

TTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTG  >  CP009273/1909197‑1909488
                                                                                                                                                 |                                                                                                                                                  
tttGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCcatctt                                                                                                                                                <  2:3081155/150‑1 (MQ=255)
    aCCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTgcg                                                                                                                                            <  1:1030219/150‑1 (MQ=255)
    aCCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTgcg                                                                                                                                            >  2:2694268/1‑150 (MQ=255)
     cccGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTgcgg                                                                                                                                           <  1:2501573/150‑1 (MQ=255)
     cccGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTgcgg                                                                                                                                           <  1:2441532/150‑1 (MQ=255)
     cccGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTgcgg                                                                                                                                           >  1:2381309/1‑150 (MQ=255)
      ccGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTgcggc                                                                                                                                          <  2:2521668/150‑1 (MQ=255)
      ccGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTgcggc                                                                                                                                          >  2:3931399/1‑150 (MQ=255)
         cTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCggcgtg                                                                                                                                       <  1:4715323/150‑1 (MQ=255)
            tCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCggcgtggcg                                                                                                                                    >  1:2234898/1‑150 (MQ=255)
            tCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCggcgtggcg                                                                                                                                    <  1:2587612/150‑1 (MQ=255)
             cACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCggcgtggcgt                                                                                                                                   >  1:803584/1‑150 (MQ=255)
              aCCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTa                                                                                                                                  >  2:2445687/1‑150 (MQ=255)
               ccTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAg                                                                                                                                 >  2:3728655/1‑150 (MQ=255)
               ccTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTGTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCggcgtggcgtcg                                                                                                                                 >  2:4663650/1‑148 (MQ=255)
                  tCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCg                                                                                                                              >  1:612456/1‑150 (MQ=255)
                   cAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCgg                                                                                                                             <  2:693659/150‑1 (MQ=255)
                    aGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGa                                                                                                                            >  2:4811596/1‑150 (MQ=255)
                     gggCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGAc                                                                                                                           <  2:2220424/150‑1 (MQ=255)
                     gggCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGAc                                                                                                                           <  2:3750629/150‑1 (MQ=255)
                      ggCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACg                                                                                                                          >  2:3413481/1‑150 (MQ=255)
                         ttGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTgg                                                                                                                       >  2:2271719/1‑150 (MQ=255)
                           gTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTATTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTCGTCGGACGTGGCt                                                                                                                     >  1:3716728/1‑150 (MQ=255)
                             ccGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTg                                                                                                                   >  2:1203905/1‑150 (MQ=255)
                                    tCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgc                                                                                                            >  2:4202171/1‑150 (MQ=255)
                                    tCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgc                                                                                                            >  2:3239107/1‑150 (MQ=255)
                                    tCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgc                                                                                                            <  2:4449367/150‑1 (MQ=255)
                                         gCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcgc                                                                                                       >  1:3860867/1‑150 (MQ=255)
                                          cTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcgcg                                                                                                      >  1:3342840/1‑150 (MQ=255)
                                           tGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcgcgg                                                                                                     >  2:2952007/1‑150 (MQ=255)
                                              aaTGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGc                                                                                                  <  2:3519670/150‑1 (MQ=255)
                                              aaTGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGc                                                                                                  <  1:992159/150‑1 (MQ=255)
                                                tGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGt                                                                                                <  1:4776783/150‑1 (MQ=255)
                                                tGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGt                                                                                                <  1:2146512/150‑1 (MQ=255)
                                                  tCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtct                                                                                              >  2:4346152/1‑150 (MQ=255)
                                                    aGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTATGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtcttc                                                                                            >  1:4992977/1‑150 (MQ=255)
                                                     gAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtcttct                                                                                           <  1:1262111/150‑1 (MQ=255)
                                                      aaaCCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTcttctt                                                                                          >  1:2559774/1‑150 (MQ=255)
                                                        aCCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGTCGCGGTGCGTCttctttc                                                                                        <  2:4627958/150‑1 (MQ=255)
                                                        aCCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCttctttc                                                                                        <  2:2501265/150‑1 (MQ=255)
                                                          cGGGGTGAGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTtctc                                                                                      >  1:4123620/1‑150 (MQ=255)
                                                           ggggTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCa                                                                                     <  1:1889517/150‑1 (MQ=255)
                                                             ggTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTGTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCAcc                                                                                   >  1:4236108/1‑150 (MQ=255)
                                                              ctcagCTGATCTACGCGAGGTGCTTATACAGTTATTCGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCAACTTAGCGGCGAGGCGTAGACGCACGACCCTTGCCATCGCCGCGCGGTGCGTCATCATACTCACca                                                                                  >  2:1259606/5‑150 (MQ=255)
                                                               tgtgCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACcag                                                                                 >  2:4893564/1‑150 (MQ=255)
                                                                   cTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGt                                                                             >  1:555303/1‑150 (MQ=255)
                                                                    tGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGtt                                                                            >  1:687637/1‑150 (MQ=255)
                                                                    tGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGtt                                                                            >  2:1492185/1‑150 (MQ=255)
                                                                     gTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGATGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTg                                                                           >  1:3391808/1‑150 (MQ=255)
                                                                      ttcttcGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGc                                                                          <  2:2352732/150‑1 (MQ=255)
                                                                      tncttcGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGc                                                                          >  1:116704/1‑150 (MQ=255)
                                                                               aGGTGCTTTTACTGATTTTGGCGCTTTCTCTACCCGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGAGGCGTGGCGTAGTCGGACGTGGCTTGCGTACGCGGCGCGTTGCGGCTTCTTTCTCACCAGCAGTTGCGGCAGCTTc                                                                 >  2:3210868/1‑150 (MQ=255)
                                                                                 gTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTcgc                                                               >  1:2002595/1‑150 (MQ=255)
                                                                                 gTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTcgc                                                               <  1:4555372/150‑1 (MQ=255)
                                                                                  tGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTcgcg                                                              >  2:3005987/1‑150 (MQ=255)
                                                                                   gCTTTGACTGTTTTTGGCGCTTTCTCGACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGt                                                             <  2:537464/150‑1 (MQ=255)
                                                                                    cTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGtt                                                            >  2:1323338/1‑150 (MQ=255)
                                                                                     ttttACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGAGGCGTGGCGTAGACGGACGTGGCTTGCGTTCGCGGCGCGTTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGttt                                                           >  1:3218682/1‑150 (MQ=255)
                                                                                        tACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGtttttt                                                        >  2:2121511/1‑150 (MQ=255)
                                                                                             tttttGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCtt                                                   <  2:1725995/150‑1 (MQ=255)
                                                                                             tttttGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCtt                                                   <  1:3941408/150‑1 (MQ=255)
                                                                                              ttttGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTg                                                  >  1:1494320/1‑150 (MQ=255)
                                                                                              ttttGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTg                                                  >  1:962513/1‑150 (MQ=255)
                                                                                                 tGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTg                                               >  2:3320209/1‑150 (MQ=255)
                                                                                                 tGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTg                                               >  2:2407327/1‑150 (MQ=255)
                                                                                                  ggCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGt                                              >  2:1363216/1‑150 (MQ=255)
                                                                                                    cgcTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTc                                            >  2:2615063/1‑150 (MQ=255)
                                                                                                     gcgttCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCa                                           <  1:3618993/147‑1 (MQ=255)
                                                                                                     gcTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCa                                           >  2:192881/1‑150 (MQ=255)
                                                                                                      cTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAg                                          >  2:2763191/1‑150 (MQ=255)
                                                                                                      cTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAg                                          >  1:4098604/1‑150 (MQ=255)
                                                                                                       tttCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGc                                         <  1:2959800/150‑1 (MQ=255)
                                                                                                       tttCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGc                                         <  1:2959598/150‑1 (MQ=255)
                                                                                                        ttCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCa                                        >  2:1816939/1‑150 (MQ=255)
                                                                                                        ttCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCa                                        <  1:3685057/150‑1 (MQ=255)
                                                                                                         tctcTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCAc                                       >  1:3539765/1‑150 (MQ=255)
                                                                                                         tctcTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCAc                                       >  2:3582273/1‑150 (MQ=255)
                                                                                                         tctcTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCAc                                       >  2:1005883/1‑150 (MQ=255)
                                                                                                         tctcTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCAc                                       >  2:1540883/1‑150 (MQ=255)
                                                                                                          ctctACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACg                                      <  2:3994455/150‑1 (MQ=255)
                                                                                                           tctACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGc                                     >  2:3964174/1‑150 (MQ=255)
                                                                                                            ctACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCt                                    >  2:642577/1‑150 (MQ=255)
                                                                                                              aCCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCtgt                                  >  1:1976107/1‑150 (MQ=255)
                                                                                                              aCCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCtgt                                  <  1:2717629/150‑1 (MQ=255)
                                                                                                                     tttGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGaa                           <  1:156200/150‑1 (MQ=255)
                                                                                                                     tttGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGaa                           <  2:248423/150‑1 (MQ=255)
                                                                                                                      ttGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAac                          >  2:4954982/1‑150 (MQ=255)
                                                                                                                       tGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAaca                         >  2:4743049/1‑150 (MQ=255)
                                                                                                                        gcgcACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAacac                        <  2:4375052/150‑1 (MQ=255)
                                                                                                                         cgcACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACAcg                       >  2:5116662/1‑150 (MQ=255)
                                                                                                                            aCGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACgcgc                    >  1:279871/1‑150 (MQ=255)
                                                                                                                            aCGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACgcgc                    >  2:3860134/1‑150 (MQ=255)
                                                                                                                             cGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCt                   >  2:4480737/1‑150 (MQ=255)
                                                                                                                                ggTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTc                <  2:4588925/150‑1 (MQ=255)
                                                                                                                                 gtgtACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCg               <  2:1876303/147‑1 (MQ=255)
                                                                                                                                   ttACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGct             >  1:1601914/1‑150 (MQ=255)
                                                                                                                                    tACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGctt            >  1:2183683/1‑150 (MQ=255)
                                                                                                                                    tACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGctt            <  2:4123620/150‑1 (MQ=255)
                                                                                                                                     aCGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGcttc           >  1:838045/1‑150 (MQ=255)
                                                                                                                                        ttCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCttcttc        <  2:2448522/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTG  >  CP009273/1909197‑1909488

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: