Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 417,622 C→G 5.2% P71A (CCG→GCG)  waaP → lipopolysaccharide core heptose (I) kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947417,6220CG 5.2% 256.1 / 3.0 77P71A (CCG→GCG) waaPlipopolysaccharide core heptose (I) kinase
Reads supporting (aligned to +/- strand):  ref base C (37/36);  new base G (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGGCGCCGA  >  NC_002947/417481‑417769
                                                                                                                                             |                                                                                                                                                   
gcgcTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCt                                                                                                                                              >  1:144689/1‑149 (MQ=255)
 cgcTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGGCAACTTGGCGGTgaac                                                                                                                                             <  8:127516/149‑4 (MQ=255)
  gcTGCAAGGCGAGGTCTATCGCGATCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCg                                                                                                                                            <  7:77978/149‑1 (MQ=255)
        aGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTg                                                                                                                                      >  2:219826/1‑149 (MQ=255)
         ggCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTgg                                                                                                                                     <  7:110105/149‑1 (MQ=255)
         ggCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTgg                                                                                                                                     <  2:79580/149‑1 (MQ=255)
          gCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGc                                                                                                                                    <  4:169322/149‑1 (MQ=255)
               gTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGa                                                                                                                               <  5:148569/149‑1 (MQ=255)
                    tCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggc                                                                                                                          <  6:58187/149‑1 (MQ=255)
                      gcgAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggcgg                                                                                                                        >  8:303072/1‑149 (MQ=255)
                            tGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggcgg                                                                                                                        >  1:134640/1‑143 (MQ=255)
                            tGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggcgg                                                                                                                        <  2:134640/143‑1 (MQ=255)
                             ggAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATc                                                                                                                  >  5:311739/1‑148 (MQ=255)
                                aGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAAc                                                                                                              <  3:148975/149‑1 (MQ=255)
                                       cgcACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCa                                                                                                       >  1:217832/1‑149 (MQ=255)
                                          aCGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCggg                                                                                                                                          >  7:164252/1‑111 (MQ=255)
                                          aCGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATATTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCggg                                                                                                                                          <  8:164252/111‑1 (MQ=255)
                                                 gcACCGAAGTGGCTGGCGAGGGCGTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGc                                                                                             <  1:248336/149‑1 (MQ=255)
                                                      gAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGcc                                                                                        >  6:204350/1‑149 (MQ=255)
                                                        aGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCa                                                                                                                                  >  8:268935/1‑105 (MQ=255)
                                                        aGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCa                                                                                                                                  <  7:268935/105‑1 (MQ=255)
                                                          tggctggcGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGAcc                                                                                    <  5:283693/149‑1 (MQ=255)
                                                              tggcGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGa                                                                                >  7:39530/1‑149 (MQ=255)
                                                                gcGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggcgg                                                                                                                        >  1:180985/1‑107 (MQ=255)
                                                                gcGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggcgg                                                                                                                        <  2:180985/107‑1 (MQ=255)
                                                                gcGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGAcc                                                                              <  3:116486/149‑1 (MQ=255)
                                                                gcGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGAc                                                                               <  8:39530/148‑1 (MQ=255)
                                                                  gAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCa                                                                                   >  3:148540/1‑142 (MQ=255)
                                                                  gAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCa                                                                                   <  4:148540/142‑1 (MQ=255)
                                                                     ggCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGAc                                                                                     >  1:36979/1‑137 (MQ=255)
                                                                     ggCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGAc                                                                                     <  2:36979/137‑1 (MQ=255)
                                                                     ggCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGAc                                                                                     >  1:32200/1‑137 (MQ=255)
                                                                     ggCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGAc                                                                                     <  2:32200/137‑1 (MQ=255)
                                                                        tttttCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGc                                                                      >  5:75371/1‑149 (MQ=255)
                                                                              gTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACgg                                                                <  2:144689/149‑1 (MQ=255)
                                                                               tCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGc                                                               <  6:311739/149‑1 (MQ=255)
                                                                                cAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGc                                                               >  8:289744/1‑148 (MQ=255)
                                                                                 aaGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCg                                                              <  2:217832/148‑1 (MQ=255)
                                                                                    aTCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCg                                                          >  7:71305/1‑149 (MQ=255)
                                                                                       cACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTgg                                                       <  7:289744/149‑1 (MQ=255)
                                                                                         ccGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGGTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGgca                                                     <  4:253958/149‑1 (MQ=255)
                                                                                              gCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTTGCAGCAAc                                                >  8:133287/1‑149 (MQ=255)
                                                                                              gCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAAc                                                >  8:128985/1‑149 (MQ=255)
                                                                                                aTCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAAccc                                              >  5:313496/1‑149 (MQ=255)
                                                                                                 tCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCg                                             >  1:211852/1‑149 (MQ=255)
                                                                                                    gCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGcc                                          >  1:55072/1‑149 (MQ=255)
                                                                                                        gggCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCcgcg                                       >  8:109660/1‑148 (MQ=255)
                                                                                                        gggCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCcgcg                                       >  3:82489/1‑148 (MQ=255)
                                                                                                         ggCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCgca                                     <  1:31014/149‑1 (MQ=255)
                                                                                                          gCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGgcag                                                    >  8:175616/1‑133 (MQ=255)
                                                                                                          gCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGgcag                                                    <  7:175616/133‑1 (MQ=255)
                                                                                                          gCGAGAGCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCgcag                                    <  6:75371/149‑1 (MQ=255)
                                                                                                            gagaTCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCgcagca                                  <  8:71305/149‑1 (MQ=255)
                                                                                                                  ttCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGt                                                                                           >  8:18699/1‑86 (MQ=255)
                                                                                                                  ttCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGt                                                                                           <  7:18699/86‑1 (MQ=255)
                                                                                                                    cAAGAACCTGTTCACCGCTAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATc                                                                                                                  >  5:229137/1‑61 (MQ=255)
                                                                                                                    cAAGAACCTGTTCACCGCTAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATc                                                                                                                  <  6:229137/61‑1 (MQ=255)
                                                                                                                      aGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTtcatc                        <  4:82489/149‑1 (MQ=255)
                                                                                                                        aaCCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTtcatca                       >  6:146011/1‑148 (MQ=255)
                                                                                                                            tGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCAcc                  >  3:120349/1‑149 (MQ=255)
                                                                                                                            cgagcACCGGCAACTTGGCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCAcc                  >  7:84624/5‑149 (MQ=255)
                                                                                                                             gagcACCGGCAACTTGGCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCgcgc                                      >  7:17722/4‑128 (MQ=255)
                                                                                                                             gagcACCGGCAACTTGGCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCgcgc                                      <  8:17722/125‑1 (MQ=255)
                                                                                                                              ttCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGt                                                         <  7:63919/108‑1 (MQ=255)
                                                                                                                              ttCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCAcc                  <  5:204350/147‑1 (MQ=255)
                                                                                                                              ttCACCCCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGt                                                         >  8:63919/1‑108 (MQ=255)
                                                                                                                                  ccGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAagag            >  6:288320/1‑149 (MQ=255)
                                                                                                                                   cGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAAccc                                              >  6:251293/1‑114 (MQ=255)
                                                                                                                                   cGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAAccc                                              <  3:65710/114‑1 (MQ=255)
                                                                                                                                   cGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAAccc                                              >  4:65710/1‑114 (MQ=255)
                                                                                                                                   cGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAAccc                                              <  5:251293/114‑1 (MQ=255)
                                                                                                                                    gCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGc                                                               >  2:185817/1‑96 (MQ=255)
                                                                                                                                       aaGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGgc       <  7:109660/149‑1 (MQ=255)
                                                                                                                                         gTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGgcgc     <  4:120349/149‑1 (MQ=255)
                                                                                                                                          ttGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGGCGcc    >  7:45706/1‑149 (MQ=255)
                                                                                                                                            gCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGCGCGTGGCAGCAACCCGGCCGCCCAGCACTCGTGCATCATCACCGAAGAGCTGGCGCCGa  >  4:204653/1‑149 (MQ=255)
                                                                                                                                            gCCGGAGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGGCGCCGa  >  4:37902/1‑149 (MQ=255)
                                                                                                                                             ccGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGTCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGGCGCCGa  <  7:128985/148‑1 (MQ=255)
                                                                                                                                             ccGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGGCGCCGa  <  7:133287/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGGCGCCGA  >  NC_002947/417481‑417769

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