Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,799,613	C→A	6.7%	A136E (GCG→GAG)	rnhB →	ribonuclease HII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,799,613	0	C	A	6.7%	186.4 / 3.3	60	A136E (GCG→GAG)	rnhB	ribonuclease HII
Reads supporting (aligned to +/- strand):  ref base C (25/31);  new base A (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

CAAGCCACCATGCTGGCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGG  >  NC_002947/1799465‑1799746                                                                                                                                                     |                                                                                                                                       cAAGCCACCATGCTGGCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTgcggc                                                                                                                                        <  6:189966/149‑1 (MQ=255)  aaGCCACCATGCTGGCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTgcggcg                                                                                                                                       >  6:147748/1‑149 (MQ=255)    gCCACCATGCTGGCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCggcg‑gc                                                                                                                                    <  7:222392/149‑1 (MQ=255)      caccaTGCTGGCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCAt                                                                                                                                  >  8:38359/1‑149 (MQ=255)                gCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATc                                                                                                                                 >  8:84524/1‑140 (MQ=255)                gCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATc                                                                                                                                 <  7:84524/140‑1 (MQ=255)                gCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTg                                                                                                                          <  5:147748/147‑1 (MQ=255)                gCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGc                                                                                                                        >  8:151388/1‑149 (MQ=255)                gCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGc                                                                                                                        <  6:200172/149‑1 (MQ=255)                  cATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGcc                                                                                                                       >  7:349047/1‑148 (MQ=255)                   aTGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGcc                                                                                                                       >  2:101858/1‑147 (MQ=255)                      cAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGt                                                                                                                  <  3:241200/149‑1 (MQ=255)                         cGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGAc                                                                                                               <  3:207947/149‑1 (MQ=255)                         cGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGAc                                                                                                               <  7:151388/149‑1 (MQ=255)                          gggCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATAGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGCTGAcc                                                                                                              >  6:240901/1‑149 (MQ=255)                             cAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGt                                                                                                           <  8:349047/149‑1 (MQ=255)                              aGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCGGAAGCTGGCGGTGCCGGCAGCTCGGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTg                                                                                                          <  7:312507/149‑1 (MQ=255)                                  gAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGTCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATc                                                                                                       <  5:240901/148‑1 (MQ=255)                                    ggggTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGt                                                                                                           <  6:16668/142‑1 (MQ=255)                                    ggggTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGt                                                                                                           >  5:16668/1‑142 (MQ=255)                                      ggTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGCCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGgag                                                                                                  <  3:327765/149‑1 (MQ=255)                                           cacaTTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCgg                                                                                                     >  4:203177/1‑141 (MQ=255)                                           cacaTTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCgg                                                                                                     <  3:203177/141‑1 (MQ=255)                                           cacaTTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGa                                                                                              <  5:243391/148‑1 (MQ=255)                                                 aCGCCAACGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGt                                                                                        <  4:38184/148‑1 (MQ=255)                                                   gCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATc                                                                                                                             >  2:8936/1‑109 (MQ=255)                                                   gCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATc                                                                                                                             <  1:8936/109‑1 (MQ=255)                                                   gCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTc                                                                                      >  3:312263/1‑148 (MQ=255)                                                     cAACGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGaa                                                                                   <  4:287960/149‑1 (MQ=255)                                                      aaaGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAAc                                                                                  <  8:118156/149‑1 (MQ=255)                                                      aaaGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAAc                                                                                  >  3:109518/1‑149 (MQ=255)                                                         gCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGa                                                                               >  2:112191/1‑149 (MQ=255)                                                             gCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCGCGCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGAt                                                                              <  1:101858/147‑1 (MQ=255)                                                                   aTCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCTGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCgg                                                                      <  3:162288/148‑1 (MQ=255)                                                                    tCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCgg                                                                                                     >  8:211628/1‑116 (MQ=255)                                                                    tCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCgg                                                                                                     <  7:211628/116‑1 (MQ=255)                                                                       aTGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTAt                                                                 >  2:203399/1‑149 (MQ=255)                                                                        tGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCAGGTTATg                                                                >  1:293142/1‑149 (MQ=255)                                                                            aaTCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGAG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCAt                                                            >  1:125371/1‑149 (MQ=255)                                                                            aaTCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGAG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCAt                                                            >  1:129471/1‑149 (MQ=255)                                                                             aTCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCAtt                                                           <  2:118274/149‑1 (MQ=255)                                                                                  tGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCgg                                                      <  4:312263/149‑1 (MQ=255)                                                                                    cccGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGc                                                    <  7:38359/149‑1 (MQ=255)                                                                                    cccGAAGCTGGCGGGGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCgg                                                      <  6:167864/147‑1 (MQ=255)                                                                                      cGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGCTCTACCCGGGTTATGGCATTGGCGGGca                                                   >  2:132449/1‑148 (MQ=255)                                                                                         aGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAg                                               <  7:26590/149‑1 (MQ=255)                                                                                                ggTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGa                                                                               >  6:61363/1‑110 (MQ=255)                                                                                                ggTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGa                                                                               <  5:61363/110‑1 (MQ=255)                                                                                                       gcagcaCCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTg                                                                                                                          >  7:66185/1‑60 (MQ=255)                                                                                                       gcagcaCCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTg                                                                                                                          <  8:66185/60‑1 (MQ=255)                                                                                                           caCCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTAc                                         >  1:146502/1‑137 (MQ=255)                                                                                                            aCCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTg                             >  3:187788/1‑148 (MQ=255)                                                                                                              cGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCAt                          >  5:314288/1‑149 (MQ=255)                                                                                                                gtggtCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCAtct                        >  7:126934/1‑149 (MQ=255)                                                                                                                      aaGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCATCTCGAGGc                  >  4:208467/1‑149 (MQ=255)                                                                                                                              tAGCCAGGTGCCGGCCATTGCGGAG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCAc          <  2:129471/149‑1 (MQ=255)                                                                                                                              tAGCCAGGTGCCGGCCATTGCGGAG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCAc          <  2:125371/149‑1 (MQ=255)                                                                                                                                 ccAGGTGCCGGCCATTGCGGCG‑GCATCCATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCGCAAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTc       >  4:22034/1‑149 (MQ=255)                                                                                                                                  cAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCt      >  8:77144/1‑149 (MQ=255)                                                                                                                                    tgTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTg    <  7:245557/148‑1 (MQ=255)                                                                                                                                    ggTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGTTGCATCTCGAGGCGTTGGCACGTCTTg    <  4:109518/149‑1 (MQ=255)                                                                                                                                     gTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTgg   <  1:112191/149‑1 (MQ=255)                                                                                                                                      tGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTggg  <  8:126934/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                       CAAGCCACCATGCTGGCCATGCAGCGGGCAGTGGAGGGGTTGCACATTACGCCAAAGCTGGCCCTGATCGATGGCAATCGCTGCCCGAAGCTGGCGGTGCCGGCAGCACCGGTGGTCAAGGGGGATAGCCAGGTGCCGGCCATTGCGGCG‑GCATCGATCCTGGCCAAGGTGACCCGTGATCGGGAGATGAGTGCGTTCGAACTGATCTACCCGGGTTATGGCATTGGCGGGCACAAGGGCTACCCGACCCCGGTGCATCTCGAGGCGTTGGCACGTCTTGGG  >  NC_002947/1799465‑1799746

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