Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,080,920	A→T	6.1%	intergenic (‑144/+19)	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,920	0	A	T	6.1%	199.4 / 2.9	66	intergenic (‑144/+19)	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 = 8.54e-01

ACACAGCGTGCAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGC  >  NC_002947/4080773‑4081068                                                                                                                                                      |                                                                                                                                                     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                                                                                                              >  1:95385/1‑139 (MQ=255)                                                   aGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAgg                                                                                                              <  2:95385/139‑1 (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                                                                                            <  7:267702/149‑1 (MQ=255)                                                           cGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                            <  3:24591/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                                <  2:137228/133‑1 (MQ=255)                                                                                                                                        ctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGc                                >  1:137228/1‑133 (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)                                                                                                                                                      aaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGc  <  2:273506/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ACACAGCGTGCAGCTCCTTGGATCAGGACGGT‑ACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTT‑GTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGC  >  NC_002947/4080773‑4081068

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