Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 654,403 T→G 10.1% K114Q (AAG→CAG)  ctpL ← methyl‑accepting chemotaxis protein CtpL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947654,4030TG10.1% 174.8 / 8.0 58K114Q (AAG→CAG) ctpLmethyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base G (3/3);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCG  >  NC_002947/654263‑654548
                                                                                                                                            |                                                                                                                                                  
ggTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCGTCGGCGTTGCCACTGTCTCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCGTGCAGCAGCAGGGCTTGTGGGTCGCCGGCAAGCTggcgggcga                                                                                                                                            <  8:163407/149‑2 (MQ=255)
ggTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCcggcgg                                                                                                                                            >  2:260961/1‑149 (MQ=255)
 gTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGCCGCCGGCAAGCTGGCcggcgac                                                                                                                                           <  3:183975/149‑3 (MQ=255)
      gTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcagca                                                                                                                                      <  1:161930/149‑1 (MQ=255)
      gTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcagca                                                                                                                                      <  7:240656/149‑1 (MQ=255)
      gTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcag                                                                                                                                        >  4:119155/1‑147 (MQ=255)
       tCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcagcag                                                                                                                                     >  4:93864/1‑149 (MQ=255)
          aGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTc                                                                                                                                  >  8:120330/1‑149 (MQ=255)
              aGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTg                                                                                                                             <  7:120330/149‑1 (MQ=255)
               gggCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTgg                                                                                                                            <  1:260961/149‑1 (MQ=255)
                   agcaAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGt                                                                                                                        <  8:177339/149‑1 (MQ=255)
                        gTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAg                                                                                                                    <  6:22105/148‑1 (MQ=255)
                        gTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAg                                                                                                                    >  5:22105/1‑148 (MQ=255)
                           gCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCt                                                                                                                <  3:93864/149‑1 (MQ=255)
                             gATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGc                                                                                                              <  3:119155/149‑1 (MQ=255)
                                   ccTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTa                                                                                                        >  1:139979/1‑149 (MQ=255)
                                      cGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATcc                                                                                                     <  3:47262/149‑1 (MQ=255)
                                          gTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGc                                                                                                 >  6:195607/1‑149 (MQ=255)
                                              ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgg                                                                                              <  5:195607/148‑1 (MQ=255)
                                               cACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                            >  8:110672/1‑149 (MQ=255)
                                                aCTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgc                                                                                           <  2:96838/149‑1 (MQ=255)
                                                    tCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCaa                                                                                       >  5:147853/1‑149 (MQ=255)
                                                      gcgGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACt                                                                                     >  5:204120/1‑149 (MQ=255)
                                                          gCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTtctc                                                                                 >  2:212283/1‑149 (MQ=255)
                                                             tAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGc                                                                              <  2:182167/149‑1 (MQ=255)
                                                             tAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGc                                                                              <  3:316158/149‑1 (MQ=255)
                                                                     cATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcag                                                                                                                                        >  5:13142/1‑84 (MQ=255)
                                                                     cATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcag                                                                                                                                        <  6:13142/84‑1 (MQ=255)
                                                                           cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGGCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCg                                                                   <  3:135578/146‑1 (MQ=255)
                                                                            tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTg                                                                                                               <  4:346729/101‑1 (MQ=255)
                                                                            tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTg                                                                                                               >  3:346729/1‑101 (MQ=255)
                                                                              gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgc                                                                                           >  3:93429/1‑119 (MQ=255)
                                                                              gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgc                                                                                           <  4:93429/119‑1 (MQ=255)
                                                                                aGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTcc                                                           >  7:232097/1‑149 (MQ=255)
                                                                                      cccccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGt                                                     >  2:100281/1‑149 (MQ=255)
                                                                                       ccccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTg                                                    >  4:47182/1‑149 (MQ=255)
                                                                                         ccAGCTCCCGCTCGGCTTGCAGCAGGAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCGTGCAGGCTATCCAGGGTGGGGCGCAACTTCTCGGCCAGCTCGACCGGCCAGTCCGGCAGTGCg                                                  >  8:244250/1‑149 (MQ=255)
                                                                                           aGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGc                                                <  7:244250/149‑1 (MQ=255)
                                                                                              tCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTg                                             >  5:18717/1‑149 (MQ=255)
                                                                                                           gcaccagaaGGGCTTGTGGGTCGCCCGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTTCAGGATATCCAGGCTGGGGCGCAACTTCTATGCCAGCTATACAGGCCAGTCCGTCAGTGCGGACTGCAACAGTTGGTTc                                >  8:150663/1‑149 (MQ=255)
                                                                                                            cagcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTc                                <  2:139979/148‑1 (MQ=255)
                                                                                                            cagcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTc                                <  8:232097/148‑1 (MQ=255)
                                                                                                            cagcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCg                               >  6:198724/1‑149 (MQ=255)
                                                                                                             agcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGc                              >  4:225745/1‑149 (MQ=255)
                                                                                                                     ggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                                 >  5:343421/1‑106 (MQ=255)
                                                                                                                     ggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                                 <  6:343421/106‑1 (MQ=255)
                                                                                                                          gctggccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggcg                 >  7:197311/7‑149 (MQ=255)
                                                                                                                           ctggccGCCGGCAAGCTGGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggcgg                >  4:117659/6‑149 (MQ=255)
                                                                                                                            gggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                        >  7:207457/1‑140 (MQ=255)
                                                                                                                            gggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                        <  8:207457/140‑1 (MQ=255)
                                                                                                                             ggccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGcc              >  6:100467/4‑149 (MQ=255)
                                                                                                                              gccGCCGGCAAGCTGGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCt             <  3:117659/147‑1 (MQ=255)
                                                                                                                               tCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTg            >  7:179963/1‑149 (MQ=255)
                                                                                                                               tCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTg            <  6:18717/149‑1 (MQ=255)
                                                                                                                                     gCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAg      >  8:259402/1‑149 (MQ=255)
                                                                                                                                      ccAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg     <  1:212283/149‑1 (MQ=255)
                                                                                                                                      ccAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg     <  6:204120/149‑1 (MQ=255)
                                                                                                                                      ccAGCAGGCAGGCGGCCAGCATTTCG‑TTCGAGGTGAAGGCCTGCAGGCTATCCATGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGTCCTTCAACTGTTGGTTCTCATCCATGGCGGCGGCCTTGCGCAgg     <  7:150663/149‑1 (MQ=255)
                                                                                                                                          cTTGCCGGCGGCCAGCAGTTCGTTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCg  >  7:52409/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCG‑TTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCG  >  NC_002947/654263‑654548

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: ‑