Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,802,047 G→A 100% intergenic (+100/‑6) ydiZ → / → yniA uncharacterized protein/fructosamine kinase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,802,0470GA96.4% 1097.5 / 21.0 309intergenic (+100/‑6)ydiZ/yniAuncharacterized protein/fructosamine kinase family protein
Reads supporting (aligned to +/- strand):  ref base G (10/1);  new base A (173/125);  total (183/126)
Fisher's exact test for biased strand distribution p-value = 3.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 315 total aligned reads displayed.

CTGAAATCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAGACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGAT  >  CP009273/1801898‑1802193
                                                                                                                                                     |                                                                                                                                                  
cTGAAATCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGgagag                                                                                                                                                    >  1:1993330/1‑150 (MQ=255)
   aaaTCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAGAca                                                                                                                                                 >  2:2260550/1‑150 (MQ=255)
   aaaTCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAAca                                                                                                                                                 >  1:864071/1‑150 (MQ=255)
    aaTCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTTGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAAAATGGGAGAAAcat                                                                                                                                                >  1:1519501/1‑150 (MQ=255)
    aaTCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAAcat                                                                                                                                                >  1:3988168/1‑150 (MQ=255)
        aCCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCAtg                                                                                                                                            >  1:2206363/1‑150 (MQ=255)
         ccATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCAtgt                                                                                                                                           >  2:3535085/1‑150 (MQ=255)
              aCGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTggcag                                                                                                                                      >  2:2953400/1‑150 (MQ=255)
               cGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTggcagg                                                                                                                                     >  2:4253626/1‑150 (MQ=255)
                  ccTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCaa                                                                                                                                  <  1:1682610/150‑1 (MQ=255)
                   cTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAAt                                                                                                                                 >  2:1296500/1‑150 (MQ=255)
                      ggaggaATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAg                                                                                                                              >  2:4112009/1‑150 (MQ=255)
                       gaggaATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAgt                                                                                                                             >  1:3283867/1‑150 (MQ=255)
                        aggaATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAgtc                                                                                                                            >  1:1117303/1‑150 (MQ=255)
                               gaatAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCtttt                                                                                                                     >  1:4256466/1‑150 (MQ=255)
                                 atAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTga                                                                                                                   >  1:127188/1‑150 (MQ=255)
                                    gCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTgagcg                                                                                                                <  2:1544704/150‑1 (MQ=255)
                                     cAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTgagcga                                                                                                               <  2:3916055/150‑1 (MQ=255)
                                       aCAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTgagcgagc                                                                                                             >  2:3175203/1‑150 (MQ=255)
                                       aCAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTgagcgagc                                                                                                             >  1:240061/1‑150 (MQ=255)
                                        cAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCa                                                                                                            <  1:3891706/150‑1 (MQ=255)
                                           aaaaTAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGtt                                                                                                         >  1:3627041/1‑150 (MQ=255)
                                           aaaaTAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGtt                                                                                                         >  1:551383/1‑150 (MQ=255)
                                            aaataaatGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTa                                                                                                        <  2:1671208/150‑1 (MQ=255)
                                              ataaatGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAgg                                                                                                      <  2:1916568/150‑1 (MQ=255)
                                              ataaatGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAgg                                                                                                      >  2:1222901/1‑150 (MQ=255)
                                              ataaatGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAgg                                                                                                      >  2:3795241/1‑150 (MQ=255)
                                                aaatGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTg                                                                                                    >  1:3763548/1‑150 (MQ=255)
                                                  atGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGaa                                                                                                  >  2:576088/1‑150 (MQ=255)
                                                   tGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAg                                                                                                 >  2:1754445/1‑150 (MQ=255)
                                                    ggTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAgg                                                                                                >  1:1735724/1‑150 (MQ=255)
                                                      tttttGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCg                                                                                              <  2:864071/150‑1 (MQ=255)
                                                        tttGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGaa                                                                                            >  2:4667007/1‑150 (MQ=255)
                                                        tttGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGaa                                                                                            <  1:753073/150‑1 (MQ=255)
                                                         ttGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGaac                                                                                           >  2:5088289/1‑149 (MQ=255)
                                                         ttGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGaaa                                                                                           <  2:1358060/150‑1 (MQ=255)
                                                          tGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAAt                                                                                          >  1:2634222/1‑150 (MQ=255)
                                                          tGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAAt                                                                                          >  2:5036617/1‑150 (MQ=255)
                                                            ggCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCg                                                                                        >  2:1524472/1‑150 (MQ=255)
                                                             gCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGa                                                                                       >  2:3653503/1‑150 (MQ=255)
                                                               aataatCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAAc                                                                                     >  2:4328920/1‑150 (MQ=255)
                                                                ataatCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACt                                                                                    <  1:2873838/150‑1 (MQ=255)
                                                                ataatCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACt                                                                                    >  1:1835694/1‑150 (MQ=255)
                                                                    tCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGt                                                                                >  1:2731673/1‑150 (MQ=255)
                                                                     cAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTa                                                                               >  2:3266952/1‑150 (MQ=255)
                                                                     cAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTa                                                                               <  2:2863715/150‑1 (MQ=255)
                                                                     cAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTa                                                                               <  2:623918/150‑1 (MQ=255)
                                                                      aGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTaa                                                                              >  1:1946930/1‑150 (MQ=255)
                                                                        tCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATg                                                                            <  1:3210518/150‑1 (MQ=255)
                                                                           gtggtgTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAAc                                                                         >  2:4270871/1‑150 (MQ=255)
                                                                           gtggtgTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAAc                                                                         >  2:3024795/1‑150 (MQ=255)
                                                                            tggtgTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACt                                                                        >  1:4356449/1‑150 (MQ=255)
                                                                             ggtgTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTg                                                                       <  1:3709871/150‑1 (MQ=255)
                                                                              gtgtGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGc                                                                      <  1:2425582/150‑1 (MQ=255)
                                                                                  gCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTgg                                                                  <  1:4477715/150‑1 (MQ=255)
                                                                                   cGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGc                                                                 <  2:4729944/150‑1 (MQ=255)
                                                                                    gTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCg                                                                >  1:3346906/1‑150 (MQ=255)
                                                                                    gTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACACCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCg                                                                >  2:4491115/1‑150 (MQ=255)
                                                                                      tAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGga                                                              >  1:4685411/1‑150 (MQ=255)
                                                                                          tCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAg                                                          <  2:737789/150‑1 (MQ=255)
                                                                                            gtgtTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCCTGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTc                                                        <  1:3031573/150‑1 (MQ=255)
                                                                                            gtgtTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTc                                                        <  1:821206/150‑1 (MQ=255)
                                                                                            gtgtTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTc                                                        <  1:4874819/150‑1 (MQ=255)
                                                                                             tgtTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTcc                                                       >  1:3157357/1‑150 (MQ=255)
                                                                                               tttttACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCAt                                                     <  2:4878239/150‑1 (MQ=255)
                                                                                               tttttACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCAt                                                     <  1:3631019/150‑1 (MQ=255)
                                                                                                  ttACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAGACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGcc                                                  <  1:2260550/150‑1 (MQ=255)
                                                                                                      acCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCAt                                              >  1:105547/1‑150 (MQ=255)
                                                                                                          cATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCa                                          <  2:2124669/150‑1 (MQ=255)
                                                                                                          cATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCa                                          >  2:3851770/1‑150 (MQ=255)
                                                                                                              cTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTg                                      >  2:3406802/1‑150 (MQ=255)
                                                                                                               ttGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTgc                                     <  2:1071967/150‑1 (MQ=255)
                                                                                                               ttGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTgc                                     <  1:901370/150‑1 (MQ=255)
                                                                                                                  cgcTAACCTTATGATCTGGCAGACAACATGGGAGAGACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCt                                  >  1:4848654/1‑150 (MQ=255)
                                                                                                                    cTAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTAt                                <  1:1143355/150‑1 (MQ=255)
                                                                                                                     tAACCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATg                               >  1:5008847/1‑150 (MQ=255)
                                                                                                                      aaCCTTATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGc                              <  1:3468882/150‑1 (MQ=255)
                                                                                                                          ttATGATCTGGCAGACAACATGGGAGAGACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGa                          >  1:2129291/1‑150 (MQ=255)
                                                                                                                          tgaggctCTGGCAGACAACCGGGGAGAAACCTCATGTGGCAGGCAATCAGTCGTCTTGTGAGCGAGCAGTTAGGTGACGGCGAAATCGAACTGCGTAATGCACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGa                          <  2:1414590/144‑1 (MQ=255)
                                                                                                                           tcggagcTGGCAGACAACCTGGGAGAAACATCATGTGGCAGGCACTCAGTCGTCTTTTGAGCGAGCAGTGAGGTGCAGGCGAAATCGAACTGCGTCATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTGGCGCTATGCAGGAc                         <  2:4438821/144‑1 (MQ=255)
                                                                                                                           tATGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGAc                         >  2:3711936/1‑150 (MQ=255)
                                                                                                                             tGATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACAt                       <  2:2477288/150‑1 (MQ=255)
                                                                                                                              gATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATg                      >  1:4781494/1‑150 (MQ=255)
                                                                                                                              gATCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATg                      >  2:277549/1‑150 (MQ=255)
                                                                                                                               aTCTGGCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGa                     >  2:1254865/1‑150 (MQ=255)
                                                                                                                                   ggCAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtttt                 <  1:3429295/150‑1 (MQ=255)
                                                                                                                                     cAGACAACATGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtttttt               <  1:1006371/150‑1 (MQ=255)
                                                                                                                                        acaacaTGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGt            >  2:4562867/1‑150 (MQ=255)
                                                                                                                                           acaTGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCaa         <  1:1944949/150‑1 (MQ=255)
                                                                                                                                            caTGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGATATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCaaa        <  2:1578947/150‑1 (MQ=255)
                                                                                                                                            caTGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCaaa        <  2:2078195/150‑1 (MQ=255)
                                                                                                                                            caTGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCaaa        <  2:776637/150‑1 (MQ=255)
                                                                                                                                             aTGGGAGAGACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAAt       >  1:4655154/1‑150 (MQ=255)
                                                                                                                                             aTGGGAGAGACATCATGTGGCAGGCAATCAGTAGTCTTTTGAGAGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTAGCGCTATGCAGGACATGCTTTTTTCGTCAAAt       >  1:4762137/1‑150 (MQ=255)
                                                                                                                                             aTGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAAt       <  1:72731/150‑1 (MQ=255)
                                                                                                                                             aTGGGAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAAt       >  1:4068030/1‑150 (MQ=255)
                                                                                                                                               gggAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAAtgt     <  1:1463677/150‑1 (MQ=255)
                                                                                                                                                ggAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAAtgtg    <  1:3016475/150‑1 (MQ=255)
                                                                                                                                                ggAGAAACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAAtgtg    >  1:2477659/1‑150 (MQ=255)
                                                                                                                                                  agaaaCATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGtgat  <  1:4414116/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTGAAATCACCATCACGACCTGGGAGGAATGGAATAGCAACAGAAAATAAATGGTTTTTGGGCAATAATCAGTCTGTGGTGTGCGTTAGCTCGTGTTTTTACACCGCATTCTTGCGCTAACCTTATGATCTGGCAGACAACATGGGAGAGACATCATGTGGCAGGCAATCAGTCGTCTTTTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGAT  >  CP009273/1801898‑1802193

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: