Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,034,118 G→T 17.1% 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,1180GT17.1% 195.3 / 28.7 92intergenic (+122/+148)rrfA/mobB5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base G (51/24);  new base T (10/6);  total (62/30)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTG  >  CP009273/4033976‑4034264
                                                                                                                                              |                                                                                                                                                  
gagTAGGGAACTGCCAGGCATCAAATAAATCGAAAGGCTCAGTCGGAAGACTGGGCCTGTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTCGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGTAGCAACGGCCCg                                                                                                                                               >  2:228897/1‑148 (MQ=255)
gagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTTTTTGTCGGTTAACGCTCTAATTAGTAGGACAAAACCGCCCGGAGCGGATTTGAAAGTTGCGAAGCAAACGCCCg                                                                                                                                               >  1:2734004/1‑148 (MQ=11)
 agTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCgg                                                                                                                                              >  2:466850/1‑148 (MQ=32)
   tAGGGAACTGCCAGGCATCAATTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGTACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAg                                                                                                                                            >  2:3155202/1‑148 (MQ=17)
   tAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAg                                                                                                                                            >  2:3155582/1‑148 (MQ=32)
   tAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTGTTATCTGGTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAg                                                                                                                                            >  2:1671664/1‑148 (MQ=17)
       gAACTGCCAGGCATAAAATAAACCGAACGGCACAGTCGGAAGAGTGGGCCTTTCTTTTTATCTGTTGTTTGTCGGTTAACGCTCTCCTGAGTAGGACAAATCCGCCGGGCGCGGATTTTAACTTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                                                        <  2:3310762/148‑1 (MQ=11)
          cTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg                                                                                                                                     >  2:3246498/1‑148 (MQ=32)
          cTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg                                                                                                                                     >  1:2227801/1‑148 (MQ=32)
          cTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg                                                                                                                                     >  1:2230544/1‑148 (MQ=32)
               aGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAg                                                                                                                                >  1:659795/1‑148 (MQ=32)
                ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGcgggcggg                                                                                                                               >  1:64118/1‑148 (MQ=17)
                ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                               >  2:860657/1‑148 (MQ=32)
                 gCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGGGGggagga                                                                                                                              >  1:1035112/1‑148 (MQ=17)
                 gCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGa                                                                                                                              >  1:893579/1‑148 (MQ=32)
                    tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgc                                                                                                                           >  1:835377/1‑148 (MQ=32)
                    tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCGGAGTAGGACAAATCCGCCGGGAGCGGAGTTGAACGTTGCGACGCCACGGCGCGGAGGGTGGCGGGCAGGAcgc                                                                                                                           >  1:835362/1‑148 (MQ=11)
                    tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCGTTCGTTTTATCTGTGGTTGGGCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgc                                                                                                                           >  2:3195336/1‑148 (MQ=255)
                       aaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccg                                                                                                                        >  1:1471796/1‑148 (MQ=32)
                           aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAt                                                                                                                    >  1:827617/1‑148 (MQ=17)
                           aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgccct                                                                                                                    >  2:3387388/1‑146 (MQ=17)
                              cGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                                                                                                                 >  2:1219923/1‑148 (MQ=17)
                              cGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                                                                                                                 >  1:2053442/1‑148 (MQ=32)
                              cGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaaa                                                                                                                 >  2:1639384/1‑148 (MQ=32)
                               gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAAc                                                                                                                <  2:1469366/148‑1 (MQ=32)
                                aaaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAAGCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt                                                                                                               <  2:1708112/148‑1 (MQ=17)
                                  aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGc                                                                                                             <  2:2608440/148‑1 (MQ=32)
                                   ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGTGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGACGCCCGCCCGAAACAGcc                                                                                                            >  1:319180/1‑148 (MQ=11)
                                    gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCa                                                                                                           >  2:1050262/1‑148 (MQ=17)
                                    gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATACACTGCCa                                                                                                           >  2:3320599/1‑148 (MQ=17)
                                       cAGTGGGAAGACTGGGCCTTTGGTTTTATCTGTTGTTTGTCGGTGACCGCTCTCCTGAGTCGGACAAATCCGCGGGGAGCGTATTTTAACTTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGc                                                                                                        <  2:1008610/148‑1 (MQ=11)
                                       cAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGACGCCCGCCCGAAACAGCCAGGc                                                                                                        >  1:2736800/1‑148 (MQ=255)
                                         gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACAGCCAGGCAt                                                                                                      >  1:1960089/1‑148 (MQ=17)
                                         gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGACGCCCGCCCTAAACTGCCCGGcct                                                                                                      >  1:1752431/1‑146 (MQ=12)
                                             gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                                  >  2:1193895/1‑148 (MQ=11)
                                             gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa                                                                                                  >  2:1193819/1‑148 (MQ=11)
                                                                                   gAACGCTCTCCTGCGTAGGACAAACCCGCCGGGCGCGGATTTCAACTTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCa                                                            <  2:1534523/148‑1 (MQ=11)
                                                                                   gAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGACGCCCGCCCTAAACAGCCCGGCATCAAATTAAGCAGAACGCCATCCTGACCGATTGCCTTTTTGCa                                                            >  1:536189/1‑148 (MQ=11)
                                                                                     aCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGAAGGCCTTTTTGCAtt                                                          >  2:549866/1‑148 (MQ=17)
                                                                                      cGCTCGCCTGAGTAGGACAAATCCGCCGGGCGCGGATTTGAACGTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTg                                                         <  2:1075988/148‑1 (MQ=12)
                                                                                       gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTgg                                                        <  2:2167259/148‑1 (MQ=33)
                                                                                       gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTgg                                                        >  2:1273816/1‑148 (MQ=33)
                                                                                       gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCAACCTGACGGATGGCCTTTTTGCAGTgg                                                        >  2:1272983/1‑148 (MQ=14)
                                                                                       gCTCGCCTGCGTAGGACAAATCCGCCGGGAGCGGATTTGAACTTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTgg                                                        <  1:2307655/148‑1 (MQ=255)
                                                                                         tctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCGGACGGATGGCCTTTTTGCAGTGgcg                                                      >  1:2309664/1‑148 (MQ=12)
                                                                                         tctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgcg                                                      >  1:2356603/1‑148 (MQ=35)
                                                                                         tctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgcg                                                      >  1:2356837/1‑148 (MQ=35)
                                                                                          ctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCCGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAAGGGCATCCTGACGGGTTGCCTTTTTGCATTTgcgc                                                     >  1:785066/1‑148 (MQ=11)
                                                                                            ccTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAATCAGAAGGCCATCCAGACGGATGGGCTTTTTTCATTTGCGGAg                                                   >  1:3409546/1‑148 (MQ=11)
                                                                                                aGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaa                                               <  2:804894/148‑1 (MQ=40)
                                                                                                aGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaa                                               >  2:756216/1‑148 (MQ=40)
                                                                                                  tAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaaaa                                             <  1:2315592/148‑1 (MQ=255)
                                                                                                    ggACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAAt                                           >  1:2675357/1‑148 (MQ=255)
                                                                                                    ggACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAAt                                           >  2:3136509/1‑148 (MQ=255)
                                                                                                     gACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATg                                          <  2:736615/148‑1 (MQ=255)
                                                                                                       cAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGGGGGCAGGACGCCCGCCCAAAACTGGCAGGCATCAAATTAAGCAGAAGGCCCTCCAGTCGGCTGGGCGTTTTGCCTTGGCGCCGCACAAAATGcc                                        >  1:176978/1‑148 (MQ=11)
                                                                                                       cAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCACGCATCAAATATAGCAGAAGGCCAACCTGACGGATGGCCATTTTGCCTTGGGGCGGAAAAAAATGcc                                        >  1:3426682/1‑148 (MQ=14)
                                                                                                         aaTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCgg                                      <  2:3515586/148‑1 (MQ=255)
                                                                                                         aaTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTATTTGCATTGGCGCAGAAAAAAATGCCgg                                      <  2:3515544/148‑1 (MQ=255)
                                                                                                          aTCCGCCGGGCGCTGATTTGAACGTGGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGa                                     <  1:3107697/148‑1 (MQ=255)
                                                                                                           tCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGTAAAAAATGCCGGAt                                    >  2:2899436/1‑148 (MQ=255)
                                                                                                           tCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGAt                                    >  2:2899150/1‑148 (MQ=255)
                                                                                                             cgccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGc                                  >  2:569620/1‑148 (MQ=255)
                                                                                                              gccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCCTCCAGACGGATGGGCTTTTTTGATTGGGGCCGAAAAAAATGCCGGGTGCg                                 >  1:2165569/1‑148 (MQ=14)
                                                                                                               ccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCCGGCCGCCCGCCCTAAACGGGCCGGCCGCAAATTAAGCAGAAGGCCCGCCGGACGGCTGGGCTTTTTGGATTGGCGGCGGAAAAAATGACgggggcgc                                >  1:772528/1‑142 (MQ=2)
                                                                                                               ccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGa                                >  2:1294857/1‑148 (MQ=255)
                                                                                                               ccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGa                                >  2:513401/1‑148 (MQ=255)
                                                                                                                cgGGCGCGGATTTGACCGTTGCGCCCCCCCGGCCCGGAGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTCAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGAc                               <  1:353940/148‑1 (MQ=17)
                                                                                                                  ggAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAATATGCCGGATGCGACGc                             >  2:2663017/1‑148 (MQ=255)
                                                                                                                  ggAGCGGATTTGAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGc                             >  2:2663068/1‑148 (MQ=255)
                                                                                                                  ggAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGACGGCCATCCGGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGc                             <  2:196285/148‑1 (MQ=255)
                                                                                                                       ggATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGGGGGGAGGGCGCCCGCCCAAAACTGCCCGGCATCAAAATAAGCAGAAGGCCCACCGGAAGGCTGGGCTTTTTTGATTGGGGGCGACAAAAATTGCGGATGCGACGCggggg                        >  1:1636772/1‑146 (MQ=2)
                                                                                                                         aTTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGcgcg                      <  1:883147/148‑1 (MQ=255)
                                                                                                                            tGAAGGTTGCGAACCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCt                   <  1:976299/148‑1 (MQ=255)
                                                                                                                            tGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCt                   <  2:2828686/148‑1 (MQ=255)
                                                                                                                             gAACGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCAACCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGTTGCGACGCTGGCGCGTCtt                  >  2:1567856/1‑148 (MQ=255)
                                                                                                                              aaCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTa                 <  1:3195336/148‑1 (MQ=255)
                                                                                                                               aCGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTAt                <  2:1700060/148‑1 (MQ=255)
                                                                                                                               aCGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCCTCAAATTAAGCAGAAGGCCCTCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTAt                >  2:13940/1‑148 (MQ=255)
                                                                                                                                cGTTGCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCCGGCATCAAATTAAGCAGAAGGCCATACAGACGGATGGCCTTTTTGGATTGGGGCCGAAAAAAAAGCCGGCTGCGACGCGGGCGCGTCTTAta               >  1:3157304/1‑147 (MQ=17)
                                                                                                                                cGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATc               >  2:152118/1‑148 (MQ=255)
                                                                                                                                  ttGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCAACCTGACGGATGGGCTTTTTTCCTTGGGGCCGAAAAAAATTCCGGATGCGACGCTGGCGCGTCTTATCCa             >  1:2012449/1‑148 (MQ=18)
                                                                                                                                    gCGAAGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGc           <  1:1219923/148‑1 (MQ=255)
                                                                                                                                    gCGAACCAACGGCCCGGAGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGACGGCCATCCTCACGGATGGCCATTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGc           <  2:2579963/148‑1 (MQ=21)
                                                                                                                                       aaGCAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTa        <  2:827617/148‑1 (MQ=255)
                                                                                                                                          cAACTGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCAACCAGTCGGATGGCCTTTTTTGATTGGCGCAGAAAAAAATTCCGGATGCGACGCTGGCGCGGCTTATACACCCTAcac     >  1:1735240/1‑146 (MQ=18)
                                                                                                                                          cAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTActc     <  2:3288576/148‑1 (MQ=255)
                                                                                                                                          cAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTActc     <  2:898951/148‑1 (MQ=255)
                                                                                                                                           aaCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTActct    <  2:1752431/148‑1 (MQ=255)
                                                                                                                                           aaCGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTActct    <  1:2258203/148‑1 (MQ=255)
                                                                                                                                             cGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTg  <  2:2687070/148‑1 (MQ=255)
                                                                                                                                             cGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTg  >  2:850279/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTG  >  CP009273/4033976‑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: