Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,034,118 G→T 21.8% 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,1180GT21.8% 279.7 / 64.4 141intergenic (+122/+148)rrfA/mobB5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base G (67/43);  new base T (18/13);  total (85/56)
Fisher's exact test for biased strand distribution p-value = 8.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

TGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTG  >  CP009273/4033971‑4034264
                                                                                                                                                   |                                                                                                                                                  
tGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACg                                                                                                                                                    >  1:2471709/1‑148 (MQ=32)
 gCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACgg                                                                                                                                                   >  1:2852577/1‑148 (MQ=32)
 gCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTGGTTTGTCGGTGAACGCTCTCCGGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACgg                                                                                                                                                   >  2:2407703/1‑148 (MQ=17)
    agagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                                                                                                                >  1:3204573/1‑148 (MQ=32)
    agagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                                                                                                                >  1:3204507/1‑148 (MQ=32)
        tAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAg                                                                                                                                            >  2:2325723/1‑148 (MQ=32)
          gggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAggg                                                                                                                                          <  2:2871346/148‑1 (MQ=32)
               cTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg                                                                                                                                     >  2:3630689/1‑148 (MQ=32)
                tGCCAGGCATCAAATAAAACGAAAGGCGCAGTCGGAATACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCCCTCCTGAGCAGGACAAAGCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCgg                                                                                                                                    <  2:499626/148‑1 (MQ=11)
                  ccAGCCATCAAATAAACCGAACCGCTCAGGCGGAAGACGGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTGGCCACGCAACGGCCCGGAGGGTGGCGGGc                                                                                                                                  <  2:2575666/148‑1 (MQ=11)
                     ggCATCACATAAACCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                               <  2:2375412/148‑1 (MQ=17)
                     ggCATCAAATAAACCGAACGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAgg                                                                                                                               <  2:3728703/148‑1 (MQ=12)
                     ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGCTTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                               >  2:3806830/1‑148 (MQ=17)
                      gCATCAAATAACACGACAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGa                                                                                                                              <  1:258126/148‑1 (MQ=12)
                         tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGTGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGGGGGGGGGAcgc                                                                                                                           >  1:1314325/1‑148 (MQ=255)
                         tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCTACGGCCCGGAGGGTGGCGGGCAGGAcgc                                                                                                                           >  2:575244/1‑148 (MQ=17)
                         acAAATAAAACGAACGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGAcgc                                                                                                                           <  1:3602907/147‑1 (MQ=17)
                          cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcc                                                                                                                          <  2:240216/148‑1 (MQ=32)
                             aTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGAcgcccgc                                                                                                                       >  1:2036564/1‑148 (MQ=17)
                             aTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGAcgcccgc                                                                                                                       >  1:2036548/1‑148 (MQ=17)
                             aTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGAcgcccgc                                                                                                                       >  1:1947800/1‑148 (MQ=17)
                             aTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgc                                                                                                                       >  2:2640217/1‑148 (MQ=32)
                              tAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgcc                                                                                                                      >  1:2185066/1‑148 (MQ=32)
                               aaaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgccc                                                                                                                     >  2:3346683/1‑147 (MQ=17)
                               aaaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgccc                                                                                                                     >  1:2738606/1‑147 (MQ=17)
                               aaaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCa                                                                                                                     >  2:3618134/1‑148 (MQ=32)
                                 aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATa                                                                                                                   >  1:1580391/1‑148 (MQ=32)
                                 aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATa                                                                                                                   >  2:476121/1‑148 (MQ=32)
                                 aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTGTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGGAGGACAAATCCGCCGGGAGCGGATGTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATa                                                                                                                   >  1:1580410/1‑148 (MQ=12)
                                 aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTGTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCCCCATa                                                                                                                   >  1:2950790/1‑148 (MQ=17)
                                  aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTaa                                                                                                                  >  1:877479/1‑148 (MQ=17)
                                  aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaa                                                                                                                  >  2:3279018/1‑148 (MQ=32)
                                  aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCGGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGGAGGACGCCCGCCATaa                                                                                                                  >  1:330735/1‑148 (MQ=17)
                                  aCGAAAGGCTCAGTCGGAAGACTGGGCCTGTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCGGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGGGGGGGGGGAGCCCGCCCTaa                                                                                                                  >  1:330733/1‑148 (MQ=11)
                                   cGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                                                                                                                 >  1:1228851/1‑148 (MQ=32)
                                   cGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                                                                                                                 >  2:3570483/1‑148 (MQ=32)
                                     aaaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt                                                                                                               >  2:3345782/1‑148 (MQ=32)
                                      aaGGCTCAGTCGGAAGACTGGGCCTTTCGTGTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTg                                                                                                              >  2:1798740/1‑148 (MQ=17)
                                       aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGc                                                                                                             >  2:1121037/1‑148 (MQ=17)
                                       aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGc                                                                                                             <  2:3097081/148‑1 (MQ=32)
                                        ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc                                                                                                            >  2:1496816/1‑148 (MQ=32)
                                        ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc                                                                                                            >  2:3138883/1‑148 (MQ=32)
                                         gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGACCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATGTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTACACTGCCa                                                                                                           >  2:3092470/1‑148 (MQ=255)
                                           tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACAGCCCgg                                                                                                         >  1:3057169/1‑148 (MQ=12)
                                           tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCGAAACTGCCAgg                                                                                                         >  1:3056641/1‑148 (MQ=17)
                                           tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAgg                                                                                                         >  1:3417589/1‑148 (MQ=32)
                                            cAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGACGCCCGCCATAAAATGCCAGGc                                                                                                        >  1:1684703/1‑148 (MQ=17)
                                            cAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCGAAACGGCCCGGc                                                                                                        >  1:1684560/1‑148 (MQ=255)
                                                cGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCa                                                                                                    >  2:3397293/1‑148 (MQ=32)
                                                cGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCa                                                                                                    >  2:1821548/1‑148 (MQ=32)
                                                 ggAAGCCGGGGCCTTTCGTTTTATCTGTTGTTTGGCGGTGAACGCGCTCCTGAGTCGGACAACTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaa                                                                                                   <  2:1120143/148‑1 (MQ=11)
                                                  gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                                  >  2:2252151/1‑148 (MQ=11)
                                                  gAAGACTGGGCATTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                                  <  2:1253328/148‑1 (MQ=11)
                                                                                        gaccGCTCTCCTGAGTAGCAAAAACCCCCGGGGAGCGGATTTGAAGTTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCa                                                            <  1:2721734/145‑1 (MQ=11)
                                                                                         aaCGCTCTCCTGAGTAGGACAACTCCGCCGGGAGCGGATTGGAACGTGGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCAt                                                           <  2:1338163/148‑1 (MQ=9)
                                                                                          aCGCTCACCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCATTTTGCAtt                                                          >  2:1981759/1‑148 (MQ=17)
                                                                                            gatatcCTGAGTAGCAAAAACCCGCGGGGAGCGTATTTTAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACCGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTgg                                                        <  1:2481370/144‑1 (MQ=11)
                                                                                            gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTgg                                                        <  1:3473593/148‑1 (MQ=33)
                                                                                            gCTCGCCTCAGGAGGCCACACCCGCCGGGACCGGACTTGAACGTCGCGAACCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTgg                                                        <  1:3708644/148‑1 (MQ=11)
                                                                                             ctctCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgc                                                       >  2:2109223/1‑148 (MQ=21)
                                                                                              tctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGTCGGATGGCCTTTTTGCATTgggg                                                      >  1:457304/1‑146 (MQ=255)
                                                                                               ctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCCTTGgcgc                                                     >  2:2992567/1‑148 (MQ=18)
                                                                                               ctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCCTTGgcgc                                                     >  2:2992672/1‑148 (MQ=18)
                                                                                               ctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACGGCCAGGCATCCAATTAAGCAGAAGGCCATCCTGACGGCTGGGCTTTTTTGATTGgcgc                                                     >  1:3047663/1‑148 (MQ=11)
                                                                                               ctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGGATTGgcgc                                                     >  1:2120458/1‑148 (MQ=21)
                                                                                               cgccTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgcgc                                                     <  2:3651317/146‑1 (MQ=25)
                                                                                                tcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCa                                                    >  2:3793292/1‑148 (MQ=37)
                                                                                                tcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCa                                                    >  2:3793540/1‑148 (MQ=37)
                                                                                                 ccTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCCACCAGACGGATGGCCGTTTTGGATTGgggccg                                                   >  1:3683202/1‑146 (MQ=11)
                                                                                                 ccTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCAAAACTGCCAGGCCTCCAATTAATCAGAAGGCCATCCGGACGGCAGGCCGTTTTTGATTGgcgccg                                                   >  1:3669955/1‑146 (MQ=11)
                                                                                                  cTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCCTCCTGACGGATGGCCTTTTTGCATTGGCGCAGa                                                  >  1:2229551/1‑148 (MQ=38)
                                                                                                     aGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCAACCTGACGGCTGGCCTTTTTGGATTGGCGCAGaaaa                                               >  1:548312/1‑148 (MQ=255)
                                                                                                     aGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCAGACGGATGGCCTTTTTTCAGTGGCGCAGaaaa                                               >  1:3475282/1‑148 (MQ=255)
                                                                                                          gACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATg                                          <  2:3243927/148‑1 (MQ=255)
                                                                                                           aCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGc                                         >  1:2556544/1‑148 (MQ=255)
                                                                                                            cAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGcc                                        <  2:2471709/148‑1 (MQ=255)
                                                                                                               aTCCGCCGGGAGCGGATTTGAACTTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGa                                     <  2:3542539/148‑1 (MQ=255)
                                                                                                               aTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGa                                     <  2:2852577/148‑1 (MQ=255)
                                                                                                                tCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGAt                                    <  2:1517492/148‑1 (MQ=255)
                                                                                                                      gggAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACg                              <  2:2505627/148‑1 (MQ=255)
                                                                                                                      gggAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAACTTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACg                              <  2:1247388/148‑1 (MQ=255)
                                                                                                                      ggcaGCGTATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACg                              <  1:862046/145‑1 (MQ=255)
                                                                                                                       ggAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGc                             >  2:2715458/1‑148 (MQ=255)
                                                                                                                              aTTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGACGCCCGCCCTAAACTGCCAGGCATCAAAATAAGCAGAATGCCCTCCTGACGGATGGCCTTTTTTGATTGGGGCAGGAAAAAATTCCGGCTGCGCCGCGGGcgcg                      >  1:551361/1‑148 (MQ=14)
                                                                                                                               tttGAAGGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGt                     <  1:3329823/148‑1 (MQ=255)
                                                                                                                                ttGAACGTTGCGACCCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTc                    <  1:3346683/148‑1 (MQ=255)
                                                                                                                                ttGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGTGACGCTGGCGCGTc                    >  2:1518763/1‑148 (MQ=255)
                                                                                                                                ttGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTc                    <  2:330735/148‑1 (MQ=255)
                                                                                                                                ttGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTc                    <  2:2332886/148‑1 (MQ=255)
                                                                                                                                ttGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTc                    <  2:3429460/148‑1 (MQ=255)
                                                                                                                                ttGAACGTGGCGACCCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTc                    <  2:330733/148‑1 (MQ=255)
                                                                                                                                 tGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCt                   <  1:3570483/148‑1 (MQ=255)
                                                                                                                                  caccGTTGCCAACCCACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCtt                  <  1:1798624/145‑1 (MQ=255)
                                                                                                                                   aaCGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTTGATTGGGGCCGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTa                 >  1:1231650/1‑148 (MQ=18)
                                                                                                                                   aaCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTa                 >  2:2475209/1‑148 (MQ=255)
                                                                                                                                              ccAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGTCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTAct      <  1:476121/147‑1 (MQ=255)
                                                                                                                                               cAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTActc     <  2:2535228/148‑1 (MQ=255)
                                                                                                                                                aaCTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAACTTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGAGCGTCTTATCCAGCCTActct    <  2:777558/148‑1 (MQ=255)
                                                                                                                                                aaCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTActct    <  2:2185066/148‑1 (MQ=255)
                                                                                                                                                  cGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTg  <  2:1314325/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTG  >  CP009273/4033971‑4034264

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: