Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,909,342 T→A 83.5% 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,3420TA83.5% 778.0 / 140.5 272E150D (GAA→GATproQRNA chaperone, putative ProP translation regulator
Reads supporting (aligned to +/- strand):  ref base T (25/20);  new base A (118/109);  total (143/129)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 277 total aligned reads displayed.

TTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTG  >  CP009273/1909196‑1909488
                                                                                                                                                  |                                                                                                                                                  
ttttGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCcatc                                                                                                                                                  >  2:2356545/1‑149 (MQ=255)
ttttGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCcatc                                                                                                                                                  >  2:3739273/1‑149 (MQ=255)
    gACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTg                                                                                                                                              >  1:1228697/1‑149 (MQ=255)
    gACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTg                                                                                                                                              >  1:1228862/1‑149 (MQ=255)
        cGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTgcggc                                                                                                                                          <  2:2537653/149‑1 (MQ=255)
            ttCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCggcgtgg                                                                                                                                      <  1:2523269/149‑1 (MQ=255)
               aCCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCggcgtggcgt                                                                                                                                   >  2:93214/1‑149 (MQ=255)
               aCCTTCAGGGCGTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCggcgtggcgt                                                                                                                                   <  1:2842251/149‑1 (MQ=255)
                  ttCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGt                                                                                                                                >  2:2422616/1‑149 (MQ=255)
                  ttCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGt                                                                                                                                >  2:2453285/1‑149 (MQ=255)
                   tCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTc                                                                                                                               <  2:2245044/149‑1 (MQ=255)
                   tCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTc                                                                                                                               >  2:2206378/1‑149 (MQ=255)
                      gggCTTGTCCGACAGTCAGAGCTGAACTGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGa                                                                                                                            <  2:3628968/149‑1 (MQ=255)
                       ggCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGAc                                                                                                                           >  1:1066360/1‑149 (MQ=255)
                        gCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACg                                                                                                                          >  1:2563106/1‑149 (MQ=255)
                              ccGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCtt                                                                                                                    >  2:2830426/1‑149 (MQ=255)
                                gACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGc                                                                                                                  <  1:3189177/149‑1 (MQ=255)
                                        gagCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcgg                                                                                                          >  1:4353207/1‑149 (MQ=255)
                                          gCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcg                                                                                                        <  2:819095/149‑1 (MQ=255)
                                            tGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcgcg                                                                                                      >  2:1835563/1‑149 (MQ=255)
                                            tGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcgcg                                                                                                      >  2:1081964/1‑149 (MQ=255)
                                            tGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcgcg                                                                                                      >  2:2944603/1‑149 (MQ=255)
                                            tGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcgcg                                                                                                      >  2:3418124/1‑149 (MQ=255)
                                             gAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcgcgg                                                                                                     <  1:4407848/149‑1 (MQ=255)
                                              aaaTGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGt                                                                                                    <  1:4127639/149‑1 (MQ=255)
                                              aaaTGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGt                                                                                                    >  2:4329830/1‑149 (MQ=255)
                                               aaTGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTg                                                                                                   <  1:99968/149‑1 (MQ=255)
                                               aaTGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTg                                                                                                   >  2:494725/1‑149 (MQ=255)
                                                aTGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGc                                                                                                  >  2:310240/1‑149 (MQ=255)
                                                  gTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGt                                                                                                >  1:3604506/1‑149 (MQ=255)
                                                  gTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGt                                                                                                >  1:3604723/1‑149 (MQ=255)
                                                    cAGAAACCGGGGTGTGCTGTTCGTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtct                                                                                              <  2:351121/149‑1 (MQ=255)
                                                     aGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtctt                                                                                             >  1:216401/1‑149 (MQ=255)
                                                     aGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtctt                                                                                             <  1:1411108/149‑1 (MQ=255)
                                                     aGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtctt                                                                                             <  1:1142637/149‑1 (MQ=255)
                                                     aGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtctt                                                                                             <  1:3813041/149‑1 (MQ=255)
                                                     aGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtctt                                                                                             <  1:3227295/149‑1 (MQ=255)
                                                     aGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtctt                                                                                             >  1:1759527/1‑149 (MQ=255)
                                                          ccGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCttctttc                                                                                        <  2:3483759/149‑1 (MQ=255)
                                                            ggggTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTtctc                                                                                      >  1:3954226/1‑149 (MQ=255)
                                                            ggggTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTtctc                                                                                      >  2:187944/1‑149 (MQ=255)
                                                             gggTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCa                                                                                     <  2:3314579/149‑1 (MQ=255)
                                                              ggTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCAc                                                                                    >  2:158166/1‑149 (MQ=255)
                                                                  tgCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACcagc                                                                                >  1:1591416/1‑149 (MQ=255)
                                                                     tGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGt                                                                             <  1:1338771/149‑1 (MQ=255)
                                                                      gTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGtt                                                                            >  2:776225/1‑149 (MQ=255)
                                                                      gTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTCCTCACCAGCAGtt                                                                            >  1:1058737/1‑149 (MQ=255)
                                                                       ttcttcGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTg                                                                           >  1:210259/1‑149 (MQ=255)
                                                                        tcttcGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGc                                                                          >  1:4049783/1‑149 (MQ=255)
                                                                          ttcGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCgg                                                                        <  1:201299/149‑1 (MQ=255)
                                                                          ttcGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCgg                                                                        >  1:3729067/1‑149 (MQ=255)
                                                                             gcgAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAg                                                                     <  1:4583927/149‑1 (MQ=255)
                                                                                aGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCtt                                                                  <  2:3211832/149‑1 (MQ=255)
                                                                                aGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCtt                                                                  >  2:1148573/1‑149 (MQ=255)
                                                                                aGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCtt                                                                  >  2:2642230/1‑149 (MQ=255)
                                                                                aGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCtt                                                                  >  1:2275880/1‑149 (MQ=255)
                                                                                 ggTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTc                                                                 >  1:1956731/1‑149 (MQ=255)
                                                                                 ggTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTc                                                                 <  2:1260698/149‑1 (MQ=255)
                                                                                 ggTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTc                                                                 <  1:2311312/149‑1 (MQ=255)
                                                                                   tGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTcgc                                                               >  1:2694247/1‑149 (MQ=255)
                                                                                   tGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTcgc                                                               >  1:353097/1‑149 (MQ=255)
                                                                                     cTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGt                                                             >  1:2593928/1‑149 (MQ=255)
                                                                                      ttttACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGtt                                                            <  1:4115311/149‑1 (MQ=255)
                                                                                         tACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGttttt                                                         <  2:2080192/149‑1 (MQ=255)
                                                                                          aCTGTTTTTGGCGCTTTCTCGACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGtttttt                                                        <  1:40675/149‑1 (MQ=255)
                                                                                             gTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGc                                                     <  2:892508/149‑1 (MQ=255)
                                                                                              tttttGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCt                                                    >  2:2614018/1‑149 (MQ=255)
                                                                                              tttttGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCt                                                    >  2:1965588/1‑149 (MQ=255)
                                                                                                    gcgcTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGt                                              <  2:99365/149‑1 (MQ=255)
                                                                                                     cgcTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGtt                                             <  1:567090/149‑1 (MQ=255)
                                                                                                      gcTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTc                                            >  1:995644/1‑149 (MQ=255)
                                                                                                       cTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCa                                           >  2:2757983/1‑149 (MQ=255)
                                                                                                       cTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCa                                           >  1:3650154/1‑149 (MQ=255)
                                                                                                        tttCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAg                                          >  2:1598503/1‑149 (MQ=255)
                                                                                                         ttCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGc                                         >  1:739404/1‑149 (MQ=255)
                                                                                                           ctctACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCAc                                       <  1:916962/149‑1 (MQ=255)
                                                                                                            tctACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACg                                      >  1:3021771/1‑149 (MQ=255)
                                                                                                             ctACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGc                                     <  2:2563106/149‑1 (MQ=255)
                                                                                                             ctACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGc                                     <  2:3954226/149‑1 (MQ=255)
                                                                                                                ccGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCtgt                                  <  1:814001/149‑1 (MQ=255)
                                                                                                                ccGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCtgt                                  <  2:836495/149‑1 (MQ=255)
                                                                                                                  ggCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGc                                >  1:2335904/1‑149 (MQ=255)
                                                                                                                  ggCTTTTGCGCACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGc                                >  1:1654353/1‑149 (MQ=255)
                                                                                                                        tGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAac                          <  2:3721143/149‑1 (MQ=255)
                                                                                                                         gcgcTCGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAaca                         <  2:2035872/149‑1 (MQ=255)
                                                                                                                           gcACGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACAcg                       >  1:3310478/1‑149 (MQ=255)
                                                                                                                             aCGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACAcgcg                     >  2:654031/1‑149 (MQ=255)
                                                                                                                              cGAGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACgcgc                    <  1:816310/149‑1 (MQ=255)
                                                                                                                                aGGTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCtt                  <  2:3699053/149‑1 (MQ=255)
                                                                                                                                 ggTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCttt                 <  1:3820795/149‑1 (MQ=255)
                                                                                                                                 ggTTTACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCttt                 <  2:4548185/149‑1 (MQ=255)
                                                                                                                                     tACGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGct             <  2:2008211/149‑1 (MQ=255)
                                                                                                                                      aCGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGctt            >  1:1814796/1‑149 (MQ=255)
                                                                                                                                       cGTTCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGcttc           <  1:47698/149‑1 (MQ=255)
                                                                                                                                         ttCAGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGcttctt         <  2:2719499/149‑1 (MQ=255)
                                                                                                                                         ttCAGCGCCATCTTTGCCGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCAGAACACGCGCTATCTcttctt         >  2:2885573/1‑149 (MQ=255)
                                                                                                                                            aGCGCCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCAGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCaa      <  1:251341/149‑1 (MQ=255)
                                                                                                                                               gcCATCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAgct   >  2:4253476/1‑149 (MQ=255)
                                                                                                                                               gcCATCTTCGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAgct   <  1:1104801/149‑1 (MQ=255)
                                                                                                                                                ccTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAgctg  >  1:4552212/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTG  >  CP009273/1909196‑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: