Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,805,968 G→A 5.7% G186D (GGC→GAC)  mrcA → penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,805,9680GA 5.7% 209.9 / 2.1 70G186D (GGC→GAC) mrcApenicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (32/34);  new base A (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

TTCTCACGCAAGACCAACGAAATCCTGCTGGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGG  >  NC_002947/5805823‑5806111
                                                                                                                                                 |                                                                                                                                                
ttCTCACGCAAGACCAACGAAATCCTGCTGGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATTGATGCAGCAGCACAGGTGTATTACGGc                                                                                                                                                 >  4:313893/1‑147 (MQ=255)
       gCAAGACCAACGAAATCCTGCTGGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATc                                                                                                                                        >  3:228611/1‑149 (MQ=255)
       gCAAGACCAACGAAATCCTGCTGGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGATGTATTACGGCAAGTCCATc                                                                                                                                        <  6:59734/149‑1 (MQ=255)
        cAAGACCAACGAAATCCTGCTGGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATcc                                                                                                                                       <  6:290849/149‑1 (MQ=255)
             ccAACGAAATCCTGCTGGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCcgtgac                                                                                                                                  <  5:274694/149‑1 (MQ=255)
                aCGAAATCCTGCTGGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCcgtgactg                                                                                                                                >  5:283194/1‑146 (MQ=255)
                     aTCCTGCTGGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCt                                                                                                                          <  4:228611/149‑1 (MQ=255)
                        ctgctgGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGgc                                                                                                                       >  5:205715/1‑149 (MQ=255)
                        ctgctgGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGgc                                                                                                                       <  4:121980/149‑1 (MQ=255)
                              gCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGa                                                                                                                  <  2:222628/148‑1 (MQ=255)
                              gCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGAt                                                                                                                 >  2:236966/1‑149 (MQ=255)
                                ccTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATgg                                                                                                               >  6:116344/1‑149 (MQ=255)
                                       aTCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCCTCCGTGACGTGAGCCTGGCGCAGATGGCgaggct                                                                                                        >  8:107339/1‑145 (MQ=255)
                                                     gACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGccc                                                                                           >  4:5137/1‑148 (MQ=255)
                                                        cAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAgg                                                                                       >  1:247995/1‑149 (MQ=255)
                                                        cAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAgg                                                                                       >  6:288164/1‑149 (MQ=255)
                                                        cAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAgg                                                                                       <  5:116344/149‑1 (MQ=255)
                                                        aaaGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCAGGCCTGCCCAAgg                                                                                       >  2:99364/2‑149 (MQ=255)
                                                           ggACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGACGGCCTGCCCAAGGCAc                                                                                    >  7:12575/1‑149 (MQ=255)
                                                               gAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTc                                                                                >  6:3167/1‑149 (MQ=255)
                                                                 aaTCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTcgc                                                                              >  8:29790/1‑149 (MQ=255)
                                                                  aTCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTcgcg                                                                             <  3:58339/149‑1 (MQ=255)
                                                                          gCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAAc                                                                      <  6:205715/148‑1 (MQ=255)
                                                                              tACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGgcgc                                                                                                                     <  3:170174/97‑1 (MQ=255)
                                                                              tACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGgcgc                                                                                                                     >  4:170174/1‑97 (MQ=255)
                                                                              tACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCt                                                                <  2:247995/149‑1 (MQ=255)
                                                                               aCGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTg                                                               >  3:34703/1‑149 (MQ=255)
                                                                                cGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTgg                                                              <  5:206801/149‑1 (MQ=255)
                                                                                cGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTg                                                               <  6:14544/148‑1 (MQ=255)
                                                                                cGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTg                                                               <  1:266380/148‑1 (MQ=255)
                                                                                    aaCAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCaa                                                          >  2:17302/1‑149 (MQ=255)
                                                                                      cAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAAcc                                                        >  4:221520/1‑149 (MQ=255)
                                                                                               ccTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTgcgcg                                               >  2:212582/1‑149 (MQ=255)
                                                                                               ccTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTgc                                                  <  3:5137/146‑1 (MQ=255)
                                                                                               ccTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTg                                                   <  1:17302/145‑1 (MQ=255)
                                                                                                  gggCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGcc                                             >  7:166556/1‑148 (MQ=255)
                                                                                                   ggCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGcc                                             <  7:29790/147‑1 (MQ=255)
                                                                                                      aaCCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCaaa                                          >  7:43385/1‑147 (MQ=255)
                                                                                                      aaCCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCaaa                                          <  1:197620/147‑1 (MQ=255)
                                                                                                      aaCCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCaaa                                          >  2:197620/1‑147 (MQ=255)
                                                                                                           cgcCTACGGTATCGATGCAGCAGCACAGGTGTATTACGACAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGcgcc                                   >  2:288164/1‑149 (MQ=255)
                                                                                                           cgcCTACGGTATCGATGCAGCAGCACAGGTGTATTACGACAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGcgcc                                   >  2:234851/1‑149 (MQ=255)
                                                                                                            gcCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGcgccg                                  <  4:34703/149‑1 (MQ=255)
                                                                                                                 cGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGa                               <  5:187560/147‑1 (MQ=255)
                                                                                                                 cGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGa                               <  5:3167/147‑1 (MQ=255)
                                                                                                                 cGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACt                             >  2:305158/1‑149 (MQ=255)
                                                                                                                  ggTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTgcg                                                 >  8:82821/1‑128 (MQ=255)
                                                                                                                  ggTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTgcg                                                 <  7:82821/128‑1 (MQ=255)
                                                                                                                  ggTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTg                            >  3:130059/1‑149 (MQ=255)
                                                                                                                   gTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTg                            >  6:119100/1‑148 (MQ=255)
                                                                                                                    tATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGa                          >  2:295930/1‑149 (MQ=255)
                                                                                                                          tGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCg                    <  3:221520/149‑1 (MQ=255)
                                                                                                                            cagcagcaCAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACACCGCTGGCCAACCCGGTGCGCGCCAAGGAGCGCCGCGACTGGATCCTCg                    <  1:236966/148‑1 (MQ=255)
                                                                                                                               cagcaCAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGcgcc                                   <  2:203806/129‑1 (MQ=255)
                                                                                                                               cagcaCAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGcgcc                                   >  1:203806/1‑129 (MQ=255)
                                                                                                                               cagcaCAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGgcgc               <  7:278297/149‑1 (MQ=255)
                                                                                                                               cagcaCAGGTGTATTACGACAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGgcgc               <  1:288164/149‑1 (MQ=255)
                                                                                                                               cagcaCAGGTGTATTACGACAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGgcgc               <  1:234851/149‑1 (MQ=255)
                                                                                                                                 gcaCAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCAt             <  8:43385/149‑1 (MQ=255)
                                                                                                                                  cacaGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATg            <  1:212582/149‑1 (MQ=255)
                                                                                                                                  cacaGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATg            >  7:87508/1‑149 (MQ=255)
                                                                                                                                     aGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTa          >  6:254385/1‑148 (MQ=255)
                                                                                                                                     aGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAAGCCCCGCGACTGGATCCTCGGGCGCATGTAc         >  2:98494/1‑148 (MQ=255)
                                                                                                                                      ggTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACa        <  1:99364/149‑1 (MQ=255)
                                                                                                                                        tgtATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAg      <  8:87508/149‑1 (MQ=255)
                                                                                                                                            ttACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTTGATCCTCGGGCGCATGTACAAGCTgg  <  8:12575/149‑1 (MQ=255)
                                                                                                                                            ttACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTgg  <  3:313893/149‑1 (MQ=255)
                                                                                                                                            ttACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTgg  <  2:138141/149‑1 (MQ=255)
                                                                                                                                            ttACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTgg  >  1:138141/1‑149 (MQ=255)
                                                                                                                                            ttACGGCAAGTCCATCCGGGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTgg  <  7:107339/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCTCACGCAAGACCAACGAAATCCTGCTGGCCCTGCAAATCGAGCGTGAGCTGACCAAGGACGAAATCCTTGAGCTGTACGTGAACAAGATCTACCTGGGCAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGG  >  NC_002947/5805823‑5806111

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: