Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,978,166	T→A	26.2%	Y245F (TAT→TTT)	lysA ←	diaminopimelate decarboxylase, PLP‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,978,166	0	T	A	26.2%	72.5 / 5.7	44	Y245F (TAT→TTT)	lysA	diaminopimelate decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base T (21/11);  new base A (2/10);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACCAGAT  >  NC_000913/2977937‑2978404                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               ggCCCGGGCTCTGATAGGGTGGTAACCCACGTACCTTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCCCAAAGTGGCGGCCCCCCATTTGTTTGACGCTCCGCCCCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGtt                                                                                                                                                                                                                                       <  2:930691/238‑1 (MQ=255)           attattttgttgtAACCACCGTATATTCGCGGGCCCTCCAGTTCGTTGAACCCGGCCACCACCAGCCAAAGGGGGCGGCCCCCCCTTTGTTTGACGCTCCGCCACGGAGT‑AATTAATACGCCAGACGGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACGGGGTGGCCAAAATGGCGGGCTATTTGCTCACCCGCGCCATCCCACAGAC‑CATAATAATGTTCGGTATCAAc                                                                                                                                                                                                                             <  2:472045/228‑1 (MQ=255)            ttattgtggggtAACCACCGGTCATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTc                                                                                                                                                                                                                                      <  2:1208501/220‑1 (MQ=255)                   gtggtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGGTTTGCTCCCGCGCGGCATTCCACAGAC‑CCTAATAAAGTTCGGTATCAACCGCccct                                                                                                                                                                                                                      >  1:463565/1‑238 (MQ=255)                    tgtaacacaCCGTACATTGCGGGGCCCAGCAGTTCTTTGACCCCGCCCACACCCAACAAAAGGGGGCGTCTCCCATTTTGTTTGACGCCCCGCCCCTGATT‑AATTAATACGCAAGACTGTGCTACCAGGAAGGGCCCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCAAAAGGGCGGGCTATTTGCTCCCGCGCGGCATTCCACAGAC‑CAAAAAAATGTTGGGTACCAACCGCCTCtt                                                                                                                                                                                                                     <  2:1141459/230‑1 (MQ=255)                      gtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGGGATTTGCTCACACCCGGCATTCCCCAGAC‑ACAAATAAAGTTCGGGTTCAACCGctcgtcctt                                                                                                                                                                                                                  >  1:1208501/1‑231 (MQ=255)                      gtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAAGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAAAGTTCGGTATCAACCGCCTCTTCa                                                                                                                                                                                                                   >  1:1060647/1‑238 (MQ=255)                         aCTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCCCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTTCCAGGGTGGCCCAAATGGCGGGCGATTTGCTCTCGCGCGGCATTCCCCAGAC‑CCTAATAATGTTCGGGATCAACCGCCCCCTCAc                                                                                                                                                                                                                  >  1:234078/1‑236 (MQ=255)                             ccGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCCTTTGTTTGACGCTCCCCCCCCGAGA‑AATTAATACGCCCGACTGCGCTACCCGGAAGCGGCCCCG‑GTCCATTT‑CCCATTTTACAGGGGGGCCCAAAAAGGGGGCGGCTTTCTCTCCCGGCGCCTTTCACAAAA‑CAAAATAAAATTCTGTGTCAAACCCCTCTTctttaaaac                                                                                                                                                                                                            >  2:675618/1‑230 (MQ=255)                                                         tGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑GTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCCCAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCCGTTGATAAGGAACAGAAAGCCCAccgcccccc                                                                                                                                                                               >  1:186774/1‑236 (MQ=255)                                                              ccGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCCCCCCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCCG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCCAAAGGCGGGCGATTTGCTCACGCGCGGCGTTTCCCAGAC‑CATAATAATGTTCGGTGTCCACCCCCCCTCCACCCTCTTGTTAAGGAAAAAAAAACACACCGCCCCCCGaaaa                                                                                                                                                                          >  2:1219285/1‑238 (MQ=255)                                                                                                      tgACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAACCCACCGCCCCCAGAAATAGCCTGTAAATCCTGACCGAAATCGATGACCTGACGAAcc                                                                                                                                  >  2:1208694/1‑239 (MQ=255)                                                                                                          gCTCCGACCCTGAGT‑ATTTAATACGCCAGACTGCGCTACCAGGAGGCGACCCGG‑TTCAATTT‑CCAGTTCCACAGGGTGGCCCAAAGGGCGGGCGTTTTGCTCACGCGCGGCATCCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg                                                                                                                              <  1:1219285/239‑1 (MQ=255)                                                                                                          ccTCCGCACC‑GAGTAAATTAATCCGCAAGACTGCGCTCCAAGGAGGCGCCCCGG‑TTCAATTC‑CCAGTTCCACGGGGTGGCCAAAAGGGCGGGCAATTTGCTCACGCGCGGCATTCCACAGCC‑AATAAAAATGTCGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg                                                                                                                              <  1:1208694/238‑1 (MQ=255)                                                                                                                      aGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCt                                                                                                                  >  1:381135/1‑239 (MQ=255)                                                                                                                      aGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATCCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGACAGCCCACCGCCCGCAGAAATAGCATGTAAATCCTGACCGAATTCGAGGACCTGACGCACCATAGcacctgaccctt                                                                                                                  >  1:1336188/1‑231 (MQ=255)                                                                                                                               agccGCAAGGCTGCGCAACCAGGAAGCCACCCGG‑TTAAATTC‑CCAGTTTCCCGGGGGGGCCAAAATGGCGGGCAATTTGCTCACGCGCGGAATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGa                                                                                                          <  1:1008958/236‑1 (MQ=255)                                                                                                                                 aCGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CCTAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTAGATGAACTGACGCCCCCTACCACCACCCACCCGTTCCaacga                                                                                                        >  2:1008958/1‑235 (MQ=255)                                                                                                                                     caagctgccttACCAGGAAGGCACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCAAAAGGGCGGGCAATTTGCCCACGCGCGGCATTCCACAGAC‑AATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCa                                                                                                    <  1:1179127/229‑1 (MQ=255)                                                                                                                                            cgcTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCCGTTCCAGATGGGCATAATCa                                                                                              >  1:209642/1‑238 (MQ=255)                                                                                                                                            cgcTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAAAAGCCTGTAAATCCTCACCGAATTCGATGACCCGACCCACCACAACCACACCACCCTGTTCCCGCTGGGCATAAtaa                                                                                              >  2:938049/1‑236 (MQ=255)                                                                                                                                             gcTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACCCCTGTTCCAAATGGGCATAATCaa                                                                                             >  1:1270711/1‑238 (MQ=255)                                                                                                                                                 ccgggagggcccccGT‑TTCAATTT‑CGAGTTTCCAGGGGGGGCCCAAAGGGGGGGGTGTTTGCGCCCGGGCGGATTTCAAAAGAC‑AAAAATAATGTGCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCa                                                                                        <  2:755384/228‑1 (MQ=255)                                                                                                                                                 ccAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACCCGCGGCATTCCACCGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAAAAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGCCCGAATTCGCTGACATGACCGCCCAAAGCACCACCCACCCGTTCCAAGTGGGCCAAAACAACCCCa                                                                                        >  2:1307120/1‑239 (MQ=255)                                                                                                                                                      aaGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTACATCCTGACCGAATTCGATGACCTTACGCACCATAGCACCACACCCCCGTCCCAACTGGGCTTAATCAACGCCAAAAcc                                                                                   >  1:128334/1‑239 (MQ=255)                                                                                                                                                         gcgcccGG‑TTTTATTTCCCATTTTCAGGGGGGGCCCAAAATGGCGGGGAATTTCCCCCCGCGCGGCATCCCACAGAC‑CAAAATAATGTCCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                                  <  2:381135/235‑1 (MQ=255)                                                                                                                                                         gcccccGG‑TTAATTTT‑CCATTTCCCAAGGGGGCCCCAAAGGGGGGGGTTTTTGCCCACCCGCGGCATCCCACAACC‑AAAAAAAAGGTGCGGCATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                                  <  2:1270711/234‑1 (MQ=255)                                                                                                                                                         agACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTTCACAGAC‑CATAATAATGTTCGGTATCCACCGCCTCTTTACCCTGTTGATAGGGAACAGAAAGCCCCCCGCCCGCAGAAATAGCCAGTAAATCCTGACCGACTTCGCTGCCCTTACCCACCATTGCACCACCCACCTTTTCCAACTGGGGTTATAAAACGCCCGAACCacc                                                                                >  2:1036991/2‑237 (MQ=255)                                                                                                                                                          gACCCGGTTTCAATTT‑CCAGTTTAACAGGGTGGCCAAAATGGCGGGCGATTTGCTCACGCGGGGCATTCC‑CAGGCACATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                                  <  1:938049/237‑1 (MQ=255)                                                                                                                                                           aCCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCCTAGCACCCCACACCTGTTCCAGATGGGCATAAt                                                                                                >  1:16505/1‑221 (MQ=255)                                                                                                                                                              cGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCCCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCCTAGCACCCCACACCGGTTCCAGCTGGGCATAATAAACGCCAGAACCAAAGTGCa                                                                           >  1:552063/1‑239 (MQ=255)                                                                                                                                                                        ttcccatgtccacaGGGTGGCCAAAAGTGCTGGCGATTTTCTCACGCCACGCATCCCACAGAC‑ACAAATGATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATg                                                                   <  1:672938/229‑1 (MQ=255)                                                                                                                                                                                           gCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCacacac                                                                                                                    >  1:272719/1‑171 (MQ=255)                                                                                                                                                                                                                  cACGCGCGGCATGCCACAGCC‑CAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTgcg                         <  1:305909/239‑1 (MQ=255)                                                                                                                                                                                                                  cACGCGCGGCATGCCAAAGAC‑CATAATAAGTTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTgcg                         <  1:997651/239‑1 (MQ=255)                                                                                                                                                                                                                            ctaccACAGAC‑ATAAAGAATGTTCGGTACAAACCGCCTCTTCACCCTGTTGAAAAGGAACAGAAAGCCCCCCGCCCGAAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt               <  2:683485/236‑1 (MQ=255)                                                                                                                                                                                                                            catCCACAGAC‑CATAAAAATGTCCGGTAACAACCGCCTCTTCACCCTGTTGAAAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt               <  1:113186/237‑1 (MQ=255)                                                                                                                                                                                                                                cccaGAC‑AAAAATATTGTGCGGTATCACCCCCCTCTTCCCCCTTTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  2:16508/237‑1 (MQ=255)                                                                                                                                                                                                                                ccaaGAC‑AAAAATAATGTTCGTTAACAACCGCCTTTTCACCCTGGTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  1:1052968/236‑1 (MQ=255)                                                                                                                                                                                                                                agacagc‑cATAATAATGTTCGGTACCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  1:912065/233‑1 (MQ=255)                                                                                                                                                                                                                                acaaGAC‑AAAAATAATTTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt           <  1:750115/236‑1 (MQ=255)                                                                                                                                                                                                                                     aC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCACTCTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGGGCGGCGGGAACATCGGGGTAcc      >  1:1354119/1‑239 (MQ=255)                                                                                                                                                                                                                                     aC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGGTCCAGATGGGCATAATCAACGCCAGAACCAATTGGCATGTGAATCCCGACCAGCTGCAGCGGATTACGTTGTATAACGCCCAGGGAGGGGGGCGGagagggggagc      >  2:96286/1‑229 (MQ=255)                                                                                                                                                                                                                                        cATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCAAACACCTGTTCCCATTGGGCATAATCAACGCACGAACCAATTTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCCGATCGGGGTACCaa    >  2:1060460/1‑238 (MQ=255)                                                                                                                                                                                                                                          taataaTGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCGGTGCGGCGGGAAGATCGGTGTACCaccg  >  1:683485/1‑236 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               TGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGT‑AATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGG‑TTCAATTT‑CCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGAC‑CATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACCAGAT  >  NC_000913/2977937‑2978404

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