Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,677,305 G→T 10.8% intergenic (+36/+68) PP_4137 → / ← chrR outer membrane siderophore receptor/chromate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,677,3050GT10.8% 298.6 / 2.6 127intergenic (+36/+68)PP_4137/chrRouter membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand):  ref base G (32/73);  new base T (4/13);  total (40/87)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

CTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGA  >  NC_002947/4677158‑4677452
                                                                                                                                                      |                                                                                                                                                    
cTATGCAGCGACCGTGACGCCGGCCTACACCGACGCCGGCCAGGACACTGCCCGTTCGACACCGGGCGAAGGATTGGGCGTGTAGGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGaa                                                                                                                                                    <  2:4521770/150‑1 (MQ=255)
cTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGaa                                                                                                                                                    <  1:1347690/150‑1 (MQ=255)
cTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGaa                                                                                                                                                    <  2:3980428/150‑1 (MQ=255)
cTATGCAGCGACCGTGACGCCGGCCTACAACCACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGaa                                                                                                                                                    <  2:180227/150‑1 (MQ=255)
      aGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGcc                                                                                                                                              <  2:1421632/150‑1 (MQ=255)
       gCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCg                                                                                                                                             <  1:3761973/150‑1 (MQ=255)
         gACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGAc                                                                                                                                           <  1:4125231/150‑1 (MQ=255)
         gACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGAc                                                                                                                                           <  2:4047943/150‑1 (MQ=255)
             gTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTAACCCCG‑AATGCGTCA‑ATGAAATTTACGACGC‑At                                                                                                                                      >  1:3725041/1‑150 (MQ=255)
              tGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑Att                                                                                                                                     <  1:3273687/150‑1 (MQ=255)
               gACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑CATGGGGGA‑TTGTAGTTGCGGACGC‑AGTc                                                                                                                                    >  2:4173435/1‑150 (MQ=255)
                aCGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTcg                                                                                                                                   <  1:4204250/150‑1 (MQ=255)
                aCGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTcg                                                                                                                                   <  2:3582694/150‑1 (MQ=255)
                  gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTTACCGCG‑AATTCGTGA‑TTGAATTTTGCGACGC‑AATagcg                                                                                                                                 >  2:239106/1‑150 (MQ=255)
                  gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTTAAAGCG‑AATTGGTGA‑TTTTATTTTCCGAAGC‑ATTagcg                                                                                                                                 >  2:239203/1‑150 (MQ=255)
                   ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCgg                                                                                                                                <  1:627479/150‑1 (MQ=255)
                   ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCgg                                                                                                                                <  1:4412326/150‑1 (MQ=255)
                   ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCgg                                                                                                                                <  1:4279505/150‑1 (MQ=255)
                   ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCgg                                                                                                                                <  1:2834889/150‑1 (MQ=255)
                   ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCgg                                                                                                                                <  2:934049/150‑1 (MQ=255)
                   ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCgg                                                                                                                                <  2:3307164/150‑1 (MQ=255)
                        cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTaaa                                                                                                                           <  1:1022220/150‑1 (MQ=255)
                        cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTaaa                                                                                                                           <  2:6886/150‑1 (MQ=255)
                        cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTaaa                                                                                                                           <  2:1494806/150‑1 (MQ=255)
                        cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTaaa                                                                                                                           <  2:1095536/150‑1 (MQ=255)
                         tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTTACCC‑CG‑AATGCGGGC‑GGGGATTTGGCGCCGC‑CATCGCGGGGAAAc                                                                                                                          >  2:3908177/1‑150 (MQ=255)
                         tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTAACCGCG‑AATGCGTCA‑TTGAATTTTCCGACGC‑AATAGCGGGTAAAc                                                                                                                          >  2:705753/1‑150 (MQ=255)
                         tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTAACCGCG‑AATGCGGCA‑ATTAATTTGGCGACGC‑AATAGCGGGTaaaa                                                                                                                          >  2:2079841/1‑149 (MQ=255)
                            aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTTTCCCGCG‑CATGCGTCA‑TTGAATTTTCCGACGC‑ATTCGCGGGTGAAcccc                                                                                                                       >  2:2192838/1‑149 (MQ=255)
                              cgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTAACCCCG‑CATGCGGCC‑TTGTATTTTGCGCCGC‑AATCGCGGGGGAACCCCCt                                                                                                                     >  2:2754257/1‑150 (MQ=255)
                              cgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTTTACCGCG‑GAAGGGTCC‑TTGGATTTGCCGGCGC‑ATTAGCGGGGAAACCCGcg                                                                                                                     >  2:3143586/1‑149 (MQ=255)
                               gacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTTACCCCCG‑GATGCGGCC‑TTGGAGTTGGCGACGC‑CGTCGCGGGGGAACCCGCTc                                                                                                                    >  2:2481412/1‑150 (MQ=255)
                                         aGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTTACCCCCC‑AAT‑CGGGA‑GTGTATTTTACGACGC‑AGTAGCGGGGGAAACCCCCCCAACACAcgg                                                                                                          >  1:3222670/1‑150 (MQ=255)
                                           gACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTTGCCGGG‑TATTGGTGA‑TTTTATTTTGGGACGC‑ATTAGCGGGTAAAACCGGTGCTACACAggggg                                                                                                        >  2:2475322/1‑149 (MQ=255)
                                              acggCCCGTCCGACCCCGGGCGACGGATTGGGCGTGTATGCCGGGGTTTCCTACAGCTGCCGCTGACGGGTGGGCGC‑GGGGTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAACCCCGCTCCTACACAGGGGCGtt                                                                                                     <  2:2098233/147‑1 (MQ=38)
                                                   ccGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTTCCCGCGCAATGGGGC‑‑GTGTATTTTGCGACGC‑AATAGCCGGTAAACCCCCACCTACACAGGGGCGGTggcca                                                                                                >  2:1998070/1‑146 (MQ=255)
                                                    cGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTTACCC‑CG‑AAAGCGGGAGTTGAATTTG‑CGACCCAAAACCCGGGTAAACCC‑CACCGGCACAAGGGCGTTGGAcag                                                                                               >  1:1064466/1‑150 (MQ=255)
                                                     gTTCTACACCGGGCGACGGATTGGGCGTGTAGGCCGGGGTTTCCTCCAGCTGCCGCTGAGGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACCC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCAcagc                                                                                              <  2:3902974/150‑1 (MQ=38)
                                                     gTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTTAATTTGCCGACGC‑ATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCAcagc                                                                                              <  2:3902850/150‑1 (MQ=255)
                                                      ttCGACACCGGGCGAAGTATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTAA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCAcagca                                                                                             <  1:2975448/150‑1 (MQ=255)
                                                      ttCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCAcagca                                                                                             <  2:1091314/150‑1 (MQ=255)
                                                              cGGGCTAAGGTTTGGGCTTTTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGG‑GGTTTAACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAgg                                                                                     <  2:9573/150‑1 (MQ=38)
                                                              cGGGCGCTGGATTGGGCGTGTATGCCGGGGTGTCCTACAGCTTCCGCTGCCGGGGGGGCGG‑GGGTTTCCCCGCG‑AATGCGTCA‑TTGATTTTGCCGACGC‑TTTCGGGTGTAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAgg                                                                                     <  2:3413942/150‑1 (MQ=255)
                                                               gggCGATGGTTTGGGGGTGTTGGCGGGGGTTTCCTACCGCTCCCGCTGCGGGGTGGGCGG‑GGGTTTACCCGCT‑ATTGCGTCT‑TTTATTTTGCCGACGC‑ATTCGGGGTTAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAggg                                                                                    <  1:3095187/150‑1 (MQ=255)
                                                                   gAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTTAACCGCG‑GATGCCGAA‑TTGTATTTTCCGGCCC‑ATTAGCGGGTAAACCCGCTCCTAAACCGGGGCGTTGCCCACCAGGATTAGGGGttg                                                                                >  2:3217131/1‑149 (MQ=255)
                                                                    aaGTATTTGGCGTGTATGCCGGGGTTTCCTCCTGCTGCGGCTGAGGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑TTTCGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCa                                                                               <  1:1012971/150‑1 (MQ=38)
                                                                    aaGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGG‑GGGTTTACCCGCG‑ATTGCGTCT‑TTTAATTTGCCTACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCa                                                                               <  1:1433343/150‑1 (MQ=39)
                                                                    aaGGATTGGGCGTGTATGCCGGGGTTTCCTACCGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGACTTGGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCa                                                                               <  1:1433108/150‑1 (MQ=255)
                                                                    aaGGATGTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGCCGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCa                                                                               <  2:1941471/150‑1 (MQ=255)
                                                                          ttGGGGTTTTTGCGGGGGTTACCTACCGCTTCGGCTGTGGGGGGGGGGG‑GGTTTCACCCGCT‑ATTGCGTAT‑TTGAATTTGCCGCCGA‑TTGCGGGGTTAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGc                                                                         <  1:1189898/150‑1 (MQ=255)
                                                                           tGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTTAATTTGCCGACGC‑ATTCGCGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                                        <  2:106568/150‑1 (MQ=255)
                                                                             gCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCt                                                                      <  2:4357378/150‑1 (MQ=255)
                                                                                 gttggcgggggTTTCCTCCCGCTTCCGCTGCGGGGGGGGCGC‑GGGTTTCCCCGCG‑AATGCGTCT‑TTGATTTTGCCGACGC‑ATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                  <  1:2507398/143‑1 (MQ=255)
                                                                                  tATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCa                                                                 <  2:1765505/150‑1 (MQ=255)
                                                                                   aTGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                                <  2:3485164/150‑1 (MQ=255)
                                                                                    tGCCGGGTTTTCCTACCGCTTCCGCTGAGGGGTGGGCGC‑GGGTTTACCCGCT‑AATGCGTAA‑TTTAATTTGCCGCCGC‑TTTCGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAAc                                                               <  1:3279557/150‑1 (MQ=255)
                                                                                    tGCCGGGTTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTAA‑TTGAATTTGCCGACGC‑ATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAAc                                                               <  2:1586116/150‑1 (MQ=255)
                                                                                      ccGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCC‑CG‑AATGCGTCA‑TTGAAATTGCCGACGC‑ATTCGCGGGGAAACCCCCTCCTACACAGGGGCGTTGCACAGCAGGAGTAGGGGTTCAGACCGCCCGGTTCAACtta                                                            >  2:747951/1‑149 (MQ=39)
                                                                                               cctcccgCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGa                                                    <  2:3535020/144‑1 (MQ=255)
                                                                                               ccTACTGCTTCCGGTGCGGGGTGGGGGG‑GGGTTTACCGGGG‑AATGCTAAA‑TTTATTTTGCCGACGC‑ATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGa                                                    <  1:1958127/150‑1 (MQ=37)
                                                                                                  aCAGCTTCCGCTGACGGGTGGGCGC‑GGGGTTACCCCCG‑AATGCGTCA‑GTGAATTTGCCGACGC‑CATCGCGGGTGAACCCCCGCCGACAAAGGGGCGGTGCACAGCCGGATTAGGGGTTCAGACCGCCCGGGTCAACATCACCCAAGACGc                                                 >  2:1416263/1‑150 (MQ=255)
                                                                                                     gCTTCGGCTGTGGGGTGGGGGC‑GGGTTTACCCGCT‑ATTGCTTAT‑TTTAATTTGCCGACGA‑TTTCGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCaaa                                              <  1:3475419/150‑1 (MQ=255)
                                                                                                     gCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCTTCA‑TTGAATTTGCCGACGA‑ATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCaaa                                              <  1:2192838/150‑1 (MQ=255)
                                                                                                      cTTCCTCTGTGGGGTGGGGGG‑GGTTTTACCCGGG‑ATTGCTTAT‑TTTATTTTGCCGACGA‑TTTGGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                             <  1:14670/150‑1 (MQ=255)
                                                                                                        tCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCT‑AATGCGTCA‑TTGATTTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCt                                           <  1:1093612/150‑1 (MQ=255)
                                                                                                           ggggACGGGGGGGCGC‑GGGTTTCCCGGCG‑AATGGGTCA‑TTGATTTGGCCGACGC‑TTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTc                                        <  1:2826448/147‑1 (MQ=255)
                                                                                                           gCTGAGGGGTGGGCGC‑GGGTTTACCCGCT‑AATGCGTCA‑TTTATTTTGCCGACGA‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTc                                        <  1:3129669/150‑1 (MQ=255)
                                                                                                           gCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTAA‑TTGATTTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTc                                        <  1:2563710/150‑1 (MQ=255)
                                                                                                             tgccGGGTGGGCGC‑GGGTTTACCCGGG‑AATGCGTAT‑TTGAATTTGCCGACGC‑ATTCGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGa                                      <  1:2171628/147‑1 (MQ=255)
                                                                                                             tGAGGGGTGGGCGC‑GGTTTTACCCGCT‑ATTGCGTAT‑TTTAATTTGCCGACGA‑TTTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGa                                      <  1:2129710/150‑1 (MQ=255)
                                                                                                             tGAGGGGGGGGCGC‑GGTTTTACCCGCG‑AATGCGTCA‑TTTAATTTGCCGACGC‑ATCCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGa                                      <  1:1105274/150‑1 (MQ=255)
                                                                                                                     aggagc‑gGGTTTACCCGCG‑AATGCGTCG‑GCAAATTCAATGACGC‑ATTCGCGGGTAAACCCGCGCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGcc                              >  1:666969/5‑150 (MQ=21)
                                                                                                                            gggTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAg                        >  1:174867/1‑150 (MQ=255)
                                                                                                                            gggTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAg                        >  1:174172/1‑150 (MQ=255)
                                                                                                                              gTTTACCCGCT‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                      <  2:1149939/150‑1 (MQ=255)
                                                                                                                              gTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                      <  2:4447065/150‑1 (MQ=255)
                                                                                                                              gTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGAGTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCACACCTGTCGACAAACGCCTGCAAGaa                      >  1:154182/1‑150 (MQ=255)
                                                                                                                               tttACCCGCG‑AATGCGTCG‑GCAAATTCAATGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc                     >  2:48964/1‑150 (MQ=37)
                                                                                                                                ttACCCGCG‑AATGCGTCA‑TTGAAGTTGCCGACGC‑AGTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACg                    <  2:4402767/150‑1 (MQ=255)
                                                                                                                                   cccGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGTGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc                 <  2:1190804/150‑1 (MQ=255)
                                                                                                                                    ccGCG‑AATGCGTCG‑GCAAATTCAATGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                >  2:3377455/1‑150 (MQ=37)
                                                                                                                                    ccGCG‑AATGCGTCG‑GCAAATTCAATGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                >  1:828821/1‑150 (MQ=37)
                                                                                                                                    ccGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                >  1:2463703/1‑150 (MQ=255)
                                                                                                                                     cgcg‑AATGCGTCA‑TTGATTTTGCCGCCGA‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc               <  1:3110287/150‑1 (MQ=255)
                                                                                                                                     cgcg‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc               <  2:3848021/150‑1 (MQ=255)
                                                                                                                                     cgcg‑AATGCGTAA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc               <  2:1883756/150‑1 (MQ=255)
                                                                                                                                      gcg‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgcg              >  1:321716/1‑150 (MQ=255)
                                                                                                                                      gcg‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgcg              >  2:1033588/1‑150 (MQ=255)
                                                                                                                                           aTGCGTCT‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCtt          <  1:289659/150‑1 (MQ=255)
                                                                                                                                           aTGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCtt          <  2:492830/150‑1 (MQ=255)
                                                                                                                                               gTCA‑TTGAGTTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa      >  1:862964/1‑150 (MQ=255)
                                                                                                                                               gTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACACATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa      >  1:862947/1‑150 (MQ=255)
                                                                                                                                                tCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAt     >  1:1197738/1‑150 (MQ=255)
                                                                                                                                                 cA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt    >  2:3513131/1‑150 (MQ=255)
                                                                                                                                                 cA‑TTGAATTTGCCGACGC‑ATTCGCGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt    <  1:1418253/150‑1 (MQ=255)
                                                                                                                                                  a‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                                <  1:485177/89‑1 (MQ=255)
                                                                                                                                                    ttGAATTTGCCGCCGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGa  <  2:3222670/150‑1 (MQ=255)
                                                                                                                                                    ttGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGa  <  1:3752746/150‑1 (MQ=255)
                                                                                                                                                    ttGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGa  <  1:2481412/150‑1 (MQ=255)
                                                                                                                                                    ttGAATTTGCCGACGC‑ATTCGCGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGa  <  2:1225182/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTACCCGCG‑AATGCGTCA‑TTGAATTTGCCGACGC‑ATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGA  >  NC_002947/4677158‑4677452

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