Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,869,613	G→C	8.1%	G41R (GGG→CGG) ‡	PP_5144 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,869,613	0	G	C	8.1%	185.9 / 5.3	59	G41R (GGG→CGG) ‡	PP_5144	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/21);  new base C (2/3);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTG  >  NC_002947/5869466‑5869759                                                                                                                                                    |                                                                                                                                                    tGACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTgg                                                                                                                                                    >  4:194427/1‑149 (MQ=255)       aaaaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCCCCTGGGGTGCa                                                                                                                                              >  2:235181/1‑149 (MQ=255)       aaaaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCa                                                                                                                                              >  4:66421/1‑149 (MQ=255)         aaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCTGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCCCAGGTGCAtc                                                                                                                                            <  2:162306/149‑2 (MQ=255)                gagGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGgcg                                                                                                                                     <  8:238710/149‑1 (MQ=255)                   gCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGa                                                                                                                                  >  3:279375/1‑149 (MQ=255)                     ccAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAcc                                                                                                                                >  3:114500/1‑149 (MQ=255)                         ccTATGTGCCGGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGa                                                                                                                            <  3:35975/149‑1 (MQ=255)                            aTGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGagg                                                                                                                          >  6:9477/1‑148 (MQ=255)                                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGaggcaggc                                                                                                                     <  5:9477/148‑1 (MQ=255)                                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCg                                                                                                                    <  4:279375/149‑1 (MQ=255)                                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCg                                                                                                                    >  5:110937/1‑149 (MQ=255)                                   tGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCgg                                                                                                                   >  5:186746/1‑148 (MQ=255)                                                cTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCg                                                                                                     >  8:214417/1‑149 (MQ=255)                                                   gcggcCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAg                                                                                                  >  2:213170/1‑149 (MQ=255)                                                         ggggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCcac                                                                                            <  3:194427/149‑1 (MQ=255)                                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggc                                                                                                           <  1:146596/132‑1 (MQ=255)                                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggc                                                                                                           >  2:146596/1‑132 (MQ=255)                                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCg                                                                                                     <  1:336541/138‑1 (MQ=255)                                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCg                                                                                                     >  2:336541/1‑138 (MQ=255)                                                              aCGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATTCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAgggg                                                                                       >  1:29029/1‑149 (MQ=255)                                                                gCCCTGCCGCTGATCGTGGCGGCCAACCGCGATGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCCGGAGGTACCTGGCTGGGCGTAGGCCCACAGGGgcg                                                                                     >  4:65760/1‑149 (MQ=255)                                                                gCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGgcg                                                                                     >  2:87501/1‑149 (MQ=255)                                                                  ccggCCGCTGATCGT‑GCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCtt                                                                                  <  4:124974/146‑1 (MQ=255)                                                                  ccTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCt                                                                                   >  4:281268/1‑149 (MQ=255)                                                                    tGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCT‑GCTGGGCGTAGGCCCACAGGGGCGCTTTg                                                                                >  4:259521/1‑149 (MQ=255)                                                                            aTCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCt                                                                         <  2:238575/149‑1 (MQ=255)                                                                             tCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGTGGCGCTTTGCCGCGCt                                                                         <  6:110937/148‑1 (MQ=255)                                                                                 ggcggcCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCa                                                                    <  1:87501/149‑1 (MQ=255)                                                                                    ggcCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACa                                                                 <  1:245731/149‑1 (MQ=255)                                                                                    ggcCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACa                                                                 <  6:186746/149‑1 (MQ=255)                                                                                        aaCCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCAGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCAACAGGGGCGCTTTGCCGCGCTGACCAACATCcg                                                             >  1:207203/1‑149 (MQ=255)                                                                                            gcgACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGAc                                                         >  2:82554/1‑149 (MQ=255)                                                                                            gcgACGAATTCTAAGCCCGCCCCACTCAGCCACTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTCCCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCAGCGCTGACCAACAGCcgcgcc                                                         >  6:324584/1‑147 (MQ=255)                                                                                              gacgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGAccc                                                       <  4:114500/149‑1 (MQ=255)                                                                                               acgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGcc                                                                              >  1:81936/1‑125 (MQ=255)                                                                                               acgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGcc                                                                              <  2:81936/125‑1 (MQ=255)                                                                                                cgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑gg                                                    >  2:243245/1‑149 (MQ=255)                                                                                                    ttCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGc                                                  >  7:43303/1‑147 (MQ=255)                                                                                                    ttCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGc                                                  >  7:38224/1‑147 (MQ=255)                                                                                                     tCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCCGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAg                                               >  2:98700/1‑149 (MQ=255)                                                                                                     tCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCCGGGGTGCATGCAGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCAAg                                               >  2:94399/1‑149 (MQ=255)                                                                                                          gcccgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGgcg                                            >  1:196221/1‑147 (MQ=255)                                                                                                             cgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAggg                                       <  3:66421/149‑1 (MQ=255)                                                                                                                   aCTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCctcg                                 >  5:267167/1‑149 (MQ=255)                                                                                                                       aGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCctcgctcg                             >  2:314978/1‑149 (MQ=255)                                                                                                                        gCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCtcgctcgc                            >  2:155328/1‑149 (MQ=255)                                                                                                                          cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCgg                          >  8:34041/1‑149 (MQ=255)                                                                                                                          cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCgg                          <  7:214417/149‑1 (MQ=255)                                                                                                                          cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCgg                          >  2:236435/1‑149 (MQ=255)                                                                                                                            cTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCgggg                        >  8:223377/1‑149 (MQ=255)                                                                                                                               gcggcCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCAGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTg                      <  1:213170/149‑1 (MQ=255)                                                                                                                                       ggaggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggt               <  1:314978/147‑1 (MQ=255)                                                                                                                                       ggaggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCT‑GCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggtg              <  3:259521/147‑1 (MQ=255)                                                                                                                                         ggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                       <  6:306120/134‑1 (MQ=255)                                                                                                                                         ggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                       <  6:311168/134‑1 (MQ=255)                                                                                                                                         ggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                       >  5:311168/4‑137 (MQ=255)                                                                                                                                         ggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                       >  5:306120/4‑137 (MQ=255)                                                                                                                                           gaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGc         <  1:243245/149‑1 (MQ=255)                                                                                                                                                   tGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTGCCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTg  >  6:127827/1‑148 (MQ=255)                                                                                                                                                    ggggTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCAACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTg  >  5:336976/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTG  >  NC_002947/5869466‑5869759

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