Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,586,861 G→T 16.7% intergenic (+32/+106) yjjN → / ← opgB L‑galactonate oxidoreductase/phosphoglycerol transferases I and II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,586,8610GT16.7% 766.5 / 90.5 248intergenic (+32/+106)yjjN/opgBL‑galactonate oxidoreductase/phosphoglycerol transferases I and II
Reads supporting (aligned to +/- strand):  ref base G (117/89);  new base T (20/22);  total (137/111)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

AAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCC  >  CP009273/4586712‑4587010
                                                                                                                                                     |                                                                                                                                                     
  aaaTCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTTATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTc                                                                                                                                                     <  1:421472/150‑1 (MQ=255)
  aaaTCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTc                                                                                                                                                     >  1:2349022/1‑150 (MQ=255)
   aaTCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTcc                                                                                                                                                    <  1:1975313/150‑1 (MQ=255)
    attaCTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCg                                                                                                                                                   <  2:435212/147‑1 (MQ=255)
      cACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCggg                                                                                                                                                 >  2:1870662/1‑150 (MQ=255)
       aCTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgggg                                                                                                                                                >  1:4118904/1‑150 (MQ=255)
          ccagACATGATTTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGAGAGTTAATTAAAGGCGTAATTACTTTATGATAAGGCGAGATTCTTAAAGTTGCCATGCAGCGTCCGGGGAAg                                                                                                                                             >  2:717409/4‑150 (MQ=255)
           ctgACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAgt                                                                                                                                            <  2:1522892/150‑1 (MQ=255)
            tgACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAgtg                                                                                                                                           >  2:2969249/1‑150 (MQ=255)
             gACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCTTCCGGGGAAgtgt                                                                                                                                          <  1:1693527/150‑1 (MQ=255)
              aCATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCTTCCGGGGAAGTGtt                                                                                                                                         <  1:1767651/150‑1 (MQ=255)
                atgatgTTAGCCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTgg                                                                                                                                       <  1:4608094/150‑1 (MQ=255)
                atgatgTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTgg                                                                                                                                       <  1:3881730/150‑1 (MQ=255)
                atgatgTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTgg                                                                                                                                       >  2:4045812/1‑150 (MQ=255)
                atgatgTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTgg                                                                                                                                       >  2:4094850/1‑150 (MQ=255)
                   atgTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGgcg                                                                                                                                    <  2:728305/150‑1 (MQ=255)
                       tAACCCTTCGCTATCCGTTCGCCACGGCGGCTGAAACCTACGAGCGCGATGTGTTTAACTATCGTGAGTTAATTAACGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGt                                                                                                                                <  2:3909072/150‑1 (MQ=255)
                        aaCCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGtt                                                                                                                               >  2:3967332/1‑150 (MQ=255)
                            cATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGtttttt                                                                                                                           <  1:577540/150‑1 (MQ=255)
                              tCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGtttttttt                                                                                                                         >  2:2716655/1‑150 (MQ=255)
                               cGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTg                                                                                                                        >  1:1021830/1‑150 (MQ=255)
                                    tCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCt                                                                                                                   <  1:451158/150‑1 (MQ=255)
                                     ccGTTCGCCACGCTGGTAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCtt                                                                                                                  <  1:4124173/150‑1 (MQ=255)
                                      cGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTta                                                                                                                 >  2:4568162/1‑150 (MQ=255)
                                      cGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTta                                                                                                                 >  2:2676526/1‑150 (MQ=255)
                                      cGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTta                                                                                                                 >  2:1122665/1‑150 (MQ=255)
                                      cGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTta                                                                                                                 >  2:1234802/1‑150 (MQ=255)
                                       gTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTtaa                                                                                                                <  2:2148866/150‑1 (MQ=255)
                                       gTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGGGTTGGGCGCTGTTTTTTTTGTTTCTtaa                                                                                                                >  2:308951/1‑150 (MQ=255)
                                         tCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTtaata                                                                                                              <  1:1682566/150‑1 (MQ=255)
                                           gCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTaataat                                                                                                            >  1:543303/1‑150 (MQ=255)
                                           gCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTaataat                                                                                                            <  2:1070644/150‑1 (MQ=255)
                                             cACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAAtgt                                                                                                          >  1:246688/1‑150 (MQ=255)
                                              aCGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCTTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAAtgtg                                                                                                         <  2:3762311/150‑1 (MQ=255)
                                               cGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCTTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATgtgt                                                                                                        <  1:441010/150‑1 (MQ=255)
                                                 cTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGtgttg                                                                                                      >  2:522040/1‑150 (MQ=255)
                                                      aGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCTTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGc                                                                                                 >  2:256762/1‑150 (MQ=255)
                                                      aGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGc                                                                                                 <  1:3035945/150‑1 (MQ=255)
                                                          aCCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCTTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTa                                                                                             <  2:2573909/150‑1 (MQ=255)
                                                          aCCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTa                                                                                             >  1:1283214/1‑150 (MQ=255)
                                                           ccTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAg                                                                                            >  1:3295637/1‑150 (MQ=255)
                                                           ccTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAg                                                                                            >  1:3670717/1‑150 (MQ=255)
                                                            cTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGa                                                                                           >  2:2518013/1‑150 (MQ=255)
                                                             tACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGaa                                                                                          >  1:3329595/1‑150 (MQ=255)
                                                              aCGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAg                                                                                         <  2:866214/150‑1 (MQ=255)
                                                              aCGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAg                                                                                         >  1:3172112/1‑150 (MQ=255)
                                                                 aGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCg                                                                                      <  2:1625809/150‑1 (MQ=255)
                                                                   cgcgATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCgtg                                                                                    <  1:618782/150‑1 (MQ=255)
                                                                   cgcgATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCgtg                                                                                    <  1:1910017/150‑1 (MQ=255)
                                                                     cgATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCTTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTa                                                                                  >  2:1083590/1‑150 (MQ=255)
                                                                      gATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAg                                                                                 <  2:4402747/150‑1 (MQ=255)
                                                                        tgtgATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGACGTGTGGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGt                                                                               >  2:3918673/1‑150 (MQ=255)
                                                                           gATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTCGAAGGCGTGTAGGTCGc                                                                            >  2:1549424/1‑150 (MQ=255)
                                                                                aCAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCt                                                                       >  2:3610291/1‑150 (MQ=255)
                                                                                aCAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCt                                                                       <  2:2808408/150‑1 (MQ=255)
                                                                                   aTCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATg                                                                    >  1:362089/1‑150 (MQ=255)
                                                                                     cGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCg                                                                  >  1:269230/1‑150 (MQ=255)
                                                                                       tGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGAc                                                                >  1:2778368/1‑150 (MQ=255)
                                                                                        gAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGAcc                                                               >  1:791369/1‑150 (MQ=255)
                                                                                         aGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCt                                                              <  2:4310293/150‑1 (MQ=255)
                                                                                         aGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCt                                                              <  2:4310380/150‑1 (MQ=255)
                                                                                         aGTTAATTAAAGGCGTAATTACGTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCt                                                              <  2:851196/150‑1 (MQ=255)
                                                                                                 aaaGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtca                                                      >  1:1174539/1‑150 (MQ=255)
                                                                                                      cGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtcagtcag                                                 >  2:3908472/1‑150 (MQ=255)
                                                                                                      cGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtcagtcag                                                 >  2:1881076/1‑150 (MQ=255)
                                                                                                       gTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCTTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAgg                                                >  2:465431/1‑150 (MQ=255)
                                                                                                        tAATTACTTTCTGATAAGGCGAGATTATGAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGa                                               <  2:3566133/150‑1 (MQ=255)
                                                                                                         aaTTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGaa                                              >  2:2454218/1‑150 (MQ=255)
                                                                                                         aaTTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGaa                                              >  2:369389/1‑150 (MQ=255)
                                                                                                          aTTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAAt                                             <  1:3488846/150‑1 (MQ=255)
                                                                                                             aCTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAc                                          <  1:4025968/150‑1 (MQ=255)
                                                                                                             aCTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCCTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAc                                          >  2:3168849/1‑150 (MQ=255)
                                                                                                              cattCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTAGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAcc                                         <  2:4441601/148‑1 (MQ=255)
                                                                                                              cTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAcc                                         <  2:4576361/150‑1 (MQ=255)
                                                                                                                 tCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCtt                                      <  2:650716/150‑1 (MQ=255)
                                                                                                                 gcTGATAAGGCGCGATTATTAAAGTTCCCATGCAGCGTCCGGGGACCTGTTGGGCGCTCTTTATTTTGTTTCTTAATCATGTGTTGAAAGACGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCtt                                      <  2:2722703/149‑1 (MQ=255)
                                                                                                                  cTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTc                                     >  1:2619090/1‑150 (MQ=255)
                                                                                                                     aTAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACg                                  <  1:3602169/150‑1 (MQ=255)
                                                                                                                     aTAAGGCGAGATTATTAAAGGTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACg                                  <  1:632857/150‑1 (MQ=255)
                                                                                                                       aaGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGtt                                >  2:2243676/1‑150 (MQ=255)
                                                                                                                         ggCGAGATTATTAAAGTTGCCATGCAGCTTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCt                              >  1:3571652/1‑150 (MQ=255)
                                                                                                                         ggCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCt                              <  2:3325644/150‑1 (MQ=255)
                                                                                                                           cGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTAc                            <  2:1044193/150‑1 (MQ=255)
                                                                                                                                ttattaAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACtt                       >  1:2259806/1‑150 (MQ=255)
                                                                                                                                   ttaAAGTTGCCATGCAGCTTCCGGGGACGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTaa                    >  2:661469/1‑150 (MQ=255)
                                                                                                                                    taAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAAt                   >  1:3704763/1‑150 (MQ=255)
                                                                                                                                     aaaGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAAtt                  >  2:2841821/1‑150 (MQ=255)
                                                                                                                                         ttGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTcca              <  1:2880101/150‑1 (MQ=255)
                                                                                                                                           gCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTccacc            >  2:2384297/1‑150 (MQ=255)
                                                                                                                                           gCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTccacc            <  1:2190443/150‑1 (MQ=255)
                                                                                                                                              aTGCAGCTTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATg         <  2:1610578/150‑1 (MQ=255)
                                                                                                                                              aTGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATg         >  1:2315909/1‑150 (MQ=255)
                                                                                                                                               tGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGc        >  1:3418344/1‑150 (MQ=255)
                                                                                                                                                gCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGcc       >  1:2350255/1‑150 (MQ=255)
                                                                                                                                                  aGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGa     >  2:3532757/1‑150 (MQ=255)
                                                                                                                                                    cGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATc   <  1:1122665/150‑1 (MQ=255)
                                                                                                                                                    cGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATc   <  1:4568162/150‑1 (MQ=255)
                                                                                                                                                    cGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATc   <  1:3228696/150‑1 (MQ=255)
                                                                                                                                                     gTCCGGGGACGTGGTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGGTCTACCACGTTAATTTCCACCATGCCGATcc  >  2:3875531/1‑150 (MQ=255)
                                                                                                                                                     gTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATcc  >  2:978166/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAAATCACTGCTGACATGATGTTAACCCATCGCTATCCGTTCGCCACGCTGGCAGAAACCTACGAGCGCGATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCC  >  CP009273/4586712‑4587010

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: