Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,097,822 G→T 92.5% E115* (GAA→TAA)  pfkA → 6‑phosphofructokinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,097,8220GT92.5% 1015.8 / 72.0 332E115* (GAA→TAA) pfkA6‑phosphofructokinase I
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base T (144/163);  total (155/177)
Fisher's exact test for biased strand distribution p-value = 8.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 333 total aligned reads displayed.

GGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATT  >  CP009273/4097675‑4097971
                                                                                                                                                   |                                                                                                                                                     
ggTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTaa                                                                                                                                                     >  2:3357579/1‑150 (MQ=255)
ggTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGaa                                                                                                                                                     >  1:1329740/1‑150 (MQ=255)
  tACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAAt                                                                                                                                                   >  1:113664/1‑150 (MQ=255)
    cGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATgg                                                                                                                                                 <  2:2240857/150‑1 (MQ=255)
            ggTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCg                                                                                                                                         >  1:4285815/1‑150 (MQ=255)
            ggTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCg                                                                                                                                         >  1:930118/1‑150 (MQ=255)
            ggTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGATATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCg                                                                                                                                         >  2:199948/1‑150 (MQ=255)
             gTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGt                                                                                                                                        <  1:3331172/150‑1 (MQ=255)
              ttCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTg                                                                                                                                       >  1:1545668/1‑150 (MQ=255)
              ttCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTg                                                                                                                                       >  1:4106685/1‑150 (MQ=255)
               tCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGc                                                                                                                                      >  1:4836254/1‑150 (MQ=255)
                 tGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCAt                                                                                                                                    >  2:4990033/1‑150 (MQ=255)
                 tGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCAt                                                                                                                                    <  2:1045591/150‑1 (MQ=255)
                  gcgcGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATc                                                                                                                                   <  1:116410/150‑1 (MQ=255)
                  gcgcGGTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATc                                                                                                                                   <  2:3158912/150‑1 (MQ=255)
                    gcgTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCgg                                                                                                                                 <  1:4488964/150‑1 (MQ=255)
                      gTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTc                                                                                                                               <  1:5124555/150‑1 (MQ=255)
                      gTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTc                                                                                                                               <  1:3635330/150‑1 (MQ=255)
                        ttCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTg                                                                                                                             >  1:2064335/1‑150 (MQ=255)
                          cccGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGcc                                                                                                                           >  1:4290537/1‑150 (MQ=255)
                          cccGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGcc                                                                                                                           <  1:2586170/150‑1 (MQ=255)
                            cGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCgg                                                                                                                         <  2:4870341/150‑1 (MQ=255)
                             ggAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCggg                                                                                                                        <  1:3439582/150‑1 (MQ=255)
                              gcagtCCGCGACGAGAACATCCGCGCCGTGGCTATCGAACACCTGCAACAACGTGGGATCGACGCGCGGGTGGTTATCGGCGGGGCCGGTTCCTACCTGGGTGCAATGCGTCTGACCTAACTGGGCTTCCCGTGCATCGGTCTGCCGGGc                                                                                                                       <  2:4401161/146‑1 (MQ=255)
                               aaTTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCa                                                                                                                      >  2:2764235/1‑150 (MQ=255)
                               aaTTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCa                                                                                                                      >  1:4627936/1‑150 (MQ=255)
                                 ttCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACt                                                                                                                    <  1:1760794/150‑1 (MQ=255)
                                    cgcgACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATc                                                                                                                 <  2:2841061/150‑1 (MQ=255)
                                     gcgACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTATATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCg                                                                                                                <  1:1915605/150‑1 (MQ=255)
                                          gagaACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAAc                                                                                                           >  2:1724343/1‑150 (MQ=255)
                                            gaACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGa                                                                                                         >  2:2615201/1‑150 (MQ=255)
                                              aCATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACa                                                                                                       >  2:206425/1‑150 (MQ=255)
                                               cATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACAt                                                                                                      >  2:1118419/1‑150 (MQ=255)
                                                   cgcgCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCaaa                                                                                                  <  1:1469773/150‑1 (MQ=255)
                                                    gcgcCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAg                                                                                                 <  1:629393/150‑1 (MQ=255)
                                                     cgcCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAgg                                                                                                <  1:699328/150‑1 (MQ=255)
                                                      gcCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGc                                                                                               <  2:3737852/150‑1 (MQ=255)
                                                        cGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCAc                                                                                             >  2:3297836/1‑150 (MQ=255)
                                                          tGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTg                                                                                           <  1:5002657/150‑1 (MQ=255)
                                                           ggCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGa                                                                                          <  1:4050093/150‑1 (MQ=255)
                                                             cTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACt                                                                                        <  1:4141021/150‑1 (MQ=255)
                                                             cTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACt                                                                                        <  2:3339108/150‑1 (MQ=255)
                                                               aTCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTac                                                                                      >  2:983392/1‑150 (MQ=255)
                                                                  gAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACt                                                                                   >  1:5020699/1‑150 (MQ=255)
                                                                  gAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACt                                                                                   <  2:4074046/150‑1 (MQ=255)
                                                                   aaaaCCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACtc                                                                                  >  2:2983351/1‑149 (MQ=255)
                                                                    aaaCCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTAt                                                                                 >  2:2228556/1‑150 (MQ=255)
                                                                     aaCCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATc                                                                                >  2:3039282/1‑150 (MQ=255)
                                                                        cTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGt                                                                             >  1:4930617/1‑150 (MQ=255)
                                                                           aaaaaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTttc                                                                          >  1:689134/1‑150 (MQ=255)
                                                                              aaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTttcttc                                                                       >  1:1026156/1‑150 (MQ=255)
                                                                              aaaCGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTttcttc                                                                       <  1:3597260/150‑1 (MQ=255)
                                                                                   tGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTgc                                                                  <  2:3251631/150‑1 (MQ=255)
                                                                                    ggTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTgcg                                                                 >  1:2660695/1‑150 (MQ=255)
                                                                                    ggTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTgcg                                                                 <  1:3713032/150‑1 (MQ=255)
                                                                                    ggTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTgcg                                                                 <  1:4785560/150‑1 (MQ=255)
                                                                                     gTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTgcgc                                                                >  1:848954/1‑150 (MQ=255)
                                                                                     gTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTgcgc                                                                >  1:2192207/1‑150 (MQ=255)
                                                                                      tATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCt                                                               >  2:4013245/1‑150 (MQ=255)
                                                                                      tATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCt                                                               <  1:4356058/150‑1 (MQ=255)
                                                                                         cGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAg                                                            >  2:4307237/1‑150 (MQ=255)
                                                                                           aCGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCCAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCa                                                          <  2:3429922/150‑1 (MQ=255)
                                                                                            cgcgCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCAc                                                         >  1:318637/1‑150 (MQ=255)
                                                                                               gcTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGt                                                      <  2:3822388/150‑1 (MQ=255)
                                                                                                cTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGtt                                                     <  1:3819241/150‑1 (MQ=255)
                                                                                                cTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGtt                                                     >  1:3996528/1‑150 (MQ=255)
                                                                                                    tggtTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAg                                                 <  2:3599060/150‑1 (MQ=255)
                                                                                                    tggtTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAg                                                 <  2:4853331/150‑1 (MQ=255)
                                                                                                        tATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGc                                             <  2:2139217/150‑1 (MQ=255)
                                                                                                         aTCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCg                                            <  2:2674467/150‑1 (MQ=255)
                                                                                                          tCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGa                                           <  1:3349904/150‑1 (MQ=255)
                                                                                                            ggcggTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATc                                         <  2:3897330/150‑1 (MQ=255)
                                                                                                            ggcggTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATc                                         <  2:1329740/150‑1 (MQ=255)
                                                                                                             gcggTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCg                                        >  1:5059724/1‑150 (MQ=255)
                                                                                                              cggTGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGa                                       <  2:1186593/150‑1 (MQ=255)
                                                                                                                 tGACGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCg                                    >  1:2303286/1‑150 (MQ=255)
                                                                                                                    cGGTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCt                                 <  2:417953/150‑1 (MQ=255)
                                                                                                                     ggTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTg                                <  2:2045739/150‑1 (MQ=255)
                                                                                                                      gTTCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGc                               >  1:4502567/1‑150 (MQ=255)
                                                                                                                        tCCTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGt                             <  1:3064852/150‑1 (MQ=255)
                                                                                                                         ccTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTg                            <  1:2861642/150‑1 (MQ=255)
                                                                                                                         ccTACATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTg                            <  2:2317439/150‑1 (MQ=255)
                                                                                                                             cATGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGacac                        >  1:1392322/1‑150 (MQ=255)
                                                                                                                              aTGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACAcc                       <  1:2698952/150‑1 (MQ=255)
                                                                                                                              aTGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACAcc                       >  2:3519506/1‑150 (MQ=255)
                                                                                                                              aTGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACAcc                       >  2:2509323/1‑150 (MQ=255)
                                                                                                                              aTGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACAcc                       <  1:3148329/150‑1 (MQ=255)
                                                                                                                              aTGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACAcc                       <  1:2788257/150‑1 (MQ=255)
                                                                                                                               tGGGTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACct                      >  1:3632866/1‑150 (MQ=255)
                                                                                                                                  gTGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCtctt                   <  2:4627936/150‑1 (MQ=255)
                                                                                                                                   tGCACTGCGGCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCCCTAGCGACAACGACATCAAAGGCACTGACTACACTCTCGGTGTCTTCACTGCGCTGAGCACCGGTCTCGAAGCGATCGACCGTCTCCGTGACCCCtcttc                  >  2:4696485/1‑150 (MQ=255)
                                                                                                                                   tGCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCtcttc                  >  1:2884892/1‑150 (MQ=255)
                                                                                                                                    gCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCtcttct                 >  2:1232839/1‑150 (MQ=255)
                                                                                                                                    gCAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCtcttct                 <  1:1687673/150‑1 (MQ=255)
                                                                                                                                     cAATGCGTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCtcttctt                >  2:1877365/1‑150 (MQ=255)
                                                                                                                                           gTCTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCAcc          >  2:668166/1‑150 (MQ=255)
                                                                                                                                             cTGACCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAg        <  2:799297/150‑1 (MQ=255)
                                                                                                                                                aCCTAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGt     <  2:589718/150‑1 (MQ=255)
                                                                                                                                                aCCTAAATGGGCGTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGt     <  2:5123678/150‑1 (MQ=255)
                                                                                                                                                   gAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTAtt  >  1:4875365/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTGAAAAAACGTGGTATCGACGCGCTGGTGGTTATCGGCGGTGACGGTTCCTACATGGGTGCAATGCGTCTGACCGAAATGGGCTTCCCGTGCATCGGTCTGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATT  >  CP009273/4097675‑4097971

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: