Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,034,118 G→T 31.2% intergenic (+122/+148) rrfA → / ← mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,034,1180GT31.2% 208.2 / 118.8 157intergenic (+122/+148)rrfA/mobB5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base G (70/38);  new base T (27/22);  total (97/60)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 160 total aligned reads displayed.

ATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTT  >  CP009273/4033970‑4034267
                                                                                                                                                    |                                                                                                                                                     
aTGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATGTGAACGTTGCGAAGCACCgg                                                                                                                                                      >  2:194683/1‑150 (MQ=17)
   cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGACAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                                                                                                                   <  2:4782716/150‑1 (MQ=17)
    gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCg                                                                                                                                                  >  2:2309317/1‑150 (MQ=17)
    gagaGGAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCGTTTCGTTTTATCTGTTGTGTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCg                                                                                                                                                  >  2:1023929/1‑150 (MQ=255)
      gagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGa                                                                                                                                                >  2:476887/1‑150 (MQ=17)
      gagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGa                                                                                                                                                >  2:410551/1‑150 (MQ=17)
       agTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAg                                                                                                                                               >  1:1282134/1‑150 (MQ=17)
          aGGGAACTGCCAGGCATCAAATAAAACTAAAGGCTCAGTCGGAAGACTGGGCCTGTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATGTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                                                                                                                            >  2:3625946/1‑150 (MQ=12)
          aGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                                                                                                                            >  1:4005655/1‑150 (MQ=32)
          aGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                                                                                                                            >  1:1340772/1‑150 (MQ=32)
              aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg                                                                                                                                        >  2:4453095/1‑150 (MQ=32)
               aCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCGGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCgg                                                                                                                                       >  2:132049/1‑150 (MQ=17)
                   ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAg                                                                                                                                   >  1:4374109/1‑150 (MQ=17)
                     aGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGa                                                                                                                                 >  2:1381306/1‑150 (MQ=17)
                     aGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGa                                                                                                                                 >  1:4319046/1‑150 (MQ=32)
                      ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGAc                                                                                                                                >  1:2482654/1‑150 (MQ=17)
                      ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGAc                                                                                                                                <  1:677251/150‑1 (MQ=17)
                         aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCGGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGGTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcc                                                                                                                             >  2:4386098/1‑150 (MQ=17)
                          tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGAcgccc                                                                                                                            <  2:3807590/150‑1 (MQ=17)
                           cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGTcgcccg                                                                                                                           >  1:728976/1‑150 (MQ=17)
                           cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccg                                                                                                                           >  1:677931/1‑150 (MQ=32)
                            aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgc                                                                                                                          >  1:2169641/1‑150 (MQ=32)
                            aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCGGAGGAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgc                                                                                                                          >  2:1230133/1‑150 (MQ=17)
                               tAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAt                                                                                                                       >  2:3045105/1‑150 (MQ=17)
                                aaaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATa                                                                                                                      >  1:2968998/1‑150 (MQ=32)
                                 aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaa                                                                                                                     >  1:2808023/1‑150 (MQ=17)
                                 aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaa                                                                                                                     >  1:4121358/1‑150 (MQ=32)
                                  caCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                                                                                                                    <  2:3012944/149‑1 (MQ=17)
                                   aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAGAAAc                                                                                                                   >  2:4034257/1‑150 (MQ=17)
                                   aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAAc                                                                                                                   >  1:4229245/1‑150 (MQ=32)
                                   aCGAAAGGCTCAGTCGGAAGACTGGGCCTGTCGTGTTATCTGTTGTTTGTCGTTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAGAAAc                                                                                                                   >  1:3467606/1‑150 (MQ=255)
                                    cGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt                                                                                                                  <  2:739376/150‑1 (MQ=17)
                                     gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTg                                                                                                                 >  2:2277824/1‑150 (MQ=32)
                                      actggCTCAGACGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCACCTGAGTAGGACAAAGCCGCGGGGAGCGTATTTGAACGTGGCGAAGCCACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGc                                                                                                                <  2:5104776/147‑1 (MQ=11)
                                      accgccTCCGTCGGAAGACGGGGCCTTTCTTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGCGCGTTTTTGAACGTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGc                                                                                                                <  1:4417590/145‑1 (MQ=11)
                                       aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGTACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc                                                                                                               >  2:1799793/1‑150 (MQ=17)
                                       aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc                                                                                                               <  2:595613/150‑1 (MQ=17)
                                        aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCa                                                                                                              >  1:2674754/1‑150 (MQ=17)
                                         ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAg                                                                                                             <  2:3881760/150‑1 (MQ=17)
                                          gCTCCGTCGGACGAGGGGGCCTTTCGTTTTATCTGTTGTTGGGCGGTGACCGCTCTCCTGTGTCGGACAAAGCCGCGGGGAGCGGATTTGACCGTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAgg                                                                                                            <  2:1553941/150‑1 (MQ=11)
                                          gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAgg                                                                                                            >  1:2382107/1‑150 (MQ=17)
                                           cTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGc                                                                                                           >  2:1951757/1‑150 (MQ=32)
                                            tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa                                                                                                          >  2:3378035/1‑150 (MQ=17)
                                            tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa                                                                                                          <  1:1155870/150‑1 (MQ=17)
                                             cAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAGAAACTGCCAGGCAt                                                                                                         >  1:4755691/1‑150 (MQ=17)
                                              aGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc                                                                                                        >  2:1656026/1‑150 (MQ=32)
                                              aGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc                                                                                                        >  1:2689571/1‑150 (MQ=32)
                                              aGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc                                                                                                        >  2:1640384/1‑150 (MQ=32)
                                               gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCa                                                                                                       <  1:3258626/150‑1 (MQ=17)
                                                 cGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCAGCaaa                                                                                                     >  1:4259962/1‑150 (MQ=17)
                                                  ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGACGCTACGGCCCGGAGGGTGGCGGGCAGGACGCGCGCCATAAACTGCCAGGCATCAAAt                                                                                                    <  2:3716932/150‑1 (MQ=9)
                                                  ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAAt                                                                                                    >  2:4373693/1‑150 (MQ=14)
                                                  ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAt                                                                                                    >  2:5271814/1‑150 (MQ=31)
                                                  ggAAGACTGGGCCTTTCGTTGTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCACCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCTAAt                                                                                                    >  2:5271764/1‑150 (MQ=9)
                                                   gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGTTGAACGCTCTCCTGAGTAGGACAAATCCGCGGGGTGCGTATTTGAACTTTGCTAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAtt                                                                                                   <  1:4753924/150‑1 (MQ=11)
                                                   gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTTTGTAGGACAAATCCGCCGGGCGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAtt                                                                                                   <  2:1608579/150‑1 (MQ=2)
                                                   gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAtt                                                                                                   >  2:4147980/1‑150 (MQ=11)
                                                   gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAtt                                                                                                   <  2:2537441/150‑1 (MQ=11)
                                                                                       gTGAACGCTCTCCTTAGTAGGAAAAATCCGCCGGTAGCGGATTTGAACTTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCa                                                               <  1:3940799/150‑1 (MQ=11)
                                                                                        tGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGACGGCCATCCTGACGGATGGCCTTTTTGCAt                                                              <  2:225761/150‑1 (MQ=14)
                                                                                         taccGCTCACCTTAGTCGGAAAAATCCGCCGGTAGCGTATTTGAACGTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCAtt                                                             <  1:1795863/147‑1 (MQ=11)
                                                                                             gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCAGACGGATGGCCGTTTTGGAGTGgcg                                                         >  1:4190701/1‑150 (MQ=11)
                                                                                             gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCCTTGgcg                                                         >  2:3099311/1‑150 (MQ=17)
                                                                                                  ccTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaa                                                    >  2:3893998/1‑150 (MQ=255)
                                                                                                   cTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCCTCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaa                                                   >  1:2077454/1‑150 (MQ=255)
                                                                                                       gTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaaaaa                                               >  2:846132/1‑150 (MQ=255)
                                                                                                       gTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCAGTGGCGCAGcaaaaaa                                               >  2:3520180/1‑150 (MQ=255)
                                                                                                        tAGGACAACTCCGCCGGGAGCGGATTTGAACGTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAAt                                              <  1:4387395/150‑1 (MQ=255)
                                                                                                        tAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAGAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGAAGGATGGCCATTTTGGATTGGCGCGGACAAAAAt                                              >  1:621439/1‑150 (MQ=11)
                                                                                                              aaaTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGa                                        >  2:1602760/1‑150 (MQ=255)
                                                                                                               aaTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGAt                                       >  2:2394102/1‑150 (MQ=255)
                                                                                                                aTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATg                                      <  1:250141/150‑1 (MQ=255)
                                                                                                                 tCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGGATTGGCGCAGAAAAAAATGCCGGATGc                                     >  2:4636036/1‑150 (MQ=255)
                                                                                                                   cgccgGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGa                                   <  1:4212624/150‑1 (MQ=255)
                                                                                                                    gccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGAc                                  >  1:5292547/1‑150 (MQ=255)
                                                                                                                     ccgGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTATTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACg                                 <  2:4986826/150‑1 (MQ=255)
                                                                                                                        ggAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCCGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGGCGCTg                              >  1:2995588/1‑150 (MQ=255)
                                                                                                                             ggATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGcgcg                         >  2:2659803/1‑150 (MQ=255)
                                                                                                                             ggATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGcgcg                         >  2:3938432/1‑150 (MQ=255)
                                                                                                                              gATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCCTCAAATTAAGCAGAAGGCCATCCAGAAGGTTGGCCGTTTTGGATTGGCGCAGAAAAAAATGCCGGGTGCGACGGTGGcgcgg                        >  1:2916986/1‑149 (MQ=17)
                                                                                                                              gATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGGAAAAAATGCCGGATGCGACGCTGGCGCGt                        >  1:350779/1‑150 (MQ=255)
                                                                                                                               aTTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGGTGCGACGCTGGcgcggc                       >  2:1142043/1‑148 (MQ=255)
                                                                                                                               aTTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTc                       >  2:1142119/1‑150 (MQ=255)
                                                                                                                                 ttGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCtt                     <  1:819782/150‑1 (MQ=255)
                                                                                                                                 ttGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCAAACTGACGGATGGCCTTTTTTGATTTGCGCAGAAAAAAATTACGGATGCCACGCAGGCGCGTCtt                     >  1:1217209/1‑150 (MQ=17)
                                                                                                                                   gAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTAt                   <  2:691890/150‑1 (MQ=255)
                                                                                                                                    aaCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCCTTGGCGCCGAAAAAAATGCCGGATGCGACGCTGGCGCGTcgtttc                  >  2:5047014/1‑145 (MQ=255)
                                                                                                                                      ggctGCCAACCACCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCATTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCa                <  1:1911875/147‑1 (MQ=21)
                                                                                                                                      ctttGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCa                <  1:5094001/148‑1 (MQ=255)
                                                                                                                                      cGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCa                >  2:3269204/1‑150 (MQ=255)
                                                                                                                                       gTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATTCCGGATGCGACGCTGGCGCGTCTTATCCAg               >  1:3020628/1‑150 (MQ=255)
                                                                                                                                       gTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAg               >  2:533621/1‑150 (MQ=255)
                                                                                                                                          gCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATTCCGGATGCGACGCTGGCGCGTCTTATCCAGCCt            <  2:2133545/150‑1 (MQ=255)
                                                                                                                                            gAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCCTCAAATTAAGCAGAAGGCCATCCAGACGGATGGCCATTTTGCATTGGCGCAGAAAAAAAAGCCGGATGCGACGCTGGCGCGTCTTATACAGCCGAc          >  1:3400318/1‑150 (MQ=18)
                                                                                                                                              aGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTActc        <  2:3466542/150‑1 (MQ=255)
                                                                                                                                              aGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGAGACGCTGGCGCGTCTTATCCAGCCTActc        <  1:3378035/150‑1 (MQ=255)
                                                                                                                                              aGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCCCCATAAACTGCCAGGCCTCAAAATAATCAGAAGGCCAACCAGACCGATGGCCTTTTTTGATTGGCGCAGAAAAAAAATCCCGAAGCGACGCAGGCGCGTCATATACAAccaaaac        >  1:1270678/1‑145 (MQ=17)
                                                                                                                                                cAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCtt      >  2:1250870/1‑150 (MQ=255)
                                                                                                                                                 aaCTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTg     <  2:2808023/150‑1 (MQ=255)
                                                                                                                                                 aaCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTg     <  2:72199/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTT  >  CP009273/4033970‑4034267

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