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,9640AC99.6% 855.9 / ‑5.4 212I390L (ATC→CTC) ulaAL‑ascorbate‑specific enzyme IIC permease component of PTS
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (127/84);  total (128/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-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 213 total aligned reads displayed.

CATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTA  >  CP009273/4410822‑4411107
                                                                                                                                              |                                                                                                                                               
cATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGaa                                                                                                                                            >  1:2927043/1‑148 (MQ=255)
  tCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAAt                                                                                                                                          <  2:3463008/148‑1 (MQ=255)
   ccTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATc                                                                                                                                         >  1:3473811/1‑148 (MQ=255)
   ccTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTAGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATc                                                                                                                                         >  1:793837/1‑148 (MQ=255)
    cTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCt                                                                                                                                        <  2:1712767/148‑1 (MQ=255)
      ggTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCttt                                                                                                                                      <  1:78587/148‑1 (MQ=255)
      ggTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCttt                                                                                                                                      <  1:1398970/148‑1 (MQ=255)
           ccTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTg                                                                                                                                 <  2:2147429/148‑1 (MQ=255)
               cGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTc                                                                                                                             <  1:740542/148‑1 (MQ=255)
                ggCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCt                                                                                                                            <  2:1631596/148‑1 (MQ=255)
                 gCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTg                                                                                                                           >  1:1267080/1‑148 (MQ=255)
                 gCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTg                                                                                                                           >  1:2752678/1‑148 (MQ=255)
                 gCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTg                                                                                                                           >  2:2543577/1‑148 (MQ=255)
                  ctcctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTgg                                                                                                                          >  2:937190/1‑148 (MQ=255)
                  ctcctcGATCCTGAGTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTgg                                                                                                                          <  2:2574894/148‑1 (MQ=255)
                    cctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGc                                                                                                                        >  2:3423020/1‑148 (MQ=255)
                    cctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGc                                                                                                                        >  2:3422940/1‑148 (MQ=255)
                     ctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCg                                                                                                                       >  2:227330/1‑148 (MQ=255)
                     ctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTCGCGTGCGGCGCTGAAGATTTGACAGGTGAAGGGGATGCTCGAAATCTTGGGTTGCGTCTCGGCg                                                                                                                       >  2:2154907/1‑148 (MQ=255)
                      tcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCgg                                                                                                                      <  2:2421649/148‑1 (MQ=255)
                      tcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCgg                                                                                                                      <  1:3529305/148‑1 (MQ=255)
                          tCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGaa                                                                                                                  >  2:1990487/1‑148 (MQ=255)
                           ccTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGaaa                                                                                                                 >  1:2094622/1‑148 (MQ=255)
                           ccTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGaaa                                                                                                                 >  2:10341/1‑148 (MQ=255)
                           ccTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGaaa                                                                                                                 >  2:2732421/1‑148 (MQ=255)
                             tGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACt                                                                                                               >  2:3585294/1‑148 (MQ=255)
                               attattCCTGGCTTTATCCCGTTGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTCACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGTTGGGGATGCGCGAAATCTTTGGTGGCGTGTGGGCGGTGAAACTCa                                                                                                             >  2:2215081/1‑148 (MQ=255)
                               attattCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCa                                                                                                             <  1:270861/148‑1 (MQ=255)
                                ttattCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCAc                                                                                                            <  2:2793086/148‑1 (MQ=255)
                                   ttCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCgg                                                                                                         >  2:1532616/1‑148 (MQ=255)
                                   ttCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCgg                                                                                                         >  1:1968869/1‑148 (MQ=255)
                                    tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGt                                                                                                        >  2:1503245/1‑148 (MQ=255)
                                    tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGt                                                                                                        >  2:2195492/1‑148 (MQ=255)
                                      cTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCTCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTAt                                                                                                      <  1:2506766/148‑1 (MQ=255)
                                      cTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTAt                                                                                                      <  2:1830571/148‑1 (MQ=255)
                                        ggCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                    >  1:3657036/1‑148 (MQ=255)
                                         gCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAg                                                                                                   <  1:1742663/148‑1 (MQ=255)
                                          cTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGt                                                                                                  >  1:1510351/1‑148 (MQ=255)
                                          cTTTATCCCGATGTTCTTCTCTAACGACACCATCGGCGTGTTCGCAAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATATGTCTGGTGATGGGGATGCTCGAAAACTTTGGTTCCGTCAGGGCGGTGAAACTCACCGGTATGAGt                                                                                                  >  1:1230175/1‑148 (MQ=255)
                                            ttATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGc                                                                                                <  2:760005/148‑1 (MQ=255)
                                              aTCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCt                                                                                              >  1:897269/1‑148 (MQ=255)
                                                     tGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATggg                                                                                       >  1:588073/1‑148 (MQ=255)
                                                     tGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATggg                                                                                       >  1:2795907/1‑148 (MQ=255)
                                                      gTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGc                                                                                      >  1:3260204/1‑148 (MQ=255)
                                                       ttcttcTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCa                                                                                     >  1:3247855/1‑148 (MQ=255)
                                                       ttcttcTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCa                                                                                     <  1:2946653/148‑1 (MQ=255)
                                                       ttcttcTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCa                                                                                     >  1:2675585/1‑148 (MQ=255)
                                                          ttctCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATgg                                                                                  <  1:1008652/148‑1 (MQ=255)
                                                          ttctCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATgg                                                                                  <  1:3819634/148‑1 (MQ=255)
                                                               tAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGAGTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGAc                                                                             >  2:550134/1‑148 (MQ=255)
                                                                aaCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACt                                                                            >  1:1202783/1‑148 (MQ=255)
                                                                aaCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACt                                                                            <  2:2704346/148‑1 (MQ=255)
                                                                aaCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACt                                                                            >  1:2922004/1‑148 (MQ=255)
                                                                 aCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTg                                                                           >  1:2383212/1‑148 (MQ=255)
                                                                 aCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATGTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTg                                                                           >  2:2005584/1‑148 (MQ=255)
                                                                   gCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGt                                                                         >  2:2061524/1‑148 (MQ=255)
                                                                     caccaTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCg                                                                       >  1:2039502/1‑148 (MQ=255)
                                                                        caTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGAtt                                                                    <  2:2639895/148‑1 (MQ=255)
                                                                          tCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTTCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCt                                                                  >  1:2917659/1‑148 (MQ=255)
                                                                          tCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTTCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCt                                                                  >  1:3739635/1‑148 (MQ=255)
                                                                          tCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTACGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCt                                                                  <  2:781113/148‑1 (MQ=255)
                                                                           cGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTTTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTg                                                                 >  1:955579/1‑148 (MQ=255)
                                                                           cGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGCTGGGCATGGCGGACTGGTCGATTCTg                                                                 >  2:1077181/1‑148 (MQ=255)
                                                                           cGGCGTGTTCGCTAACCACATCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTg                                                                 >  1:953459/1‑148 (MQ=255)
                                                                            ggCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTTGATGGGCATGGCGGACTGGTCGATTCTgg                                                                <  2:2821623/148‑1 (MQ=255)
                                                                                tgtTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAcc                                                            >  2:3362275/1‑148 (MQ=255)
                                                                                tgtTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAcc                                                            >  2:1461063/1‑148 (MQ=255)
                                                                                 gtTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCCTGGCGGACTGGTCGATTCTGGCAccg                                                           >  2:3636600/1‑148 (MQ=255)
                                                                                    cGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAccgccg                                                        <  1:1773101/148‑1 (MQ=255)
                                                                                    cGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAccgccg                                                        <  1:1806477/148‑1 (MQ=255)
                                                                                    cGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAccgccg                                                        <  1:528886/148‑1 (MQ=255)
                                                                                     gCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCga                                                       <  1:1281562/148‑1 (MQ=255)
                                                                                      cTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgat                                                      >  1:3576013/1‑148 (MQ=255)
                                                                                       tAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatg                                                     <  2:167947/148‑1 (MQ=255)
                                                                                        aaCCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatga                                                    >  1:226665/1‑148 (MQ=255)
                                                                                         aCCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatgat                                                   <  2:325681/148‑1 (MQ=255)
                                                                                             cTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCaa                                               >  1:1003850/1‑148 (MQ=255)
                                                                                             cTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCaa                                               <  1:663174/148‑1 (MQ=255)
                                                                                                  gcggcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGctt                                          <  1:778410/148‑1 (MQ=255)
                                                                                                   cggcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGcttc                                         <  2:885651/148‑1 (MQ=255)
                                                                                                     gcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGcttctt                                       <  1:2601826/148‑1 (MQ=255)
                                                                                                        ggCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTtctc                                    <  2:2017379/148‑1 (MQ=255)
                                                                                                         gCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTcc                                   >  1:2441685/1‑148 (MQ=255)
                                                                                                            tGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATc                                <  2:1321998/148‑1 (MQ=255)
                                                                                                                 cgcTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTAt                           <  2:3783383/148‑1 (MQ=255)
                                                                                                                     gAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGcc                       <  1:1657415/148‑1 (MQ=255)
                                                                                                                      aaGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCt                      <  2:902982/148‑1 (MQ=255)
                                                                                                                        gATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCttt                    >  2:374302/1‑148 (MQ=255)
                                                                                                                        gATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCttt                    >  2:2294577/1‑148 (MQ=255)
                                                                                                                         aTTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTa                   >  1:3729385/1‑148 (MQ=255)
                                                                                                                           ttGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATg                 >  2:609201/1‑148 (MQ=255)
                                                                                                                                   tGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcat         <  2:3528557/148‑1 (MQ=255)
                                                                                                                                    gATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcatc        <  1:2641639/148‑1 (MQ=255)
                                                                                                                                    gATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcatc        <  1:2641906/148‑1 (MQ=255)
                                                                                                                                     aTGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAGGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcatca       <  1:1305068/148‑1 (MQ=255)
                                                                                                                                      tGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTcatcat      <  2:2039502/148‑1 (MQ=255)
                                                                                                                                      tGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTcatcat      <  2:1182648/148‑1 (MQ=255)
                                                                                                                                      tGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTcatcat      <  2:3328282/148‑1 (MQ=255)
                                                                                                                                      tGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTcatcat      <  2:3089141/148‑1 (MQ=255)
                                                                                                                                         ggATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTATACATCGGTATCGCCTTTATGGCCGTCCTCATTGt   >  1:668552/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTA  >  CP009273/4410822‑4411107

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: