Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 910,669 A→C 20.6% intergenic (‑356/+139) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273910,6690AC20.6% 91.7 / 2.8 75intergenic (‑356/+139)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base A (40/11);  new base C (10/2);  total (59/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGC  >  CP009273/910521‑910817
                                                                                                                                                    |                                                                                                                                                      
cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGa                                                                                                                                                        <  1:3288096/149‑1 (MQ=32)
 ggggCAGCCCCCATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt                                                                                                                                                       <  1:808300/149‑1 (MQ=17)
 ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt                                                                                                                                                       <  2:740558/149‑1 (MQ=32)
  gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                      <  2:3436146/149‑1 (MQ=32)
   ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                     >  2:1728201/1‑149 (MQ=32)
    gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                    >  1:650104/1‑149 (MQ=32)
    gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                    >  1:650120/1‑149 (MQ=32)
      aGCCCACATGGATGTGGGCTGAGGGCGCGGGGCCCCCCCCCCCAACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgcc                                                                                                                                                  <  2:3917540/149‑1 (MQ=11)
            caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                            <  1:1232096/149‑1 (MQ=32)
            caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                            <  1:770124/149‑1 (MQ=32)
            caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                            <  2:2474748/149‑1 (MQ=32)
             aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGccccca                                                                                                                                           >  1:3091169/1‑148 (MQ=17)
              tGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGCTTGGGCccccccag                                                                                                                                          >  1:205602/1‑147 (MQ=255)
              tGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccccg                                                                                                                                          >  1:3314436/1‑147 (MQ=17)
                 aTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑‑cgc                                                                                                                                     <  2:3855096/149‑1 (MQ=32)
                 aTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑‑cgc                                                                                                                                     >  2:3202688/1‑149 (MQ=32)
                  tgtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑‑cgcg                                                                                                                                    >  2:1706588/1‑149 (MQ=32)
                      ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTTTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGCGGCGGGTTGGCGCCCCCCg‑‑cccccccc                                                                                                                                >  1:1450843/1‑142 (MQ=11)
                          gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGGGGCGGTTTGCCGCCCCCCG‑‑CCCGCTCCATgc                                                                                                                            >  1:1524690/1‑149 (MQ=11)
                           aGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                      >  2:3617344/1‑124 (MQ=32)
                           aGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                      <  1:3617344/124‑1 (MQ=32)
                            gggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCCCCGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCAGGGGGGGGCGGCGGGTGGGCGCCCCCGG‑‑CGCGCGCCGTgcgc                                                                                                                          >  1:1825067/1‑149 (MQ=11)
                            gggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGCGGGGGGTGGGCGCCCCCCG‑‑CGCGCGCCATgcgc                                                                                                                          >  1:1825101/1‑149 (MQ=11)
                             ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGGGGCGGCGCGTGGCCCCCCCCCG‑‑CGCGCGCCGTGCGcc                                                                                                                         >  1:1818733/1‑149 (MQ=11)
                                gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGt                                                                                                                      >  2:334123/1‑149 (MQ=32)
                                gcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCCGt                                                                                                                      >  2:453863/1‑149 (MQ=17)
                                         aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCAGTGGCAATATg                                                                                                             >  2:2783809/1‑149 (MQ=255)
                                          gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGGTTGGGCGCCCCCCG‑‑CGCGCGCCGTGCGCCCGGGGCCAAATGt                                                                                                            >  1:465212/1‑149 (MQ=11)
                                            gATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGGGGCGGCGCTTGGGCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCCATATTTTg                                                                                                          >  1:1989186/1‑149 (MQ=11)
                                              tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGCGGCGGCGCGTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCt                                                                                                        >  2:3451667/1‑149 (MQ=14)
                                                 tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                     >  2:2730466/1‑149 (MQ=255)
                                                  aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCACTGGCCATATGTTCCTTAGc                                                                                                    >  2:2114091/1‑149 (MQ=255)
                                                   ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGGGGCAATATGTTGCTTAGCt                                                                                                   >  2:3168877/1‑149 (MQ=255)
                                                   ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCCGCGATTTGCCCCCCCCTG‑‑CGCGCTCCTTGCGCCAGGGGCAATATGTTGCTTATCt                                                                                                   >  2:1109860/1‑149 (MQ=14)
                                                    cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGCGGGGGTTTGCCCCCCCCCG‑‑CGCGCGCCCTTCGCCACGTGCCATATGTTGGTTAtctc                                                                                                  >  1:3870962/1‑149 (MQ=11)
                                                     tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGATTTGCCGCCCCCAG‑‑CGCGCTCCGTGCGCCACTGGCCATATGTTGCTTATCTCa                                                                                                 >  1:306278/1‑149 (MQ=14)
                                                     tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCGTGCGCCAGTGGCACTATGTTGCGTAGCTCa                                                                                                 >  2:3962610/1‑149 (MQ=255)
                                                       gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGCGGGGGGGGGGGGGGTTTGGCGCCCCCAG‑‑CGCGCGCCGTGCGCCCGTGGCAAAATTTTTGTTTGCTCAtt                                                                                               >  1:1219881/1‑148 (MQ=11)
                                                       gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                               >  2:2526383/1‑149 (MQ=255)
                                                        cgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCGTGCGCCCGTGGCAAAATGTTGCGTAGCTCATGa                                                                                              >  2:856720/1‑149 (MQ=255)
                                                            cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGGTGGCCGCCCCCCG‑‑CGCCCGCCCTGCGCCCCGGGCAAAATGTTTCTTATGTCATGGAAgg                                                                                          >  1:1478560/1‑149 (MQ=11)
                                                              aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCCTGAAAGGAg                                                                                        >  2:3607358/1‑149 (MQ=255)
                                                                  ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTG‑‑CGCGCTCCTTTCGCCAGTGGCAAAATTTTGCTTATCTAAAGAAAGGTGCGca                                                                                    >  1:2285484/1‑149 (MQ=14)
                                                                  ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGGATGGCCGCCCCCAG‑‑CGCCCACCATGCGCCCCTGGCAAAATTTTGCTTATCTCATGAAACGAGCGca                                                                                    >  1:3200464/1‑149 (MQ=14)
                                                                        cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCCCTTTGGCGCCCCCGG‑‑CGCGCGCCTTGCGCCACTGGCAAAATGTTTCTTAGCACAATAAAGGAGCGcaacaaaa                                                                              >  2:1426214/1‑147 (MQ=17)
                                                                             gggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCACTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGaagatg                                                                         >  2:3772695/1‑149 (MQ=255)
                                                                              ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCGTGGCCGCCCCCAG‑‑CGCGCTCCCTGCGCAAGTGGCAAAAATTTGCTTATCGCATGAAAGGAGCGCAACAACAagatga                                                                        >  1:979097/1‑149 (MQ=17)
                                                                                   aGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                   >  2:3437305/1‑149 (MQ=255)
                                                                                    gTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATaaaa                                                                  >  2:3934954/1‑148 (MQ=255)
                                                                                         gggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGGCCCCTG‑‑CGCGCGCCGTGCGCCACTGGCAATATGTTGGTTAGCTAATGAAAGGAGCGCAACCAGATGATGAATCAACAtata                                                             >  2:223423/1‑149 (MQ=21)
                                                                                          ggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCAG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCGTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAAATaa                                                            >  2:941573/1‑149 (MQ=255)
                                                                                           gCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCGG‑‑CGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATAAACATATaac                                                           >  2:3519203/1‑149 (MQ=255)
                                                                                               cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGGGGCGAATTGCCGCCCCCAG‑‑CGCGCTCCATTCGCCCCTGGCCAAATTTTTCTTATCTCATGAAAGGTGCGCAACAAAATGATTAAAAAACAAATaaaaaaa                                                       >  1:2068988/1‑147 (MQ=9)
                                                                                                gcgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCTG‑‑CCCGCTCCTTGCGCCAGTGGCAATAAGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATTATTCATACACATATAACAAAAt                                                      >  2:2166145/1‑149 (MQ=18)
                                                                                                           aTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCATGCGCCAGGGGCAAAATGTTGCTTAGCTCATGTAAGGAGCGCAACAAGATGTTGTATCAACATATAAAAAACTCTTAAAAAAAg                                           >  1:3711740/1‑149 (MQ=21)
                                                                                                            ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCGCCGTGCGCCAGTGGCAATATTTTTCTTAGCTCATTAAAGGTGCGCAACAAGAATATGAATCAAAAAATAAAAACAACTTAAAAAAAgg                                          >  1:616105/1‑149 (MQ=18)
                                                                                                            ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCCATTGGGCCCCCCCAG‑‑CGCGCTCCATGCGCCACTGGCAAAATTTTTCTTATCGCAATGAATGAGCGCAACAAGATGATGTATAAAAAAATAACAACAACTTaaaaaaaat                                          >  1:3566635/1‑147 (MQ=12)
                                                                                                             tttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCATGGCCGCCCCCAG‑‑CGCGCGCCTTGCGCCCGTGGCCATATTTTGCTTAGCTCATGTAAGGAGCGCAAAAACATGATGAATCAACATATAACAACAACTTAAAAAAAGGc                                         >  2:1423392/1‑149 (MQ=25)
                                                                                                              ttggccgcccgcgcccgGGGGGGCGGGGGGCGCGGGCGCGGGGCCCCCCCCCG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGcc                                        <  2:2657110/144‑1 (MQ=255)
                                                                                                                    gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTTCGCCAGGGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAACAACATATAACAACATCATAAAAAAAGGCCTGACAt                                  >  2:2847179/1‑149 (MQ=255)
                                                                                                                    gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAGATAACCACATCTTAAAAAAAGGCCAGACAt                                  >  2:2847170/1‑149 (MQ=255)
                                                                                                                     cccgcgCCGGGGGGGGCGGGGGGCGGGGGCTAGGGCCCCCCCCCCC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                                 <  2:2068988/144‑1 (MQ=37)
                                                                                                                         gggggggggggGCGGGGGGGGGGGGCTTTGGCCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                             <  2:3314436/144‑1 (MQ=255)
                                                                                                                          gCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCACCGTTCGCCAGTGGCAAAATTTTTCTTATCTAAAGAAACGAGCGCAACAAAATTATGAAAAAACAAATAACAAAAAAATAAAAAAAAGCCAGACAATaaccc                            >  2:19927/1‑145 (MQ=12)
                                                                                                                            gcGGGGGGGGGGGGGGGGGGGGGGCTTGGGCCCCCCCCGCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGcc                            <  2:3977542/148‑1 (MQ=255)
                                                                                                                               ggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                       >  2:2510622/1‑149 (MQ=255)
                                                                                                                                 gggggggggggCGGCGGCGTGGCGCCCCCCCCGC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                     <  1:1733313/143‑1 (MQ=38)
                                                                                                                                    gCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                  >  1:3560242/1‑149 (MQ=255)
                                                                                                                                       gggggCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt               >  1:584283/1‑149 (MQ=255)
                                                                                                                                        ggggggggggCGGTGGCCCCCCCCCGGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt               <  2:798193/149‑1 (MQ=255)
                                                                                                                                        ggggCGGGCGCGAGGGGCCGCCCCCGC‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt              <  2:3312138/149‑1 (MQ=255)
                                                                                                                                              gcggcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaattaa        >  2:3478553/1‑149 (MQ=255)
                                                                                                                                              gcggcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaattaa        >  2:3478547/1‑149 (MQ=255)
                                                                                                                                                 gcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCCGGCCTTCTGCGTTAATTAATCa     >  2:3882852/1‑149 (MQ=255)
                                                                                                                                                 gcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa     >  1:2694649/1‑149 (MQ=255)
                                                                                                                                                 gcgATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa     >  2:3882253/1‑149 (MQ=255)
                                                                                                                                                    aTTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc  >  1:2926979/1‑149 (MQ=255)
                                                                                                                                                    aTTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc  >  2:727920/1‑149 (MQ=255)
                                                                                                                                                    aTTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc  >  2:2789821/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGC  >  CP009273/910521‑910817

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: