Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,786,058 Δ1 bp 6.7% coding (1151/1242 nt) mhpT ← 3‑(3‑hydroxy‑phenyl)propionate transporter MhpT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,786,0580A. 6.7% 185.8 / 7.0 60coding (1151/1242 nt)mhpT3‑(3‑hydroxy‑phenyl)propionate transporter MhpT
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base . (2/2);  total (29/31)
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.98e-01

GCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCACCGCCACCGAAGCGCCCAAGCCCGTTGCGCG  >  NC_002947/3785913‑3786194
                                                                                                                                                 |                                                                                                                                         
gCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCa                                                                                                                                        <  2:31546/149‑1 (MQ=255)
 ccGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCAc                                                                                                                                       >  5:239816/1‑149 (MQ=255)
  cGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCAc                                                                                                                                       >  1:284093/1‑148 (MQ=255)
    aaCCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAg                                                                                                                                    >  4:22140/1‑149 (MQ=255)
         gAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTg                                                                                                                               <  5:306427/149‑1 (MQ=255)
          aGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTgg                                                                                                                              <  7:46064/149‑1 (MQ=255)
           gCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGa                                                                                                                             <  2:42906/149‑1 (MQ=255)
                 ggCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCcgcg                                                                                                                      >  7:160505/1‑149 (MQ=255)
                   cATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGc                                                                                                                    <  1:195277/149‑1 (MQ=255)
                     ttGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCCCCCGGCTCCGGGCGTTGCGTGCGGGGTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCgg                                                                                                                     <  7:210759/146‑1 (MQ=255)
                     ttGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCgg                                                                                                                     >  8:210759/1‑146 (MQ=255)
                       gcgcAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCAt                                                                                                                 <  8:69292/148‑1 (MQ=255)
                       gcgcAAGCAGCACTAGCCCCCCAAGGAACTGCGATGCAGGGTCCCGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCAt                                                                                                                 <  8:255701/148‑1 (MQ=255)
                         gcAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATcag                                                                                                              >  4:118375/1‑149 (MQ=255)
                           aaGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGGCTCCAGGCCTggtg                                                                                                                            <  1:128885/133‑3 (MQ=255)
                           aaGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGGCTCCAGGCCTggtg                                                                                                                            >  2:128885/1‑131 (MQ=255)
                             gcagcaCTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCagcagc                                                                                                          <  2:284093/149‑1 (MQ=255)
                                    tAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCg                                                                                                   >  1:200728/1‑149 (MQ=255)
                                         ccccAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCgcacct                                                                                              >  8:118979/1‑149 (MQ=255)
                                          cccAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCgcacctg                                                                                             >  7:89652/1‑149 (MQ=255)
                                           ccAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCgcacctg                                                                                             >  6:100717/1‑148 (MQ=255)
                                            cAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCgcacctgcac                                                                                          >  7:53904/1‑149 (MQ=255)
                                              aGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCgcacctgcacc                                                                                         >  1:38199/1‑149 (MQ=255)
                                                gAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGcacctgcacctg                                                                                       <  3:209967/149‑1 (MQ=255)
                                                    tGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCtgctgc                                                                                   <  2:38199/149‑1 (MQ=255)
                                                      cTATGCAGGGGCCTGGGCTGCAGCCGCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCaa                                                                                 <  6:251211/149‑1 (MQ=255)
                                                        aTGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAg                                                                                >  4:16414/1‑148 (MQ=255)
                                                         tGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAgg                                                                               <  3:22140/148‑1 (MQ=255)
                                                           cAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGc                                                                                                                                  <  5:23192/96‑1 (MQ=255)
                                                           cAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGc                                                                                                                                  >  6:23192/1‑96 (MQ=255)
                                                           cAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCacctgcacctgc                                                                                      <  1:281091/139‑1 (MQ=255)
                                                           cAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCacctgcacctgc                                                                                      >  2:281091/1‑139 (MQ=255)
                                                            aGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGGCTCCAGGCCTGGTGCTCGCGGc                                                                                                                    <  2:34685/108‑1 (MQ=38)
                                                            aGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGGCTCCAGGCCTGGTGCTCGCGGc                                                                                                                    >  1:34685/1‑108 (MQ=38)
                                                                tCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACctgctg                                                                                    >  2:182524/1‑136 (MQ=255)
                                                                tCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACctgctg                                                                                    <  1:182524/136‑1 (MQ=255)
                                                                  cTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCgg                                                                     >  3:153106/1‑149 (MQ=255)
                                                                   tGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGgc                                                                    >  4:114143/1‑149 (MQ=255)
                                                                      gCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGgccgc                                                                 <  6:239816/149‑1 (MQ=255)
                                                                          cAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCaga                                                             <  5:100717/149‑1 (MQ=255)
                                                                          cAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCa                                                               <  2:200728/147‑1 (MQ=255)
                                                                           aGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCagag                                                            <  8:160505/149‑1 (MQ=255)
                                                                             ccTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGa                                                          >  5:213887/1‑149 (MQ=255)
                                                                             ccTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGa                                                          >  2:238548/1‑149 (MQ=255)
                                                                                   ccGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCACCCTCACCAGCACCTCCATAAAGGATCTGCACGCCCGCCAccggcccggc                                                     >  8:292215/1‑148 (MQ=255)
                                                                                     ggCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCAt                                                  <  6:213887/149‑1 (MQ=255)
                                                                                          cGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAAc                                             <  3:16414/149‑1 (MQ=255)
                                                                                          cGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCC‑‑‑‑‑‑GCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAg                                         <  2:242260/147‑1 (MQ=255)
                                                                                           gggCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAAcc                                            <  5:103837/149‑1 (MQ=255)
                                                                                              cGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAg                                         <  3:118375/149‑1 (MQ=255)
                                                                                                          cttCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGAc                                                         <  6:36869/120‑1 (MQ=255)
                                                                                                              ggTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCccaccgccacc                         >  6:203837/1‑149 (MQ=255)
                                                                                                                       gTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGc                                                     <  7:202120/112‑1 (MQ=255)
                                                                                                                       gTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGc                                                     >  8:202120/1‑112 (MQ=255)
                                                                                                                           cgcACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCACCGCCACCGAAGCGCCCAAGc            <  7:292215/149‑1 (MQ=255)
                                                                                                                                gTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCACCGCCACCGAAGCGCCCAAGCCCGtt       >  7:30859/1‑149 (MQ=255)
                                                                                                                                gTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCACCGCCACCGAAGCGCCCAAGCCCGtt       <  8:89652/149‑1 (MQ=255)
                                                                                                                                  gccgccAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCACCGCCACCGAAGCGCCCAAGCCCGTTgc     >  4:301532/1‑149 (MQ=255)
                                                                                                                                  gccgccAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCACCGCCACCGAAGCGCCCAAGCCCGTTgc     >  1:145803/1‑149 (MQ=255)
                                                                                                                                     gccAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCACCGCCACCGAAGCGCCCAAGCCCGTTgcgcg  >  1:174249/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
GCCGAACCGGAGCAGGAGGCATTGCGCAAGCAGCACTAGCCCCCCAAGGAAGTGCTATGCAGGGTCCTGGGCTGCAGCCTCACCCGGCTCCGGGCGTTGCGTGCGGGTTCGGTGTTTCAGTAGCGCACGTGCCGCCAACGCCGCGAGCACCAGGCCTGGAGCCGCGGCCATCAGCAGCCCCCCCGCACCTGCACCTGCTGCAAGGATCTGCCCGGCCGCCAGAGGACCGGCCATGGAACCCAGGCGCCCCACCGCCACCGAAGCGCCCAAGCCCGTTGCGCG  >  NC_002947/3785913‑3786194

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: