Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,909,992 G→C 5.3% intergenic (‑1/‑281) PP_4317 ← / → PP_4318 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,909,9920GC 5.3% 601.6 / 21.0 209intergenic (‑1/‑281)PP_4317/PP_4318hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (98/100);  new base C (5/6);  total (103/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 212 total aligned reads displayed.

CTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACA  >  NC_002947/4909843‑4910141
                                                                                                                                                     |                                                                                                                                                      
cTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATg                                                                                                                                                        <  2:4009690/150‑1 (MQ=255)
cTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATg                                                                                                                                                        <  2:2565535/150‑1 (MQ=255)
  gttgttGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGa                                                                                                                                                      >  2:7350790/1‑150 (MQ=255)
   ttgttgCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAt                                                                                                                                                     <  1:6744307/150‑1 (MQ=255)
   ttgttgCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAt                                                                                                                                                     <  2:6767616/150‑1 (MQ=255)
       tgctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCGTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGccc                                                                                                                                                 <  2:8029266/150‑1 (MQ=255)
          tgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTc                                                                                                                                              >  1:7588190/1‑150 (MQ=40)
          tgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTc                                                                                                                                              >  2:3307858/1‑150 (MQ=40)
           gttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                             <  2:1127966/150‑1 (MQ=255)
           gttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                             <  1:8987788/150‑1 (MQ=255)
           attgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                             <  2:3157806/149‑1 (MQ=255)
              gctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                                          <  2:3850512/150‑1 (MQ=39)
                  ttgccgttgccgttgccgttgccgttgccgTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                                                                                                                      >  1:5179065/1‑150 (MQ=37)
                  ttgccgttgccgttgccgttgccgttgccgTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                                                                                                                      >  2:1993931/1‑150 (MQ=37)
                   tgccgttgccgttgccgttgccgttgccgtTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                    <  1:5137194/150‑1 (MQ=37)
                    gccgttgccgttgccgttgccgttgccgttGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAg                                                                                                                                   <  1:1849009/150‑1 (MQ=37)
                      cgttgccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGc                                                                                                                                 >  1:4153678/1‑150 (MQ=37)
                        ttgccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGccc                                                                                                                               >  1:5838597/1‑150 (MQ=37)
                         tgccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGcccc                                                                                                                              <  1:8321407/150‑1 (MQ=37)
                         tgccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGcccc                                                                                                                              <  1:1793998/150‑1 (MQ=37)
                          gccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCa                                                                                                                             <  1:8301372/150‑1 (MQ=37)
                          gccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCa                                                                                                                             <  2:3674805/150‑1 (MQ=37)
                            cgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGCGCATGGATGCCCGTCAAGCGCTGAGGCCCCAgg                                                                                                                           >  2:8617991/1‑150 (MQ=37)
                             gttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGa                                                                                                                          >  1:2405677/1‑150 (MQ=37)
                             gttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGa                                                                                                                          <  1:6020655/150‑1 (MQ=37)
                             cttgccgttgcccttgccgatgATCTTGATCTAGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGCCCCCAGGa                                                                                                                          <  2:5203000/149‑1 (MQ=255)
                              ttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCACGCGCTGAGGCCCCAGGAt                                                                                                                         <  2:7611329/150‑1 (MQ=37)
                               tgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATg                                                                                                                        <  1:7478307/150‑1 (MQ=37)
                                 ccgttgccgttgccgtttATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATggg                                                                                                                      <  1:6989639/150‑1 (MQ=37)
                                   gttgccgttgccgttgATCTTGATCTTGATCTTGATCTTTATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                                    <  2:285721/150‑1 (MQ=37)
                                     tgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCATCAAGCGCTGAGGCCCCAGGATGGGGCCt                                                                                                                  <  1:5178666/150‑1 (MQ=37)
                                      gccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAGGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCtt                                                                                                                 <  2:4892654/150‑1 (MQ=37)
                                       ccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCTGGGGCCtca                                                                                                                <  1:4017093/150‑3 (MQ=255)
                                       ccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCTGGGGCCtca                                                                                                                <  1:3022996/150‑3 (MQ=255)
                                       ccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTc                                                                                                                >  1:7313668/1‑150 (MQ=37)
                                       ccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTc                                                                                                                <  1:8424017/150‑1 (MQ=37)
                                       ccgttgccgttgATCTTGATCTAGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCTGGGGCCtca                                                                                                                <  1:4122441/150‑3 (MQ=255)
                                        cgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                               >  1:5356358/1‑150 (MQ=35)
                                          ttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGCGGCAGGTGCCGGGCATGGATGCCCGTCAAGCGCCGAAGGCCCCATCCTGGGGCCtcagc                                                                                                              <  1:2552360/150‑5 (MQ=11)
                                          ttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgc                                                                                                             <  1:8828520/150‑1 (MQ=34)
                                          ttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgc                                                                                                             >  1:4487913/1‑150 (MQ=34)
                                            gccgtttATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgcgc                                                                                                           <  1:4594071/150‑1 (MQ=17)
                                            gccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgcgc                                                                                                           <  2:8076911/150‑1 (MQ=32)
                                             ccgttgATCTTGGTCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGTGGCAGGTGACGGGCATGGAGGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                          <  1:393914/150‑1 (MQ=12)
                                              cgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCgg                                                                                                         >  1:4839053/1‑150 (MQ=31)
                                               gttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGt                                                                                                        >  1:3157517/1‑150 (MQ=11)
                                               gttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGt                                                                                                        >  1:2491597/1‑150 (MQ=11)
                                               gttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGt                                                                                                        >  1:1979240/1‑150 (MQ=11)
                                                             tgatcttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGa                                                                                          <  2:5775418/150‑1 (MQ=31)
                                                              gatcttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAg                                                                                         >  1:116619/1‑150 (MQ=32)
                                                                 cttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGTGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGcc                                                                                      >  1:6682159/1‑150 (MQ=17)
                                                                 cttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGTGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGcc                                                                                      >  1:6682283/1‑150 (MQ=17)
                                                                 cttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGcc                                                                                      >  1:2048131/1‑150 (MQ=32)
                                                                  ttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCg                                                                                     >  1:1679991/1‑150 (MQ=32)
                                                                    gatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGa                                                                                   >  2:6963039/1‑150 (MQ=32)
                                                                     atcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAg                                                                                  <  2:704602/150‑1 (MQ=32)
                                                                      tcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGt                                                                                 >  1:5386986/1‑150 (MQ=32)
                                                                      tcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGt                                                                                 <  1:1092110/150‑1 (MQ=32)
                                                                         tgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAg                                                                              >  1:3324275/1‑150 (MQ=32)
                                                                              ttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGAc                                                                         >  2:3372008/1‑150 (MQ=32)
                                                                                  tcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCg                                                                     >  2:6463063/1‑150 (MQ=32)
                                                                                    ttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGa                                                                   >  2:1835769/1‑150 (MQ=32)
                                                                                     tgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGag                                                                  >  2:3989755/1‑150 (MQ=32)
                                                                                     tgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGag                                                                  >  1:6441974/1‑150 (MQ=32)
                                                                                       atcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGagcg                                                                >  1:3557183/1‑150 (MQ=32)
                                                                                        tcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGagcgc                                                               >  2:1575723/1‑150 (MQ=32)
                                                                                          ttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGagcgcag                                                             >  1:3715073/1‑150 (MQ=32)
                                                                                            gatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGagcgcagcg                                                           >  2:2284051/1‑150 (MQ=32)
                                                                                                ttggcgACTTCTGGAGGCCTAACGCAGGCCTTGCGTCGGCTGGTGAGGGGCATGTATGCCCGAAAAGGGCTGAGGCCCCTGTATGGGGCTTTCGGCGGGGTCCTGCCGGGAGCGAGGCCGGAGTGTGGGGACCCCGGAGCGCAGCGCAgg                                                       <  2:4140864/147‑1 (MQ=11)
                                                                                                ttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAgg                                                       >  1:7274372/1‑150 (MQ=32)
                                                                                                    cgACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGcc                                                   >  2:4212700/1‑150 (MQ=32)
                                                                                                     gACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCg                                                  >  2:3783046/1‑150 (MQ=32)
                                                                                                        ttCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGgat                                               >  2:6887558/1‑150 (MQ=32)
                                                                                                        ttCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCGGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGgat                                               <  2:6528523/150‑1 (MQ=17)
                                                                                                             gaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATgg                                          >  1:8543205/1‑150 (MQ=32)
                                                                                                               ggcggcccGCCGGCGTGGCGGCGGCCGGGGCCGGGCAGGGCGGCCGGTCCCGCGCGGCGGCCCCGGGCGGGGGCCGTCGGCGCGGCCCTGCCGGGAGCGCGGCGGGAGGGCGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                        <  2:1189334/143‑1 (MQ=2)
                                                                                                               ggCCGAACGCCGGCCTTGCGGAGGCTGGTGCGGGGCATGGTGGCCCGTCAAGCGCTGAGGCCCCCGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                        <  2:8543205/150‑1 (MQ=11)
                                                                                                               ggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                        <  2:4487913/150‑1 (MQ=32)
                                                                                                               ggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                        <  2:8260771/150‑1 (MQ=32)
                                                                                                                gCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAg                                       >  2:5818538/1‑150 (MQ=32)
                                                                                                                  cGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcg                                     >  1:3066504/1‑150 (MQ=32)
                                                                                                                  cGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcg                                     >  1:4316678/1‑150 (MQ=32)
                                                                                                                      cGCTGGCCTTGCGGTGGCTGGTGACGGGCATGTTTGCCCGTCACGCGCTGTGGCCCCAGTTTGGGGCCTTCGGCGCGGTCCTGCGGGGAGCGTGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGc                                 <  1:8171642/150‑1 (MQ=11)
                                                                                                                         aGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtt                              >  1:613365/1‑150 (MQ=32)
                                                                                                                             cTTGCGGTGGCCGGGGACGGGCATGGATGCCCGTCAAGCGCTGTGGCCCCTGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttttt                          <  2:7852797/150‑1 (MQ=11)
                                                                                                                             cTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttttt                          >  2:5350318/1‑150 (MQ=32)
                                                                                                                              ttGCGGTGGCTGGTGACGGGCATGGATGCCCGTCACGCGCTGAGGCCCCTGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                         <  2:4649804/150‑1 (MQ=255)
                                                                                                                              ttGCGGAGGCTGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                         <  2:4448367/150‑1 (MQ=17)
                                                                                                                              ttGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCCTCGGCGCGGTCCTGCCGGGAGCGAGGCCg                                                                                     <  2:8865421/90‑1 (MQ=18)
                                                                                                                                 cGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGtt                      >  1:2124152/1‑150 (MQ=32)
                                                                                                                                   gtggCTGGTGCCGGGCTTGGATGCCGGTAACGCGGGGCGGCCCCCGTTGGGGGCTTTCGGCGGGGTCGGGCGGGGTGCGCGGCGGGTGTGCGGGGACCCCGGAGCGCAGCGCAGGGGCCGGTTGATGGGAGCGCAGCGTTTTTTGGTTAc                    <  1:8059195/148‑1 (MQ=2)
                                                                                                                                   gAGGCAGGTGACGGGCCTGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTGGTTGGTTAc                    >  2:468036/1‑150 (MQ=12)
                                                                                                                                     ggCTGGTGACGGGCTTGGTGGCCCGTCACGCGCTGTGCCCCCCGTCGGGGGCCTTCGGGGCGGTCCTGCGGGTTGCGCGGCGGGTGTGCGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACtt                  <  1:6531684/150‑1 (MQ=11)
                                                                                                                                        gcttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt               >  2:7579356/4‑150 (MQ=12)
                                                                                                                                        aGGTGACGGGCATGGTTGCCCGTCATGCGCTGTGGCCCCTGGTTGGGGCCTTGGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt               <  2:4153678/150‑1 (MQ=11)
                                                                                                                                            gACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcg           >  2:7620502/1‑150 (MQ=32)
                                                                                                                                             aCGGGCTTGGAGGCCCGTCACGCGCTGAGGCCCCAGGTTGGGGCCTTCGGCGGGGTCCTGCCGGGAGCGAGGCCGGTGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcgc          <  2:2405677/150‑1 (MQ=11)
                                                                                                                                               gggCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGACGCGt        >  2:7960498/1‑150 (MQ=17)
                                                                                                                                                ggCTTGGTTGCCCGTAATGCGCTGAGGCCCCTGGATGGGGCCTTCGGCGCGGTCCTGCCGGTTGCGAGGCCGGTGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGtt       <  2:4092940/150‑1 (MQ=11)
                                                                                                                                                     ggATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACa  >  1:1446768/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                      
CTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACA  >  NC_002947/4909843‑4910141

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: