Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	443,910	T→C	5.9%	L24L (TTG→CTG)	bioD →	dethiobiotin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	443,910	0	T	C	5.9%	257.6 / 4.3	84	L24L (TTG→CTG)	bioD	dethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base T (40/39);  new base C (3/2);  total (43/41)
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.85e-01

CGCCAACCCGCGGGGCTGCCAGCCACCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCAGGTCAAT  >  NC_002947/443762‑444059                                                                                                                                                       |                                                                                                                                                      cGCCAACCCGCGGGGCTGCCAGCCACCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGtt                                                                                                                                                      <  2:229161/150‑1 (MQ=255) cGCCAACCCGCGGGGCTGCCAGCCACCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGtt                                                                                                                                                      <  1:1240726/150‑1 (MQ=255)  gCCAACCCGCGGGGCTGCCAGCCACCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGATCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTg                                                                                                                                                     >  2:1293500/1‑150 (MQ=255)                 tGCCAGCCACCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGtt                                                                                                                                                      <  2:103504/134‑1 (MQ=255)                 tGCCAGCCACCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGtt                                                                                                                                                      >  1:103504/1‑134 (MQ=255)                 tGCCAGCCACCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTg                                                                                                                                      >  2:23590/1‑150 (MQ=255)                  gccagccaCCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGc                                                                                                                                     >  1:1452149/1‑150 (MQ=255)                      gccaCCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAggg                                                                                                                                 <  2:78059/150‑1 (MQ=255)                         cccTCTCAAGTGGTCTATAGTGTGTTGAGCAAGCCACTGGCGTAACGGCAGAAAGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCGGGC‑‑CTGTTGCATGCGGCGCGTTTGCAggctcc                                                                                                                              <  1:1056565/149‑5 (MQ=255)                         aCCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCAt                                                                                                                              <  2:576082/150‑1 (MQ=255)                         aCCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCAt                                                                                                                              >  1:1570387/1‑150 (MQ=255)                           cTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGa                                                                                                                            <  2:1570387/150‑1 (MQ=255)                             aTCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGc                                                                                                                          >  2:83381/1‑150 (MQ=255)                              tCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCa                                                                                                                         >  2:1430944/1‑150 (MQ=255)                                 aGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGc                                                                                                                      >  2:1538072/1‑150 (MQ=255)                                   tGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTg                                                                                                                    <  1:1217203/150‑1 (MQ=255)                                    ggTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTgg                                                                                                                   <  2:1452149/150‑1 (MQ=255)                                      tCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGt                                                                                                                 >  2:127698/1‑150 (MQ=255)                                      tCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGt                                                                                                                 >  1:447266/1‑150 (MQ=255)                                        tATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGGACGCTGGGTGc                                                                                                               >  1:801378/1‑150 (MQ=255)                                        tATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCCCGCTGGGTGc                                                                                                               >  1:69460/1‑150 (MQ=255)                                            gtgtgtTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAg                                                                                                           >  2:1363747/1‑150 (MQ=255)                                               tgtTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAg                                                                                                                                   >  2:500780/1‑123 (MQ=255)                                               tgtTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAg                                                                                                                                   <  1:500780/123‑1 (MQ=255)                                               tgtTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCg                                                                                                        <  2:1448318/150‑1 (MQ=255)                                                gtTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGACGGCATGCCgg                                                                                                       >  1:887943/1‑150 (MQ=255)                                                  tGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGtt                                                                                                     <  2:515880/150‑1 (MQ=255)                                                     gcAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGcc                                                                                                  <  1:127698/150‑1 (MQ=255)                                                       aaGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTg                                                                                                                    <  1:800154/130‑1 (MQ=255)                                                       aaGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTg                                                                                                                    >  2:800154/1‑130 (MQ=255)                                                       aaGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTc                                                                                                >  1:1341732/1‑150 (MQ=255)                                                       aaGCCACTGGCGTAAGGGGAGCACGATCAGGGACGCCTTTTTCATTGCGGGGGCCGATGTGAATGTCGG‑‑ACACCGACAACGCCGGTGCCATCCACTGCGTGCTGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCTGTTGCCTc                                                                                                >  1:780481/1‑150 (MQ=255)                                                          attgcgcaggcctttcggcgtcatggtgcagcccgaggcaaccggcttggcacccagcgtgctcatgccctgcaaacgCGCCGCTGGC‑‑CTGTTGCATGCGGcgcg                                                                                                                                          <  2:1477879‑M1/27‑1 (MQ=255)                                                          attgcgcaggcctttcggcgtcatggtgcagcccgaggcaaccggcttggcacccagcgtgctcatgccctgcaaacgCGCCGCTGGC‑‑CTGTTGCATGCGGcgcg                                                                                                                                          >  1:1477879‑M1/79‑105 (MQ=255)                                                            aCTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCt                                                                                           >  2:1238418/1‑150 (MQ=255)                                                                  gTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTg                                                                                                                    >  1:974696/1‑119 (MQ=255)                                                                  gTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTg                                                                                                                    <  2:974696/119‑1 (MQ=255)                                                                  gTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCg                                                                                                        >  1:434115/1‑131 (MQ=255)                                                                  gTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCg                                                                                                        <  2:434115/131‑1 (MQ=255)                                                                  gTAAGCGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCa                                                                                     <  1:23590/150‑1 (MQ=255)                                                                   tAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCAt                                                                                    <  2:1058418/150‑1 (MQ=255)                                                                     ctttcggcgtcatggtgcagcccgaggcaaccggcttggcacccagcgtgctcatgccctgcaaacgCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGa                                                                                                                            <  1:725755‑M1/41‑1 (MQ=255)                                                                     ctttcggcgtcatggtgcagcccgaggcaaccggcttggcacccagcgtgctcatgccctgcaaacgCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGa                                                                                                                            >  2:725755‑M1/68‑108 (MQ=255)                                                                     aGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAggg                                                                                                                                 >  2:1236773/1‑103 (MQ=255)                                                                     aGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAggg                                                                                                                                 <  1:1236773/103‑1 (MQ=255)                                                                     aGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGa                                                                                  <  1:83381/150‑1 (MQ=255)                                                                     aGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGa                                                                                  >  2:977583/1‑150 (MQ=255)                                                                      ggggAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGAc                                                                                 <  1:1238418/150‑1 (MQ=255)                                                                       gggAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACg                                                                                >  2:1415905/1‑150 (MQ=255)                                                                       gggAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACg                                                                                <  2:899785/150‑1 (MQ=255)                                                                           gCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGa                                                                            >  2:265613/1‑150 (MQ=255)                                                                              cGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAg                                                                         <  2:69460/150‑1 (MQ=255)                                                                                aTGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGCTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGcc                                                                                                         <  2:981119/116‑1 (MQ=255)                                                                                aTGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGCTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGcc                                                                                                         >  1:981119/1‑116 (MQ=255)                                                                                   aGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTg                                                                    <  2:1341732/150‑1 (MQ=255)                                                                                         gCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGCTGCATgcggc                                                                                                                                             >  2:424598/1‑71 (MQ=255)                                                                                         gCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGCTGCATgcggc                                                                                                                                             <  1:424598/71‑1 (MQ=255)                                                                                           cTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAg                                                            <  1:726053/150‑1 (MQ=255)                                                                                            tttttCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGc                                                           >  1:783784/1‑150 (MQ=255)                                                                                             ttttCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCg                                                          >  2:926639/1‑150 (MQ=255)                                                                                                 cATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGc                                                      <  2:768300/150‑1 (MQ=255)                                                                                                        ggTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGgcg                                               >  2:20978/1‑150 (MQ=255)                                                                                                          taccgataccgaTGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCt                                             >  2:1572490/1‑150 (MQ=255)                                                                                                            ccgataccgatGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGa                                           <  1:1415905/150‑1 (MQ=255)                                                                                                                 accgatGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATcgac                                      <  2:1019210/150‑1 (MQ=255)                                                                                                                 accgatGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATcgac                                      <  1:1363747/150‑1 (MQ=255)                                                                                                                   cgatGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATcgacga                                    >  1:1579456/1‑150 (MQ=255)                                                                                                                     atGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGaaa                                  <  1:1293500/150‑1 (MQ=255)                                                                                                                      tGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAg                                 >  2:1116230/1‑150 (MQ=255)                                                                                                                        tCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCa                               <  2:1579456/150‑1 (MQ=255)                                                                                                                         cGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCAc                              <  1:1430944/150‑1 (MQ=255)                                                                                                                               gACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCaa                        >  1:1102737/1‑150 (MQ=255)                                                                                                                               gACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCaa                        <  1:265613/150‑1 (MQ=255)                                                                                                                                       tCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCGCGGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCtt                >  1:500050/1‑150 (MQ=255)                                                                                                                                       tCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCtt                >  2:334819/1‑150 (MQ=255)                                                                                                                                           CGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGa            >  2:780481‑M1/1‑150 (MQ=255)                                                                                                                                            GCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATa                                                             <  2:314912‑M1/100‑1 (MQ=255)                                                                                                                                            GCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATa                                                             >  1:314912‑M1/1‑100 (MQ=255)                                                                                                                                             CTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGc          <  2:916480‑M1/150‑1 (MQ=255)                                                                                                                                              tGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACTAAAGCACgg                            >  1:1208889/1‑131 (MQ=255)                                                                                                                                              tGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACgg                            <  2:1208889/131‑1 (MQ=255)                                                                                                                                              tGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCa         <  1:20978/150‑1 (MQ=255)                                                                                                                                              tGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCa         >  2:604205/1‑150 (MQ=255)                                                                                                                                              tGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCa         <  1:926639/150‑1 (MQ=255)                                                                                                                                               ggC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCAg        <  1:916480/150‑1 (MQ=255)                                                                                                                                               ggC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCAg        <  1:1116230/150‑1 (MQ=255)                                                                                                                                                    cTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCAGGTc     >  1:554322/1‑150 (MQ=255)                                                                                                                                                      gTTGGATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCAGGTCaa   >  1:612912/1‑150 (MQ=255)                                                                                                                                                       ttGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCAGGTCAAt  >  1:1037830/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CGCCAACCCGCGGGGCTGCCAGCCACCTATCAAGTGGTCTATGGTGTGTTGCGCAAGCCACTGGCGTAAGGGGAGCACGATGAGCCAGGCCTTTTTCATTGCGGGTACCGATACCGATGTCGGCAAGACCACCATCGCCGCTGGC‑‑CTGTTGCATGCGGCGCGTTTGCAGGGCATGAGCACGCTGGGTGCCAAGCCGGTTGCCTCGGGCTGCACCATGACGCCGAAAGGCCTGCGCAATAGCGATGCCCTGGCGCTGATCGACGAAAGCACGGTCAAGCTGCCTTATGAGCAGGTCAAT  >  NC_002947/443762‑444059

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