Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,688,371 T→G 5.8% L163R (CTG→CGG)  PP_3253 → carboxylate‑amine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,688,3710TG 5.8% 219.2 / 4.5 69L163R (CTG→CGG) PP_3253carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base G (2/2);  total (32/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.16e-01

TACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAACCGGTGCACTGGCTGC  >  NC_002947/3688229‑3688518
                                                                                                                                              |                                                                                                                                                    
tACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCGGT‑AGGc                                                                                                                                               >  3:346666/1‑149 (MQ=255)
  ccGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTg                                                                                                                                             >  1:182716/1‑149 (MQ=255)
    gCCATGTGGCAAGGCGTAGCCTGCTCAATAGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGc                                                                                                                                           <  4:256432/149‑1 (MQ=255)
      cATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGTTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCa                                                                                                                                         <  5:348767/149‑1 (MQ=255)
             cAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTAcc                                                                                                                                  <  1:46469/149‑1 (MQ=255)
                ggCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCggg                                                                                                                               <  1:161873/149‑1 (MQ=255)
                ggCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCGGT‑AGGCTGGCCAGCCCACAggg                                                                                                                               <  4:346666/149‑1 (MQ=255)
                  cGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTa                                                                                                                             <  4:337454/149‑1 (MQ=255)
                    tAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACt                                                                                                                           >  5:135887/1‑149 (MQ=255)
                      gCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACtt                                                                                                                          <  3:58877/148‑1 (MQ=255)
                       ccTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCtgt‑g                                                                                                                                                  <  8:304207/123‑1 (MQ=255)
                       ccTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCtgt‑g                                                                                                                                                  <  8:160057/123‑1 (MQ=255)
                       ccTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCtgt‑g                                                                                                                                                  >  7:304207/1‑123 (MQ=255)
                       ccTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCtgt‑g                                                                                                                                                  >  7:160057/1‑123 (MQ=255)
                       ccTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGa                                                                                                                        <  4:30653/149‑1 (MQ=255)
                           cTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCtgt‑g                                                                                                                                                  <  3:78625/119‑1 (MQ=255)
                           cTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCtgt‑g                                                                                                                                                  >  4:78625/1‑119 (MQ=255)
                              aaTGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCg                                                                                                                 >  7:69747/1‑149 (MQ=255)
                              aaTGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCg                                                                                                                 >  7:61355/1‑149 (MQ=255)
                                    cTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCGGT‑AGGCTGGCCAGCCCACAggggc                                                                                                                             <  5:209000/131‑3 (MQ=255)
                                    cTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCGGT‑AGGCTGGCCAGCCCACAggggc                                                                                                                             >  6:209000/1‑129 (MQ=255)
                                       cATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGa                                                                                                         <  2:182716/148‑1 (MQ=255)
                                       cATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGAt                                                                                                        <  5:345908/149‑1 (MQ=255)
                                           tgCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGTGGGGGAGGCCAGCCGACCGGGTACTTGGGGTATAGCCGAGTGATCTg                                                                                                     >  1:307878/1‑149 (MQ=255)
                                                    gCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTg                                                                                            <  1:18533/148‑1 (MQ=255)
                                                    gCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTg                                                                                            <  4:295001/148‑1 (MQ=255)
                                                     cGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGc                                                                                          <  8:61355/149‑1 (MQ=255)
                                                     cGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGc                                                                                          <  8:69747/149‑1 (MQ=255)
                                                      gtaccgccAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGcc                                                                                         <  7:11383/146‑1 (MQ=255)
                                                          cgccAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCcaca                                                                                      >  4:107158/1‑148 (MQ=255)
                                                              aGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATgg                                                                                 >  6:315810/1‑149 (MQ=255)
                                                              aGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATgg                                                                                 <  8:6554/149‑1 (MQ=255)
                                                                 ttGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCt                                                                              >  8:219823/1‑149 (MQ=255)
                                                                 ttGTGACCGAATGCAACTCATCAACCGACTTCTCAACTGGTTACAGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGACTACCGGGTAATTGAGATATCGCCGAGAGATATGCTGTGAGTTGCCAAACATCGGCt                                                                              >  4:137164/1‑149 (MQ=255)
                                                                      aCCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGccc                                                                         >  8:40679/1‑149 (MQ=255)
                                                                        cGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGa                                                                       <  8:249928/149‑1 (MQ=255)
                                                                              cAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACAt                                                                                   <  8:314348/131‑1 (MQ=255)
                                                                              cAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACAt                                                                                   >  7:314348/1‑131 (MQ=255)
                                                                              cAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCAtt                                                                 >  8:319484/1‑149 (MQ=255)
                                                                               aaCTGATCAACCGCCTTTTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGATTGGCCACACATgg                                                                                 >  7:53349/1‑132 (MQ=255)
                                                                               aaCTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTg                                                                >  5:342045/1‑149 (MQ=255)
                                                                                  tGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCGACACATGGGCTTGCCCGAAGCATTGccc                                                             >  2:182916/1‑149 (MQ=255)
                                                                                      cAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACt                                                         >  6:369569/1‑149 (MQ=255)
                                                                                        aCCGCCTTCTGGCCTGGTTACCGTTGTTACTGGCTTTGAGCACTTCGTCGCCCCTGT‑GTGCTGGCCAGCCTACCTGGTACTTGAGCTATCGCCTAGTGATCTGCGGTGATGGGCCACACATGGTCTTGCCCGAATCATTGCCCGACTgg                                                       <  3:137164/149‑1 (MQ=255)
                                                                                          cGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCAtt                                                                 >  8:295393/1‑137 (MQ=255)
                                                                                          cGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCAtt                                                                 <  7:295393/137‑1 (MQ=255)
                                                                                            ccTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGt                                                                                             <  1:73537/107‑1 (MQ=255)
                                                                                            ccTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGt                                                                                             >  2:73537/1‑107 (MQ=255)
                                                                                            ccTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGgccg                                                   >  7:132704/1‑149 (MQ=255)
                                                                                                    ccttggtACCGTTGTTACTGGCGTTGAGCACTTGGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCTAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGa                                            <  8:265343/142‑1 (MQ=255)
                                                                                                    ccTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAg                                           <  7:319484/149‑1 (MQ=255)
                                                                                                     cTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAgc                                          <  7:40679/149‑1 (MQ=255)
                                                                                                       ggTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAgcgc                                        >  4:287674/1‑149 (MQ=255)
                                                                                                         gtACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTa                                      <  4:261477/148‑1 (MQ=255)
                                                                                                               ttgttACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGc                                <  3:107158/149‑1 (MQ=255)
                                                                                                                  ttACTGGCGTTGAGCACTTCGTCGCCCCTGT‑CGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGccc                              <  6:313571/148‑1 (MQ=255)
                                                                                                                      tGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTg                         >  4:179996/1‑149 (MQ=255)
                                                                                                                        gCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCa                       <  5:315810/149‑1 (MQ=255)
                                                                                                                          gTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGc                     <  3:179996/149‑1 (MQ=255)
                                                                                                                              aGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAAc                 >  4:228194/1‑149 (MQ=255)
                                                                                                                              aGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAAc                 <  3:287674/149‑1 (MQ=255)
                                                                                                                                cACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAACCg               >  2:208866/1‑149 (MQ=255)
                                                                                                                                   ttCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAg                      <  2:62713/139‑1 (MQ=255)
                                                                                                                                   ttCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAg                      >  1:62713/1‑139 (MQ=255)
                                                                                                                                   ttCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAACCGGTg            >  5:163962/1‑149 (MQ=255)
                                                                                                                                     cgtcgCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAACCGGTGc           <  6:135887/148‑1 (MQ=255)
                                                                                                                                      gtcgCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGa                                            >  6:75921/1‑114 (MQ=255)
                                                                                                                                      gtcgCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGa                                            <  5:75921/114‑1 (MQ=255)
                                                                                                                                        cgCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAACCGGTGCACTg       >  3:273191/1‑149 (MQ=255)
                                                                                                                                         gCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGc                                <  1:1542/123‑1 (MQ=255)
                                                                                                                                         gCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGc                                >  2:1542/1‑123 (MQ=255)
                                                                                                                                           cccTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAACCGGTGCACTGGCt    >  2:73087/1‑149 (MQ=255)
                                                                                                                                            ccTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGt                                  >  6:90633/1‑118 (MQ=255)
                                                                                                                                            ccTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGt                                  <  5:90633/118‑1 (MQ=255)
                                                                                                                                             cTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAACCGGTGCACTGGCTGc  <  3:228194/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTTACCGTTGTTACTGGCGTTGAGCACTTCGTCGCCCCTGT‑GGGCTGGCCAGCCTACCGGGTACTTGAGCTATCGCCGAGTGATCTGCGGTGAGTGGCCACACATGGGCTTGCCCGAAGCATTGCCCGACTGGGCCGCCTATGAGCGCTACCGTGCCCTGTTGCAGCGAACCGGTGCACTGGCTGC  >  NC_002947/3688229‑3688518

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