Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,141,207 T→G 8.4% N119H (AAC→CAC)  PP_4525 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,141,2070TG 8.4% 169.3 / 4.9 58N119H (AAC→CAC) PP_4525LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base G (3/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGATCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCGGACAT  >  NC_002947/5141060‑5141351
                                                                                                                                                    |                                                                                                                                                  
tGGTGATCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGtt                                                                                                                                                    >  2:261171/1‑148 (MQ=255)
      tCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGc                                                                                                                                             <  1:261171/149‑1 (MQ=255)
        aGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCAGAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCgg                                                                                                                                           <  4:189418/149‑1 (MQ=255)
          tACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGt                                                                                                                                         <  6:236240/149‑1 (MQ=255)
            ccAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTaa                                                                                                                                       <  8:217845/149‑1 (MQ=255)
            ccAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTGAACCTGcaagt                                                                                                                                         <  8:256664/147‑5 (MQ=255)
            ccAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTGAACCTGcaagt                                                                                                                                         >  7:256664/1‑143 (MQ=255)
              aGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGa                                                                                                                                     >  1:63753/1‑149 (MQ=255)
              aGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGa                                                                                                                                     <  2:301970/149‑1 (MQ=255)
                  gCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTgg                                                                                                                                 >  8:67339/1‑149 (MQ=255)
                  gCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTg                                                                                                                                  <  1:9520/148‑1 (MQ=255)
                  ccTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATACAGCGTGTGGTTCAAAAGCCCTTCAAGCAACTCGCCGCTGGGTCCTGCACGCACTTGCAGATTCACATGCGGGTAAGACAgg                                                                                                                                 >  4:264825/2‑149 (MQ=255)
                   cTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTTGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGt                                                                                                                                <  3:45421/149‑1 (MQ=255)
                                 gCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAg                                                                                                                   >  4:64359/1‑148 (MQ=255)
                                             gCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAgcgcgc                                                                                                                      <  8:209646/133‑1 (MQ=255)
                                             gCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAgcgcgc                                                                                                                      >  7:209646/1‑133 (MQ=255)
                                              cccGGCAAGGCTCGGTGGCACGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAACCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGATCAACTGCGGGTAAGACAGGTGATAGCGCGCCAGCAACGCCGGCAAGt                                                                                                     >  5:112596/1‑149 (MQ=255)
                                                 ggCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTggg                                                                                                  <  8:268377/149‑1 (MQ=255)
                                                   cAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTc                                                                                                <  3:64359/149‑1 (MQ=255)
                                                   cAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTGCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTc                                                                                                <  5:185575/149‑1 (MQ=255)
                                                               gCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCt                                                                                     <  5:87174/148‑1 (MQ=255)
                                                                 ccGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCg                                                                                                           <  8:249864/124‑1 (MQ=255)
                                                                 ccGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCg                                                                                                           <  8:245304/124‑1 (MQ=255)
                                                                 ccGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCg                                                                                                           >  7:249864/1‑124 (MQ=255)
                                                                 ccGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCg                                                                                                           >  7:245304/1‑124 (MQ=255)
                                                                  cGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGAGTGGGTGCTGCACGCACTTGCAGGTACACCTGCGGGTAACACTGGTGAAAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGtcatga                                                                                   >  1:83840/1‑143 (MQ=255)
                                                                      cacaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATgg                                                                           >  8:133641/1‑149 (MQ=255)
                                                                      cacaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATgg                                                                           >  2:26239/1‑149 (MQ=255)
                                                                       acaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGt                                                                          >  8:1229/1‑149 (MQ=255)
                                                                       acaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGCGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCCCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGt                                                                          >  5:298050/1‑149 (MQ=255)
                                                                           gagCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTcccc                                                                      <  2:101472/149‑1 (MQ=255)
                                                                                gcaTCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCAc                                                                 <  2:63753/149‑1 (MQ=255)
                                                                                 caTCGAGCGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCAGTTGCAGTTTCAGCTGCGTGTAAGACTGTTGATAGCGCGCCAGCAACGCCGGCAAGTTGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCTGCACg                                                                <  6:112596/149‑1 (MQ=255)
                                                                                  aTCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACg                                                                >  8:188149/1‑148 (MQ=255)
                                                                                          tgtgGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCa                                                       <  7:133641/149‑1 (MQ=255)
                                                                                           gtgGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAg                                                      >  2:44174/1‑149 (MQ=255)
                                                                                              gTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCg                                                   >  4:68667/1‑149 (MQ=255)
                                                                                               ttCAACAGCGCTTCAAGCAAATCGCCGCAGGGTGCTGCACGCACTTGCAGGTGAACCTGCga                                                                                                                                           >  2:178639/1‑61 (MQ=255)
                                                                                                 cAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGc                                                 <  1:44174/148‑1 (MQ=255)
                                                                                                   aCAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGa                                              >  7:285788/1‑149 (MQ=255)
                                                                                                    cAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGAc                                             >  2:249822/1‑149 (MQ=255)
                                                                                                         cTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCca                                        >  3:154247/1‑149 (MQ=255)
                                                                                                          ttCAAGGAACTCGCCGCTGGGTGCTGCTCGCAGTTGCAGGTTCACCTGCGGGTAAGAGTGGTGTTTGCGCGGCAGCAACGCCGGCAAGTGGGTCGCCGGCGTGCTTTA‑‑CATGGTccc                                                                       <  8:185517/117‑1 (MQ=255)
                                                                                                          ttCAAGCAACTCGCCGCTGGGTGCTCCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCCCCAGCAACGCCCGCAAGTGGCTCCCCGCCGTGCTGTA‑‑CATCGAccc                                                                       >  7:185517/1‑117 (MQ=255)
                                                                                                            cAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCcacg                                      >  4:28361/1‑148 (MQ=255)
                                                                                                                      gccgcTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCa                           >  2:273460/1‑149 (MQ=255)
                                                                                                                      gccgcTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGAGCTGTA‑‑CATGGTCCCCAGCACGAAGTccc                                                         <  8:65246/119‑1 (MQ=255)
                                                                                                                      gccgcTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGAGCTGTA‑‑CATGGTCCCCAGCACGAAGTccc                                                         >  7:65246/1‑119 (MQ=255)
                                                                                                                       ccgcTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAg                          <  3:28361/149‑1 (MQ=255)
                                                                                                                        cgcTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAgc                         <  8:285788/149‑1 (MQ=255)
                                                                                                                                     cacgcacTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTcgc               <  6:290813/147‑1 (MQ=255)
                                                                                                                                          acTTGCAGGTTCACCTGCGGGTAAGACTTGTGATAGCGCGCCAGCAACGCCGGCAATTGGGGCGCCGCCGTGGTGTACTCCAGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAg          <  1:249822/149‑1 (MQ=255)
                                                                                                                                            ttGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCg       >  3:52280/1‑148 (MQ=255)
                                                                                                                                            ccgCAGGTGAACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTtcgtcgt               <  1:64083/138‑2 (MQ=255)
                                                                                                                                            ccgCAGGTGAACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTtcgtcgt               >  2:64083/3‑139 (MQ=255)
                                                                                                                                              gCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCGGa     <  1:273460/148‑1 (MQ=255)
                                                                                                                                              gCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCGGa     <  5:249804/148‑1 (MQ=255)
                                                                                                                                                aGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCGGACa   >  6:112059/1‑148 (MQ=255)
                                                                                                                                                aGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCGGACAt  <  4:154247/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGTGATCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTA‑‑CATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCGGACAT  >  NC_002947/5141060‑5141351

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: