Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 654,403 T→G 5.8% 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.8% 208.8 / 2.4 68K114Q (AAG→CAG) ctpLmethyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base G (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

AGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCG  >  NC_002947/654261‑654548
                                                                                                                                               |                                                                                                                                                 
aGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCcggc                                                                                                                                             <  7:262249/149‑1 (MQ=255)
 gggTGGCGTCTAGCAAGGGCAGCAAGTAGCGAATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGAAGCAGCAGGGCTTGTGGGTCGACGGCCAGCTTGCcggc                                                                                                                                             >  1:40962/1‑149 (MQ=255)
  ggTGGCGTCTAGCAAGGGCAGCAAGTAGCG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGCCGCCGGCAAGCTGGCcggcga                                                                                                                                           <  8:300973/149‑2 (MQ=255)
      gCGTCTAGCAAGGGCAGCAAGTAGCG‑ATTGGCCTCGGCGTTGCCACTTTCGCGGGCGTCGCTGGCCATCTGATCGAAGTTAGCCCCCAGCTCCCGCTCGGCTTGCAGCAGAAGGGCTTGTGGGTCGCCGGCCAGATTTCCGGCGGCcag                                                                                                                                       >  6:303995/1‑149 (MQ=255)
         tCTAGCAAGGGCAGCAAGTAGCG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcagcag                                                                                                                                    <  5:51306/149‑1 (MQ=255)
             gCAAGGGCAGCAAGTAGCG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCg                                                                                                                                >  5:106651/1‑149 (MQ=255)
                   gcagcaAGTAGCG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCg                                                                                                                          <  8:14048/149‑1 (MQ=255)
                   gcagcaAGTAGCG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCg                                                                                                                          <  3:102095/149‑1 (MQ=255)
                           tAGCG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGc                                                                                                                  <  2:69425/149‑1 (MQ=255)
                              cG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                               <  4:56895/149‑1 (MQ=255)
                               g‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                               >  3:119929/1‑148 (MQ=255)
                                 aTTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCa                                                                                                             <  6:180248/149‑1 (MQ=255)
                                 aTTGGCCTCGGCGTGGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCa                                                                                                             <  8:6084/149‑1 (MQ=255)
                                    ggCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAgg                                                                                                           <  3:179895/148‑1 (MQ=255)
                                       cTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTAt                                                                                                       >  1:267405/1‑149 (MQ=255)
                                        tCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGGCGCGGGCCAGCTTGCGGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATc                                                                                                      <  2:380483/149‑1 (MQ=255)
                                             gTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATcc                                                                                                     <  3:100105/145‑1 (MQ=255)
                                             gTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATcc                                                                                                     >  4:100105/1‑145 (MQ=255)
                                                gCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                               <  7:237395/132‑1 (MQ=255)
                                                gCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                               >  8:237395/1‑132 (MQ=255)
                                                 ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTggg                                                                                             <  7:333521/149‑1 (MQ=255)
                                                 ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTggg                                                                                             >  5:302283/1‑149 (MQ=255)
                                                 ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTggg                                                                                             <  5:263303/149‑1 (MQ=255)
                                                  cACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                            >  2:172149/1‑149 (MQ=255)
                                                  cACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                            >  1:143018/1‑149 (MQ=255)
                                                       tCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCCGCTTGCAGCAGCAGTGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGATGTCGGTGAAGGCCTGCAGCCTATCCAGGCTGGGGCGCaa                                                                                       >  3:252706/1‑149 (MQ=255)
                                                          cgGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAg                                                                                                                    >  4:188150/1‑117 (MQ=255)
                                                          cgGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAg                                                                                                                    <  3:188150/117‑1 (MQ=255)
                                                             gCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCGTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAActcctc                                                                                 <  5:304116/149‑1 (MQ=255)
                                                              cGTAGCTGGCCATCTGCTCGAAGTTCGCGCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCg                                                                                <  3:156203/149‑1 (MQ=255)
                                                              cGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCg                                                                                <  7:153481/149‑1 (MQ=255)
                                                                tAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGc                                                                              <  1:172149/149‑1 (MQ=255)
                                                                              cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTAcc                                                                    <  6:230862/145‑1 (MQ=255)
                                                                              cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                                 <  6:106651/148‑1 (MQ=255)
                                                                                cGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                                 >  2:90330/1‑146 (MQ=255)
                                                                                cGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCa                                                               <  1:90330/148‑1 (MQ=255)
                                                                                 gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCa                                                               <  5:90407/147‑1 (MQ=255)
                                                                                 gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGt                                                             >  2:69381/1‑149 (MQ=255)
                                                                                 gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGGTCGGTGGCGGTGAAGGGCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGGTCTACCGGCCAGt                                                             >  3:39332/1‑149 (MQ=255)
                                                                                         cccccAGCTCCCGCTCGGCTTGCAGCAACAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTAc                                                                     >  4:292213/1‑133 (MQ=255)
                                                                                         cccccAGCTCCCGCTCGGCTTGCAGCAACAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTAc                                                                     <  3:292213/133‑1 (MQ=255)
                                                                                            ccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGGCGCCTGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGCCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCg                                                  <  1:285444/149‑1 (MQ=255)
                                                                                             cAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCgg                                                 >  3:291459/1‑149 (MQ=255)
                                                                                              aGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGc                                                >  5:141352/1‑149 (MQ=255)
                                                                                                   gcGCTCGGCTTGCAGCAGCATGGCTTTTGGGTCGCGGGCCAGCTTGCCGGCGTGCAGCAGTTCTTTGGGTGTGAAGGCCTGCATGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCTGCCAGTCCGTCAGTGCGGCCttc                                            <  4:252706/147‑3 (MQ=255)
                                                                                                   ccGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCTGCCAGTCCGTCAGTGCGGCCTGc                                            >  5:97245/1‑148 (MQ=255)
                                                                                                   ccGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGc                                            >  7:334883/1‑148 (MQ=255)
                                                                                                   ccGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGc                                            <  1:76399/148‑1 (MQ=255)
                                                                                                     gCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAAc                                         >  7:51724/1‑149 (MQ=255)
                                                                                                         ggCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTTGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGtt                                     >  1:305755/1‑149 (MQ=255)
                                                                                                           cTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTgg                                   >  8:152207/1‑149 (MQ=255)
                                                                                                               cagcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTc                                <  3:174414/148‑1 (MQ=255)
                                                                                                                agcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGc                              >  6:224794/1‑149 (MQ=255)
                                                                                                                   agcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTc                           <  1:69381/149‑1 (MQ=255)
                                                                                                                            tgtgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGGTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggc                  <  6:302283/149‑1 (MQ=255)
                                                                                                                              tgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCTGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCTGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAAATGTTGGTTCGCCTCCATggcggcg                 <  2:40962/148‑1 (MQ=255)
                                                                                                                              ctggccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggcgg                >  6:174132/6‑149 (MQ=255)
                                                                                                                                ggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGGCTGCAACTGTTGGTTCGCCTCCATggcggcggc               <  2:143018/148‑1 (MQ=255)
                                                                                                                                 gccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCt             >  8:96274/3‑149 (MQ=255)
                                                                                                                                  tCGCCGGCCAGCTTGCCGGCGGGCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCGAGGCTGGGGCGAAACTTCTCGGCCAGCTCTCCCGGCCAGTCCGTCAGTGCGGCCTGCACCTGTTGGTTCGCCTCCATGGCGGCGGCCTg            >  7:59310/1‑149 (MQ=255)
                                                                                                                                  tCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTg            >  5:314486/1‑149 (MQ=255)
                                                                                                                                    gCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGgc          >  8:44916/1‑149 (MQ=255)
                                                                                                                                      cGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGgcg         >  3:137538/1‑148 (MQ=255)
                                                                                                                                        gCCAGCTTCCCGGCGGCCAGCAGGTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCATCTTCTCGGCCAGCTCTCCCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAg      <  4:39332/149‑1 (MQ=255)
                                                                                                                                           aGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGc   <  6:141352/149‑1 (MQ=255)
                                                                                                                                           aGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGc   <  1:181895/149‑1 (MQ=255)
                                                                                                                                           aGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGc   <  7:96274/149‑1 (MQ=255)
                                                                                                                                            gCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCg  >  8:78392/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
AGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCG‑ATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCG  >  NC_002947/654261‑654548

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