Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 654,403 T→G 5.3% 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,4030TG 5.3% 244.8 / 4.2 75K114Q (AAG→CAG) ctpLmethyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base G (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

TGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGT  >  NC_002947/654258‑654549
                                                                                                                                                 |                                                                                                                                                  
tGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGcc                                                                                                                                                 <  1:337056/149‑1 (MQ=255)
   aGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCcggc                                                                                                                                              >  5:374368/1‑149 (MQ=255)
     ggTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCAAGCTGGCcggcga                                                                                                                                            <  6:192193/149‑2 (MQ=255)
         gCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCTGCGGTGCCACTGTCGCGGGCGTAGCTGGCCATCGGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCca                                                                                                                                         <  8:316590/148‑1 (MQ=255)
         gCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcag                                                                                                                                        <  7:3700/149‑1 (MQ=255)
         gCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCca                                                                                                                                         <  1:183447/148‑1 (MQ=255)
           gTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcag                                                                                                                                        >  1:167130/1‑147 (MQ=255)
              tAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGtt                                                                                                                                   >  2:121811/1‑149 (MQ=255)
                gCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCg                                                                                                                                 >  1:263671/1‑149 (MQ=255)
                gCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCg                                                                                                                                 >  3:203790/1‑149 (MQ=255)
                gCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCg                                                                                                                                 <  3:122650/149‑1 (MQ=255)
                gCAACGGCAGCAAGTATCGAGTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGt                                                                                                                                    <  6:82362/146‑1 (MQ=255)
                    gggCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTgg                                                                                                                             <  3:250100/149‑1 (MQ=255)
                     ggCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGc                                                                                                                            <  6:374368/149‑1 (MQ=255)
                        agcaAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGt                                                                                                                         <  4:226710/149‑1 (MQ=255)
                        agcaAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGt                                                                                                                         >  3:285135/1‑149 (MQ=255)
                                  gATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                                >  5:67513/1‑148 (MQ=255)
                                  gATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                                >  1:205872/1‑148 (MQ=255)
                                   aTTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                                >  6:250599/1‑147 (MQ=255)
                                   aTTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCa                                                                                                              <  2:287921/149‑1 (MQ=255)
                                      ggCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCGGGCGGCCAGCAGTTCGTTGGCGGTGAAGTCCTGCAgg                                                                                                            <  6:257053/148‑1 (MQ=255)
                                       gCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCAGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGATGGCGGTGAAGGCCTGCAGGc                                                                                                           >  7:268221/1‑148 (MQ=255)
                                            ggCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTg                                                                                                                        >  6:250631/1‑130 (MQ=255)
                                            ggCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTg                                                                                                                        <  5:250631/130‑1 (MQ=255)
                                            ggCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATcc                                                                                                      <  6:242476/148‑1 (MQ=255)
                                            ggCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATc                                                                                                       <  5:250599/147‑1 (MQ=255)
                                                  gCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTg                                                                                                <  8:274404/148‑1 (MQ=255)
                                                   ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTggg                                                                                              <  1:258065/149‑1 (MQ=255)
                                                   ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTggg                                                                                              <  7:132617/149‑1 (MQ=255)
                                                   ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgg                                                                                               <  2:263671/148‑1 (MQ=255)
                                                   ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgg                                                                                               <  6:168270/148‑1 (MQ=255)
                                                         tCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgcgc                                                                                          <  2:89501/147‑1 (MQ=255)
                                                         tCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGGCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgcgc                                                                                          <  2:84499/147‑1 (MQ=255)
                                                           gcgGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACt                                                                                      >  6:176928/1‑149 (MQ=255)
                                                              ggCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGACGCCGGCCAGCTTGCCGGCCCCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTtct                                                                                   >  4:34170/1‑149 (MQ=255)
                                                                cGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCg                                                                                 <  4:203790/149‑1 (MQ=255)
                                                                   aGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGc                                                                               >  2:38805/1‑148 (MQ=255)
                                                                   aGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCAAGCTGGCcggcga                                                                                                                                            <  3:117773/87‑2 (MQ=255)
                                                                   aGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCAAGCTGGCcggcga                                                                                                                                            >  4:117773/1‑86 (MQ=255)
                                                                    gCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGc                                                                               >  7:17063/1‑147 (MQ=255)
                                                                            tCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGc                                                                                                               >  2:14224/1‑107 (MQ=255)
                                                                            tCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGc                                                                                                               <  1:14224/107‑1 (MQ=255)
                                                                            tCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                             >  5:333264/1‑125 (MQ=255)
                                                                            tCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                             <  6:333264/125‑1 (MQ=255)
                                                                            tCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTAcc                                                                     >  5:29054/1‑149 (MQ=255)
                                                                               gCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                                  >  7:259105/1‑149 (MQ=255)
                                                                                ctagAAGTTCGCCCCCAGCTCCCGCTAGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                                  <  2:305580/145‑1 (MQ=255)
                                                                                cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                                  >  5:209631/1‑148 (MQ=255)
                                                                                       ttCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAgg                                                                                                            <  6:39416/99‑1 (MQ=255)
                                                                                       ttCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAgg                                                                                                            >  5:39416/1‑99 (MQ=255)
                                                                                       ttCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGt                                                          <  4:285135/149‑1 (MQ=255)
                                                                                         cGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTc                                                         <  6:67513/148‑1 (MQ=255)
                                                                                          gCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGACAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGATGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCa                                                        <  1:293095/148‑1 (MQ=255)
                                                                                               cAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCgg                                                  >  3:187042/1‑149 (MQ=255)
                                                                                                       gCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAAc                                          >  8:322373/1‑149 (MQ=255)
                                                                                                           ggCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCCGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGtt                                      >  6:123965/1‑149 (MQ=255)
                                                                                                              ttGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGt                                   >  4:18181/1‑149 (MQ=255)
                                                                                                                      gcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTcc                           <  6:209631/149‑1 (MQ=255)
                                                                                                                        agGGCTTGTGGGTCGCCGACCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCCATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                         >  5:142564/1‑149 (MQ=255)
                                                                                                                          ggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATgg                       >  4:26062/1‑149 (MQ=255)
                                                                                                                          ggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATgg                       <  4:289119/149‑1 (MQ=255)
                                                                                                                            cTTGTGGGTCGCCGACCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCCATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcg                     <  6:142564/149‑1 (MQ=255)
                                                                                                                               gtgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTCAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCCACTGTTGGTTCGCCTCAATggcggcg                  >  6:365931/1‑149 (MQ=255)
                                                                                                                                tgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                             <  6:61067/73‑1 (MQ=255)
                                                                                                                                tgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                             >  5:61067/1‑73 (MQ=255)
                                                                                                                                tgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCa                                                        >  2:29595/1‑110 (MQ=255)
                                                                                                                                tgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCa                                                        <  1:29595/110‑1 (MQ=255)
                                                                                                                                  ggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGcc               >  5:23089/1‑149 (MQ=255)
                                                                                                                                   gTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCt              >  4:81959/1‑149 (MQ=255)
                                                                                                                                     cGCCGGCAAGCTGGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTgg            >  8:181364/1‑149 (MQ=255)
                                                                                                                                          gCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAg       <  6:29054/149‑1 (MQ=255)
                                                                                                                                           ccAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg      <  5:176928/149‑1 (MQ=255)
                                                                                                                                           ccAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg      <  2:167130/149‑1 (MQ=255)
                                                                                                                                            cAGCTTGCCGGAGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAT‑‑‑GGCGGCCTGGCGCAGGGCGt  >  6:121931/1‑149 (MQ=255)
                                                                                                                                             aGCGTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGc    <  8:268221/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGT  >  NC_002947/654258‑654549

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: