Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,273,500 G→T 7.4% A126D (GCT→GAT)  glcC ← transcriptional regulator GlcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,273,5000GT 7.4% 215.2 / 5.4 67A126D (GCT→GAT) glcCtranscriptional regulator GlcC
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base T (3/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

ATGCGACGCCTCACTGATGGCCCGGTGAAAGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCGGTGAGATCGGCGACGAACGAACCCC  >  NC_002947/4273361‑4273647
                                                                                                                                           |                                                                                                                                                   
aTGCGACGCCTCACTGATGGCCCGGGGAAAGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAAc                                                                                                                                             <  1:360399/148‑1 (MQ=255)
   cGACGCCTCACTGATGGCCCGGTGAAAGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAAcgccg                                                                                                                                         >  5:403936/1‑149 (MQ=255)
                aTGGCCCGGTGAAAGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGa                                                                                                                            >  1:156599/1‑149 (MQ=255)
                     ccGGTGAAAGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAg                                                                                                                       >  3:26370/1‑149 (MQ=255)
                      cGGTGAAAGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGc                                                                                                                      <  1:14361/149‑1 (MQ=255)
                         tGAAAGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTc                                                                                                                   >  7:72014/1‑149 (MQ=255)
                           aaaGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCa                                                                                                                 <  7:110473/149‑1 (MQ=255)
                            aaGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAg                                                                                                                <  8:72014/149‑1 (MQ=255)
                               gCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCaa                                                                                                             >  2:333639/1‑149 (MQ=255)
                                  tGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCCGCCTCCATCAACGc                                                                                                          >  5:252125/1‑149 (MQ=255)
                                     gcACGGCGGGCATGCTCACGGCCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCg                                                                                                       >  8:352482/2‑149 (MQ=255)
                                       aCGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGa                                                                                                     >  8:119839/1‑149 (MQ=255)
                                           cggGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTc                                                                                                 <  3:232283/149‑1 (MQ=255)
                                                 tGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCaagcag                                                                                           >  7:304647/1‑149 (MQ=255)
                                                 tGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTAAGCCTCCAGCAACGCCCGGACCTCcagcag                                                                                           <  1:246019/149‑1 (MQ=255)
                                                   ctgaCGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGt                                                                                         <  7:352482/146‑1 (MQ=255)
                                                    tCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTc                                                                                        >  4:373581/1‑149 (MQ=255)
                                                        gggCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAAc                                                                                    <  4:26370/149‑1 (MQ=255)
                                                         ggCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACa                                                                                   <  2:156599/149‑1 (MQ=255)
                                                         ggCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACa                                                                                   <  1:333639/149‑1 (MQ=255)
                                                            aTCGAGCGCCTCCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCTCGGACCTCCAGCAGGTCGAACagag                                                                                >  4:83955/1‑149 (MQ=255)
                                                               gAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAAc                                                                                    >  6:224961/1‑142 (MQ=255)
                                                               gAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAAc                                                                                    <  5:224961/142‑1 (MQ=255)
                                                                   gcCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTTGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGTAAGGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCgg                                                                         <  7:61875/149‑1 (MQ=255)
                                                                   gcCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTATCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCcagca                                                                                            >  7:229912/1‑130 (MQ=255)
                                                                   gcCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTATCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCcagca                                                                                            >  7:223826/1‑130 (MQ=255)
                                                                   gcCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTATCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCcagca                                                                                            <  8:223826/130‑1 (MQ=255)
                                                                   gcCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTATCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCcagca                                                                                            <  8:229912/130‑1 (MQ=255)
                                                                        cgcAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGc                                                                     <  8:126707/148‑1 (MQ=255)
                                                                               tCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGaa                                                             >  7:367439/1‑149 (MQ=255)
                                                                               tCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGaa                                                             >  1:157743/1‑149 (MQ=255)
                                                                               tCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGaa                                                             >  4:324341/1‑149 (MQ=255)
                                                                                cGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGcc                                                                                                         <  3:390253/104‑1 (MQ=255)
                                                                                cGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGcc                                                                                                         >  4:390253/1‑104 (MQ=255)
                                                                                  tGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAg                                                          <  3:324341/149‑1 (MQ=255)
                                                                                    cgcAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGt                                                        <  3:373581/149‑1 (MQ=255)
                                                                                      cagccagcATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTg                                                       >  4:179072/1‑148 (MQ=255)
                                                                                       agccagcATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCa                                                     <  3:260015/149‑1 (MQ=255)
                                                                                        gccagcATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCAt                                                    >  2:25981/1‑149 (MQ=255)
                                                                                         ccagcATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCAt                                                    >  6:53802/1‑148 (MQ=255)
                                                                                           agcATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCaa                                                 <  8:224402/149‑1 (MQ=255)
                                                                                           agcATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCaa                                                 <  5:53802/149‑1 (MQ=255)
                                                                                           agcATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCa                                                  >  7:224402/1‑148 (MQ=255)
                                                                                            gcATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCaaa                                                <  3:129128/149‑1 (MQ=255)
                                                                                              aTTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAAc                                                            <  6:387244/135‑1 (MQ=255)
                                                                                              aTTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAAc                                                            >  5:387244/1‑135 (MQ=255)
                                                                                              aTTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAgg                                              >  1:89608/1‑149 (MQ=255)
                                                                                                  cttcGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCg                                            >  1:173107/1‑147 (MQ=255)
                                                                                                   ttcGTAGCGCCGGCGGATTAGCAGGAGGTCGACGACGGTATCCCGCAACGCCGCCAGACGCGCAGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGAGAACTGAACAGGTGCATCAAAGGCGAAg                                         >  5:108734/1‑149 (MQ=255)
                                                                                                     cGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCaa                                                                                                                                              >  5:137605/1‑46 (MQ=255)
                                                                                                     cGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCaa                                                                                                                                              <  6:137605/46‑1 (MQ=255)
                                                                                                     cGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTa                                          <  6:386055/146‑1 (MQ=255)
                                                                                                        tgggcCGGCGGATTAGCTGGCGGTCGACGTCGGTAGCCCGCAACGGCGCGAGCCGCGCCGACTCATCCTCGAGCTACGCCGGGACCTCCAGCAGGTCGAACAGAGTTCGAGGTTGCGATGTGAACAGGTGCATCAAAGGCGTATCAc                                      <  6:108734/144‑1 (MQ=255)
                                                                                                          cgcCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGc                                  <  3:179072/149‑1 (MQ=255)
                                                                                                             cggcggATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCg                                                                    >  1:248462/1‑112 (MQ=255)
                                                                                                             cggcggATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCg                                                                    <  2:248462/112‑1 (MQ=255)
                                                                                                               gcggATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGc                                                                                                               <  6:160849/67‑1 (MQ=255)
                                                                                                               gcggATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGcc                             >  5:43243/1‑149 (MQ=255)
                                                                                                                 ggATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCg                           >  7:235984/1‑149 (MQ=255)
                                                                                                                   aTTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCGGt                         <  4:25896/149‑1 (MQ=255)
                                                                                                                   aTTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCGGt                         <  6:43243/149‑1 (MQ=255)
                                                                                                                   aTTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCGGt                         <  8:320281/149‑1 (MQ=255)
                                                                                                                    ttAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCGGTg                        <  7:119839/149‑1 (MQ=255)
                                                                                                                    ttAGCACGCGTTCGACGGCGCGAGCCCGCACCGCCGCCCGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCCGCCGGTCGTACAGAGGGCGCGGTTGCGAACTCAACAGGTGCATCAAACGCGTAGCACCCGCCTGCCCGGTg                        <  8:304647/149‑1 (MQ=255)
                                                                                                                       gCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAAGAGAGTGCGCGGTTGCGATCTGAACAGGt                                                        <  3:251576/114‑1 (MQ=255)
                                                                                                                                cGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCGGTGAGATCGGcgac             >  3:187048/1‑148 (MQ=255)
                                                                                                                                         tAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCGGTGAGATCGGCGACGAACGAAcc    <  2:157743/148‑1 (MQ=255)
                                                                                                                                           gCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCGGTGAGATCGGCGACGAACGAAcccc  <  6:403936/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
ATGCGACGCCTCACTGATGGCCCGGTGAAAGGCATGGTCACGGCGGGCATGCTCACGGGCATCGAGCGCCTGCGCAGCCTCGTGCGCAGCCAGCATTTCTTCGTAGCGCCGGCGGATTAGCAGGCGGTCGACGTCGGTAGCCCGCAACGCCGCCAGCCGCGCCGACTCAGCCTCCAGCAACGCCCGGACCTCCAGCAGGTCGAACAGAGTGCGCGGTTGCGAACTGAACAGGTGCATCAAAGGCGTAGCACCCGCCTGCCCGGTGAGATCGGCGACGAACGAACCCC  >  NC_002947/4273361‑4273647

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