Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,984,509	G→T	6.5%	A204E (GCG→GAG)	ycaC‑II ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,984,509	0	G	T	6.5%	222.7 / 3.0	71	A204E (GCG→GAG)	ycaC‑II	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (37/29);  new base T (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGC  >  NC_002947/5984364‑5984656                                                                                                                                                   |                                                                                                                                                    gATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTg                                                                                                                                                  >  1:175014/1‑149 (MQ=255)     ccGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTGCTTCGGTCGGCTTTAGGGGGGGGTTTAGCCGCGAAGAGGCGAGTAC‑AGGTTTACATCAGCTCGAGGCCAGCGCTGCGTa                                                                                                                                              <  6:235318/149‑1 (MQ=38)     ccGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGt                                                                                                                                               >  1:205045/1‑148 (MQ=255)     ccGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGt                                                                                                                                               <  2:205045/148‑1 (MQ=255)     ccGGATTCGTGTCAGGCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCg                                                                                                                                                <  6:115707/147‑1 (MQ=255)             gtggCAGTCTGCCATTGAATCAGGCATAGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGGTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtc                                                                                                                                       <  8:287058/148‑1 (MQ=38)              tgtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtcat                                                                                                                                     >  2:19187/1‑149 (MQ=255)                tCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCTCTGCGTAGCTGTtcat                                                                                                                                     >  7:40531/1‑147 (MQ=255)                tCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCTCTGCGTAGCTGTtcat                                                                                                                                     >  7:35722/1‑147 (MQ=255)                    tCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCCGTAC‑AGGTTTACATCAGCCCGCGGCCCGCGCTGCGTAGCAGTTCATCACGtt                                                                                                                               >  5:270000/1‑149 (MQ=37)                    tCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTTAATAATCGCGAGGCCAGCGCTGCGTAGCTGTTCATCATGtt                                                                                                                               >  5:265220/1‑149 (MQ=37)                    tCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTAAGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCCGCGCTGGGTAGCTGTTCATGAGGtt                                                                                                                               >  5:164804/1‑149 (MQ=38)                      tGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCGGTAGGGGCGGGGTTACCCGCGCAGAGGCCAGTAC‑AGGTTTTCAACAGCCCGAGGCCAGCGCTGCGTAGCTGTTCCTAATGTTGc                                                                                                                             >  8:145500/1‑149 (MQ=37)                        ccTTTGAATCACGCATAGCCTGTGCCGGCCTGTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCgg                                                                                                                           <  2:191882/149‑1 (MQ=255)                        ccTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGTCAGCGCTGCGTAGCTGTTCATCAGGTTGCgg                                                                                                                           >  7:278012/1‑149 (MQ=255)                          tttGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCg                                                                                                                            >  5:157970/1‑146 (MQ=255)                          tttGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCg                                                                                                                            <  6:157970/146‑1 (MQ=255)                            tGAATCACGCATAGCCTGGGCCGGCCTCTTCGCGGTTAACCCCGCTCCTACGTGCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTa                                                                                                                         <  3:125324/147‑1 (MQ=38)                             gAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGt                                                                                                                       <  8:272116/148‑1 (MQ=255)                                          ccTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCTCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGc                                                                                                         <  8:40531/149‑1 (MQ=255)                                          ccTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCGGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCTCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGc                                                                                                         <  8:35722/149‑1 (MQ=255)                                                 cGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGGAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcag                                                                                                    >  6:280877/1‑147 (MQ=255)                                                 cGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCGGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagca                                                                                                  <  2:175014/149‑1 (MQ=255)                                                  ggCCTCTTCGCGGGTAATCGCGGTCCTAGGGTCGGGTGTTGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagcag                                                                                                 <  6:270000/149‑1 (MQ=255)                                                  ggCCTATTCGCTGTTAAACCCGATCATAGGGTCGGGTGTAGTGGCGGGTTTACCCGCGAAGATGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCTCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagcag                                                                                                 <  6:265220/149‑1 (MQ=255)                                                   gCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAgg                                                                                                >  6:121475/1‑149 (MQ=255)                                                         tCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTAcc                                                                                           >  5:289240/1‑148 (MQ=255)                                                                tAACCCCGCTCCTCCGGGCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCtt                                                                                   <  7:145500/149‑1 (MQ=255)                                                                  aaCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCg                                                                                 <  7:29833/149‑1 (MQ=255)                                                                    cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt                                                                               >  1:135205/1‑149 (MQ=255)                                                                     ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATATGCTGCGACAGCAGGTTACCCAGCCCTTCGATg                                                                              >  2:233150/1‑149 (MQ=255)                                                                     ccGCCCCTACAGCCGACCGTAGGAGCGGGATTACCCGCGAAGAGGCCGGCACTGGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt                                                                               >  1:234509/1‑149 (MQ=25)                                                                           cTAGGGTCGGGTGTTGGGGCGGTTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTg                                                                        <  6:164804/149‑1 (MQ=255)                                                                                      tGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTcgcg                                                             >  8:165400/1‑149 (MQ=255)                                                                                       gTAGGAGCGGGATTACCCGCGAAGAGGCCGGCAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCgg                                                            <  2:234509/149‑1 (MQ=255)                                                                                         aGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTg                                                          >  2:273741/1‑149 (MQ=255)                                                                                           gggCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATGAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTGCGATGTCGTTGCGCCAGTAGCGGTg                                                          >  7:251470/1‑147 (MQ=255)                                                                                           gagCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTg                                                          >  3:88225/3‑147 (MQ=255)                                                                                             gCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACGCAGCGCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGt                                                      >  4:237422/1‑149 (MQ=255)                                                                                             gCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGt                                                      >  8:20110/1‑149 (MQ=255)                                                                                               gggTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCtt                                                                                   >  3:68374/1‑118 (MQ=255)                                                                                               gggTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCtt                                                                                   <  4:68374/118‑1 (MQ=255)                                                                                               gggTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTc                                                    >  6:66584/1‑149 (MQ=255)                                                                                                       ccGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTA‑CCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCa                                                 <  3:36369/143‑1 (MQ=255)                                                                                                       ccGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTA‑CCAGCACTTCGATGTCGTTGCGCAAGTCGCGGTGCAGTTCTCAgagatc                                           >  4:36369/1‑144 (MQ=255)                                                                                                         gcgAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCAcc                                          <  5:66584/149‑1 (MQ=255)                                                                                                               agGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAAc                                    <  1:273741/149‑1 (MQ=255)                                                                                                                ggCCAGTAC‑AGGTTTACATCAGCGCGAGGTCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCATTCGCGGTGCAGTTCGCAGGCCACCGAGAAc                                    <  8:278012/148‑1 (MQ=255)                                                                                                                 gCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCAcc                                          >  3:150858/1‑141 (MQ=255)                                                                                                                 gCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCAcc                                          <  4:150858/141‑1 (MQ=255)                                                                                                                 gCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAAc                                    >  8:157114/1‑147 (MQ=255)                                                                                                                  ccAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATATGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAg                                 <  1:233150/149‑1 (MQ=255)                                                                                                                     gTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGtt                               >  6:154918/1‑148 (MQ=255)                                                                                                                     gTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTtc                              <  5:121475/149‑1 (MQ=255)                                                                                                                       aC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTtca                             >  3:267187/1‑148 (MQ=255)                                                                                                                        c‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTtcatc                           <  6:289240/149‑1 (MQ=255)                                                                                                                        c‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTtcat                            >  3:234098/1‑148 (MQ=255)                                                                                                                            gTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGt                                                           >  5:120209/1‑114 (MQ=255)                                                                                                                            gTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGt                                                           <  6:120209/114‑1 (MQ=255)                                                                                                                            gTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCatcatc                        >  2:141530/1‑149 (MQ=255)                                                                                                                             tttACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCt                       <  5:280877/149‑1 (MQ=255)                                                                                                                                   tCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCt                                                                                                           <  1:3404/59‑1 (MQ=255)                                                                                                                                   tCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCt                                                                                                           >  2:3404/1‑59 (MQ=255)                                                                                                                                   tCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCAc                 >  6:269040/1‑149 (MQ=255)                                                                                                                                     aGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCAccc               >  2:189228/1‑149 (MQ=255)                                                                                                                                       cgcgAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGc             <  2:135205/149‑1 (MQ=255)                                                                                                                                          gAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTg          >  7:257272/1‑149 (MQ=255)                                                                                                                                          gAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTg          >  2:247586/1‑149 (MQ=255)                                                                                                                                            ggCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGTCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTggg        <  4:88225/149‑1 (MQ=255)                                                                                                                                               cAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTgg         <  1:9413/145‑1 (MQ=255)                                                                                                                                               cAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTgg         >  2:9413/1‑145 (MQ=255)                                                                                                                                                  cgcTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGc  <  4:267187/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTAC‑AGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGC  >  NC_002947/5984364‑5984656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑