Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,080,919	A→T	6.1%	intergenic (‑143/+20)	amaC ← / ← dpkA	D‑lysine aminotransferase/delta 1‑piperideine‑2‑carboxylate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,080,919	0	A	T	6.1%	198.8 / 2.9	66	intergenic (‑143/+20)	amaC/dpkA	D‑lysine aminotransferase/delta 1‑piperideine‑2‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base T (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

AACACAGCGTGCAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGC  >  NC_002947/4080772‑4081062                                                                                                                                                      |                                                                                                                                                aaCACAGCGTGCAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTa                                                                                                                                                 >  4:24591/1‑148 (MQ=255)  acacAGCGTGCAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTaaa                                                                                                                                               <  5:298532/149‑1 (MQ=255)   cacaGCGTGCAGCTCCTTGGATCAGGACGGTGACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTaa                                                                                                                                                >  1:196407/1‑148 (MQ=255)   cacaGCGTGCAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAc                                                                                                                                              >  4:305538/1‑149 (MQ=255)      aGCGTGCAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCcg                                                                                                                                           >  1:173057/1‑149 (MQ=255)         gTGCAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCcgcg                                                                                                                                         >  4:35262/1‑148 (MQ=255)            cAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAgg                                                                                                                                     >  2:252453/1‑149 (MQ=255)             aGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                    >  1:216157/1‑149 (MQ=255)             aGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                    <  3:305538/149‑1 (MQ=255)             aGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAggg                                                                                                                                    <  2:173057/149‑1 (MQ=255)                  cTTGGATCAGTACGGT‑ACGGAAAGTCAGTGTTGGGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGc                                                                                                                                  <  8:241547/146‑1 (MQ=255)                  cTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTTTACCCGCGaaaaa                                                                                                                                    <  4:271084/144‑4 (MQ=255)                  cTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTTTACCCGCGaaaaa                                                                                                                                    >  3:271084/1‑141 (MQ=255)                   ttGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCaaa                                                                                                                              >  2:186999/1‑149 (MQ=255)                   ttGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCaaa                                                                                                                              >  3:170786/1‑149 (MQ=255)                      gATCAGGACGGA‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCa                                                                                                                           <  3:43999/149‑1 (MQ=255)                           ggACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCaa                                                                                                                      <  2:42759/149‑1 (MQ=255)                                       aaaGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTC‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTc                                                                                                           <  6:187681/149‑1 (MQ=255)                                           tCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGc                                                                                                       <  1:186999/149‑1 (MQ=255)                                            cAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGTGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGc                                                                                                       <  1:252453/148‑1 (MQ=255)                                             aGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGACTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCt                                                                                                     >  6:150036/1‑149 (MQ=255)                                                   gAGCATACCCGGAGTGCGGGTTCAGGTACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAAc                                                                                               >  5:229546/1‑149 (MQ=255)                                                   gAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTc                                                                                                  <  1:35526/146‑1 (MQ=255)                                                   gAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTc                                                                                                  >  2:35526/1‑146 (MQ=255)                                                    aGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAgg                                                                                                        <  2:95385/139‑1 (MQ=255)                                                    aGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAgg                                                                                                        >  1:95385/1‑139 (MQ=255)                                                     gCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                             >  5:29726/1‑149 (MQ=255)                                                     gCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                             >  8:132685/1‑149 (MQ=255)                                                        tACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTg                                                                                          <  6:29726/149‑1 (MQ=255)                                                         aCCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                         >  5:65897/1‑133 (MQ=255)                                                         aCCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                         <  6:65897/133‑1 (MQ=255)                                                          cccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGTGTAAACCCGCGAAGGGGACAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCt                                                                                        <  6:164861/149‑1 (MQ=255)                                                            cGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                      <  3:24591/149‑1 (MQ=255)                                                            cGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                      <  7:267702/149‑1 (MQ=255)                                                                   cGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                         <  8:157126/123‑1 (MQ=255)                                                                   cGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                         >  7:157126/1‑123 (MQ=255)                                                                   cGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCAcct                                                                               <  1:30724/149‑2 (MQ=255)                                                                          aGGTACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAcccc                                                                        >  6:17676/1‑149 (MQ=255)                                                                          aGGAACATCCACTATTGTGCAAGGTTAGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAccc                                                                         >  1:2951/1‑149 (MQ=255)                                                                           ggAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAccccc                                                                       <  5:213520/149‑1 (MQ=255)                                                                            gAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCTACTCCTGCTCGGTCACCGCCACCCcct                                                                      >  4:177661/1‑149 (MQ=255)                                                                            gAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcct                                                                      >  1:273506/1‑149 (MQ=255)                                                                                   caatATTGAGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGCGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCACATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAACCCctcctcctc                                                               >  4:189592/4‑149 (MQ=255)                                                                                      tATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAcc                                                                          >  2:81890/1‑135 (MQ=255)                                                                                       aTTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGcc                                                           <  3:35262/149‑1 (MQ=255)                                                                                           tgCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCAcc                                                        <  7:132685/148‑1 (MQ=255)                                                                                           tgCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCt                                                       >  7:7695/1‑149 (MQ=255)                                                                                           tgCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCt                                                       >  5:330334/1‑149 (MQ=255)                                                                                             cAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTc                                                      <  8:154527/148‑1 (MQ=255)                                                                                              aaGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgc                                                    <  5:17676/149‑1 (MQ=255)                                                                                              aaGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACATcgc                                                    >  3:253660/1‑149 (MQ=255)                                                                                                      tCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCa                                                                                                          <  7:181494/88‑1 (MQ=255)                                                                                                      tCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCa                                                                                                          >  8:181494/1‑88 (MQ=255)                                                                                                         aTGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACgg                                            >  3:267188/1‑147 (MQ=255)                                                                                                            cAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGAGGGTAAACCCGCGAAGGGGCCAAATCATCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCAACGCCACACCCTCCTACTCGGCCACCTCTCGCTCACGGTAgcg                                       >  2:188173/1‑149 (MQ=255)                                                                                                                aTTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGcgccgc                                   >  7:232775/1‑149 (MQ=255)                                                                                                                   gCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTc                                 <  8:232775/148‑1 (MQ=255)                                                                                                                           gggATGAGGGCTGACTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCAt                        <  5:150036/149‑1 (MQ=255)                                                                                                                             gATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgc                      <  6:229546/149‑1 (MQ=255)                                                                                                                                gAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcgc                    <  4:267188/148‑1 (MQ=255)                                                                                                                                  gggCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctcctcc                                                                 >  2:270218/1‑101 (MQ=255)                                                                                                                                  gggCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctcctcc                                                                 <  1:270218/101‑1 (MQ=255)                                                                                                                                        gctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGcc             <  2:2951/147‑1 (MQ=255)                                                                                                                                         ctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGc                          >  1:137228/1‑133 (MQ=255)                                                                                                                                         ctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGc                          <  2:137228/133‑1 (MQ=255)                                                                                                                                          tttttCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTATCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCAc         >  2:81012/1‑149 (MQ=255)                                                                                                                                          cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGc                               >  4:96165/4‑127 (MQ=255)                                                                                                                                          cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGc                               <  3:96165/124‑1 (MQ=255)                                                                                                                                                  gggTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGc  >  8:219864/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                AACACAGCGTGCAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGC  >  NC_002947/4080772‑4081062

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