Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,876,562 Δ1 bp 7.0% coding (1010/1212 nt) xcpS → type II secretion pathway protein XcpS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,876,5620T. 7.0% 348.2 / 20.4 115coding (1010/1212 nt)xcpStype II secretion pathway protein XcpS
Reads supporting (aligned to +/- strand):  ref base T (41/66);  new base . (4/4);  total (45/70)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

Only 100 of 123 total aligned reads displayed.

CTCACCGCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCCTGCTGATCGGCCTGA  >  NC_002947/3876414‑3876711
                                                                                                                                                    |                                                                                                                                                      
cTCACCGCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCctgctg                                                                                                                                                       >  1:2743131/1‑150 (MQ=255)
 tCACCGCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCtgctgc                                                                                                                                                      >  1:458310/1‑150 (MQ=255)
 tCACCGCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCtgctgc                                                                                                                                                      >  1:983350/1‑150 (MQ=255)
 gctccGCCATTGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCtgctgc                                                                                                                                                      >  2:43762/4‑150 (MQ=255)
    ccGCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTgcgcg                                                                                                                                                  >  1:162209/1‑150 (MQ=255)
         aTGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGc                                                                                                                                             >  1:772631/1‑150 (MQ=255)
         aTGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGc                                                                                                                                             <  2:2130340/150‑1 (MQ=255)
         aTGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGc                                                                                                                                             <  2:2857279/150‑1 (MQ=255)
         aTGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGATCGTACAAGGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGc                                                                                                                                             <  1:2065494/150‑1 (MQ=255)
         aTGGGCATGGCCCGCGGCCTGCTCGGCAACACAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCACTGGCGGCCGGGCACCGGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAggg                                                                                                                                             >  2:541925/1‑149 (MQ=255)
           gggCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGa                                                                                                                                           <  1:1271026/150‑1 (MQ=255)
           gggCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGa                                                                                                                                           <  1:774362/150‑1 (MQ=255)
           gggCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGa                                                                                                                                           <  1:2474909/150‑1 (MQ=255)
           gggCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGa                                                                                                                                           <  1:278062/150‑1 (MQ=255)
           gggCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGAATTGTTCTCGCAGGCGa                                                                                                                                           <  2:908357/150‑1 (MQ=255)
             gCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGc                                                                                                                                         >  2:679485/1‑150 (MQ=255)
             gCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGc                                                                                                                                         >  2:1083530/1‑150 (MQ=255)
              cATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCa                                                                                                                                        <  1:2813389/150‑1 (MQ=255)
              cATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCa                                                                                                                                        <  2:2873257/150‑1 (MQ=255)
               aTGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGACGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAg                                                                                                                                       >  2:416999/1‑150 (MQ=255)
                   cccGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTg                                                                                                                                   <  2:1430457/150‑1 (MQ=255)
                   cccGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTg                                                                                                                                   >  1:2783409/1‑150 (MQ=255)
                         gCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAAcc                                                                                                                             <  1:1741026/150‑1 (MQ=255)
                           cTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTg                                                                                                                           >  1:1770498/1‑150 (MQ=255)
                            tacTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTgg                                                                                                                          <  2:1107170/148‑1 (MQ=255)
                            tGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGTGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCGCGCAGCAGGCgccgg                                                                                                                                        <  1:3067430/136‑5 (MQ=255)
                            tGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGTGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCGCGCAGCAGGCgcagg                                                                                                                                        >  2:3067430/1‑132 (MQ=255)
                            tGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTgg                                                                                                                          <  2:2743131/150‑1 (MQ=255)
                            tGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTgg                                                                                                                          <  1:1206704/150‑1 (MQ=255)
                               tCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCg                                                                                                                       >  2:719253/1‑150 (MQ=255)
                                 ggCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCCGCTGCCCGCAGCCCCGCACTCTGGCCCCCTgcccgcc                                                                                                                     >  2:455413/1‑144 (MQ=255)
                                  gCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCgaga                                                                                                                    >  1:268423/1‑150 (MQ=255)
                                       gcagcaGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTg                                                                                                               <  2:458310/150‑1 (MQ=255)
                                       gcagcaGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTg                                                                                                               >  1:3012378/1‑150 (MQ=255)
                                        cagcagCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTgg                                                                                                              <  2:772631/150‑1 (MQ=255)
                                        cagcagCGCAAGGCCTGGAGCAGGCCAGCCAGCGTGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCTGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTgg                                                                                                              >  2:1817954/1‑150 (MQ=255)
                                              cgcAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTAGCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCt                                                                                                        <  1:1928395/150‑1 (MQ=255)
                                                   ggCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCt                                                                                                                                    >  1:242974/1‑117 (MQ=255)
                                                   ggCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCt                                                                                                                                    <  2:242974/117‑1 (MQ=255)
                                                     ccTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGa                                                                                                 <  1:1100676/150‑1 (MQ=255)
                                                        ggAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCGCGCAGCAGGCgcagg                                                                                                                                        >  1:1664860/1‑104 (MQ=255)
                                                        ggAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCGCGCAGCAGGCgcagg                                                                                                                                        <  2:1664860/108‑5 (MQ=255)
                                                          aGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTcc                                                                                            >  2:2435410/1‑150 (MQ=255)
                                                            cAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCAt                                                                                          <  2:162209/150‑1 (MQ=255)
                                                              ggCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGa                                                                                        >  1:2541359/1‑150 (MQ=255)
                                                               gccagccaGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCAt                                                                                          <  1:1825808/147‑1 (MQ=255)
                                                               gccagccaGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCAt                                                                                          >  2:1825808/1‑147 (MQ=255)
                                                               gccagccaGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGAc                                                                                       <  2:671347/150‑1 (MQ=255)
                                                                 cagccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCa                                                                                     >  1:2455358/1‑150 (MQ=255)
                                                                  agccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGTCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAg                                                                                    <  1:1083530/150‑1 (MQ=255)
                                                                  agccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAg                                                                                    >  2:1707531/1‑150 (MQ=255)
                                                                    ccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTgg                                                                                                                          <  1:1821493/110‑1 (MQ=255)
                                                                    ccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTgg                                                                                                                          >  2:1821493/1‑110 (MQ=255)
                                                                    ccagcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGa                                                                                  >  2:1657566/1‑150 (MQ=255)
                                                                      agcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAAc                                                                                                                              >  2:2977497/1‑104 (MQ=255)
                                                                       gcGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAAt                                                                               <  1:719253/150‑1 (MQ=255)
                                                                         gggTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCg                                                                             <  2:614681/150‑1 (MQ=255)
                                                                               gCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGtt                                                                       <  1:1624875/150‑1 (MQ=255)
                                                                                 gAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGa                                                                                                                                           >  1:1315731/1‑80 (MQ=255)
                                                                                 gAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGTGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTgg                                                                     <  2:2281054/150‑1 (MQ=255)
                                                                                 gAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTgg                                                                     <  1:2435410/150‑1 (MQ=255)
                                                                                      gCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGa                                                                >  1:1434218/1‑150 (MQ=255)
                                                                                       cTGCCGCTTTCCGATGCTCTGGCGGCCGGGCCCCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGACGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGaa                                                               <  2:1323557/150‑1 (MQ=255)
                                                                                        tGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAAt                                                              <  2:1434218/150‑1 (MQ=255)
                                                                                          ccGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGa                                                                                        <  2:2131207/122‑1 (MQ=255)
                                                                                          ccGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGa                                                                                        >  1:2131207/1‑122 (MQ=255)
                                                                                           cGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGAGCTGCGCCGACTTCCATGACCAGGAAGTCGGCCGTTGGGTGGAATGGt                                                           >  2:2061999/1‑150 (MQ=255)
                                                                                              tttCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCg                                                        <  2:193879/150‑1 (MQ=255)
                                                                                                tCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTc                                                      <  1:679485/150‑1 (MQ=255)
                                                                                                  cGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGaa                                                               >  2:2494528/1‑139 (MQ=255)
                                                                                                   gATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCaaa                                                   >  2:3033656/1‑150 (MQ=255)
                                                                                                        tctAGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCATGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGtt                                              <  2:576731/150‑1 (MQ=255)
                                                                                                         ctGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTc                                             <  1:3033656/150‑1 (MQ=255)
                                                                                                          tGGCGGCCGGGCACCTGAGGACGCCGGTATACCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCg                                            >  1:884949/1‑150 (MQ=255)
                                                                                                                ccGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAACCTGTTCGCACAgc                                      <  2:2352010/150‑1 (MQ=255)
                                                                                                                  gggCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGAc                                                                                       <  2:1469617/99‑1 (MQ=255)
                                                                                                                  gggCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGAc                                                                                       >  1:1469617/1‑99 (MQ=255)
                                                                                                                  gggCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCg                                                                             <  1:436784/109‑1 (MQ=255)
                                                                                                                      aCCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCtga                                <  1:1707531/150‑1 (MQ=255)
                                                                                                                        cTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGtt                                              <  2:2172587/134‑1 (MQ=255)
                                                                                                                        cTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTGCGTCAAACTGtt                                              >  1:2172587/1‑134 (MQ=255)
                                                                                                                         tGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCtgatga                             <  1:1501148/150‑1 (MQ=255)
                                                                                                                          ggTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGAc                            <  1:3035600/150‑1 (MQ=255)
                                                                                                                           gTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGAcc                           <  2:2783409/150‑1 (MQ=255)
                                                                                                                           gTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGAcc                           >  2:2107544/1‑150 (MQ=255)
                                                                                                                             gaggCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGTAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCtt                         <  1:43762/147‑1 (MQ=255)
                                                                                                                                gCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCtgatga                             >  2:269227/1‑143 (MQ=255)
                                                                                                                                gCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCtgatga                             <  1:269227/143‑1 (MQ=255)
                                                                                                                                    gTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGc                  <  2:980305/150‑1 (MQ=255)
                                                                                                                                        cccTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCctgc              <  2:3012378/150‑1 (MQ=255)
                                                                                                                                         ccTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCctgct             >  1:2398462/1‑150 (MQ=255)
                                                                                                                                         ccTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCctgct             >  1:3004343/1‑150 (MQ=255)
                                                                                                                                         ccTGCACCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGCTCGAACCGCTGCTGATGACCTTCATCGGCctgct             >  1:1178197/1‑150 (MQ=255)
                                                                                                                                          cTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATgg                                                            <  1:2952244/102‑1 (MQ=255)
                                                                                                                                          cTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATgg                                                            >  2:2952244/1‑102 (MQ=255)
                                                                                                                                            tcgcCTGCGCGCAGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGtt                                                          <  2:1203935/101‑1 (MQ=255)
                                                                                                                                            tcgcCTGCGCGCAGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGtt                                                          >  1:1203935/2‑102 (MQ=255)
                                                                                                                                            tcgcCTGCGCGCAGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATc                     <  2:800060/138‑1 (MQ=255)
                                                                                                                                             cgcCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCCTGCTGATc         <  2:2541359/150‑1 (MQ=255)
                                                                                                                                                ctgctgCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCCTGCTGATCGGc      <  1:369376/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTCACCGCCATGGGCATGGCCCGCGGCCTGCTCGGCAACGCAGCAGCGCAGGGCCTGGAGCAGGCCAGCCAGCGGGTGGGCGAAGGGCTGCCGCTTTCCGATGCTCTGGCGGCCGGGCACCTGGTGACGCCGGTATCCCTGCGCCTGCTGCGCGCAGGCGAGCAGTCTGGCAACCTGGGCGAGATGCTGGAGCGCTGCGCCGACTTCCATGACCAGGAAATCGGCCGTTGGGTGGAATGGTTCGTCAAACTGTTCGAACCGCTGCTGATGACCTTCATCGGCCTGCTGATCGGCCTGA  >  NC_002947/3876414‑3876711

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