Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,909,993 G→C 5.7% intergenic (‑2/‑280) 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,9930GC 5.7% 592.2 / 21.7 209intergenic (‑2/‑280)PP_4317/PP_4318hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (98/96);  new base C (5/7);  total (103/106)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

GTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAA  >  NC_002947/4909845‑4910142
                                                                                                                                                    |                                                                                                                                                      
gttgttGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGa                                                                                                                                                       >  2:7350790/1‑150 (MQ=255)
 ttgttgCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAt                                                                                                                                                      <  2:6767616/150‑1 (MQ=255)
   gttgctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGc                                                                                                                                                    >  1:6196299/1‑150 (MQ=255)
     tgctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGccc                                                                                                                                                  <  1:8314633/150‑1 (MQ=255)
     tgctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCGTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGccc                                                                                                                                                  <  2:8029266/150‑1 (MQ=255)
        tgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTc                                                                                                                                               >  1:7588190/1‑150 (MQ=40)
         gttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                              <  2:1127966/150‑1 (MQ=255)
         gttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                              <  1:8987788/150‑1 (MQ=255)
            gctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                                           <  2:3850512/150‑1 (MQ=39)
            gctgttgcCGTTGCCGTTGCCGTTGCCGGTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCTAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCacct                                                                                                                                           <  1:7292000/150‑4 (MQ=14)
              tgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcgc                                                                                                                                         >  2:4544000/1‑150 (MQ=38)
               gttgccgttgccgttgccgttgccgttgccGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                                                                                                                        <  2:6512578/150‑1 (MQ=37)
               gttgccgttgccgttgccgttgccgttgccGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                                                                                                                        >  2:4464981/1‑150 (MQ=37)
               gttgccgttgccgttgccgttgccgttgccGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                                                                                                                        <  2:3310896/150‑1 (MQ=37)
                ttgccgttgccgttgccgttgccgttgccgTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                                                                                                                       >  2:1993725/1‑150 (MQ=37)
                 tgccgttgccgttgccgttgccgttgccgtTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                     <  1:5137194/150‑1 (MQ=37)
                 tgccgttgccgttgccgttgccgttgccgtTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGACGCCCGTCAAGCGCTGa                                                                                                                                     <  1:4192877/150‑1 (MQ=37)
                  gccgttgccgttgccgttgccgttgccgttGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAg                                                                                                                                    <  1:1849009/150‑1 (MQ=37)
                       tgccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGcccc                                                                                                                               <  1:1793998/150‑1 (MQ=37)
                        gccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCa                                                                                                                              <  2:3674805/150‑1 (MQ=37)
                          cgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGCGCATGGATGCCCGTCAAGCGCTGAGGCCCCAgg                                                                                                                            >  2:8617991/1‑150 (MQ=37)
                           cttgccgttgcccttgccgatgATCTTGATCTAGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGCCCCCAGGa                                                                                                                           <  2:5203000/149‑1 (MQ=255)
                            ttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCACGCGCTGAGGCCCCAGGAt                                                                                                                          <  2:7611329/150‑1 (MQ=37)
                             tgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATg                                                                                                                         <  1:7478307/150‑1 (MQ=37)
                               ccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATggg                                                                                                                       <  1:6820083/150‑1 (MQ=37)
                                   tgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCATCAAGCGCTGAGGCCCCAGGATGGGGCCt                                                                                                                   <  1:5178671/150‑1 (MQ=37)
                                    gccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAGGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCtt                                                                                                                  <  2:4892654/150‑1 (MQ=37)
                                     ccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCTGGGGCCtca                                                                                                                 <  1:5274756/150‑3 (MQ=255)
                                     ccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCTGGGGCCtca                                                                                                                 <  1:3022996/150‑3 (MQ=255)
                                      cgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                                <  2:2043462/150‑1 (MQ=35)
                                       gttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCgg                                                                                                               >  1:4092940/1‑150 (MQ=35)
                                        ttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGCGGCAGGTGCCGGGCATGGATGCCCGTCAAGCGCCGAAGGCCCCATCCTGGGGCCtcagc                                                                                                               <  1:2552360/150‑5 (MQ=11)
                                        ttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgc                                                                                                              <  1:8828520/150‑1 (MQ=34)
                                        ttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCATGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGTTGACGGGCATGGATGCCCGTCAAGCGCCGAAGGCCCCATCCTGGGGCCtcagc                                                                                                               <  1:8378068/150‑5 (MQ=11)
                                          gccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgcgc                                                                                                            <  2:8076911/150‑1 (MQ=32)
                                           ccgttgATCTTGGTCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCTGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                           <  1:3826497/150‑1 (MQ=17)
                                           ccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                           <  1:4367416/150‑1 (MQ=32)
                                           ccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                           >  2:3230997/1‑150 (MQ=32)
                                            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)
                                             gttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGt                                                                                                         >  1:5371620/1‑150 (MQ=11)
                                                          ttgatcttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGCGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCg                                                                                            <  1:231565/150‑1 (MQ=11)
                                                           tgatcttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGa                                                                                           <  2:5775418/150‑1 (MQ=31)
                                                            gctcttgatgttgttctggctcttgatcttgttcttCGCGCCTTGTGGAGGCCGAACGCTGGCCTTGCGGTGGCAGGGGACGGGCTTGGATGCCCGTCAAGCGCGGAGGCCCCTGGTTGGGGGCGTCGGCGCGGTCCTGCCGGGAGCGAg                                                                                          <  2:3149185/148‑1 (MQ=11)
                                                            gatcttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAg                                                                                          >  1:116619/1‑150 (MQ=32)
                                                               cttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGTGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGcc                                                                                       >  1:6682160/1‑150 (MQ=17)
                                                               cttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGcc                                                                                       >  1:2048131/1‑150 (MQ=32)
                                                               cttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGcc                                                                                       >  1:414786/1‑150 (MQ=32)
                                                                ttgatcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCg                                                                                      >  1:6735697/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:1092110/150‑1 (MQ=32)
                                                                       tgatcttgatctttatcttgatcttCGCGACTTCTGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCATGCGCTTAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAg                                                                               <  1:8617991/150‑1 (MQ=255)
                                                                       tgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAg                                                                               >  1:3324275/1‑150 (MQ=32)
                                                                        gatcttgatcttgatcttgatcttCGCGACTTCCGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAgg                                                                              <  1:8457842/150‑1 (MQ=17)
                                                                          tcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAgggg                                                                            >  1:5994609/1‑150 (MQ=32)
                                                                            ttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGAc                                                                          >  1:738913/1‑150 (MQ=32)
                                                                                  ttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGa                                                                    >  1:452751/1‑150 (MQ=32)
                                                                                   tgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGag                                                                   >  1:6441974/1‑150 (MQ=32)
                                                                                   tgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGag                                                                   >  2:3989909/1‑150 (MQ=32)
                                                                                     atcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGagcg                                                                 >  1:3557183/1‑150 (MQ=32)
                                                                                              ttggcgACTTCTGGAGGCCTAACGCAGGCCTTGCGTCGGCTGGTGAGGGGCATGTATGCCCGAAAAGGGCTGAGGCCCCTGTATGGGGCTTTCGGCGGGGTCCTGCCGGGAGCGAGGCCGGAGTGTGGGGACCCCGGAGCGCAGCGCAgg                                                        <  2:4140864/147‑1 (MQ=11)
                                                                                              ttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAgg                                                        >  1:4801220/1‑150 (MQ=32)
                                                                                                  ggACTTCCGGCGGCCGACCGCCGGCCTTGCGGCGGCTGGTGACGGGCTTGGCTGCCCGTCAAGGGCTGTGGCCCCCGGCTGGGGCCTTCGGCGGGGTCCTGCGGGGAGCGCGGCCGGCGTGCGGGGACCCCGGAGCGCAGCGCAGGGGcc                                                    <  1:3701690/149‑1 (MQ=11)
                                                                                                  cgACTTCAGGTGGCCGACCGCAGGCCTTGCGGTGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCTGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGTGGCCGGCGTGAGGGGACCCCGGAGCGCAGCGCAGGGGcc                                                    <  1:8296802/150‑1 (MQ=11)
                                                                                                  cgACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGcc                                                    >  2:4212700/1‑150 (MQ=32)
                                                                                                   gACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCg                                                   >  2:3783046/1‑150 (MQ=32)
                                                                                                      ttCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGgat                                                >  2:6887558/1‑150 (MQ=32)
                                                                                                      ttCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGgat                                                >  2:2688580/1‑150 (MQ=32)
                                                                                                          ggaggCCGACCGCAGGCCTTGCGGAGGCAGGTGCGGGGCTTGGAGGCCCGACAAGGGCTGCGGCCCCAGGATGGGGCCTTGGGCGCGGTCCTGCGGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGatgatg                                            <  2:3952489/150‑1 (MQ=11)
                                                                                                           gaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATgg                                           >  1:5666830/1‑150 (MQ=32)
                                                                                                             ggcggcccGCCGGCGTGGCGGCGGCCGGGGCCGGGCAGGGCGGCCGGTCCCGCGCGGCGGCCCCGGGCGGGGGCCGTCGGCGCGGCCCTGCCGGGAGCGCGGCGGGAGGGCGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                         <  2:1189334/143‑1 (MQ=2)
                                                                                                             ggCCGAACGCCGGCCTTGCGGAGGCTGGTGCGGGGCATGGTGGCCCGTCAAGCGCTGAGGCCCCCGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                         <  2:8543205/150‑1 (MQ=11)
                                                                                                             ggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                         <  2:8260771/150‑1 (MQ=32)
                                                                                                              gCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAg                                        >  2:5818538/1‑150 (MQ=32)
                                                                                                                cGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcg                                      >  1:3066504/1‑150 (MQ=32)
                                                                                                                 gAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcgc                                     >  2:2441044/1‑150 (MQ=32)
                                                                                                                  aaCGCAGGCCTTGCGGCGGCCGGGGCCGGGCCGGGATCCCCGCCACGCGCTGCGGCCCCCGGCTGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGCGGCCGGCGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCa                                    <  2:2363941/150‑1 (MQ=11)
                                                                                                                         gCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttt                             >  1:6129618/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)
                                                                                                                            ttGCGGTGGCAGGTGCCGGGCATGCCTGCCCGTCACGCGCTGCGGCCCCTGGTGGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGTGGCCGGCGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                          <  1:708762/150‑1 (MQ=11)
                                                                                                                            ttGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCCTCGGCGCGGTCCTGCCGGGAGCGAGGCCg                                                                                      <  2:8865421/90‑1 (MQ=18)
                                                                                                                            ttGCGGAGGCAGGTGACGGGCAGGGATGCCCGTCACGCGCTGAGGCCCCCGGAGGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGCGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                          <  2:1679991/150‑1 (MQ=11)
                                                                                                                              gCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGt                        >  2:7660500/1‑150 (MQ=32)
                                                                                                                               cGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGtt                       >  1:2124152/1‑150 (MQ=32)
                                                                                                                                 gtggCTGGTGCCGGGCTTGGATGCCGGTAACGCGGGGCGGCCCCCGTTGGGGGCTTTCGGCGGGGTCGGGCGGGGTGCGCGGCGGGTGTGCGGGGACCCCGGAGCGCAGCGCAGGGGCCGGTTGATGGGAGCGCAGCGTTTTTTGGTTAc                     <  1:8059195/148‑1 (MQ=2)
                                                                                                                                 gAGGCAGGTGACGGGCCTGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTGGTTGGTTAc                     >  2:468036/1‑150 (MQ=12)
                                                                                                                                 gAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTAc                     <  2:6579022/150‑1 (MQ=32)
                                                                                                                                   ggCTGGTGACGGGCTTGGTGGCCCGTCACGCGCTGTGCCCCCCGTCGGGGGCCTTCGGGGCGGTCCTGCGGGTTGCGCGGCGGGTGTGCGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACtt                   <  1:6531684/150‑1 (MQ=11)
                                                                                                                                      gcttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt                >  2:7579356/4‑150 (MQ=12)
                                                                                                                                      aGGTGACGGGCATGGTTGCCCGTCATGCGCTGTGGCCCCTGGTTGGGGCCTTGGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt                <  2:4153678/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)
                                                                                                                                                    gATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACaa  >  1:398953/1‑150 (MQ=32)
                                                                                                                                                    |                                                                                                                                                      
GTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAA  >  NC_002947/4909845‑4910142

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: