Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,410,964 A→C 88.7% 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,9640AC88.7% 659.7 / 38.6 204I390L (ATC→CTC) ulaAL‑ascorbate‑specific enzyme IIC permease component of PTS
Reads supporting (aligned to +/- strand):  ref base A (13/11);  new base C (103/77);  total (116/88)
Fisher's exact test for biased strand distribution p-value = 8.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

GCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAA  >  CP009273/4410821‑4411108
                                                                                                                                               |                                                                                                                                                
gCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGa                                                                                                                                              >  1:462785/1‑148 (MQ=255)
 cATCCTGGTCGCCTGCGGCTCCTCGATCCTGTTTATTCCTGGCTTTATCCCGTTGTTCTTCTCTAACGCCACCTTCGGCGTGTTCGCTAACCACTTAGGCGGCTGGCGTGCGGCGCTGTAGATTTGTCTGGTGATGGGGATGATCGaa                                                                                                                                             <  1:1092187/148‑1 (MQ=255)
 cATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGaa                                                                                                                                             <  2:1028734/148‑1 (MQ=255)
 cATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGaa                                                                                                                                             <  1:3189130/148‑1 (MQ=255)
   tCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAAt                                                                                                                                           <  1:413828/148‑1 (MQ=255)
      tGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCtt                                                                                                                                        >  2:3150942/1‑148 (MQ=255)
      tGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCtt                                                                                                                                        >  2:3151065/1‑148 (MQ=255)
        gTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTg                                                                                                                                      >  2:902272/1‑148 (MQ=255)
         tCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTgg                                                                                                                                     <  1:154907/148‑1 (MQ=255)
              tGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCg                                                                                                                                <  2:3060147/148‑1 (MQ=255)
               gCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGt                                                                                                                               <  1:2858071/148‑1 (MQ=255)
               gCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGt                                                                                                                               >  1:116212/1‑148 (MQ=255)
                cGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGGTGGGGATGCTCGAAATCTTTGGTTGCGTc                                                                                                                              <  2:613744/148‑1 (MQ=255)
                  gCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTg                                                                                                                            <  2:2944943/148‑1 (MQ=255)
                    tcctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTggg                                                                                                                          <  2:153787/148‑1 (MQ=255)
                    tcctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTggg                                                                                                                          <  1:1747087/148‑1 (MQ=255)
                    tcctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTggg                                                                                                                          >  1:59531/1‑148 (MQ=255)
                     cctcGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGc                                                                                                                         >  1:815694/1‑148 (MQ=255)
                          aTCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGa                                                                                                                    >  1:3216994/1‑148 (MQ=255)
                            ccTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGaaa                                                                                                                  >  2:452558/1‑148 (MQ=255)
                            ccTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGaaa                                                                                                                  >  2:630190/1‑148 (MQ=255)
                            ccTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGaaa                                                                                                                  >  2:198856/1‑148 (MQ=255)
                              tGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACt                                                                                                                <  1:1618570/148‑1 (MQ=255)
                               gATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTc                                                                                                               >  2:1303843/1‑148 (MQ=255)
                                 ttattCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCAc                                                                                                             >  2:17784/1‑148 (MQ=255)
                                  tattCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCAcc                                                                                                            >  1:2760626/1‑148 (MQ=255)
                                   attCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCg                                                                                                           <  1:2096/148‑1 (MQ=255)
                                   attCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGNGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCg                                                                                                           <  1:1679799/148‑1 (MQ=255)
                                     tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGCTGCTCGAAATCTGTGGTTGCGTCTGGGCGGTGAAACTCACCGGt                                                                                                         >  2:2536320/1‑148 (MQ=255)
                                     tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGAGTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGt                                                                                                         >  1:442273/1‑148 (MQ=255)
                                        tGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATg                                                                                                      <  2:2241925/148‑1 (MQ=255)
                                         ggCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                     <  2:2246892/148‑1 (MQ=255)
                                         ggCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                     <  2:92102/148‑1 (MQ=255)
                                          gCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAg                                                                                                    >  2:1385189/1‑148 (MQ=255)
                                             ttATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTAGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGc                                                                                                 >  1:3147600/1‑148 (MQ=255)
                                                tCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTg                                                                                              >  2:2305368/1‑148 (MQ=255)
                                                tCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTg                                                                                              >  2:1036143/1‑148 (MQ=255)
                                                tCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTg                                                                                              >  1:642094/1‑148 (MQ=255)
                                                tCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTg                                                                                              >  1:2629588/1‑148 (MQ=255)
                                                tCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTg                                                                                              >  2:2207993/1‑148 (MQ=255)
                                                     aTGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATgg                                                                                         >  1:1408977/1‑148 (MQ=255)
                                                        ttcttcTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCa                                                                                      <  1:2966756/148‑1 (MQ=255)
                                                         tcttctCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCAt                                                                                     >  1:508961/1‑148 (MQ=255)
                                                          cttctCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATg                                                                                    >  1:3269542/1‑148 (MQ=255)
                                                           ttctCTAACGCCACCATCGGCGTGTTCGCTAAACACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTGTGGTGATTGGGATGCTCGAAATATTTTGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATgg                                                                                   >  1:1517689/1‑148 (MQ=255)
                                                            tctcTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGc                                                                                  >  2:1696865/1‑148 (MQ=255)
                                                            tctcTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGc                                                                                  >  1:1715794/1‑148 (MQ=255)
                                                               ctAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGa                                                                               <  2:719851/148‑1 (MQ=255)
                                                                 aaCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACt                                                                             <  2:202808/148‑1 (MQ=255)
                                                                 aaCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACt                                                                             <  2:1905499/148‑1 (MQ=255)
                                                                  aCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTg                                                                            <  1:1446721/148‑1 (MQ=255)
                                                                   cGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTCTGAGTGCCTGGATGGGCCTGGCGGagggg                                                                           >  2:2076624/1‑144 (MQ=255)
                                                                     ccaccaTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTc                                                                         >  2:2365538/1‑148 (MQ=255)
                                                                     ccaccaTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTc                                                                         >  2:794069/1‑148 (MQ=255)
                                                                     ccaccaTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTc                                                                         >  1:2943888/1‑148 (MQ=255)
                                                                     ccaccaTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTc                                                                         >  1:2066898/1‑148 (MQ=255)
                                                                      caccaTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGAc                                                                              >  2:2722428/1‑142 (MQ=255)
                                                                      caccaTCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGAc                                                                              <  1:2721931/142‑1 (MQ=255)
                                                                           tCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCt                                                                   >  2:693559/1‑148 (MQ=255)
                                                                           tCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCt                                                                   >  1:2221160/1‑148 (MQ=255)
                                                                            cGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTg                                                                  >  2:67937/1‑148 (MQ=255)
                                                                             ggCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTgg                                                                 <  2:1748958/148‑1 (MQ=255)
                                                                                 tgtTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCAcc                                                             >  2:1767036/1‑148 (MQ=255)
                                                                                       cTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgat                                                       >  1:1969536/1‑148 (MQ=255)
                                                                                         aaCCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatga                                                     >  1:1092434/1‑148 (MQ=255)
                                                                                          aCCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCGTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatgat                                                    >  2:1157893/1‑148 (MQ=255)
                                                                                           ccACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGatgatg                                                   <  1:2243081/148‑1 (MQ=255)
                                                                                           ccACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGatgatg                                                   >  2:1428507/1‑148 (MQ=255)
                                                                                            cACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGc                                                  >  2:887043/1‑148 (MQ=255)
                                                                                             aCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCa                                                 >  1:2457084/1‑148 (MQ=255)
                                                                                             aCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCa                                                 <  1:103000/148‑1 (MQ=255)
                                                                                                tCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAgg                                              >  1:173477/1‑148 (MQ=255)
                                                                                                 cggcggCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGc                                             <  1:1222145/148‑1 (MQ=255)
                                                                                                   gcggcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGctt                                           <  2:1083743/148‑1 (MQ=255)
                                                                                                      gcTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGcttctt                                        <  1:3150957/148‑1 (MQ=255)
                                                                                                         ggCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTtctc                                     <  1:887043/148‑1 (MQ=255)
                                                                                                           cGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCa                                   >  2:3227452/1‑148 (MQ=255)
                                                                                                           cGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCa                                   <  1:3305234/148‑1 (MQ=255)
                                                                                                             tGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATc                                 >  2:2274923/1‑148 (MQ=255)
                                                                                                                ggcgCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGt                              <  2:1103379/148‑1 (MQ=255)
                                                                                                                  cgcTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCCTCGGTAt                            >  1:404835/1‑148 (MQ=255)
                                                                                                                  cgcTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTAt                            <  1:2670454/148‑1 (MQ=255)
                                                                                                                  cgcTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTAt                            <  1:296213/148‑1 (MQ=255)
                                                                                                                    cTGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCg                          >  1:48135/1‑148 (MQ=255)
                                                                                                                     tGAAGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGc                         <  1:1485115/148‑1 (MQ=255)
                                                                                                                       aaGATTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTAACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGTCACCGCCGATGATGCAATGCTTCTTCTCCATCGGTATCGCCt                       <  2:248775/148‑1 (MQ=255)
                                                                                                                        aGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCtt                      >  2:2715775/1‑148 (MQ=255)
                                                                                                                            ttGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATg                  <  2:2221160/148‑1 (MQ=255)
                                                                                                                            ttGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATg                  >  1:1275397/1‑148 (MQ=255)
                                                                                                                            ttGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATg                  <  1:2228004/148‑1 (MQ=255)
                                                                                                                            ttGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATg                  >  2:3038710/1‑148 (MQ=255)
                                                                                                                                  ggTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCTCCTTTATGGCCGtc            <  2:371023/148‑1 (MQ=255)
                                                                                                                                    tGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcat          >  1:2318414/1‑148 (MQ=255)
                                                                                                                                    tGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcat          >  1:176198/1‑148 (MQ=255)
                                                                                                                                    tGATGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcat          >  1:1334550/1‑148 (MQ=255)
                                                                                                                                      aTGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcatca        <  2:508961/148‑1 (MQ=255)
                                                                                                                                      aTGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcatca        <  2:1211633/148‑1 (MQ=255)
                                                                                                                                      aTGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGtcatca        <  2:2171080/148‑1 (MQ=255)
                                                                                                                                       tGGGGATGCTCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTcatcat       >  2:1427819/1‑148 (MQ=255)
                                                                                                                                            atgctCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTaa  >  2:349144/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAA  >  CP009273/4410821‑4411108

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: