Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,410,964 A→C 100% I390L (ATC→CTC)  ulaA → L‑ascorbate‑specific enzyme IIC permease component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,410,9640AC98.7% 1156.9 / ‑0.9 297I390L (ATC→CTC) ulaAL‑ascorbate‑specific enzyme IIC permease component of PTS
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base C (170/123);  total (173/124)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 302 total aligned reads displayed.

GTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATT  >  CP009273/4410817‑4411110
                                                                                                                                                   |                                                                                                                                                  
        ccTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCtt                                                                                                                                          >  2:3143749/1‑150 (MQ=255)
         cTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCttt                                                                                                                                         >  1:2111915/1‑150 (MQ=255)
          tGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTg                                                                                                                                        >  2:4953515/1‑150 (MQ=255)
            gTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGt                                                                                                                                      <  1:2007607/150‑1 (MQ=255)
            gTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGt                                                                                                                                      <  1:4483432/150‑1 (MQ=255)
              cGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTg                                                                                                                                    >  1:3453884/1‑150 (MQ=255)
                tcTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCg                                                                                                                                  >  1:4740063/2‑150 (MQ=255)
                  tGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTc                                                                                                                                <  1:4953839/150‑1 (MQ=255)
                  tGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTc                                                                                                                                >  1:1886314/1‑150 (MQ=255)
                     ggCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTgg                                                                                                                             <  1:3513545/150‑1 (MQ=255)
                       ctcctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGc                                                                                                                           >  1:2543746/1‑150 (MQ=255)
                       ctcctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGc                                                                                                                           >  1:644230/1‑150 (MQ=255)
                        tcctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTTGGCg                                                                                                                          <  2:214067/150‑1 (MQ=255)
                        tcctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCg                                                                                                                          >  1:212122/1‑150 (MQ=255)
                        tcctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCg                                                                                                                          >  2:1179024/1‑150 (MQ=255)
                         cctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCgg                                                                                                                         >  1:2787008/1‑150 (MQ=255)
                          ctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGt                                                                                                                        >  1:3216349/1‑150 (MQ=255)
                          ctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGt                                                                                                                        <  1:4414047/150‑1 (MQ=255)
                                  tGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCa                                                                                                                >  1:3141485/1‑150 (MQ=255)
                                   gATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCAc                                                                                                               >  1:2376141/1‑150 (MQ=255)
                                    attattCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCAcc                                                                                                              <  2:4342675/150‑1 (MQ=255)
                                      tattCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCgg                                                                                                            >  1:3790028/1‑150 (MQ=255)
                                      tattCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCgg                                                                                                            <  1:23187/150‑1 (MQ=255)
                                        ttCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTa                                                                                                          <  1:2388820/150‑1 (MQ=255)
                                        ttCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTa                                                                                                          <  2:2486277/150‑1 (MQ=255)
                                          ccTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGAGTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATg                                                                                                        >  2:2815685/1‑150 (MQ=255)
                                               cTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGc                                                                                                   >  2:889501/1‑150 (MQ=255)
                                                tttATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGcc                                                                                                  >  1:5151392/1‑150 (MQ=255)
                                                tttATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGcc                                                                                                  >  2:180370/1‑150 (MQ=255)
                                                 ttATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCt                                                                                                 >  1:112495/1‑150 (MQ=255)
                                                  tATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTg                                                                                                >  2:1204715/1‑150 (MQ=255)
                                                  tATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTg                                                                                                <  2:3109196/150‑1 (MQ=255)
                                                   cgcccGATGTTCTTCTCTCACGCCACCATCGGCGGGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTgg                                                                                               <  2:2743516/148‑1 (MQ=255)
                                                      ccGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATg                                                                                            <  2:2983258/150‑1 (MQ=255)
                                                       cGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATgg                                                                                           >  2:580833/1‑150 (MQ=255)
                                                        gaagTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGCAATCTATGGTTGCGTCTGGGCGGTGAAACACACCGGTATGAGTGCCTGGATggg                                                                                          >  2:162893/4‑150 (MQ=255)
                                                        gATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATggg                                                                                          <  2:4867636/150‑1 (MQ=255)
                                                        gATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATggg                                                                                          >  2:801010/1‑150 (MQ=255)
                                                          tGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCa                                                                                        <  2:473931/150‑1 (MQ=255)
                                                          tGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCa                                                                                        <  2:4112635/150‑1 (MQ=255)
                                                           gTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCAt                                                                                       >  1:98073/1‑150 (MQ=255)
                                                            ttcttcTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCTGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATg                                                                                      <  1:2212230/150‑1 (MQ=255)
                                                            ttcttcTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATg                                                                                      <  1:5002588/150‑1 (MQ=255)
                                                            ttcttcTCTAACGCCACCATCGGCGTGTTCGCTAACCACCTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCCCCGGTATGAGTGCCTGGATGGGCATg                                                                                      >  1:2736249/1‑150 (MQ=255)
                                                                 ctctAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGa                                                                                 >  2:1920095/1‑150 (MQ=255)
                                                                   ctAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACt                                                                               >  2:505870/1‑150 (MQ=255)
                                                                   ctAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACt                                                                               >  2:3842102/1‑150 (MQ=255)
                                                                   cgaaCGCCCCCCTCGGCGTGTTCGCGAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGACACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACt                                                                               <  2:4969568/148‑1 (MQ=255)
                                                                      aCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGt                                                                            >  1:3978818/1‑150 (MQ=255)
                                                                              aTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTg                                                                    <  2:2932131/150‑1 (MQ=255)
                                                                               tCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTgg                                                                   <  2:2655518/150‑1 (MQ=255)
                                                                                 ggCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCa                                                                 <  2:1338711/150‑1 (MQ=255)
                                                                                  gCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAc                                                                <  2:502239/150‑1 (MQ=255)
                                                                                   cGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGTAGATTTGTCTGGTGATGGGGATGCTCGAAAACTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAcc                                                               <  2:4169964/150‑1 (MQ=255)
                                                                                      gtTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGCTTCTGGCAccgcc                                                            >  2:4761909/1‑150 (MQ=255)
                                                                                      gtTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAccgcc                                                            <  2:3750663/150‑1 (MQ=255)
                                                                                      gtTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAccgcc                                                            <  1:3850311/150‑1 (MQ=255)
                                                                                          gcgaaCCACTTCGGCGGCTGGCGTGCGGCGCGGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAACCTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatg                                                        <  2:4631166/147‑1 (MQ=255)
                                                                                          gCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatg                                                        >  1:3650865/1‑150 (MQ=255)
                                                                                          gCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatg                                                        >  2:1758092/1‑150 (MQ=255)
                                                                                           cTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatga                                                       >  2:2590755/1‑150 (MQ=255)
                                                                                           cTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatga                                                       >  1:2624058/1‑150 (MQ=255)
                                                                                            tAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatgat                                                      <  2:552311/150‑1 (MQ=255)
                                                                                               ctaCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCa                                                   >  1:1761731/3‑150 (MQ=255)
                                                                                               ccACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCa                                                   >  2:232568/1‑150 (MQ=255)
                                                                                                 aCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCaaa                                                 <  1:194970/150‑2 (MQ=255)
                                                                                                  cTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAgg                                                >  1:3482390/1‑150 (MQ=255)
                                                                                                   ttCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGc                                               >  1:1649410/1‑150 (MQ=255)
                                                                                                   ttCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGc                                               >  1:3680198/1‑150 (MQ=255)
                                                                                                    tCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGct                                              >  1:1247890/1‑150 (MQ=255)
                                                                                                     cggcggCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCATGGctt                                             <  1:248272/150‑1 (MQ=255)
                                                                                                       gcggcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGcttct                                           <  1:2848365/150‑1 (MQ=255)
                                                                                                         ggcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCttcttc                                         >  1:4909309/1‑150 (MQ=255)
                                                                                                         ggcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCttcttc                                         <  2:2367695/150‑1 (MQ=255)
                                                                                                         ggcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCttcttc                                         <  2:2242883/150‑1 (MQ=255)
                                                                                                          gcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTtcttct                                        <  2:4936489/150‑1 (MQ=255)
                                                                                                            tGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTcc                                      <  1:888016/150‑1 (MQ=255)
                                                                                                            tGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTcc                                      >  2:3724401/1‑150 (MQ=255)
                                                                                                             ggCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCa                                     <  2:150463/150‑1 (MQ=255)
                                                                                                                gTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCg                                  <  2:2428483/150‑1 (MQ=255)
                                                                                                                gTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCg                                  <  1:2923598/150‑1 (MQ=255)
                                                                                                                 tGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCgg                                 <  2:4467333/150‑1 (MQ=255)
                                                                                                                  gcggcgCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGt                                <  1:1059278/150‑1 (MQ=255)
                                                                                                                  gcggcgCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGt                                >  2:4703587/1‑150 (MQ=255)
                                                                                                                   cggcgCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTa                               <  1:67214/150‑1 (MQ=255)
                                                                                                                     gcgcTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATc                             <  2:2215536/150‑1 (MQ=255)
                                                                                                                      cgcTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCg                            >  2:1960209/1‑150 (MQ=255)
                                                                                                                       gcTGAAGATTTGTCTGGTGATGGGGATGCTCTAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGc                           >  1:4797155/1‑150 (MQ=255)
                                                                                                                         tGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACAGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCt                         <  1:2695262/150‑1 (MQ=255)
                                                                                                                          gAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCtt                        <  2:2068775/150‑1 (MQ=255)
                                                                                                                           aaGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGAGGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCttt                       <  2:3069459/150‑1 (MQ=255)
                                                                                                                              aTTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATg                    >  2:405597/1‑150 (MQ=255)
                                                                                                                               tttGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATgg                   <  1:24378/150‑1 (MQ=255)
                                                                                                                                    cTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtc              >  1:3067122/1‑150 (MQ=255)
                                                                                                                                           tGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTg       <  2:3970324/150‑1 (MQ=255)
                                                                                                                                           tGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTg       <  1:4567731/150‑1 (MQ=255)
                                                                                                                                             gggATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTa     <  2:2322608/150‑1 (MQ=255)
                                                                                                                                             gggATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTa     >  1:791282/1‑150 (MQ=255)
                                                                                                                                               gatgctCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGTGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAAt   >  2:1876717/1‑150 (MQ=255)
                                                                                                                                               gatgctCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCAGTGTAAt   >  2:3495902/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATT  >  CP009273/4410817‑4411110

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: