Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 910,654 C→G 11.4% intergenic (‑341/+154) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273910,6540CG11.4% 238.0 / 3.8 96intergenic (‑341/+154)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base C (63/21);  new base G (8/4);  total (71/25)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCG  >  CP009273/910509‑910803
                                                                                                                                                 |                                                                                                                                                        
ttCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggg                                                                                                                                                      <  2:687593/150‑1 (MQ=32)
ttCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggg                                                                                                                                                      <  2:2245495/150‑1 (MQ=32)
  cccTGACGGCCGGGGCAGCCCACATGGTTGTGGGCTGCGGGCGCGTTTTACAGGGATGTTCCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                                                                                                                                                    <  1:821055/150‑1 (MQ=255)
  cccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATACGTCGCGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCggggg                                                                                                                                                    >  1:234242/1‑150 (MQ=12)
  cccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                                                                                                                                                    <  1:821034/150‑1 (MQ=32)
  cccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                                                                                                                                                    <  1:1809407/150‑1 (MQ=32)
  cccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                                                                                                                                                    >  1:234248/1‑150 (MQ=32)
  cccTGAAGGCCGGGGCAGCCCACAGGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                                                                                                                                                    <  1:4087004/150‑1 (MQ=17)
  cccGGATGGCCGGGGCAGCCCACTTGGTTGTGGGCTGCGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCGGTAAGTGATAAGTCGTGGGCACGGCGCAGCGGCGATTTTGTTGGCCAGAGCCCGGGGGTGCAggggg                                                                                                                                                    <  1:821048/150‑1 (MQ=11)
   ccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggc                                                                                                                                                   <  2:964977/150‑1 (MQ=32)
     tGAAGGCCGGGGCAGCCCACCTGGATGTGGGCGGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                                                 <  2:117327/150‑1 (MQ=17)
     tGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                                                 >  2:1507420/1‑150 (MQ=32)
     tGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                                                 >  1:1799617/1‑150 (MQ=32)
     tGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                                                 <  1:4001386/150‑1 (MQ=32)
       aaGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcggcg                                                                                                                                               >  2:2556716/1‑150 (MQ=32)
       aaGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcggcg                                                                                                                                               >  2:291738/1‑150 (MQ=32)
            cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGgcggccgcgga                                                                                                                                          >  1:909309/1‑148 (MQ=12)
            cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGGCCGGGggggcggcgga                                                                                                                                          >  1:3045075/1‑148 (MQ=11)
            cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTGTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCCGAGCCCGGGGGTGGCGgggggggggggggg                                                                                                                                          >  1:7368/1‑146 (MQ=11)
              gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                        <  2:1609124/150‑1 (MQ=32)
              gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                        <  2:1603298/150‑1 (MQ=32)
              gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                        <  2:3881246/150‑1 (MQ=32)
               ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                      <  2:3366260/150‑1 (MQ=32)
               ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                      >  2:1474703/1‑150 (MQ=32)
                  aGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGGGGGGCTTGgccg                                                                                                                                   >  1:894546/1‑150 (MQ=11)
                    cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                 <  2:1392139/150‑1 (MQ=32)
                     ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGGGGCGGCGGGTGGGCGccc                                                                                                                                >  1:847632/1‑150 (MQ=11)
                     ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGGTTGGCCGccc                                                                                                                                >  1:847609/1‑150 (MQ=17)
                      cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                               <  2:3700746/150‑1 (MQ=32)
                      cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                               >  2:2903534/1‑150 (MQ=32)
                      cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                               >  2:2905131/1‑150 (MQ=32)
                        caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccca                                                                                                                             >  2:1377553/1‑149 (MQ=17)
                        caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCt                                                                                                                             <  2:395415/150‑1 (MQ=32)
                         agggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                            <  1:1570686/148‑1 (MQ=17)
                          tGGTTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTATGGGATAAGTCGAGGGCACCGCGCAGGGGCGATTTTGTTCGCCAGAGCCGGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTgc                                                                                                                          <  2:2100961/150‑1 (MQ=11)
                          tGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATCAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgc                                                                                                                          <  2:2192096/150‑1 (MQ=17)
                              tgtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgcgc                                                                                                                      >  2:2952415/1‑150 (MQ=32)
                                tgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGGGGTTTGGCGCCCCCAgcgcgcgc                                                                                                                    >  1:3177161/1‑148 (MQ=11)
                                 gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGGGGCGGCGGTTGGCCCCCCCCCGCGCGCGcc                                                                                                                   >  1:2300484/1‑150 (MQ=11)
                                     tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTgc                                                                                                               >  2:3996534/1‑150 (MQ=17)
                                      gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTgcg                                                                                                              >  1:1499978/1‑150 (MQ=17)
                                      gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATgcg                                                                                                              >  1:1499914/1‑150 (MQ=12)
                                      gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgcg                                                                                                              >  2:1613590/1‑150 (MQ=32)
                                      gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATGTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGCGGCGGCGGTTGGGCGCCCCCCGCGCGCGCCTTgcg                                                                                                              >  1:4262207/1‑150 (MQ=11)
                                       aGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgcgc                                                                                                             >  2:2009300/1‑150 (MQ=32)
                                             cgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTgg                                                                                                       >  2:2345099/1‑150 (MQ=17)
                                              gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGGGGCGGCCCTTGGGCGCCCCCACCGCGCGCCATGCGCCacaggc                                                                                                      >  1:3072047/1‑145 (MQ=11)
                                              gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCCGCGCTTGGGCGCCCCCCCCGCGCGCCATGCGccccgggc                                                                                                      >  1:3908300/1‑144 (MQ=11)
                                              gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCCGGGGCGGCGGGGCGTGGGCGCCCCCAGCCCCCGCCCTGCGccccaggc                                                                                                      >  1:3908290/1‑144 (MQ=11)
                                                 ttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGGGCTTGGCCGCCCCCGGCGCGCGCCGTGCGCCCCTGGccaa                                                                                                   >  1:3315298/1‑147 (MQ=11)
                                                    cAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGGGGCGGCGGTTGGCCGCCCCCAGCCCGCTCCATTCGCCACTGGccaaaat                                                                                                >  1:1331492/1‑144 (MQ=11)
                                                      gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                             >  2:2562063/1‑150 (MQ=40)
                                                       ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGGGGCGAGTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCACAATGTTg                                                                                            >  2:2744927/1‑150 (MQ=255)
                                                         gtggtACCTCGCGCCCGACCGGGTTGCCGTTCGGGATTCGTGGGGGGCACCGCGCCGGGGTTTTTTTGTTCGCCAGTGCCCGGGGGTGGCGGGGGCGGCGGCGATTGGCCGCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGc                                                                                           <  2:2063373/146‑1 (MQ=9)
                                                         aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGGGGCGCGTGGCCGCCCCCAGCGCGCGCCATGCGCCCGTGGCAATATGTTGCt                                                                                          >  1:2137917/1‑150 (MQ=14)
                                                          tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCGTGCGCCAGGGGCAATATGTTGCtt                                                                                         >  2:525918/1‑150 (MQ=255)
                                                             tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCGCCATGCGCCCGGGGCACAATGTTTCGTAGc                                                                                      >  1:4425895/1‑150 (MQ=11)
                                                                cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCGTGCGCCACTGGCCATATGTTGCGTAGCTCa                                                                                   >  1:2550508/1‑150 (MQ=14)
                                                                cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCAGCGCGCTCCATGCGCCACTGGCAATATGTTTCTTAGCtaa                                                                                   >  1:4158146/1‑148 (MQ=14)
                                                                   gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTTCGCCAGTGGCAATATGTTGCTTAGCTCCTGa                                                                                >  2:2930443/1‑150 (MQ=255)
                                                                   gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCGCATGa                                                                                >  2:2931671/1‑150 (MQ=255)
                                                                     gcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTATCTCATGaaa                                                                              >  2:1324183/1‑150 (MQ=255)
                                                                              ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaa                                                                     >  2:2971047/1‑150 (MQ=255)
                                                                              ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGATCGcaa                                                                     >  2:2970134/1‑150 (MQ=255)
                                                                                 aGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCAGCCCGCACCATGCGCCCCGGGCAATAATTTTCTTAGCTCAATAAAGGAGCGcaacaa                                                                  >  1:3333993/1‑150 (MQ=14)
                                                                                    cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATTTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                               >  2:347599/1‑150 (MQ=255)
                                                                                    cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                               >  2:348371/1‑150 (MQ=255)
                                                                                    cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                               >  2:348270/1‑150 (MQ=255)
                                                                                    cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                               >  2:347826/1‑150 (MQ=255)
                                                                                    cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCCCCCTCCATGCGCCAGTGGCAAAATTTTTCTTATCTCAATAAAGGAGCGCAACAAcat                                                               >  2:2342874/1‑150 (MQ=14)
                                                                                    cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGGCCGCCCCCAGCGCCCGCCTTGCGCCCGTGGCCATCTGTTGCTGTGCGCCTGACAGGGGGGCAACAAgct                                                               >  2:2754196/1‑148 (MQ=14)
                                                                                    cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGATTGGCCGCCCCCAGCGCGCTCCATGCCCCAGTGGCAAAATGTTGCGTAGCTCAAGAAAGGAGCCCAACAAcat                                                               >  2:2754869/1‑150 (MQ=17)
                                                                                      tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTTGCCGCCCCCAGCGCGCTCCATGCGCCAGTGGCCAAATTTTTCATAGCACAATTAAGGAGCGCAACAAaatga                                                             >  1:1080772/1‑150 (MQ=14)
                                                                                            aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCACTGGCAATATTTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATc                                                       >  2:1430991/1‑150 (MQ=255)
                                                                                            aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCACCTTGCCCCAGTGGCAATAAGTTGCTTAGCACATGAAAGGAGAGCAACAAGATGATGAAta                                                       >  2:3196821/1‑149 (MQ=18)
                                                                                                gTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGGCCGCCCCCGGCGCGCGCCATGCGCCCGTGGCAATATTTTTCTTAGCTCATGTAAGGGGCGCAAACAGATGATGAATCAACa                                                   >  1:2193923/1‑150 (MQ=18)
                                                                                                gTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGGCCGCCCCCAGCGCGCTCCATGCGCCCCGGGCAAAAAGTTTATTAGCTCAAGAAAGGAGCGCAAAAACATTATGAATCAACa                                                   >  1:4029761/1‑150 (MQ=17)
                                                                                                   gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGCATGGCCGCCCCCAGCGCGCGCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCAAGAAAGGACCGCAACAACATGAAGAATCaaactat                                                >  2:4417427/1‑145 (MQ=18)
                                                                                                       gCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTTGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATTATGAATAAAAAAATaacc                                            >  2:4383541/1‑149 (MQ=18)
                                                                                                       gCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCGCCTTGCGCCCGTGGCCAAATGTTGCTTAGCTCATGAAAGGGGCGCAACAAGATGATGAATCAACAAATaaca                                            >  2:650513/1‑150 (MQ=18)
                                                                                                       gCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATAAAAATATaaca                                            >  2:4383249/1‑150 (MQ=18)
                                                                                                       gCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATTCGCCAGTGGCAATATTTTTCTTATCTCATTAAATGAGCGCAAAAAGAATTTTAATAAAAATATaaaa                                            >  1:2626015/1‑148 (MQ=14)
                                                                                                            gcgcAGCGGTTTTTTTTTTCGCCAGTGGGGGGGGTGGCTGGGGGCGGCGTCATTGCCCCCCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATc                                       <  1:2151867/150‑1 (MQ=21)
                                                                                                               ccggggctATTTTGTCCGCGTGCGCGGGGGGGGGGGGGGGGGGGCGGCGTGCGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTa                                    <  2:376387/146‑1 (MQ=255)
                                                                                                               cAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAATGAGCGCAACAAGATGATTTAACAACATATAACAAAATATTa                                    >  2:2259951/1‑150 (MQ=18)
                                                                                                                      gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCCCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTAAAGAAAGGAGCGCAACAACATGATGAATCAAAAAATAAAAAAATAATaaaaaaaa                             >  2:3861021/1‑149 (MQ=18)
                                                                                                                            gTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGa                       >  2:1367942/1‑150 (MQ=255)
                                                                                                                              tCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTTGCCCCCCCCGGCGCGCGCCGTTCGCCCGTGGCAAAAATGTGGGTTGCTCCAGTAAGGAGCGCAACAAGAAGATGAAACAACAAATAACAACAACTTAAAAAAAGGCCGGaaa                     >  2:1311968/1‑148 (MQ=16)
                                                                                                                              tCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCCGTGGCCATATGTTGCTTAGCTCATGAAAGGAGCGCAACAACATGATGAATCAACATATAACAAACTCATAAAAAAAGGCCTGACa                     >  2:3650054/1‑150 (MQ=25)
                                                                                                                               cGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCCTGAAAGGAGCGCAACAAGAAGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACaa                    >  2:3345558/1‑149 (MQ=255)
                                                                                                                                     agCCCGGGGGTGCAGGGGGGGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAAAAATTTGCTTAGGTCATGAAAGGAGCGCAACAAAATGATGAATCAACATAAAACAACATCTTAAAAAAAAGCCTGACAATAAGcc              >  1:2205160/1‑150 (MQ=25)
                                                                                                                                        tggggggTGGGGGGGGGGGGTGTGTTTGCCCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg            <  2:3092636/148‑1 (MQ=255)
                                                                                                                                         cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc          >  2:1404603/1‑150 (MQ=255)
                                                                                                                                         cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGc          >  2:1814879/1‑150 (MQ=255)
                                                                                                                                          gggggTGGTGGGGGGGGGGGGTTTTGGCCGCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc         <  2:2181556/150‑1 (MQ=255)
                                                                                                                                          gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc         >  2:2385677/1‑150 (MQ=255)
                                                                                                                                                 cAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg  >  2:789607/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                        
TTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCG  >  CP009273/910509‑910803

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