Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 654,398 C→A 10.9% L115L (CTG→CTT)  ctpL ← methyl‑accepting chemotaxis protein CtpL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947654,3980CA10.9% 165.5 / 10.0 55L115L (CTG→CTT) ctpLmethyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base A (3/3);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGG  >  NC_002947/654263‑654545
                                                                                                                                       |                                                                                                                                                   
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)
              aGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTg                                                                                                                          <  7:120330/149‑1 (MQ=255)
               gggCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTgg                                                                                                                         <  1:260961/149‑1 (MQ=255)
                   agcaAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGt                                                                                                                     <  8:177339/149‑1 (MQ=255)
                        gTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAg                                                                                                                 <  6:22105/148‑1 (MQ=255)
                        gTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAg                                                                                                                 >  5:22105/1‑148 (MQ=255)
                           gCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCt                                                                                                             <  3:93864/149‑1 (MQ=255)
                             gATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGc                                                                                                           <  3:119155/149‑1 (MQ=255)
                                   ccTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTa                                                                                                     >  1:139979/1‑149 (MQ=255)
                                      cGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATcc                                                                                                  <  3:47262/149‑1 (MQ=255)
                                          gTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGc                                                                                              >  6:195607/1‑149 (MQ=255)
                                              ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgg                                                                                           <  5:195607/148‑1 (MQ=255)
                                               cACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                         >  8:110672/1‑149 (MQ=255)
                                                aCTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgc                                                                                        <  2:96838/149‑1 (MQ=255)
                                                    tCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCaa                                                                                    >  5:147853/1‑149 (MQ=255)
                                                      gcgGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACt                                                                                  >  5:204120/1‑149 (MQ=255)
                                                          gCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTtctc                                                                              >  2:212283/1‑149 (MQ=255)
                                                             tAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGc                                                                           <  2:182167/149‑1 (MQ=255)
                                                             tAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGc                                                                           <  3:316158/149‑1 (MQ=255)
                                                                     cATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcag                                                                                                                                    >  5:13142/1‑84 (MQ=255)
                                                                     cATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcag                                                                                                                                    <  6:13142/84‑1 (MQ=255)
                                                                           cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGGCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCg                                                                <  3:135578/146‑1 (MQ=255)
                                                                            tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                            <  4:346729/101‑1 (MQ=255)
                                                                            tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                            >  3:346729/1‑101 (MQ=255)
                                                                              gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgc                                                                                        >  3:93429/1‑119 (MQ=255)
                                                                              gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgc                                                                                        <  4:93429/119‑1 (MQ=255)
                                                                                aGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTcc                                                        >  7:232097/1‑149 (MQ=255)
                                                                                      cccccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGt                                                  >  2:100281/1‑149 (MQ=255)
                                                                                       ccccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTg                                                 >  4:47182/1‑149 (MQ=255)
                                                                                         ccAGCTCCCGCTCGGCTTGCAGCAGGAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCGTGCAGGCTATCCAGGGTGGGGCGCAACTTCTCGGCCAGCTCGACCGGCCAGTCCGGCAGTGCg                                               >  8:244250/1‑149 (MQ=255)
                                                                                           aGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGc                                             <  7:244250/149‑1 (MQ=255)
                                                                                              tCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTg                                          >  5:18717/1‑149 (MQ=255)
                                                                                                           gcaccagaaGGGCTTGTGGGTCGCCCGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTTCAGGATATCCAGGCTGGGGCGCAACTTCTATGCCAGCTATACAGGCCAGTCCGTCAGTGCGGACTGCAACAGTTGGTTc                             >  8:150663/1‑149 (MQ=255)
                                                                                                            cagcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTc                             <  2:139979/148‑1 (MQ=255)
                                                                                                            cagcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTc                             <  8:232097/148‑1 (MQ=255)
                                                                                                            cagcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCg                            >  6:198724/1‑149 (MQ=255)
                                                                                                             agcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGc                           >  4:225745/1‑149 (MQ=255)
                                                                                                                     ggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                              <  6:343421/106‑1 (MQ=255)
                                                                                                                     ggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                              >  5:343421/1‑106 (MQ=255)
                                                                                                                          gctggccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggcg              >  7:197311/7‑149 (MQ=255)
                                                                                                                           ctggccGCCGGCAAGCTGGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggcgg             >  4:117659/6‑149 (MQ=255)
                                                                                                                            gggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                     >  7:207457/1‑140 (MQ=255)
                                                                                                                            gggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                     <  8:207457/140‑1 (MQ=255)
                                                                                                                             ggccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGcc           >  6:100467/4‑149 (MQ=255)
                                                                                                                              gccGCCGGCAAGCTGGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCt          <  3:117659/147‑1 (MQ=255)
                                                                                                                               tCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTg         >  7:179963/1‑149 (MQ=255)
                                                                                                                               tCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTg         <  6:18717/149‑1 (MQ=255)
                                                                                                                                     gCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAg   >  8:259402/1‑149 (MQ=255)
                                                                                                                                      ccAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg  <  1:212283/149‑1 (MQ=255)
                                                                                                                                      ccAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg  <  6:204120/149‑1 (MQ=255)
                                                                                                                                      ccAGCAGGCAGGCGGCCAGCATTTCGTTCGAGGTGAAGGCCTGCAGGCTATCCATGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGTCCTTCAACTGTTGGTTCTCATCCATGGCGGCGGCCTTGCGCAgg  <  7:150663/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                   
GGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGG  >  NC_002947/654263‑654545

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