Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,329,743	G→T	D261Y (GAT→TAT)	dacB →	peptidoglycan DD‑endopeptidase DacB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,329,743	0	G	T	85.7%	131.5 / 18.4	49	D261Y (GAT→TAT)	dacB	peptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base T (22/20);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

AAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCG  >  NC_000913/3329508‑3329984                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   aaaGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAgg                                                                                                                                                                                                                                     <  2:322523/251‑1 (MQ=255)   aGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCg                                                                                                                                                                                                                                   <  1:97692/251‑1 (MQ=255)   aGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCg                                                                                                                                                                                                                                   <  1:69771/251‑1 (MQ=255)            atatGGCTTTTATTCGCGTGGCATCTTATTCCCCCGTTACGAGTTTCAGCCAGGTACGGCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGGTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAcc                                                                                                                                                                                                                          <  1:296638/251‑1 (MQ=255)                            cgTGGCATCTGATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGGCGAAGCGCAATACTACGAACTGGATGGGGTGCCAGGCGACCTGAACCGATTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATCCGGGTGCATTTCTGAAATATGAGTTAAAACa                                                                                                                                                                                                                                       >  1:107085/1‑222 (MQ=255)                             gTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgct                                                                                                                                                                                                         <  1:113260/251‑1 (MQ=255)                               ggCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGCTCACCAGGCGCGGAACACtgctgc                                                                                                                                                                                                       >  2:194111/1‑251 (MQ=255)                                      tattaCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGa                                                                                                                                                                                                 >  1:9271/1‑250 (MQ=255)                                          aCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAg                                                                                                                                                                                           <  1:167759/251‑1 (MQ=255)                                           ccccGTTACGATGTTCAGCCCGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTCCTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGAAGCCTGCAACACCGTCCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATCCAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTCAAACAGGCGGGGATCAACCGGAGCGGGACActtctggccccgactcacg                                                                                                                                                                                           >  2:275379/1‑238 (MQ=255)                                               gTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTaa                                                                                                                                                                                      <  2:379741/251‑1 (MQ=255)                                               gTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAAAACTGCTGCGCCAGAC‑TCAGGTTa                                                                                                                                                                                       >  1:381200/1‑250 (MQ=255)                                                  aCGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGAAAATGCTGCGGCAGAA‑TCAGGTTAACg                                                                                                                                                                                    >  1:12006/1‑250 (MQ=255)                                                          cAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGaa                                                                                                                                                                           <  2:9271/251‑1 (MQ=255)                                                                  tACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACCTC‑CGTTAACGAACCTggcaagggtgtt                                                                                                                                                                   >  2:60005/1‑241 (MQ=255)                                                                   aCGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACAct                                                                                                                                                                                                            <  1:98459/210‑1 (MQ=255)                                                                   aCGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTg                                                                                                                                                                  <  1:455137/251‑1 (MQ=255)                                                                   aCGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCCGGAGCGGAACAct                                                                                                                                                                                                            >  2:98459/1‑210 (MQ=255)                                                                     gCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGGACCTGGaaagagataggg                                                                                                                                                                 >  1:442208/1‑240 (MQ=255)                                                                     gCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACg                                                                                                                                                                         <  2:442208/242‑1 (MQ=255)                                                                     gCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGcc                                                                                                                                                                >  1:257983/1‑251 (MQ=255)                                                                             ccccGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTCGAGCGGGAAACTGCTGCGCCAGAC‑CCAGGGTAACGAAACTGGAACGGTAGCTGCCAGCAAAc                                                                                                                                                         >  2:317262/1‑250 (MQ=255)                                                                                cGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTc                                                                                                                                                     <  1:416002/251‑1 (MQ=255)                                                                                     ttCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGGTAAAGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGccc                                                                                                                                                >  1:181621/1‑251 (MQ=255)                                                                                         gCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGcccc                                                                                                                                               <  1:429424/248‑1 (MQ=255)                                                                                                                 gATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGGCGGTATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACGCTGCCGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTAGGCCCCGCTGCACGATCTGCTTAAGATTATGCt                                                                                                                    <  1:274770/251‑1 (MQ=255)                                                                                                                   tgtgGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCATATGGAGCGGAAAACTGCTGAGCCAGAC‑TCAGGGTAACGAAACTGGAACGGGAGTTGCCAGCAAAAAGTCGGCCCCGCTGCCTGATCCGCTTAAGATTATTCTg                                                                                                                   >  1:89102/1‑250 (MQ=255)                                                                                                                   tgtgGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGcc                                                                                                                                                 <  1:29105/220‑1 (MQ=255)                                                                                                                   tgtgGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAA‑TCCGGGTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGcc                                                                                                                                                 >  2:29105/1‑220 (MQ=255)                                                                                                                         gCCAGGCGACCGGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACTTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAg                                                                                                             <  1:125069/250‑1 (MQ=255)                                                                                                                          ccAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCCCGATCTGCTTACGATTATTCTGAAAAAGTc                                                                                                           >  2:105681/1‑251 (MQ=255)                                                                                                                           cAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAAAC‑TCAGGGTAACGAACCTGGAACGGGAGTTGACAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAggcg                                                                                                          >  2:181099/1‑248 (MQ=255)                                                                                                                                       aCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGGTAAAGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATc                                                                                              >  1:406161/1‑251 (MQ=255)                                                                                                                                        cccatttACGCGGACGGGATGCCTGCCAAAACGTGCTGAGCCGCTTCCGTAGGCTTTTGCCGGGCAGGATGGAGCCAGCGATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCTCCTGGAGCGGAACACTGCTGCGCCAGAC‑TCCGCTGAACGACCCTGGAACTTTAGTTGCCAGTTACCAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTTGGACAACATGATCg                                                                                             <  2:321361/247‑1 (MQ=255)                                                                                                                                         cGCTTTACGCTGACGGGATGCCTGACACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAACGGAACAATCATGGGCCAGAA‑TCAGGTTAAAGCACCCGGAGCGGGCGTTGCCAGTAAACAGTCGGCCCCACTGCCCGCACTGATTCAAATTATGATGAAAAAATCGGccaaacagatcgc                                                                                            >  2:489722/1‑242 (MQ=255)                                                                                                                                          ggtttACGCTGACGTGGTGCCTGCCACAACGTTCTGAGCCGCTCCCGTAGGCCTTTGTCGTGCAGGATGGGGCCAGCTAGGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGATCGGAACACTGCTGTGCCAGAC‑TGAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCGCCGCTGCGCGAGGTGCTTAAGATTATGCTGAAAAAGTCGCGCAACATGATCGcc                                                                                           <  2:204789/249‑1 (MQ=255)                                                                                                                                                        ggATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGAt                                                                                               <  2:238713/233‑1 (MQ=255)                                                                                                                                                        ggATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGGCGGAAAACATGAt                                                                                               >  1:238713/1‑233 (MQ=255)                                                                                                                                                            gcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGTGGCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCag                                                                          <  2:19698/250‑2 (MQ=255)                                                                                                                                                                                 cGCTCCCGTTGGCTTTAGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGTCAGAC‑TCAGGTTAACGAACCTGGAACGGTACTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGAGACGGTTTTCCGCATGATTGGCCATGCGCGCTTCAAt                                                    <  1:124428/251‑1 (MQ=255)                                                                                                                                                                                            gCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCTCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGATTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAAc                                         <  1:275379/251‑1 (MQ=255)                                                                                                                                                                                                       gCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCg                              <  2:257983/251‑1 (MQ=255)                                                                                                                                                                                                                  ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGcc                   <  1:447788/251‑1 (MQ=255)                                                                                                                                                                                                                  ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAAAACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGGTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGGTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATggcggggcggggcgggcggc                   >  1:403628/1‑246 (MQ=255)                                                                                                                                                                                                                    aGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGCCAGAC‑TCAGGGTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTTAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGGCGGGGCGGGCGcccg                 >  2:228394/1‑249 (MQ=255)                                                                                                                                                                                                                       tATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCCCGATATGCTTCAGATTATGCTGAAAAAGTCGGACAAAATGATCGCCGATACGGGTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGGCGGGTCgggcccccttcg              >  1:40046/1‑241 (MQ=255)                                                                                                                                                                                                                          gCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCa           <  1:60005/251‑1 (MQ=255)                                                                                                                                                                                                                                 cAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGAACACTGCTGCGCCAGAA‑TCAGGGTAACGAAACTGGAACGGTAGTTGCCCGTAAAAAGTCCGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAAAATGGCGGGCCGGGTCCGGCGCCGTGCGGCACATcccc    >  2:282100/1‑249 (MQ=255)                                                                                                                                                                                                                                   aTTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGCACACTGCAGCGGCAGAC‑TCAGGGTCACGAACCTGGAACGGGAGTTGCCAGGAAACAGTCGGCCCCGCTGCACGATCTGCTTAATATTTTGCTGAAAAAGTCGGACAACATGATCGCCGATACGGGTTTCCGCATGATAGGCCATGCGCGCTTCAATGTTCCTGGAACATggcggggcggggcgggcacctcgcgtccgattccccg  >  2:405661/1‑229 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   AAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGAC‑TCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCG  >  NC_000913/3329508‑3329984

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑