Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,848 A→T 6.7% intergenic (‑137/+43) aroQ‑I ← / ← dsbD‑I type II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947650,8480AT 6.7% 177.5 / 2.4 60intergenic (‑137/+43)aroQ‑I/dsbD‑Itype II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base T (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCAGC  >  NC_002947/650702‑650995
                                                                                                                                                  |                                                                                                                                                   
gCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAAt                                                                                                                                                   <  2:313369/149‑1 (MQ=255)
   ttGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCg                                                                                                                                                <  1:204048/149‑1 (MQ=255)
      ccATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTc                                                                                                                                             <  8:59712/149‑1 (MQ=255)
      ccATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGt                                                                                                                                              <  6:378388/148‑1 (MQ=255)
           tGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTTAATACGGTCACtaa                                                                                                                                        >  8:202012/1‑147 (MQ=255)
             cGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCAtt                                                                                                                                      <  8:195102/149‑1 (MQ=255)
                   ccTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGGTCGCGCGAAGTATGTGTCCGGATATTTAAATCCGGTCACTCATTTCAg                                                                                                                                  <  1:35689/147‑1 (MQ=255)
                    cTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGa                                                                                                                                >  8:197809/1‑148 (MQ=255)
                    cTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAt                                                                                                                               >  7:52065/1‑149 (MQ=255)
                    cTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAt                                                                                                                               >  2:43964/1‑149 (MQ=255)
                    cTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGCCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGa                                                                                                                                <  7:197809/148‑1 (MQ=255)
                       ggATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAtttt                                                                                                                            >  5:19979/1‑149 (MQ=255)
                                gggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTc                                                                                                                   >  8:311969/1‑149 (MQ=255)
                                 ggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcg                                                                                                                  >  5:105092/1‑149 (MQ=255)
                                 ggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcg                                                                                                                  >  8:281814/1‑149 (MQ=255)
                                      ttCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGa                                                                                                             <  7:202012/149‑1 (MQ=255)
                                        cGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCTTTTCAGGATTTTTCCGAGCTCGCGGGaa                                                                                                            <  1:209667/148‑1 (MQ=255)
                                               gcgcGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGa                                                                                                    >  7:15484/1‑149 (MQ=255)
                                                cgcgGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGAc                                                                                                   <  6:105092/149‑1 (MQ=255)
                                                   ggcgAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGa                                                                                                >  7:108289/1‑149 (MQ=255)
                                                       aGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACttt                                                                                            >  1:83932/1‑149 (MQ=255)
                                                           tgCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCcg                                                                                        >  4:354354/1‑149 (MQ=255)
                                                              cGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCAGTTCAGGAGTTTTCGGAGCTCGCGGGATGGTATTGACCGCCTTTGCcgag                                                                                      >  8:347159/1‑146 (MQ=255)
                                                                  aaGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgc                                                                                 <  2:83932/149‑1 (MQ=255)
                                                                   aGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgcg                                                                                >  2:2426/1‑149 (MQ=255)
                                                                    tttaCGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgcgc                                                                               <  4:188153/146‑1 (MQ=255)
                                                                      aaCGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGt                                                                             >  8:199898/1‑149 (MQ=255)
                                                                           cAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGaa                                                                        >  8:318001/1‑149 (MQ=255)
                                                                                 tCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgc                                                                   <  8:52065/148‑1 (MQ=255)
                                                                                   cAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgcgc                                                                 <  7:281814/148‑1 (MQ=255)
                                                                                      ttCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgc                                                                   >  7:177187/1‑143 (MQ=255)
                                                                                      ttCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgc                                                                   <  8:177187/143‑1 (MQ=255)
                                                                                         ccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCa                                                          >  8:66588/1‑149 (MQ=255)
                                                                                           ggCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAcc                                                        <  1:2426/149‑1 (MQ=255)
                                                                                            gCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCg                                                       <  7:199898/149‑1 (MQ=255)
                                                                                             cAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTcc                                                                                                                         >  1:241161/1‑82 (MQ=255)
                                                                                             cAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTcc                                                                                                                         <  2:241161/82‑1 (MQ=255)
                                                                                                 aGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCAtt                                                  <  8:15484/149‑1 (MQ=255)
                                                                                                     aGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtc                                              >  7:322211/1‑149 (MQ=255)
                                                                                                      gCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtct                                             <  8:29246/149‑1 (MQ=255)
                                                                                                       cACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtctc                                            <  6:153304/149‑1 (MQ=255)
                                                                                                         cGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGc                                          <  1:43964/149‑1 (MQ=255)
                                                                                                            tGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCaa                                       <  6:19979/149‑1 (MQ=255)
                                                                                                             gTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAc                                                         <  8:36014/130‑1 (MQ=255)
                                                                                                             gTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAc                                                         >  7:36014/1‑130 (MQ=255)
                                                                                                              tCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAAt                                      >  1:75955/1‑148 (MQ=255)
                                                                                                              tCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAAt                                      >  5:320333/1‑148 (MQ=255)
                                                                                                                         cgcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAcc                                  <  4:314120/141‑1 (MQ=255)
                                                                                                                         cgcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAcc                                  >  3:314120/1‑141 (MQ=255)
                                                                                                                         cgcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGc                          >  2:177067/1‑149 (MQ=255)
                                                                                                                             aaGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTa                      <  8:108289/149‑1 (MQ=255)
                                                                                                                                  tgtgACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTc                 >  3:3649/1‑149 (MQ=255)
                                                                                                                                  tgtgACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTc                 >  6:99186/1‑149 (MQ=255)
                                                                                                                                      aCCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCAt                                                   <  2:136575/111‑1 (MQ=255)
                                                                                                                                      aCCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCAt                                                   >  1:136575/1‑111 (MQ=255)
                                                                                                                                      aCCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTg             >  8:108498/1‑149 (MQ=255)
                                                                                                                                        cGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCg                                           <  1:137038/117‑1 (MQ=255)
                                                                                                                                        cGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCg                                           >  2:137038/1‑117 (MQ=255)
                                                                                                                                          gATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGt         >  6:158768/1‑149 (MQ=255)
                                                                                                                                               ttAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCg      >  3:313929/1‑147 (MQ=255)
                                                                                                                                                  aaTCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgcagc  >  5:246190/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCAGC  >  NC_002947/650702‑650995

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: