Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,914,111	G→A	10.3%	intergenic (+3/+150)	spuI → / ← argA	glutamylpolyamine synthetase/amino acid N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,914,111	0	G	A	10.3%	184.0 / 9.4	68	intergenic (+3/+150)	spuI/argA	glutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base A (3/4);  total (33/35)
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.54e-01

GCGAGCTGGACGACAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGG  >  NC_002947/5913963‑5914259                                                                                                                                                      |                                                                                                                                                     gcgAGCTGGACGACAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACg                                                                                                                                                      <  6:316401/149‑1 (MQ=255)     gCTGGACGACAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACGTCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAaaaa                                                                                                                                                  >  6:345280/1‑149 (MQ=255)      cTGGACGACAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAAc                                                                                                                                                 >  3:328880/1‑149 (MQ=255)           cgacAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAAcgccgc                                                                                                                                            <  2:139399/149‑1 (MQ=255)            gacAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACgccgcc                                                                                                                                           >  1:172854/1‑149 (MQ=255)             acAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGccc                                                                                                                                          <  5:57679/149‑1 (MQ=255)                 cGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTc                                                                                                                                       >  7:47389/1‑148 (MQ=255)                       ccTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGggcg                                                                                                                                >  3:66779/1‑149 (MQ=255)                           aaCAAGTACATCGATCCGAAGTACGTCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGggcggc                                                                                                                              <  5:345280/147‑1 (MQ=255)                            aCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtt                                                                                                                           >  3:366103/1‑149 (MQ=255)                            aCAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtt                                                                                                                           >  8:196863/1‑149 (MQ=255)                             cAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAATGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACATGCATACCTTGTAAACGAAAACGCCGCCCATCAAGGGCGGCGtt                                                                                                                           >  1:261735/1‑149 (MQ=255)                              aaGTACATCGATCCGATGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGTAGAGTTCGAGCAGTCGATCTCCGACCTCGAGTACATCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGtttt                                                                                                                         <  7:261502/149‑1 (MQ=255)                               aGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttttt                                                                                                                        <  6:203701/149‑1 (MQ=255)                                gTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttg                                                                                                                       >  6:243548/1‑149 (MQ=255)                                gTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttg                                                                                                                       <  7:338927/149‑1 (MQ=255)                                   cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                    <  2:172854/149‑1 (MQ=255)                                   cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgtt                                                                                                                     <  7:151356/148‑1 (MQ=255)                                   cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgtt                                                                                                                     <  4:66779/148‑1 (MQ=255)                                           cGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCgg                                                                                                            >  3:99113/1‑149 (MQ=255)                                           cGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCgg                                                                                                            >  3:262878/1‑149 (MQ=255)                                             aaGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGgtg                                                                                                          <  5:243548/149‑1 (MQ=255)                                                  cATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGgt                                                                                                           <  3:3677/143‑1 (MQ=255)                                                  cATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGgt                                                                                                           >  4:3677/1‑143 (MQ=255)                                                  cATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCC                                                                                                     <  2:261735‑M1/149‑1 (MQ=255)                                                   aTCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCC                                                                                                    >  3:280557‑M1/1‑149 (MQ=255)                                                      gACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCT                                                                                                 >  7:204067‑M1/1‑149 (MQ=255)                                                       aCATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGCAAACAAAAACGCCGCCCTTGAGGGGCGGCGTTttcgtttacacg                                                                                                               >  6:57838/1‑129 (MQ=37)                                                       aCATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGCAAACAAAAACGCCGCCCTTGAGGGGCGGCGTTttcgtttacacg                                                                                                               <  5:57838/134‑6 (MQ=37)                                                         aTCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAACACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCG                                                                                              >  2:295512‑M1/1‑149 (MQ=255)                                                          tCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGC                                                                                             >  3:322580‑M1/1‑149 (MQ=255)                                                           cTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCG                                                                                            >  5:4600‑M1/1‑149 (MQ=255)                                                                tcgCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCAC                                                                                       >  7:149094‑M1/1‑149 (MQ=255)                                                                   cgTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCC                                                                                    <  4:328880‑M1/149‑1 (MQ=255)                                                                      gCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGC                                                                                 >  3:203537‑M1/1‑149 (MQ=255)                                                                        aa‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTC                                                                               >  4:28966‑M1/1‑149 (MQ=255)                                                                                gCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                        <  6:7424/149‑1 (MQ=255)                                                                                gCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                        <  8:47389/149‑1 (MQ=255)                                                                                           gagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCg                                                               >  3:317000/1‑147 (MQ=255)                                                                                           gagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                             <  8:149094/149‑1 (MQ=255)                                                                                           gagTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGaa                                                             <  3:28966/149‑1 (MQ=255)                                                                                                     aCTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGcc                                                   >  8:340574/1‑149 (MQ=255)                                                                                                       tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagcct                                                    >  4:229650‑M1/1‑121 (MQ=255)                                                                                                       tCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGGAATTTCGAAGGCTGAAGCCtg                                                 >  6:170940/1‑149 (MQ=255)                                                                                                             tcCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaatta                                           >  8:44495‑M1/1‑115 (MQ=255)                                                                                                               cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaa                                              <  3:229650‑M1/144‑32 (MQ=255)                                                                                                               cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGAc                                             <  4:317000/145‑1 (MQ=255)                                                                                                               cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                         <  4:262878/149‑1 (MQ=255)                                                                                                               cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                         <  6:4600/149‑1 (MQ=255)                                                                                                               cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                         <  4:203537/149‑1 (MQ=255)                                                                                                               cGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCTTGAGGGGCGGCGtt                                                                                                                           <  4:97107/67‑1 (MQ=21)                                                                                                                gACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccc                                        <  4:322580‑M1/149‑38 (MQ=255)                                                                                                                   cTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                     <  1:84540/149‑1 (MQ=255)                                                                                                                    tCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTg                                    <  8:86341‑M2/149‑2 (MQ=255)                                                                                                                     cGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCCTCAAGGGCGGCGTTtttgt                                                                                                                      <  5:105120/66‑1 (MQ=255)                                                                                                                     cGAGTACAACTGGTACCTGCATACCGTGTAAACAAAAACGCCGCCCCTCAAGGGCGGCGTTtttgt                                                                                                                      >  6:105120/1‑66 (MQ=255)                                                                                                                      gAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgg                                   <  7:44495‑M1/148‑43 (MQ=255)                                                                                                                           cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgg                             <  4:366103‑M1/149‑49 (MQ=255)                                                                                                                           cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACACCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgg                             >  4:280557‑M2/1‑141 (MQ=255)                                                                                                                           cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcg                              <  7:196863‑M1/148‑48 (MQ=255)                                                                                                                            aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                <  3:265719/129‑1 (MQ=255)                                                                                                                            aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                >  4:265719/1‑129 (MQ=255)                                                                                                                            aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACCCCCGCGCCCGCAGGGTAATTTCGAATGCTGAAGCCGGTGACTTACCAGTgggagcggg                            <  3:220530‑M2/149‑10 (MQ=255)                                                                                                                            aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCACTTCGCGCGCAACCCCCCACCCGCACGGTAAATTAGAAGGCTGAAACCTGTGACATAACACTGGGacccgg                            >  3:69193/1‑144 (MQ=255)                                                                                                                                 gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgggagcgggcgt                         >  7:340574‑M2/1‑135 (MQ=255)                                                                                                                                 gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGCCGCCCTCTTCGCGGGCACGCCCGCTCCCacaggtaagtcacaggcttcagccttcgaaattaccctgcgggagcgggcgtgc                       <  3:141469‑M1/149‑55 (MQ=255)                                                                                                                                     cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTggga                                 >  7:41795‑M2/1‑131 (MQ=255)                                                                                                                                     cTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTggga                                 <  8:41795‑M2/135‑5 (MQ=255)                                                                                                                                                    aCGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTAGCGGGCACCCCCCCGCCCGCACGGGATTTTGGGAGGCTGGAAGCGGGGGGGTAACGGTGGGGGGGGGGGTGGCCGCGGAGAGGGGGGgccc    >  4:351761/1‑148 (MQ=255)                                                                                                                                                      gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCCCGCCCGCGGCCCGAGGGGTATTTGGAAGGCGGGAGGCGGGGGCTTAGCTGTGGGGGCGGGGGGGGCCGGGCAGCGGGCGGGCCCgg  >  8:29504/1‑149 (MQ=255)                                                                                                                                                      gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGTGTGCCCGCGAAGAGGGCGGCCCAgg  >  6:267355/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GCGAGCTGGACGACAGCGAGATCCTGAACAAGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAA‑GGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTACAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGG  >  NC_002947/5913963‑5914259

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