Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 384,794 C→G 5.3% A423A (GCC→GCG)  PP_0320 → methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947384,7940CG 5.3% 469.3 / 14.4 151A423A (GCC→GCG) PP_0320methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (66/76);  new base G (4/4);  total (71/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

Only 100 of 157 total aligned reads displayed.

TCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGAAGTCG  >  NC_002947/384645‑384942
                                                                                                                                                      |                                                                                                                                                    
tCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGc‑c                                                                                                                                                      >  1:1887654/1‑150 (MQ=255)
tCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGc‑c                                                                                                                                                      >  1:1132408/1‑150 (MQ=255)
tCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGc‑c                                                                                                                                                      >  1:1346050/1‑150 (MQ=255)
 ccGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑Cg                                                                                                                                                     <  1:1365777/150‑1 (MQ=255)
  cGAGAGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑Cgg                                                                                                                                                    <  1:297179/150‑1 (MQ=255)
   gAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGc                                                                                                                                                   <  2:2274143/150‑1 (MQ=255)
      cGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCGGCATTGC‑GGGCGAt                                                                                                                                                <  1:499776/150‑1 (MQ=255)
         cTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCCCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATgtc                                                                                                                                             <  2:1498250/150‑3 (MQ=255)
          ggcggcgGGCCGTGGAGCAGGTAACGCAGGTGGTGGACAAGACCGCGGAGGGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCgg                                                                                                                                            <  2:1942411/149‑1 (MQ=255)
              aCGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCa                                                                                                                                        >  1:459416/1‑150 (MQ=255)
                 ggCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGNCTGTCCNTGCAGGACATCGNCCGCAATGC‑CGGCGATGCGGCCCAgtc                                                                                                                                     <  1:2238500/150‑3 (MQ=255)
                  gCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCg                                                                                                                                    >  2:300536/1‑150 (MQ=255)
                  gCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCg                                                                                                                                    >  2:2008154/1‑150 (MQ=255)
                  gCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCg                                                                                                                                    >  2:2617280/1‑150 (MQ=255)
                   ccGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGt                                                                                                                                   <  1:1931324/150‑1 (MQ=255)
                    cGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGATCGCTCCGGGCCCCAGCAGAAGATGAACGAAATGGTCGCCACTGCTGTACACGAGCTGGCCCTGACCGTCCAGGACATCGCACGCAATGC‑CGGCGATGCGGCCCAGGCGTc                                                                                                                                  <  1:861302/150‑1 (MQ=255)
                          gCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTc                                                                                                                            <  1:2822751/150‑1 (MQ=255)
                                       gtggtgGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgagacgag                                                                                                               <  1:670353/150‑1 (MQ=255)
                                        tggtggACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAgg                                                                                                              >  1:2424219/1‑150 (MQ=255)
                                           tggACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGt                                                                                                           <  2:2343216/150‑1 (MQ=255)
                                            ggACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGtt                                                                                                          >  2:877881/1‑150 (MQ=255)
                                             gACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCac                                                                                                                       <  2:815809/136‑1 (MQ=255)
                                                aacaCCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTTCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAg                                                                                                      <  1:2767299/150‑1 (MQ=255)
                                                  cacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCTACCGCTGTACACGAGATGGGCCTGCCCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTAGCAGTCGGGAAGAGCGGCGGCGTTGGAggg                                                                                                    >  2:642070/1‑149 (MQ=255)
                                                  cacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGc                                                                                                    >  1:1714768/1‑150 (MQ=255)
                                                  cacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGc                                                                                                    >  2:1522519/1‑150 (MQ=255)
                                                  cacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGc                                                                                                    >  1:1269939/1‑150 (MQ=255)
                                                   acCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCgt                                                                                                                          >  1:2669586/1‑126 (MQ=255)
                                                   acCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCgt                                                                                                                          <  2:2669586/127‑2 (MQ=255)
                                                   acCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGcc                                                                                                   >  1:1208137/1‑150 (MQ=255)
                                                     cgccgAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCcg                                                                                                 >  2:2519519/1‑150 (MQ=255)
                                                       ccgAGCGCTCCGGGCGCCAGCAGGAGCTGCCGGCAATGGTCGCCACGGCTGGACACGAAATAGGCCTGACAGTAC‑GCCGAACACAAGTGGGGCTCCGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTGGCAGTCCcgcg                                                                                               >  1:2995837/1‑150 (MQ=255)
                                                       ccgAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCcgcg                                                                                               >  2:2565800/1‑150 (MQ=255)
                                                        cgAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGa                                                                                              <  1:33637/150‑1 (MQ=255)
                                                         gAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAg                                                                                             >  1:496512/1‑150 (MQ=255)
                                                          ggcgcTCCGCGCGCCAGCAGGAGATGACCGAAATTCTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAgg                                                                                            <  1:642070/149‑1 (MQ=255)
                                                              cTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggtggt                                                                                        >  1:928217/1‑150 (MQ=255)
                                                              cTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggtggt                                                                                        >  2:1717630/1‑150 (MQ=255)
                                                              cTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggtggt                                                                                        <  2:928217/150‑1 (MQ=255)
                                                               tCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGc‑c                                                                                                                                                      >  2:2977935/1‑87 (MQ=255)
                                                                     cgcCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGa                                                                                 <  2:459416/150‑1 (MQ=255)
                                                                      gcCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCa                                                                                                                                        <  2:660349/94‑1 (MQ=255)
                                                                      gcCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCa                                                                                                                                        >  1:660349/1‑94 (MQ=255)
                                                                      gcCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGAt                                                                                <  2:2513823/150‑1 (MQ=255)
                                                                      gcCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGAt                                                                                <  1:2273639/150‑1 (MQ=255)
                                                                          gcagGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCAt                                                                            >  1:1006293/1‑150 (MQ=255)
                                                                            agGGGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGc                                                                                                                         >  2:1230693/1‑103 (MQ=255)
                                                                            agGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTc                                                                          >  2:3037208/1‑150 (MQ=255)
                                                                              gagaTGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGt                                                                        <  1:363771/150‑1 (MQ=255)
                                                                              gagaTGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGt                                                                        >  2:1905436/1‑150 (MQ=255)
                                                                               agaTGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTg                                                                       <  1:198987/150‑1 (MQ=255)
                                                                                   gACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGGCCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgaga                                                                                                                   >  2:2144621/1‑102 (MQ=255)
                                                                                   gACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGGCCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgaga                                                                                                                   <  1:2144621/102‑1 (MQ=255)
                                                                                   gACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCAt                                                                   <  1:2028672/150‑1 (MQ=255)
                                                                                     ccGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCg                                                                 >  1:1224966/1‑150 (MQ=255)
                                                                                      cGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGa                                                                >  2:511449/1‑150 (MQ=255)
                                                                                      cGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGa                                                                >  1:2425601/1‑150 (MQ=255)
                                                                                        aaaTGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAgg                                                              <  2:1006293/150‑1 (MQ=255)
                                                                                        aaaTGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAgg                                                              <  1:2553625/150‑1 (MQ=255)
                                                                                                  cACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggt                                                                                           <  2:1188601/111‑1 (MQ=255)
                                                                                                  cACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTTGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCggg                                                    >  2:1175480/1‑150 (MQ=255)
                                                                                                     tgctgTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAGCGATCCATTCGTGGCATCGAGGGCATGTCGGGCGa                                                 >  1:2174044/1‑150 (MQ=255)
                                                                                                     tgctgTACACGAGATGGGCCGGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGa                                                 <  1:1175440/150‑1 (MQ=255)
                                                                                                          tACACGAGATGGGCGTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCg                                            <  1:1002370/150‑1 (MQ=255)
                                                                                                            cacGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAg                                                                                                      >  1:1134130/1‑90 (MQ=255)
                                                                                                            cacGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAg                                                                                                      <  2:1134130/90‑1 (MQ=255)
                                                                                                            cacGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGc                                          >  1:1948972/1‑150 (MQ=255)
                                                                                                             acGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCa                                         >  1:2874055/1‑150 (MQ=255)
                                                                                                              cGAGATGGGCCTGACCGTGCAGGACTTCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCaa                                        <  2:2241950/150‑1 (MQ=255)
                                                                                                                agaTGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAgg                                      <  2:496512/150‑1 (MQ=255)
                                                                                                                     ggCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGcc                                 <  1:2960857/150‑1 (MQ=255)
                                                                                                                     ggCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGcc                                 <  1:2526162/150‑1 (MQ=255)
                                                                                                                      gCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggtgg                                                                                         >  2:2685663/1‑93 (MQ=255)
                                                                                                                      gCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggtgg                                                                                         <  1:2685663/93‑1 (MQ=255)
                                                                                                                       ccTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGAt                                                                                >  1:376512/1‑101 (MQ=255)
                                                                                                                       ccTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGAt                                                                                <  2:376512/101‑1 (MQ=255)
                                                                                                                        ctacccGTGCAGGGTGTCGCCGTGAATAC‑CGGCGATGGGGTCGAGGCGGCGCAGTCGGCACGAGACGAGGCGTGGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGAt                              <  2:330823/146‑1 (MQ=255)
                                                                                                                        cTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGAt                              >  1:1686535/1‑150 (MQ=255)
                                                                                                                           aCCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCg                           <  1:3037208/150‑1 (MQ=255)
                                                                                                                            ccGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCAGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCgg                          <  1:1415377/150‑1 (MQ=255)
                                                                                                                             cGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGt                         >  2:1918717/1‑150 (MQ=255)
                                                                                                                               tGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTca                       <  1:1522519/150‑1 (MQ=255)
                                                                                                                               tGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTca                       <  2:1224966/150‑1 (MQ=255)
                                                                                                                                 cAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcagc                     <  2:2129674/150‑1 (MQ=255)
                                                                                                                                 cAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcagc                     >  1:2129674/1‑150 (MQ=255)
                                                                                                                                   ggACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGt                                                        >  2:1353050/1‑113 (MQ=255)
                                                                                                                                   ggACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGt                                                        <  1:1353050/113‑1 (MQ=255)
                                                                                                                                   ccgcATCGCCGGCATTGC‑GGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcagcca                   >  1:415550/4‑150 (MQ=255)
                                                                                                                                    gACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgagacgag                                                                                                               <  2:1957940/57‑1 (MQ=255)
                                                                                                                                     gcATCGCCGGCATTGC‑GGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCaa                                                                                    >  1:115057/2‑83 (MQ=38)
                                                                                                                                     gcATCGCCGGCATTGC‑GGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCaa                                                                                    <  2:115057/82‑1 (MQ=38)
                                                                                                                                       aTCGCCGGCATTGC‑GGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTg               >  1:1420302/1‑150 (MQ=255)
                                                                                                                                       aTCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTg               <  1:877881/150‑1 (MQ=255)
                                                                                                                                       aTCGCCCGAAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTg               >  2:1061017/1‑150 (MQ=255)
                                                                                                                                          gCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGcc            >  1:2020740/1‑150 (MQ=255)
                                                                                                                                           ccggCATTGC‑GGGCGATGTGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggtggt                                                                                        <  2:412451/70‑1 (MQ=255)
                                                                                                                                           cccGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGTCATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCcg           >  2:3001526/1‑150 (MQ=255)
                                                                                                                                            ccGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAg                                       <  1:1760639/121‑1 (MQ=255)
                                                                                                                                            ccGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAg                                       >  2:1760639/1‑121 (MQ=255)
                                                                                                                                                    c‑cGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGAAGTCg  >  1:1337987/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGC‑CGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGAAGTCG  >  NC_002947/384645‑384942

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