Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,032,423	G→A	V219V (GTC→GTT)	yedJ ←	putative HD superfamily phosphohydrolase YedJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,032,423	0	G	A	97.4%	110.9 / ‑5.7	39	V219V (GTC→GTT)	yedJ	putative HD superfamily phosphohydrolase YedJ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (26/11);  minor base C (0/1);  total (27/12)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTG  >  NC_000913/2032177‑2032649                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                    cgcgCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAacac                                                                                                                                                                                                                                 >  2:208127/1‑251 (MQ=255)       ccAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCAt                                                                                                                                                                                                                          <  1:61022/251‑1 (MQ=255)          gTGATTCTCACCTACCCCATTAAGCTGCGCCACCAACGTGTTAACATCATAAATGTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAAGCGGGTACTGATATATGTTCCGGCATAGAGTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATTATTACAGCCATTTAGGGGCACGCGCTCAGCCAGCGGATGAAATCGCATCGATCGCCTGGTGATCAACACC‑GTCATtct                                                                                                                                                                                                                       <  2:170070/251‑1 (MQ=255)                  cACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGcc                                                                                                                                                                                                               >  1:17619/1‑251 (MQ=255)                            atgacgccactccaccaTCGTTTTAACACCATAAATTTGCAGCAGTGTCTCCATCTGTGCCTGTGCGGCATCCCCGTGGCTGTATGAATCGTTTACTTTTTTATTTCCCTGCATAGTTTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTCCAGCCATTTAGGGCCACGCCCCCAGCCACCGGATGAAACCGCATCTACCCCCTGGTGATAACCCCC‑TCCATTCCCCCCCGCCAGGTCGGCAc                                                                                                                                                                                                     <  2:408198/240‑1 (MQ=255)                              taatatGCGCCACCAACGTTTTAACATCATAAATTTGCAGCCGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGACTCGGTTACTGCTATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTCCCATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGCTGAAAACGCATCTATCACCTTGTGATCAACACA‑TTCATTCTCCCCCGCCAGtccgcaactg                                                                                                                                                                                                   >  2:251647/6‑242 (MQ=255)                               aaGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTg                                                                                                                                                                                                   >  1:263116/1‑250 (MQ=255)                                   tGCGCCACCAACGTTTTAACATCATAACTTTGCAGCAGTTTCTCCAGCCTTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCCCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCCCTGAGCtt                                                                                                                                                                                              >  2:126120/1‑251 (MQ=255)                                                    aaCATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGATCAGCCAGCGGATGATAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGACACCTAGCTTCGCCCTAAACTCGACt                                                                                                                                                                              >  2:310915/1‑250 (MQ=255)                                                                              ccAGCATTGCCTGTGCGGCATTCCCGGTGCGGTATGAATCGGTTACTGATATATTTGCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCa                                                                                                                                                     <  1:256652/249‑1 (MQ=255)                                                                                cgCATGGCCTGTGCGGCATTCCCGGGGCTGTATGAACCGGTTACTGATATATTTTCCGGAATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGAGGAAAACGCATCTATCGCCTGTTGATCAACCCC‑TCCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAgctg                                                                                                                                                 <  2:378349/250‑1 (MQ=255)                                                                                    ttGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGCTATATATTCCTGCATAGAATTCACCCGCCATCTGAGATGGCCGCCAGGTTACCATAATTACAGCCATTTAGCGCCACGCGCTCAGCCAGCCGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCCGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAACATGCGCGATGTGCTGCGCCAGCTGCTTg                                                                                                                                             >  1:217724/1‑251 (MQ=255)                                                                                          gtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCtt                                                                                                                                       >  1:416825/1‑251 (MQ=255)                                                                                                 aTTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGttt                                                                                                                                >  1:144781/1‑251 (MQ=255)                                                                                                   ccccGGTGCTGTATAAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCGGAGATGGCCGACAGGTTACCATAATTCCCGCCATTTAGGGCCACGCCCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACCCC‑TTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCTTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGttt                                                                                                                                <  2:186349/248‑1 (MQ=255)                                                                                                             gTATGAATCGGTTACTAATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACGGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTg                                                                                                                      <  1:262036/249‑1 (MQ=255)                                                                                                              tATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCgg                                                                                                                   >  2:8852/1‑251 (MQ=255)                                                                                                               aTGAATCGGTTACTGATATATTTTCATGTATAGCTTTCACCGGCACTCTGAGATGGCCCACAGGTTCACATAATTACAGCCATTTAGGGCCACGCGCCAAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCCGTTCGGCACTGCGCTTCGCCATAAACTCGACTACAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGt                                                                                                                  >  2:205079/1‑251 (MQ=255)                                                                                                                          aCTGATATATTTTACTGACTACATTTCAACGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATAACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTCAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCa                                                                                                        >  2:448643/1‑250 (MQ=255)                                                                                                                           cttatatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCAc                                                                                                                                                                                                     <  1:199593/153‑1 (MQ=255)                                                                                                                           cTGATATATTTTCCTGCATAGATTTCAACGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCAc                                                                                                                                                                                                     >  2:199593/1‑156 (MQ=255)                                                                                                                              atatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTAGCATGGTTTGCGGTAGCTTCAGCGACtt                                                                                                    >  2:67274/1‑250 (MQ=255)                                                                                                                              atatatTTTCCTGCATAGATTTCACCGCCAACCTGAGAGGGCCTCCAGGTAAACATATTTACAGCCATTTAGGGCCACGCGCCCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACCCC‑TCCTTTCTCCCCCGCCAGTTCGGCACTGAGCGTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCGGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACtt                                                                                                    <  2:416825/250‑1 (MQ=255)                                                                                                                                      tCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCGTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGa                                                                                            <  1:448643/250‑1 (MQ=255)                                                                                                                                                    ccccGGCCATCTGAGAGGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCTCTGTGTGATCACCACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGTTTTTCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGcgcg                                                                              <  1:182587/248‑1 (MQ=255)                                                                                                                                                        ggCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCg                                                                         >  1:143044/1‑251 (MQ=255)                                                                                                                                                        ggCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCg                                                                         >  1:278359/1‑251 (MQ=255)                                                                                                                                                             tCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTat                                                                    >  1:251824/1‑251 (MQ=255)                                                                                                                                                               tGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCACCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGCAATGATCCAGCGCGTAGCGTTTatc                                                                   >  1:453445/1‑250 (MQ=255)                                                                                                                                                                         gACAGGATAACATAATTCCAGGCATTTAGCGACACGAGCTCAGCCAGCAGCTGAAACAGCATCTATAACCATGGGAACAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACCCGACTAAGAACTGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTATGCATGGTTTGAGGTAGCCTCAGCGACTTGGCCTGGAAATGATCCTGCGTGGAGCGTTTATCATCAAGCGGGc                                                        >  2:211899/1‑251 (MQ=255)                                                                                                                                                                            aGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGAt                                                     >  1:370595/1‑251 (MQ=255)                                                                                                                                                                                  aCATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCtgtg                                               >  1:340935/1‑251 (MQ=255)                                                                                                                                                                                             gCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCt                                    >  1:176461/1‑251 (MQ=255)                                                                                                                                                                                                    aGGGCCAAGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCAt                             >  1:159500/1‑251 (MQ=255)                                                                                                                                                                                                     gggCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATc                            >  1:211390/1‑251 (MQ=255)                                                                                                                                                                                                      ggCCACGCGCTCAGCCAGCGGATGAAAACGCATCCATCACCTTGTGATCAACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATc                            >  1:210257/1‑250 (MQ=255)                                                                                                                                                                                                                     cagcGTATGAAAACGCATATACCACCTTGT‑ACCCACCCCCTTCTTTCCCCCCCTCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGccc            <  2:370595/251‑1 (MQ=255)                                                                                                                                                                                                                           aTGAAAACGCATCTATCACCATGTGATCAACGCA‑TTCAGTCACCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTGAAAAATGAGCGTTGAGCTGCGCCCGATGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGACTGGAAATGATCCAGCGCGTAGCGTTTGTCACCAAGCGGGCGATGCTTTGCGAACGGATCTTCACCATCAAACAGCGCAACGCCCAATGcc      >  2:323770/1‑251 (MQ=255)                                                                                                                                                                                                                               aaaCGCATCTATCACATTGTGATCAACACC‑TTCATTCTCCCCCGCCAGATCGGCACTGAGCTTCGCCATAAACTCGACTAAAAACTGCGCGTTGTGCTGCGCAAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTCTCATAAAGCGGGCGATGCTGTGCGAACGGCTCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTg  >  2:407081/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                    CGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACC‑TTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTG  >  NC_000913/2032177‑2032649

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