Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,966,454	C→T	6.1%	V413V (GTG→GTA)	amtB ←	ammonium transporter AmtB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,966,454	0	C	T	6.1%	190.5 / 2.2	65	V413V (GTG→GTA)	amtB	ammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base C (25/36);  new base T (2/2);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

AAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGTGAGGATGGCAC  >  NC_002947/5966307‑5966602                                                                                                                                                     |                                                                                                                                                     aaaaaCTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGc                                                                                                                                                      <  6:322663/148‑1 (MQ=255)      cTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCAGCAGGCCCATCACCACATCCAGCACCtt                                                                                                                                                 <  1:22888/149‑1 (MQ=255)      cTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTg                                                                                                                                                <  1:177310/149‑1 (MQ=255)         gCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTga                                                                                                                                               <  3:139856/147‑1 (MQ=255)         gCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTga                                                                                                                                               <  8:61590/147‑1 (MQ=255)           aaaaaaaTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGAGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAAt                                                                                                                                           >  1:19816/1‑149 (MQ=255)                 aTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTa                                                                                                                                     <  2:375573/149‑1 (MQ=255)                  tGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg                                                                                                                                    <  1:48342/149‑1 (MQ=255)                  tGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg                                                                                                                                    <  6:152115/149‑1 (MQ=255)                    cccGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGAGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGt                                                                                                                                  <  2:19816/149‑1 (MQ=255)                       ggCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGAc                                                                                                                               >  6:78031/1‑149 (MQ=255)                          aaGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGACTCACGTAGGTGACGAt                                                                                                                            >  6:86812/1‑149 (MQ=255)                            gCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGCCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATgg                                                                                                                          <  1:313115/149‑1 (MQ=255)                                ggCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGgt                                                                                                                      <  2:210872/149‑1 (MQ=255)                                ggCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGgt                                                                                                                      <  4:110895/149‑1 (MQ=255)                                   aaGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTa                                                                                                                   <  7:124850/149‑1 (MQ=255)                                   aaGTTTTTACCGCGTCGCAGTTAAAGGTTGTAGCCGCGCTCGTTGTGATGAGCCAGGTCGAGTCCGACCGACTCTTCTTCTTGGTTGACCC‑GCAGGCCCATCACCACATCCATCACCTTGAGAATCAAGTAGGTGACGATGGCGGTGTa                                                                                                                   <  4:351598/149‑1 (MQ=255)                                       ttttACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGAc                                                                                                                <  2:27799/148‑1 (MQ=255)                                           aCCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACgg                                                                                                           >  7:12081/1‑149 (MQ=255)                                           aCCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACg                                                                                                            <  5:78031/148‑1 (MQ=255)                                            ccGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACgg                                                                                                           <  2:365869/148‑1 (MQ=255)                                            ccGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACgg                                                                                                           <  7:357965/148‑1 (MQ=255)                                             cgcgTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTg                                                                                                         <  6:146315/149‑1 (MQ=255)                                                     aGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATcac                                                                                                 <  5:127020/149‑1 (MQ=255)                                                          cAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATc                                                                                                                                          <  4:250406/103‑1 (MQ=255)                                                          cAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATc                                                                                                                                          >  3:250406/1‑103 (MQ=255)                                                            ggTTGTAGCCGAGCTCGTTGTGGTGAGACAGGTCGAGGCCGACCGAATCTTCTTCTTCGTTGAGCC‑GCAGGCCCATCACCAGATCCATCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTg                                                                                          <  5:86812/149‑1 (MQ=255)                                                                    cggcgcTCGTTGTGCTGAGCAAGGTCGAGGCCTACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCa                                                                                                              <  2:315033/119‑1 (MQ=255)                                                                     cgcgCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATc                                                                                 >  3:201051/1‑149 (MQ=255)                                                                             ttGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTg                                                                         >  1:182150/1‑149 (MQ=255)                                                                                 gCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTg                                                                                     <  4:369910/133‑1 (MQ=255)                                                                                  cTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGTACCTTGAGAATCACGTAGGTGACGATGGc                                                                                                                         <  5:253338/96‑1 (MQ=255)                                                                                  cTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGTACCTTGAGAATCACGTAGGTGACGATGGc                                                                                                                         >  6:253338/1‑96 (MQ=255)                                                                                  cTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTg                                                                                     <  5:11896/132‑1 (MQ=255)                                                                                  cTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTg                                                                                     >  6:11896/1‑132 (MQ=255)                                                                                         aGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAAt                                                                                                    <  8:26783/110‑1 (MQ=255)                                                                                         aGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAAt                                                                                                    >  7:26783/1‑110 (MQ=255)                                                                                          ggTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCg                                                            >  3:131005/1‑149 (MQ=255)                                                                                            tCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGt                                                          <  8:71297/149‑1 (MQ=255)                                                                                            tCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGt                                                          >  4:105345/1‑149 (MQ=255)                                                                                             cGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGt                                                          >  5:148483/1‑148 (MQ=255)                                                                                                ggCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGaaa                                                                                                      <  8:271162/101‑1 (MQ=255)                                                                                                ggCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGaaa                                                                                                      >  7:271162/1‑101 (MQ=255)                                                                                                    gaccgacTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTAc                                                  >  8:206332/1‑149 (MQ=255)                                                                                                         acTCTTCTTCTTCGTTGACCC‑GGAGGGCCATCACCTCATCCAGCACCTTGAGAATCACTTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACg                                                      <  5:203935/140‑1 (MQ=255)                                                                                                         acTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGAACACGGTGAAAAACAAAACATTTAAAATAAAACAGACTTGCATGCAGATGTCGGTGACggagatcgc                                              >  6:203935/1‑141 (MQ=255)                                                                                                           tcttcttctTCGTTGACCC‑GCAGGCCCATCACCACATCCAGTACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAgcc                                           >  3:11712/1‑149 (MQ=255)                                                                                                             ttcttcttcGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAgccgc                                         >  1:344590/1‑149 (MQ=255)                                                                                                               cttcttcGTTGACCC‑GCAGGCCCATCACCACATCCAGTACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAgccgc                                         <  4:11712/147‑1 (MQ=255)                                                                                                               cttcttcGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAgccgc                                         <  2:344590/147‑1 (MQ=255)                                                                                                                ttcttcGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAgccgcc                                        >  7:232111/1‑147 (MQ=255)                                                                                                                  cttcGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGAAGATGGCGGTGTACACAACGGTCAAAATCACACCTTTGGCCTGGATCCAGACTTGCATTCAGATGTCGGTGACGGTACCGAAGCCGCCCAgc                                    >  5:40466/1‑149 (MQ=255)                                                                                                                       ttGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCg                                >  3:346440/1‑148 (MQ=255)                                                                                                                         gACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGa                              <  8:12081/148‑1 (MQ=255)                                                                                                                           ccc‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGACGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGcc                           >  5:73326/1‑149 (MQ=255)                                                                                                                                  ggCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGaa                      >  2:291476/1‑148 (MQ=255)                                                                                                                                   gCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATa                    >  5:260926/1‑149 (MQ=255)                                                                                                                                          accacATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGt             >  2:126839/1‑149 (MQ=255)                                                                                                                                            cacaTCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATacac                 >  1:217454/1‑143 (MQ=255)                                                                                                                                            cacaTCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATacac                 <  2:217454/143‑1 (MQ=255)                                                                                                                                            cacaTCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGTGa           <  2:220431/149‑1 (MQ=255)                                                                                                                                            cacaTCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGTGa           <  6:260926/149‑1 (MQ=255)                                                                                                                                               aTCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGTGAGGa        <  3:105345/149‑1 (MQ=255)                                                                                                                                               aTCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACAACTTTTGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGTGAGGa        >  4:56577/1‑149 (MQ=255)                                                                                                                                                  cAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGTGAGGATgg     <  4:201051/149‑1 (MQ=255)                                                                                                                                                    gCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGTGAGGATGGCa   <  1:371597/149‑1 (MQ=255)                                                                                                                                                     cACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGTGAGGATGGCAc  <  6:148483/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCC‑GCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGACGATGGCGGTGTAGACCACGGTGAAAATCACACCTTTGGCCTGGATCCAGACTTGCATGCCGATGTCGGTGACGGTACCGAAGCCGCCCAGCGCCGGAGCCGCGAATACACCCGTGAGGATGGCAC  >  NC_002947/5966307‑5966602

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