Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,077,170	C→T	7.7%	R228Q (CGG→CAG)	PP_2687 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,077,170	0	C	T	7.7%	151.2 / 2.1	52	R228Q (CGG→CAG)	PP_2687	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base T (2/2);  total (24/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 = 9.87e-01

AGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGAACATTTCATAGCGTGAATC  >  NC_002947/3077023‑3077315                                                                                                                                                    |                                                                                                                                                    aGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACAccggccg                                                                                                                                                    <  7:50562/149‑1 (MQ=255)  ggCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACAccggccgg                                                                                                                                                   <  5:122980/149‑1 (MQ=255)  ggCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACAccggccgg                                                                                                                                                   <  5:237045/149‑1 (MQ=255)          ttCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTggg                                                                                                                                           >  2:265437/1‑149 (MQ=255)            cAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGt                                                                                                                                          <  7:205556/148‑1 (MQ=255)                  ggcggcCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGgcgc                                                                                                                                   <  1:265437/149‑1 (MQ=255)                       ccATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATag                                                                                                                              <  1:18784/149‑1 (MQ=255)                            cgcTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTc                                                                                                                         <  6:189248/149‑1 (MQ=255)                                gggCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGtt                                                                                                                      >  8:203723/1‑148 (MQ=255)                                   cggcACTCCCGGTACTGCCCCGCAGGCCTGCCAGGATCAGGTCCAGATCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACACGGACACCGGCCGGGCGGTGGGTGCGGCGCTATAGAGCTCGTTGAAg                                                                                                                  <  3:98548/149‑1 (MQ=255)                                    ggcACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTgg                                                                                                                                            >  8:53500/1‑122 (MQ=255)                                    ggcACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTgg                                                                                                                                            <  7:53500/122‑1 (MQ=255)                                      cACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCt                                                                                                               <  7:203723/149‑1 (MQ=255)                                       aCGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCtt                                                                                                              >  8:138931/1‑149 (MQ=255)                                              gTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑Ga                                                                                                     <  6:299777/149‑1 (MQ=255)                                                  tGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGt                                                                                                  >  6:170494/1‑148 (MQ=255)                                                  tGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTc                                                                                                 >  7:172511/1‑149 (MQ=255)                                                         cAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAgcgcgc                                                                                          <  4:249835/149‑1 (MQ=255)                                                             ccTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCCCGCgcag                                                                                      <  3:269221/149‑1 (MQ=255)                                                                 ccAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGcagcag                                                                                   <  8:172511/148‑1 (MQ=255)                                                                           gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCAtt                                                                        >  2:93288/1‑149 (MQ=255)                                                                           gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCAtt                                                                        >  2:49468/1‑149 (MQ=255)                                                                           gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCAtt                                                                        >  1:21203/1‑149 (MQ=255)                                                                           gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCAtt                                                                        <  1:11528/149‑1 (MQ=255)                                                                            tCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTg                                                                       <  4:200284/149‑1 (MQ=255)                                                                             ccAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGt                                                                                                                                          <  6:52062/83‑1 (MQ=255)                                                                             ccAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTgg                                                                      >  7:8712/1‑149 (MQ=255)                                                                              cAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGt                                                                     <  7:138931/149‑1 (MQ=255)                                                                                gTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGtaa                                                                   >  8:263249/1‑149 (MQ=255)                                                                                 ttCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGtaa                                                                   >  6:140003/1‑148 (MQ=255)                                                                                   ccTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGtaata                                                                 <  5:94992/148‑1 (MQ=255)                                                                                    cTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAAc                                                               >  7:251743/1‑149 (MQ=255)                                                                                        aTGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACcgcg                                                           >  2:97632/1‑149 (MQ=255)                                                                                           cAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCgcgggc                                                        <  5:121147/149‑1 (MQ=255)                                                                                              gCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCgcgggcg                                                       >  6:238669/1‑147 (MQ=255)                                                                                              gCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGa                                                     <  7:263249/149‑1 (MQ=255)                                                                                                aTGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAgt                                                   >  7:272651/1‑149 (MQ=255)                                                                                                 tGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAgtc                                                  <  5:52550/149‑1 (MQ=255)                                                                                                    agcgagcgCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCTGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAgtcgtc                                               >  6:206242/1‑149 (MQ=255)                                                                                                    agcgagcgCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCTGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAgtcgtc                                               >  6:203306/1‑149 (MQ=255)                                                                                                      cgagcgCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑ACTTAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGtcgtcg                                              <  1:93288/149‑1 (MQ=255)                                                                                                        agcgCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTCAC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCt                                            <  1:97632/149‑1 (MQ=255)                                                                                                        agcgCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTg                                           <  5:170494/149‑1 (MQ=255)                                                                                                         gcgcGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTg                                                                       <  1:284422/120‑1 (MQ=255)                                                                                                         gcgcGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTg                                                                       >  2:284422/1‑120 (MQ=255)                                                                                                                     cACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAgg                              >  8:273760/1‑149 (MQ=255)                                                                                                                     cACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAg                               <  5:140003/148‑1 (MQ=255)                                                                                                                        ccGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCt                           >  1:19392/1‑149 (MQ=255)                                                                                                                            aaCAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTTACCTCTTCGGAGGCCTCGcc                       >  8:146287/1‑149 (MQ=255)                                                                                                                              cagcagCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGa                     <  8:251743/149‑1 (MQ=255)                                                                                                                                   gcaACTCGGACACCGGCTGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGAACAtt                <  5:203306/149‑1 (MQ=255)                                                                                                                                   gcaACTCGGACACCGGCTGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGAACAtt                <  5:206242/149‑1 (MQ=255)                                                                                                                                        tCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACctc                                      <  2:9033/122‑1 (MQ=255)                                                                                                                                        tCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGAACATTTCa             >  2:320538/1‑147 (MQ=255)                                                                                                                                        tCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGAACATTTCa             >  6:74414/1‑147 (MQ=255)                                                                                                                                               ccggccggGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGAACATTTCATAGCGTGaa    <  5:74414/149‑1 (MQ=255)                                                                                                                                                 ggccggGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGAACATTTCATAGCGTGAATc  <  5:238669/149‑1 (MQ=255)                                                                                                                                                 ggccggGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGAACATTTCATAGCGTGAATc  <  8:8712/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGGCACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGT‑AC‑GAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTCGTCGCTGACCTCTTCGGAGGCCTCGCCGAACATTTCATAGCGTGAATC  >  NC_002947/3077023‑3077315

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑