Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,837,490	C→T	5.5%	A97A (GCC→GCT)	PP_5116 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,837,490	0	C	T	5.5%	216.0 / 2.1	72	A97A (GCC→GCT)	PP_5116	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base T (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

CGTCAAATTCGTTGAAGCACGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCT  >  NC_002947/5837343‑5837637                                                                                                                                                    |                                                                                                                                                        cGTCAAATTCGTTGAAGCACGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGccc                                                                                                                                                        <  1:200917/149‑1 (MQ=255)    cAAATTCGTTGAAGCACGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACACTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCcag                                                                                                                                                      >  2:28445/1‑148 (MQ=255)     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                                                                                                                                <  3:254013/147‑1 (MQ=255)                           tGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGcc                                                                                                                                >  4:254013/1‑147 (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                                                                                                                           <  2:25765/144‑1 (MQ=255)                                   gTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTgggg                                                                                                                           >  1:25765/1‑144 (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                                                               >  6:240882/1‑132 (MQ=255)                                                                                                       gACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCt                                                               <  5:240882/132‑1 (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                               <  4:102115/148‑1 (MQ=255)                                                                                                                       gATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTgg                               <  7:207900/148‑1 (MQ=255)                                                                                                                       gATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCAACCCGGCCTTg                                <  8:44454/147‑1 (MQ=255)                                                                                                                        aTGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCg                             <  7:330185/149‑1 (MQ=255)                                                                                                                        aTGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCg                             <  5:55745/149‑1 (MQ=255)                                                                                                                        aTGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCg                             >  6:55745/1‑149 (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)                                                                                                                                                    |                                                                                                                                                        CGTCAAATTCGTTGAAGCACGGGAACTGCCAATCGTTGACGTGATACTACGGTTCAAGGCCGGCACCACGCAGGACACTCTGTACCCTGGCCTGGCAGCCTTGACGCTCTCCATGCTCGATGAGGGTAGCCAAGCCTACACGGCTGCCCAGCAGGCTGAGCACCTGGAGCGCCTGGGG‑‑GCGGTCATGGAAAAACAAGTACGTCTGGAACATGCGACGCTGCG‑CCTTCGCAG‑CCTTAGCCCCCCTTCTCTGCTCGACCCGGCCTTGGCGGTGTTCACCGACCTCGTGGCCCACCCT  >  NC_002947/5837343‑5837637

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