Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,414,217 G→A 5.8% S86S (AGC→AGT) ‡ PP_1236 ← glycine cleavage system transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,414,2170GA 5.8% 589.6 / 15.7 189S86S (AGC→AGT) ‡PP_1236glycine cleavage system transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (99/79);  new base A (7/4);  total (106/83)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

GCGCCAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGACCACCGC  >  NC_002947/1414069‑1414366
                                                                                                                                                    |                                                                                                                                                     
 cgcCAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCt                                                                                                                                                     <  1:2070522/150‑1 (MQ=255)
    cAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCg                                                                                                                                                  >  1:1229776/1‑150 (MQ=255)
     aggtaggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCgg                                                                                                                                                 <  1:3443048/150‑1 (MQ=255)
      ggtaggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGa                                                                                                                                                <  2:4787138/150‑1 (MQ=255)
       gtaggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGAc                                                                                                                                               >  1:4976427/1‑150 (MQ=255)
        taggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGAcc                                                                                                                                              <  2:859323/150‑1 (MQ=255)
          ggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCAc                                                                                                                                            >  2:4643210/1‑150 (MQ=255)
            tgtCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCagga                                                                                                                                          <  2:2176034/150‑4 (MQ=255)
             gtCGCAGGTCATGGCTTCCAGCTCGACGTGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCaggaa                                                                                                                                         <  1:1001541/150‑5 (MQ=255)
              tCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCg                                                                                                                                        >  2:1349925/1‑150 (MQ=255)
                gCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAg                                                                                                                                      >  2:1183689/1‑150 (MQ=255)
                gCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAg                                                                                                                                      >  1:4689384/1‑150 (MQ=255)
                gCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAg                                                                                                                                      >  1:3122416/1‑150 (MQ=255)
                 cAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAgg                                                                                                                                     >  1:2045707/1‑150 (MQ=255)
                     tCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCa                                                                                                                                 >  2:567000/1‑150 (MQ=255)
                       aTGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAgg                                                                                                                               >  2:3379265/1‑150 (MQ=255)
                           cTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGt                                                                                                                           <  1:4084521/150‑1 (MQ=255)
                             tCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGc                                                                                                                         <  1:3597760/150‑1 (MQ=255)
                             gccAGCTCGCCGCGGTGGTCACGGACGAACGGGCCCAGCTCGTTGATGCGGTCCGGGCGGTAGGCAGCGCTCACGTCGGCCACATCCGGACGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGc                                                                                                                         <  2:1252856/149‑1 (MQ=255)
                                  cTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGcttctt                                                                                                                    >  1:242681/1‑150 (MQ=255)
                                   tCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTg                                                                                                                   >  2:2865901/1‑150 (MQ=255)
                                    cGACGGGGTGGTCAAGGAATAACTGGCACAGCTCTTTTTTGATGTCCGGGCGGTAGGCAGCGCTCACTTCGGCCACATCCGGAAGGGCCTGGGGGCGCACTTCGTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGc                                                                                                                  <  1:4223793/150‑1 (MQ=255)
                                      aCGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGccc                                                                                                                >  2:3591628/1‑150 (MQ=255)
                                                cAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCa                                                                                                      >  2:3474354/1‑150 (MQ=255)
                                                 aaGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCACGCCTGGCAg                                                                                                     >  2:4812502/1‑150 (MQ=255)
                                                 aaGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAg                                                                                                     >  1:5051118/1‑150 (MQ=255)
                                                 aaGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAg                                                                                                     >  2:3119559/1‑150 (MQ=255)
                                                 aaGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAg                                                                                                     >  1:5067559/1‑150 (MQ=255)
                                                 aaGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAg                                                                                                     >  1:4511036/1‑150 (MQ=255)
                                                  aGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAgg                                                                                                    <  2:2942307/150‑1 (MQ=255)
                                                  aGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAgg                                                                                                    >  1:1411313/1‑150 (MQ=255)
                                                  aGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAgg                                                                                                    >  2:2523789/1‑150 (MQ=255)
                                                    gaataaCTGGCACAGCTGGTTGTTGTTGTCGGGGGGGTAGGCAGGGCTCACGTAGGCCACATAGGGACGGGCCGGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGt                                                                                                  <  1:3276581/150‑1 (MQ=255)
                                                    gaagaaCTGGCACAGCTCGTTGATGATGCCCGGGCGGTAGGCTGCGCTCACGTAGGCCACATACGTAAGGGCCTGCGGGCGCAATTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGt                                                                                                  <  1:3972844/150‑1 (MQ=255)
                                                       gaaCTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGa                                                                                               >  1:456465/1‑150 (MQ=255)
                                                       gaaCGGGCACAGCTGGTTTTTGAGGTCCGGGGGGTAGGCAGCGCTCACGTCGGCCACATACGGACGGGCCTGGGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGa                                                                                               <  1:2367303/150‑1 (MQ=255)
                                                           aggCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGACGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCg                                                                                           <  1:3416581/149‑1 (MQ=255)
                                                            ggCACTGCTCGTTGATGATGTCGGGGCGGTCGGCAGCGCTCACGTAGGCCACATCCGGATGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCga                                                                                          <  1:4853117/150‑1 (MQ=255)
                                                            ggCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCga                                                                                          <  2:5051118/150‑1 (MQ=255)
                                                             gCACAGCTGGTTGATGTTGTCCGGGCGGTAGGCAGCGCTCACGTTGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCgag                                                                                         <  1:2826205/150‑1 (MQ=255)
                                                              cacaGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCgaga                                                                                        <  2:3499347/150‑1 (MQ=255)
                                                                  gCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTg                                                                                    >  1:121874/1‑150 (MQ=255)
                                                                   cTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGc                                                                                   <  1:3379265/150‑1 (MQ=255)
                                                                    tCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGcc                                                                                  >  2:252048/1‑150 (MQ=255)
                                                                      gTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAg                                                                                >  2:3040453/1‑150 (MQ=255)
                                                                       ttGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                               <  2:970609/150‑1 (MQ=255)
                                                                       ttGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGCCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                               <  2:3331003/150‑1 (MQ=255)
                                                                               gTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCa                                                                       >  2:5141589/1‑150 (MQ=255)
                                                                                     gCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTgcc                                                                 <  2:4130259/150‑1 (MQ=255)
                                                                                              aGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTg                                                        >  2:1014265/1‑150 (MQ=255)
                                                                                                cgcTCCCGTCGGCCACATCCGGACGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgca                                                      <  1:2307267/150‑1 (MQ=255)
                                                                                                cgcTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgca                                                      <  1:3749225/150‑1 (MQ=255)
                                                                                                 gcTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTg                                                                                                 >  2:4546613/1‑106 (MQ=255)
                                                                                                 gcTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcag                                                     >  1:2524301/1‑150 (MQ=255)
                                                                                                  cTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcagc                                                    <  2:4024169/150‑1 (MQ=255)
                                                                                                   tCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcagca                                                   >  1:1206912/1‑150 (MQ=255)
                                                                                                   tCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcagca                                                   >  2:5173672/1‑150 (MQ=255)
                                                                                                       gTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAg                                               >  1:3571044/1‑150 (MQ=255)
                                                                                                         aGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCCAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAggg                                             >  2:1338843/1‑150 (MQ=255)
                                                                                                           gCCACATACGGAAGGGCGTGGGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGTCGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGgcg                                           <  2:2524301/150‑1 (MQ=255)
                                                                                                            ccACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGgcgc                                          >  1:2931097/1‑150 (MQ=255)
                                                                                                              acaTACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTg                                        >  1:3234818/1‑150 (MQ=255)
                                                                                                              acaTACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTg                                        >  1:962764/1‑150 (MQ=255)
                                                                                                                   cGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGtctc                                   <  2:273028/150‑1 (MQ=255)
                                                                                                                   cGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGtctc                                   <  1:4521680/150‑1 (MQ=255)
                                                                                                                    ggACGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCg                                  <  1:2384333/150‑1 (MQ=255)
                                                                                                                    ggAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCg                                  >  2:5022090/1‑150 (MQ=255)
                                                                                                                     gAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGc                                 >  1:2075033/1‑150 (MQ=255)
                                                                                                                      aaGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGcc                                >  1:1097676/1‑150 (MQ=255)
                                                                                                                        gggCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGt                              <  2:5214030/150‑1 (MQ=255)
                                                                                                                        gggCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGt                              <  2:1411313/150‑1 (MQ=255)
                                                                                                                         ggCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCgcg                             >  2:158014/1‑148 (MQ=255)
                                                                                                                           ccTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggc                           <  2:1794107/150‑1 (MQ=255)
                                                                                                                            cTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggcg                          >  1:1631283/1‑150 (MQ=255)
                                                                                                                            cTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggcg                          >  2:902729/1‑150 (MQ=255)
                                                                                                                             tGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggcgg                         <  2:1463693/150‑1 (MQ=255)
                                                                                                                                  gcgcACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAg                    <  2:3319290/150‑1 (MQ=255)
                                                                                                                                  gcgcACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAg                    >  2:4813969/1‑150 (MQ=255)
                                                                                                                                   cgcACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAgg                   <  2:1875428/150‑1 (MQ=255)
                                                                                                                                   cgcACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAgg                   <  2:767351/150‑1 (MQ=255)
                                                                                                                                      aCTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGgcgc                >  1:5240748/1‑150 (MQ=255)
                                                                                                                                      aCTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGgcgc                <  2:5349941/150‑1 (MQ=255)
                                                                                                                                        ttCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGTAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGa              >  2:3054446/1‑150 (MQ=255)
                                                                                                                                        ttCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGa              <  1:822987/150‑1 (MQ=255)
                                                                                                                                        ttCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGa              >  2:354919/1‑150 (MQ=255)
                                                                                                                                        ttCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGa              >  2:634434/1‑150 (MQ=255)
                                                                                                                                        ttCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGa              >  1:1724672/1‑150 (MQ=255)
                                                                                                                                        ttCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGa              <  1:4050621/150‑1 (MQ=255)
                                                                                                                                         acgtGGTCCGCAGCGCGGACCAGGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAg             >  1:307593/4‑150 (MQ=255)
                                                                                                                                         acgtGGTCCGCAGCGCGGACCACGTCGAGGGTCAGGCCGTGGTTCTTGGCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCAGCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAg             >  2:4229550/4‑150 (MQ=255)
                                                                                                                                          ccTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAgg            >  1:2483740/1‑150 (MQ=255)
                                                                                                                                           cTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGt           <  2:4418409/150‑1 (MQ=255)
                                                                                                                                            tGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTg          >  1:5403187/1‑150 (MQ=255)
                                                                                                                                            tGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTg          >  1:1680545/1‑150 (MQ=255)
                                                                                                                                            tGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTg          >  1:1957300/1‑150 (MQ=255)
                                                                                                                                            tGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTg          >  2:1309968/1‑150 (MQ=255)
                                                                                                                                            tGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTg          <  2:3038785/150‑1 (MQ=255)
                                                                                                                                             ggTCCGCAGCGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGa         >  1:3677112/1‑150 (MQ=255)
                                                                                                                                              gTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGac        >  2:799482/1‑150 (MQ=255)
                                                                                                                                                 cgcgCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGaccac     <  2:3234818/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCCAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGACCACCGC  >  NC_002947/1414069‑1414366

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