Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,823,288	C→T	G161D (GGC→GAC)	recA ←	DNA recombination/repair protein RecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,823,288	0	C	T	97.7%	144.1 / ‑6.1	44	G161D (GGC→GAC)	recA	DNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/21);  minor base G (0/1);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TTCGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG  >  NC_000913/2823044‑2823525                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               ttCGCTTCCCACCCCGTTTTCGCCCTCTTTCACCGCTCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCTCCGGTTGGGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTTATGAAGAACAGCAGCGGGGTGGTCTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGGGAGAGTCGGCGAtt                                                                                                                                                                                                                                           <  2:203962/250‑1 (MQ=255)                                     gcCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCa                                                                                                                                                                                                      <  1:191001/251‑1 (MQ=255)                                             gacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑c                                                                                                                                                                                              >  2:299799/1‑250 (MQ=255)                                              acgGATGTCGAGACGAACAGAGGCGTAAAATTTCAGCCCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCCCCCCAATTTTCATACGGATCTGGTTGAGGAAGATCAGCGGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGGCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCTTCAGTGCCGCCACGGAGT‑Caa                                                                                                                                                                                            <  2:43563/251‑1 (MQ=255)                                                  aTGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAacgac                                                                                                                                                                                        <  1:181118/251‑1 (MQ=255)                                                      cGAGACGAACAGAGGCGTAGAATGTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCAAACATCACACCAATTTTCATACGGTTCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATa                                                                                                                                                                                    <  2:181045/251‑1 (MQ=255)                                                                     cGTAGAATTTCAGCGCGTTACCACCGGTCGTGTTTTCCGGGTTACCGCACATCACACCAATTTTCCTACGGATCTGGTTGATGAAGATCAGCAGCGTGTTAGACTGCTTCAGGTTACCCGCCAGCTTCCGCATCGCCTGGCTCATCATACGTGCCGCACGGCCCATGTGCGCGTCTTCGCTTTCGCCTTCGCTTTCAGCTTTCGGCGTCCGTGCCGCCACGGCGT‑CAccgaccgtacctcctcctgcgcct                                                                                                                                                                     >  2:13404/1‑232 (MQ=255)                                                                             ttCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAAcgcg                                                                                                                                                              <  1:242541/250‑1 (MQ=255)                                                                              taaGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGGGTTGGACTGCTTCAGGTTACCCGCCAGCTTAAGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGCTTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGCAcgcgtc                                                                                                                                                            >  2:201/3‑249 (MQ=255)                                                                                aGCGCCTTACCACCCGTAGCGGTTTACGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGcc                                                                                                                                                            >  2:129458/1‑249 (MQ=255)                                                                                              ggTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAttt                                                                                                                                            >  1:93451/1‑251 (MQ=255)                                                                                                tAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTcc                                                                                                                                          >  1:289029/1‑251 (MQ=255)                                                                                                          cGGGTTACCGACCATCACACCATGTTTCATACGGATTTGGTGGATGAAGATCAGCAGCGTGTGGGACGGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCt                                                                                                                                <  2:65725/251‑1 (MQ=255)                                                                                                          cGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGCAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGATTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCt                                                                                                                                >  1:35741/1‑251 (MQ=255)                                                                                                           gggTTACCGACCCTCCCACCACTTTTCATACGGATCTGGTTGATGAAGCTCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCCGCTTACGCATCGCCTGGCTCATCATACGTGCCGCACGGCCCATGTGAGCGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCAGCCCCGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTc                                                                                                                               >  2:259006/1‑251 (MQ=255)                                                                                                           gggTTACCAAACATCCCACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACGGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTc                                                                                                                               <  2:254898/251‑1 (MQ=255)                                                                                                                aCCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCt                                                                                                                                >  1:207625/1‑245 (MQ=255)                                                                                                                aCCGAACACCACACCAATTTTCAGACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGCCTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCt                                                                                                                                <  2:207625/245‑1 (MQ=255)                                                                                                                       aTCACACCAATTTTCATACGGATCTGGTTGATGACGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCgg                                                                                                                   >  2:200248/1‑251 (MQ=255)                                                                                                                          acacCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTg                                                                                                                >  1:142554/1‑251 (MQ=255)                                                                                                                                aTTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGca                                                                                                          <  1:200248/251‑1 (MQ=255)                                                                                                                                aTTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGGGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGca                                                                                                          <  1:14688/251‑1 (MQ=255)                                                                                                                                  tttCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGcac                                                                                                         <  1:209527/250‑1 (MQ=255)                                                                                                                                      aTACGGATCTGGTTGATGAAGAACAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCAcagcag                                                                                                    <  2:289029/251‑1 (MQ=255)                                                                                                                                          ggATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTg                                                                                                <  1:25942/251‑1 (MQ=255)                                                                                                                                                     atgaAGCTCAGCCGCGTGTTGGACCGCTTCAGGTAACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCCTGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGAc                                                                                     >  2:54827/1‑251 (MQ=255)                                                                                                                                                          gATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCTGGCTGGGAGCACAGCAGGTTGTCGATATCGACGccc                                                                                 >  1:165042/1‑250 (MQ=255)                                                                                                                                                          gATCAGCAGCGTGTTGGACGGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCa                                                                                <  2:168139/251‑1 (MQ=255)                                                                                                                                                             gcgcagcggttgtgACTGTTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTTATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGtt                                                                             <  2:166745/238‑1 (MQ=255)                                                                                                                                                                cagcGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTAc                                                                          >  1:98566/1‑251 (MQ=255)                                                                                                                                                                 agcGTGTTGGACTGCTTCAGGTTCCCCTCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTAca                                                                         <  1:104543/251‑2 (MQ=255)                                                                                                                                                                 agcGTGTTGGACTGCGTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACg                                                                         <  1:104528/251‑1 (MQ=255)                                                                                                                                                                                  cAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTcc                                                         >  2:239698/1‑250 (MQ=255)                                                                                                                                                                                      ttACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAgcgc                                                    >  2:9682/1‑251 (MQ=255)                                                                                                                                                                                         cccGCCAGCTTACGCATCGCCTGGCTCACAATACGCACCACACGGCCCATGTGACAGTAGTCAACTTCGCCATCAACTTCCGATTTCGGCGCCAGCACAGCAACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGCACGCGCCAGGACGTAACAGATTTCAAGAGCCTGCGAGCCGGTGACCGGCTGGAAGCACCGCAGGGTGTCGATATCGACGCCCAGTTTACGTGCGTAGGTTTGGACCAGCGCgtg                                                 >  2:267538/1‑251 (MQ=255)                                                                                                                                                                                                 cTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCac                                         >  2:171385/1‑250 (MQ=255)                                                                                                                                                                                                               gtcTCATCATAAGGGCCGCCAGGCCCATGTGAGAGTCGACGATTTCGACGTCGATTTCCGCTTTCGACGGCAGTTCCGCCACGGAGT‑CAACGACGATACCGTCTACTGCGCCAGTTCGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGAACAGCAGGTGTTCGATATCGACGTCCAGTTTTCGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATaatcgaaa                             <  2:253981/247‑7 (MQ=255)                                                                                                                                                                                                                   catcatACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGt                         <  1:92903/249‑1 (MQ=255)                                                                                                                                                                                                                         aCGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCtt                 >  1:166936/1‑251 (MQ=255)                                                                                                                                                                                                                         aCGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCt                  <  1:220180/250‑1 (MQ=255)                                                                                                                                                                                                                          cGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGATCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTc                >  1:262064/1‑251 (MQ=255)                                                                                                                                                                                                                             gCCGCAAAGCCCATGTTAGAGTTGTTGAGTTCGCCTGCGATTGCCGCTTTCGGCGTCAGTGCCGCCACGG‑GTGCAACGACGATAACGGCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGACTGTTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTGGTCGATATCGACGCCCAGGTTACGTGCGTAGATTGGGTCCAGCGCGTGTGCAGCATCGATAAACGCACAGGTTTTACCTTctc              <  2:2540/250‑3 (MQ=255)                                                                                                                                                                                                                                      cccATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACgctgcgctgc    >  1:143258/1‑251 (MQ=255)                                                                                                                                                                                                                                        cATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTCTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAGTTCCAGTGCCGGCTCGCCGGTGTCCGGCTGGGAGTACAGCAGGTTGTCGATATCTACGCTCAGTTTACGTGCGTATATTGGGGCCAGCGCGTGTGCAGCATCTATAAACGCACATGTTTTACCTTCTCGCTGCGTTgcgg  <  2:290709/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               TTCGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGT‑CAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGG  >  NC_000913/2823044‑2823525

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