Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,837,495	A→C	5.4%	Q99P (CAG→CCG) ‡	PP_5116 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,837,495	0	A	C	5.4%	236.3 / 2.7	73	Q99P (CAG→CCG) ‡	PP_5116	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/32);  new base C (2/2);  total (39/34)
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.65e-01

AAATTCGTTGAAGCACGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCTGTCTT  >  NC_002947/5837347‑5837642                                                                                                                                                     |                                                                                                                                                        aaaTTCGTTGAAGCACGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCcagca                                                                                                                                                         <  4:328892/149‑1 (MQ=255)         tGAAGCACGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAg                                                                                                                                                  >  3:71368/1‑148 (MQ=255)            aGCACGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCAcc                                                                                                                                              <  5:224655/149‑1 (MQ=255)               aCGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCTCAGCCTGCTGGGCagccgt                                                                                                                                           >  5:345140/1‑144 (MQ=255)               aCGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCTCAGCCTGCTGGGCagccgt                                                                                                                                           <  6:345140/149‑6 (MQ=255)               aCGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTg                                                                                                                                            >  3:102115/1‑148 (MQ=255)                       tGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGcc                                                                                                                                     >  4:254013/1‑147 (MQ=255)                       tGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGcc                                                                                                                                     <  3:254013/147‑1 (MQ=255)                          cAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTgggg                                                                                                                                >  8:207900/1‑149 (MQ=255)                           aaTCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTgggg‑‑g                                                                                                                             <  7:345171/149‑1 (MQ=255)                            aTCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑Gc                                                                                                                            <  6:301831/149‑1 (MQ=255)                             tCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCg                                                                                                                           >  1:358253/1‑149 (MQ=255)                               gTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTgggg                                                                                                                                >  1:25765/1‑144 (MQ=255)                               gTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTgggg                                                                                                                                <  2:25765/144‑1 (MQ=255)                                ttGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGt                                                                                                                         >  4:170654/1‑148 (MQ=255)                                ttGACGTGATACTACGCTTCAAGGCAGGCACCACGCAGGACACTCTGTACACTGGCCTGGCAGCCTTGACGCTCTCAATGCTCGATGAGGGTAGCCAAGCCAACACGGCTGCCCAGCAGGCTGACCACCTGCAGCGCCTGGGG‑‑GCGGt                                                                                                                         >  8:177682/1‑148 (MQ=255)                                   aCGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATg                                                                                                                     <  5:277128/149‑1 (MQ=255)                                        aTACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGt                                                                                                                         <  4:352082/140‑1 (MQ=255)                                        aTACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGt                                                                                                                         >  3:352082/1‑140 (MQ=255)                                         tactacGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGaaaaa                                                                                                               >  2:11188/1‑149 (MQ=255)                                            tacGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAAc                                                                                                              >  7:44454/1‑147 (MQ=255)                                                  tCAAGTCCGGCACCACGCAGTACACTCTGTAGCGTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCAc                                                                                                                                               <  8:238549/110‑1 (MQ=255)                                                  tCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGt                                                                                                      >  7:33422/1‑149 (MQ=255)                                                  tCAAGCCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCACCAGGCTGACCAc                                                                                                                                               >  7:238549/1‑110 (MQ=255)                                                    aaGGCCGGCACCACGCAGGACTCTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCGAAGCCTACACGGCTGCCCAGGAGGGTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTc                                                                                                     <  7:177682/148‑1 (MQ=255)                                                      ggCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTg                                                                                                   >  5:223307/1‑148 (MQ=255)                                                          ggCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACa                                                                                              <  1:28445/149‑1 (MQ=255)                                                          ggCACCACGAAGGACAATCTGTACCCTGGACTGGCAGCCTTGACGCTCTCCATGCACGATGAGGGTAGCCAAGACAACACGGCTGCCCAGCAAGCTGAGCACCTGGAGCGCCTAGGC‑‑GCGGTCATGGAAAAACAAGTACGTCTGCAACa                                                                                              >  8:339887/1‑149 (MQ=255)                                                            caccacGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATg                                                                                            >  8:159778/1‑149 (MQ=255)                                                             accacGCAGGACACTCTGTACCCTGGCCTGCCATCCTTGACGCTCTCAATGCTCGATGAGGCTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTTGAACAtcc                                                                                           >  2:253777/1‑147 (MQ=255)                                                                acGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGACCATGCGAc                                                                                        >  3:262769/1‑149 (MQ=255)                                                                    aGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTg                                                                                    >  1:213862/1‑149 (MQ=255)                                                                     ggACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTgc                                                                                   <  4:71368/149‑1 (MQ=255)                                                                        cacTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑cc                                                                               >  1:205724/1‑149 (MQ=255)                                                                           tctGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGc                                                                                                                                                 >  4:141069/1‑83 (MQ=255)                                                                           tctGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTgc                                                                                   <  2:267183/143‑1 (MQ=255)                                                                           tctGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTgc                                                                                   >  1:267183/1‑143 (MQ=255)                                                                           tctGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTc                                                                            >  2:242036/1‑149 (MQ=255)                                                                           tctGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTc                                                                            >  3:349537/1‑149 (MQ=255)                                                                               tACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAg                                                                        <  7:159778/149‑1 (MQ=255)                                                                                      gcctggcAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGc                                                                <  8:122289/149‑1 (MQ=255)                                                                                           gcAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGcccccc                                                           >  4:173725/1‑149 (MQ=255)                                                                                           gcAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGcccccc                                                           <  2:205724/149‑1 (MQ=255)                                                                                               ccTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTtct                                                       <  1:11188/149‑1 (MQ=255)                                                                                                   gACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCt                                                                    <  5:240882/132‑1 (MQ=255)                                                                                                   gACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCt                                                                    >  6:240882/1‑132 (MQ=255)                                                                                                        tctcCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGAc                                              >  5:66621/1‑149 (MQ=255)                                                                                                            cATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGATCACCTGGATCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGGTGCGCCCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCg                                           <  1:253777/149‑1 (MQ=255)                                                                                                             aTGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGc                                         <  6:223307/149‑1 (MQ=255)                                                                                                               gCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCt                                       >  8:330185/1‑149 (MQ=255)                                                                                                                cTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATTGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGTAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCtt                                      <  1:242036/149‑1 (MQ=255)                                                                                                                 tCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGGGTGGGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCt                                       >  1:111195/1‑149 (MQ=255)                                                                                                                 tCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTg                                     >  8:156586/1‑149 (MQ=255)                                                                                                                   gATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTgg                                    <  7:207900/148‑1 (MQ=255)                                                                                                                   gATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTgg                                    <  4:102115/148‑1 (MQ=255)                                                                                                                   gATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCAACCCGGCCTTg                                     <  8:44454/147‑1 (MQ=255)                                                                                                                    aTGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCg                                  >  6:55745/1‑149 (MQ=255)                                                                                                                    aTGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCg                                  <  7:330185/149‑1 (MQ=255)                                                                                                                    aTGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCg                                  <  5:55745/149‑1 (MQ=255)                                                                                                                    aTGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGACCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGc                                   <  5:317759/148‑1 (MQ=255)                                                                                                                    aTGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGACCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGc                                   >  6:317759/1‑148 (MQ=255)                                                                                                                          gTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCTGTCATGGAAAAACAAGTCCGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTc                            <  7:339887/149‑1 (MQ=255)                                                                                                                          gTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAGACCTTAGCCCCCCTTCTCTGCTCGACCCGGCCtt                                      <  1:48659/140‑1 (MQ=255)                                                                                                                          gTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCtt                                      >  2:48659/1‑139 (MQ=255)                                                                                                                              ccAAGCCTACACGTCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCa                           <  3:7496/146‑1 (MQ=255)                                                                                                                              ccAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCa                           >  4:7496/1‑146 (MQ=255)                                                                                                                              ccAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTACCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCAcc                         >  5:319042/1‑148 (MQ=255)                                                                                                                                       cacGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTgg               <  2:358253/149‑1 (MQ=255)                                                                                                                                          ggCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGccc            >  8:117102/1‑149 (MQ=255)                                                                                                                                            gtGCTCAGCCTGCTGGGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCt                                                  <  2:31791/108‑1 (MQ=255)                                                                                                                                            gtGCTCAGCCTGCTGGGCACCTGGAGCGCCAGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCt                                                  >  1:31791/2‑109 (MQ=255)                                                                                                                                               cccAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCt       <  4:349537/149‑1 (MQ=255)                                                                                                                                                ccAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCt       <  5:143287/148‑1 (MQ=255)                                                                                                                                                ccAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCt       >  6:143287/1‑148 (MQ=255)                                                                                                                                                 cagcagGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCTGt     <  3:170654/149‑1 (MQ=255)                                                                                                                                                   gcagGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCTGTCt   >  5:262788/1‑149 (MQ=255)                                                                                                                                                   gcagGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCTGTCt   >  2:30913/1‑149 (MQ=255)                                                                                                                                                    cagGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCTGTCtt  <  2:111195/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                        AAATTCGTTGAAGCACGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCTGTCTT  >  NC_002947/5837347‑5837642

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