Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,869,614	G→A	8.5%	G41E (GGG→GAG) ‡	PP_5144 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,869,614	0	G	A	8.5%	165.2 / 4.8	59	G41E (GGG→GAG) ‡	PP_5144	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/21);  new base A (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 = 6.36e-01

TGACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTGCA  >  NC_002947/5869466‑5869761                                                                                                                                                     |                                                                                                                                                     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                                                                                                                      >  5:110937/1‑149 (MQ=255)                                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCg                                                                                                                      <  4:279375/149‑1 (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                                                                                                             >  2:146596/1‑132 (MQ=255)                                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggc                                                                                                             <  1:146596/132‑1 (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                                                                                <  2:81936/125‑1 (MQ=255)                                                                                               acgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGcc                                                                                >  1:81936/1‑125 (MQ=255)                                                                                                cgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑gg                                                      >  2:243245/1‑149 (MQ=255)                                                                                                    ttCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGc                                                    >  7:38224/1‑147 (MQ=255)                                                                                                    ttCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGc                                                    >  7:43303/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)                                                                                                                                                     gggTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTGCCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTGCa  <  5:127827/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TGACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTGCA  >  NC_002947/5869466‑5869761

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