Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 541,418 T→C 10.3% S1265P (TCC→CCC) ‡ rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947541,4180TC10.3% 523.4 / 36.7 184S1265P (TCC→CCC) ‡rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (77/88);  new base C (11/8);  total (88/96)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 187 total aligned reads displayed.

AGCTGGCCGACCTGCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAA  >  NC_002947/541269‑541566
                                                                                                                                                     |                                                                                                                                                    
   tGGCCGACCTGCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCt                                                                                                                                                   >  1:5096689/1‑150 (MQ=255)
    ggCCGACCTGCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTa                                                                                                                                                  >  1:1223030/1‑150 (MQ=255)
           ccgCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTACGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACTACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTg                                                                                                                                           <  1:4310034/148‑1 (MQ=255)
            tGCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTgg                                                                                                                                          >  1:4253220/1‑150 (MQ=255)
                   gagCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCca                                                                                                                                   <  1:773299/150‑1 (MQ=255)
                      cGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCcagca                                                                                                                                <  2:296192/150‑1 (MQ=255)
                        gCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCagcagc                                                                                                                              <  2:1381317/150‑1 (MQ=255)
                              tGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTgg                                                                                                                        <  2:3972358/150‑1 (MQ=255)
                              tGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTgg                                                                                                                        <  2:4368360/150‑1 (MQ=255)
                              tGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTgg                                                                                                                        <  2:1682223/150‑1 (MQ=255)
                               ggTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTggg                                                                                                                       <  1:4796967/150‑1 (MQ=255)
                                 tGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGggtg                                                                                                                     >  1:4966894/1‑150 (MQ=255)
                                  gCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggtaa                                                                                                                    <  2:4001217/150‑3 (MQ=255)
                                  gCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGggtgg                                                                                                                    <  2:468795/150‑1 (MQ=255)
                                   cTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggtaac                                                                                                                   <  1:882298/150‑4 (MQ=255)
                                    tGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTa                                                                                                                  >  1:2462606/1‑150 (MQ=255)
                                       tCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAgg                                                                                                               <  2:701609/150‑1 (MQ=255)
                                         gATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcg                                                                                                             >  2:2326284/1‑150 (MQ=255)
                                         gATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcg                                                                                                             >  1:547765/1‑150 (MQ=255)
                                               cGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTc                                                                                                       <  2:1394509/150‑1 (MQ=255)
                                                  aCCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggt                                                                                                    <  1:4005143/150‑1 (MQ=255)
                                                   ccGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggtg                                                                                                   <  2:3274547/150‑1 (MQ=255)
                                                   ccGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggtg                                                                                                   <  1:1891681/150‑1 (MQ=255)
                                                   ccGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggtg                                                                                                   >  1:4618335/1‑150 (MQ=255)
                                                         acaaGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACACGATGCTCGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGc                                                                                             <  2:4790390/150‑1 (MQ=255)
                                                         acaaGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGc                                                                                             >  2:4275730/1‑150 (MQ=255)
                                                          caaGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCg                                                                                            >  2:3115610/1‑150 (MQ=255)
                                                          caaGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCg                                                                                            >  2:3115631/1‑150 (MQ=255)
                                                              ttCGAGCGTCCTGTGACCGTTGGTTACATGTACCTGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTc                                                                                        <  1:670501/150‑1 (MQ=255)
                                                                cGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                      <  2:1910039/150‑1 (MQ=255)
                                                                cGACCGTCAGGTGACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCgg                                                                                      >  1:2379402/1‑150 (MQ=255)
                                                                agatcGTCCTGTGACAGTTGTTTACATGTACATGCTCAAGCTTAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                      <  1:4706619/146‑1 (MQ=255)
                                                                 gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                     <  1:3582312/150‑1 (MQ=255)
                                                                 gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                     >  1:2321056/1‑150 (MQ=255)
                                                                 gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                     <  2:2273755/150‑1 (MQ=255)
                                                                 gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                     >  2:1377603/1‑150 (MQ=255)
                                                                 gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                     <  1:815241/150‑1 (MQ=255)
                                                                   gCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGga                                                                                   <  1:3081372/150‑1 (MQ=255)
                                                                    cgtcaggtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCgaa                                                                                  >  2:4122274/7‑149 (MQ=255)
                                                                     gtcaggtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCgaaa                                                                                 >  1:4135894/6‑148 (MQ=255)
                                                                     gtcaggtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCgaaa                                                                                 >  1:1358858/6‑148 (MQ=255)
                                                                     gTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGgagt                                                                                 >  2:2533232/1‑149 (MQ=255)
                                                                      tCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGAt                                                                                <  1:3023591/150‑1 (MQ=255)
                                                                        cTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATgg                                                                              >  2:1242060/1‑150 (MQ=255)
                                                                         tgtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGa                                                                             >  2:4695903/1‑150 (MQ=255)
                                                                         tgtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGa                                                                             >  1:1248805/1‑150 (MQ=255)
                                                                          gtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGaa                                                                            <  1:3760505/150‑1 (MQ=255)
                                                                          gtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGaa                                                                            >  2:2285071/1‑150 (MQ=255)
                                                                          gtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGag                                                                            <  2:2119363/150‑2 (MQ=255)
                                                                          gtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGag                                                                            <  1:939919/150‑2 (MQ=255)
                                                                           tgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAg                                                                           >  2:3023082/1‑150 (MQ=255)
                                                                              ccGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAgtgt                                                                        >  2:2835390/1‑150 (MQ=255)
                                                                               cGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGtgtg                                                                       >  2:2256031/1‑150 (MQ=255)
                                                                                gTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTgg                                                                      >  1:2060428/1‑150 (MQ=255)
                                                                                 ttGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTggg                                                                     >  2:948561/1‑150 (MQ=255)
                                                                                    gTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGgcgc                                                                  >  2:5007926/1‑150 (MQ=255)
                                                                                     ttACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCt                                                                 <  2:387088/150‑1 (MQ=255)
                                                                                      tacatctatatgTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTg                                                                <  2:242477/150‑1 (MQ=255)
                                                                                       acatgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTgg                                                               <  1:288845/150‑1 (MQ=255)
                                                                                         atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                             <  2:1149519/150‑1 (MQ=255)
                                                                                         atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                             <  2:3465683/150‑1 (MQ=255)
                                                                                         atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                             <  1:2550726/150‑1 (MQ=255)
                                                                                         atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                             <  1:4495470/150‑1 (MQ=255)
                                                                                          tgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAg                                                                                                                <  1:2920105/98‑1 (MQ=255)
                                                                                          tgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAg                                                            >  1:1047128/1‑150 (MQ=255)
                                                                                             acatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCAt                                                         <  2:3081585/150‑1 (MQ=255)
                                                                                                 gCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACgg                                                     <  1:2444006/150‑1 (MQ=255)
                                                                                                 gCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACgg                                                     >  2:831023/1‑150 (MQ=255)
                                                                                                     aaGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGtt                                                                                                        <  2:2397931/95‑1 (MQ=255)
                                                                                                      agctccaccaccTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCgg                                                >  1:1759016/7‑150 (MQ=255)
                                                                                                      aGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCgg                                                <  2:2927816/150‑1 (MQ=255)
                                                                                                       gCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGc                                               >  1:1792976/1‑150 (MQ=255)
                                                                                                          gAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATa                                            >  2:2837723/1‑150 (MQ=255)
                                                                                                           caccaccTGGTCGACGACAAGATCCATGCTCGTTAGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTac                                           >  2:3786283/2‑150 (MQ=255)
                                                                                                           aaCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATac                                           >  2:983391/1‑150 (MQ=255)
                                                                                                             ccACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATacac                                         >  2:2912443/1‑150 (MQ=255)
                                                                                                               aCTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACAccc                                       >  1:1688751/1‑150 (MQ=255)
                                                                                                                  tggtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGc                                    <  1:2285071/150‑1 (MQ=255)
                                                                                                                  tggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGc                                    <  1:5028228/150‑1 (MQ=255)
                                                                                                                    gtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                  <  1:2244274/150‑1 (MQ=255)
                                                                                                                    gtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                  <  1:4815179/150‑1 (MQ=255)
                                                                                                                    gtggACGACAAGATGCACGCGCGGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                  <  2:2657645/150‑1 (MQ=255)
                                                                                                                      ggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGACGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGa                                >  1:717675/1‑150 (MQ=255)
                                                                                                                       gacgacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCa                                                                                               >  1:1192521/1‑86 (MQ=255)
                                                                                                                       gacgacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCa                                                                                               <  2:1192533/86‑1 (MQ=255)
                                                                                                                          gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTGGGGGAGATGGATGTGTGGGCGCTGGTAGCATACGGCGCGGCATAGACCCTGGAAGAAATg                            >  2:4419005/1‑150 (MQ=255)
                                                                                                                          gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATg                            <  1:417624/150‑1 (MQ=255)
                                                                                                                          gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATg                            >  1:216758/1‑150 (MQ=255)
                                                                                                                                aTGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTcaca                      <  2:1769069/150‑1 (MQ=255)
                                                                                                                                     cgcgcgGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa                 <  1:4445926/150‑1 (MQ=255)
                                                                                                                                      gcgcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAg                >  2:5015831/1‑150 (MQ=255)
                                                                                                                                        gcgGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTc              <  2:5096689/150‑1 (MQ=255)
                                                                                                                                          gTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCgg            >  1:3269723/1‑150 (MQ=255)
                                                                                                                                            tCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgac          >  1:2100529/1‑150 (MQ=255)
                                                                                                                                              cACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacga        >  1:3057273/1‑150 (MQ=255)
                                                                                                                                               aCTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacgac       >  2:2797303/1‑150 (MQ=255)
                                                                                                                                                cTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGacgacg      >  2:587947/1‑150 (MQ=255)
                                                                                                                                                   gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGa   >  1:3066885/1‑150 (MQ=255)
                                                                                                                                                   gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGAGGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGa   >  1:1108134/1‑150 (MQ=255)
                                                                                                                                                    ttCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGaa  <  2:4333/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCTGGCCGACCTGCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAA  >  NC_002947/541269‑541566

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