Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	232,740	A→G	6.5%	Y280C (TAC→TGC)	PP_0180 →	cytochrome c family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	232,740	0	A	G	6.5%	177.3 / 2.1	62	Y280C (TAC→TGC)	PP_0180	cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (2/2);  total (31/31)
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.21e-01

GACCTGGAACCTCGCCGACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGACAACATCGATACCCAGGCGCGCAAG  >  NC_002947/232600‑232885                                                                                                                                             |                                                                                                                                                   gACCTGGAACCTCGCCGACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGc                                                                                                                                              >  5:305002/1‑147 (MQ=255) gACCTGGAACCTCGCCGACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGc                                                                                                                                              >  7:274051/1‑147 (MQ=255)  aCCTGGAACCTCGCCGACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAg                                                                                                                                           <  2:36/149‑1 (MQ=255)    cTGGAACCTCGCCGACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGca                                                                                                                                         >  1:97950/1‑149 (MQ=255)     tGGAACCTCGCCGACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGca                                                                                                                                         >  5:366895/1‑148 (MQ=255)         aCCTCGCCGACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTg                                                                                                                                    <  2:383516/149‑1 (MQ=255)              gCCGACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGgag                                                                                                                                 >  2:63573/1‑147 (MQ=255)                cGACCTGGCTCGCCAAACCCCGGCGGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCGAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGagaaga                                                                                                                             <  8:302425/149‑1 (MQ=255)                 gACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAagag                                                                                                                            >  3:81606/1‑149 (MQ=255)                        cTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCCCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑ggcg                                                                                                                    <  3:5533/149‑1 (MQ=255)                         tCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑ggcgg                                                                                                                   >  4:232417/1‑149 (MQ=255)                           gCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑ggcggc                                                                                                                  <  5:97580/148‑1 (MQ=255)                           gCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑ggcggc                                                                                                                  <  2:365092/148‑1 (MQ=255)                           gCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGcc                                                                                                                 >  3:72036/1‑149 (MQ=255)                             cAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTa                                                                                                               >  3:154293/1‑149 (MQ=255)                                aCCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGc                                                                                                                                              <  3:210680/116‑1 (MQ=255)                                aCCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGc                                                                                                                                              >  4:210680/1‑116 (MQ=255)                                 ccccGGCCGAAGTTGCTGCTAACGAAGTCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACcgc                                                                                                           <  1:343471/149‑1 (MQ=255)                                  cccGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATGCACTGCCAGca                                                                                                                                         >  5:11468/1‑119 (MQ=255)                                  cccGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATGCACTGCCAGca                                                                                                                                         >  5:16962/1‑119 (MQ=255)                                  cccGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATGCACTGCCAGca                                                                                                                                         <  6:11468/119‑1 (MQ=255)                                  cccGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATGCACTGCCAGca                                                                                                                                         <  6:16962/119‑1 (MQ=255)                                         gAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTg                                                                                                                                    <  8:61872/117‑1 (MQ=255)                                         gAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTg                                                                                                                                    >  7:61872/1‑117 (MQ=255)                                          aaGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGAc                                                                                                  >  7:41664/1‑149 (MQ=255)                                                 tgctAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGAt                                                                                            >  5:124387/1‑148 (MQ=255)                                                 tgctAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGAt                                                                                            >  8:228864/1‑148 (MQ=255)                                                 tgctAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGAt                                                                                            >  7:41018/1‑148 (MQ=255)                                                      aCGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGc                                                                                       <  7:228864/148‑1 (MQ=255)                                                         aaGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCTCGGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTa                                                                                    <  1:63573/149‑1 (MQ=255)                                                            gCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGAc                                                                                <  4:81606/149‑1 (MQ=255)                                                                   ggCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCgg                                                                         >  4:221761/1‑149 (MQ=255)                                                                          gAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTggcgg                                                                   <  3:232417/148‑1 (MQ=255)                                                                               aTTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTAcc                                                             >  5:333452/1‑149 (MQ=255)                                                                                  ccGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTg                                                           >  7:329906/1‑148 (MQ=255)                                                                                  ccGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTg                                                           >  4:260658/1‑148 (MQ=255)                                                                                   cgcgCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTgg                                                          >  3:376474/1‑148 (MQ=255)                                                                                     cgcCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAg                                                       <  8:41664/149‑1 (MQ=255)                                                                                        ccAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGtt                                                    >  1:337121/1‑149 (MQ=255)                                                                                               gCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCt                                             <  3:221761/149‑1 (MQ=255)                                                                                               gCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCt                                             >  2:137646/1‑149 (MQ=255)                                                                                               gCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCt                                             >  6:73106/1‑149 (MQ=255)                                                                                                 ccAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCtgg                                           >  3:187466/1‑149 (MQ=255)                                                                                                   aGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCtggtg                                         <  6:366895/149‑1 (MQ=255)                                                                                                   aGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCtggtg                                         <  4:69660/149‑1 (MQ=255)                                                                                                   aGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCtggtg                                         <  6:305002/149‑1 (MQ=255)                                                                                                   aGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCtggtg                                         >  3:69660/1‑149 (MQ=255)                                                                                                       aCCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGaaa                                     <  8:39194/149‑1 (MQ=255)                                                                                                        ccTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAg                                    <  4:154293/149‑1 (MQ=255)                                                                                                          tCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCt                                  >  5:1575/1‑149 (MQ=255)                                                                                                            aGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCg                                <  1:137646/149‑1 (MQ=255)                                                                                                             ggTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGc                               <  8:329906/149‑1 (MQ=255)                                                                                                               tAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTg                             <  4:376474/149‑1 (MQ=255)                                                                                                                  atcgcgGCCCGGCCCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGTCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGac                          <  6:231571/147‑1 (MQ=255)                                                                                                                     gcgGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGa                           >  2:319033/1‑145 (MQ=255)                                                                                                                     gcgGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGa                           <  1:319033/145‑1 (MQ=255)                                                                                                                         cccGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGACAACATCg                   <  2:97950/149‑1 (MQ=255)                                                                                                                          ccGGCGCAATTGCTCGAATCCACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGACAACATCGa                  <  8:328001/149‑1 (MQ=255)                                                                                                                                       tCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGACAACATCGATACCCAGgcgcg      >  5:69676/1‑148 (MQ=255)                                                                                                                                       tCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGACAACATCGATACCCAGgcgcg      >  1:347130/1‑148 (MQ=255)                                                                                                                                          aaTACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGACAACATCGATACCCAGGCGCGCAAg  <  6:124387/149‑1 (MQ=255)                                                                                                                                           aTACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGACAACATCGATACCCAGGCGCGCAAg  <  4:187466/148‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                   GACCTGGAACCTCGCCGACCTGGCTCGCCAAACCCCGGCCGAAGTTGCTGCTAACGAAGGCACTGCGGCCGTCGAGGCATTCCGCGCCCAGCGTGCCCAGCCACCTCAGGTAAAGCGCGGCCCGGCGCAATTGCTCGAATACACTGCCAGCACACTGGAGAAGAGCCT‑GGCGGCCTACCGCGCAGGCGACCACGATCAGGCCTATGACCTGTCGGTGGCGGCCTACCTGGAAGGTTTCGAGCTGGTGGAAAGCTCGCTGGACAACATCGATACCCAGGCGCGCAAG  >  NC_002947/232600‑232885

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