Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	674,259	G→A	C175C (TGC→TGT)	leuS ←	leucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	674,259	0	G	A	98.0%	155.9 / NA	50	C175C (TGC→TGT)	leuS	leucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/22);  minor base T (0/1);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTACAACCCATAAAGGTGTCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAGCAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGT  >  NC_000913/674014‑674503                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     gtaCAACCCATAAAGGTGTCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCCGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGcc                                                                                                                                                                                                                                                 >  1:618820/1‑251 (MQ=255)    cAACCCATAAAGGTGTCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTc                                                                                                                                                                                                                                              >  1:363080/1‑251 (MQ=255)                gtgtCCTGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTgttc                                                                                                                                                                                                                                  <  1:16395/251‑1 (MQ=255)                gtggccGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTGCCGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCATCTCGTCAGCGTAAGCAGTGATTTTGATAAACCTCGTCGTGATCTTTTTACGTTCAACTTTGGTATCGCAGAGCCATCAACAGGCGTCGATAACTTgttc                                                                                                                                                                                                                                  <  2:247389/247‑1 (MQ=255)                   tCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAAAGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTgttcgtt                                                                                                                                                                                                                               >  1:61787/1‑251 (MQ=255)                       ggcgggTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCGCCACGCCTTCGGAACTACCGATCCAGTTACGCTGCATGGTTTTAACGTTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCTTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCATCAGCCGTCGATAACTTGTTCGTTCGcc                                                                                                                                                                                                                           <  2:165163/251‑1 (MQ=255)                                 gtAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCt                                                                                                                                                                                                                 >  2:488725/1‑251 (MQ=255)                                        gTCAGCGGTTTGTCATAGTCGTTAACGTTGAAGGTGACCTCCACGCCTTCGGAACGAACGATCCAGTTACGCTGCATGGTTTTACCGGTGGCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGTTCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCTATAACTTGCTCGTTCGCCAGTACGGTCTGGTCg                                                                                                                                                                                                            <  2:412744/249‑1 (MQ=255)                                            gcgtgtgttCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCggg                                                                                                                                                                                                      <  2:536081/243‑1 (MQ=255)                                             cGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCCGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGCCGTTCGGGc                                                                                                                                                                                                     >  1:99488/1‑251 (MQ=255)                                             cGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCGTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGc                                                                                                                                                                                                     <  1:123507/251‑1 (MQ=255)                                                 ttgtCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCAcc                                                                                                                                                                                                  <  2:460854/250‑1 (MQ=255)                                                     cATAGTCGTTAACGTTGAAGGTGAGCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTg                                                                                                                                                                                             <  1:381552/251‑1 (MQ=255)                                                            gTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGAGCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTAGCGCAGCGCCAGCAACAGCCGTCGATAACTTGGTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAg                                                                                                                                                                                      <  2:557725/251‑1 (MQ=255)                                                                    tGAAGGTGATCTCCTTGACTGGGGAACGTCCGATCCAGTTACTCTGCATGGTTTTAACGGTGTCTGGCCAGTCATCCAGTTTATCCAGCGCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCTCTGCAGGATCTCTTTACTTTAAACTTTGGTATCGCAGGGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGACTGTACGGTCTGGACGTTCGGGCACCAGGTGACCGCAGAAGtcgt                                                                                                                                                                               <  2:613542/250‑3 (MQ=255)                                                                       aGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTcttctt                                                                                                                                                                            >  1:284325/1‑250 (MQ=255)                                                                                         cGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCtttttt                                                                                                                                                          <  1:430835/250‑1 (MQ=255)                                                                                         cGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCtttttt                                                                                                                                                          >  1:18901/1‑250 (MQ=255)                                                                                          ggAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTtat                                                                                                                                                        >  2:467412/1‑251 (MQ=255)                                                                                               gACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGCTCTCTTTACGTTCACCTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGc                                                                                                                                                   >  1:528915/1‑251 (MQ=255)                                                                                                aCCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCt                                                                                                                                                  >  1:601572/1‑251 (MQ=255)                                                                                                aCCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCt                                                                                                                                                  <  1:340353/251‑1 (MQ=255)                                                                                                   gATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGCTTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGACGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCgg                                                                                                                                               >  2:42840/1‑251 (MQ=255)                                                                                                                    tGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCgg                                                                                                                                                                                                       >  1:301396/1‑178 (MQ=255)                                                                                                                    tGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGGCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTCCGTTCGCCAGTACGGTCTGGTCGTTCgg                                                                                                                                                                                                       <  2:301396/178‑1 (MQ=255)                                                                                                                     ggTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTccc                                                                                                                             >  1:173284/1‑251 (MQ=255)                                                                                                                          tAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACg                                                                                                                         <  1:492329/250‑1 (MQ=255)                                                                                                                           aaCGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACgg                                                                                                                        <  1:433490/250‑1 (MQ=255)                                                                                                                                            gATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACa                                                                                                       <  1:39039/250‑1 (MQ=255)                                                                                                                                              tCCAGTTTATCCAGATCGTTGCGCAGCTCGGCAGCGTAAGCAGTGATTTTGGTAAACCACTGCGGGATCTCTTTACGTTCTACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGATGACCGCAGAAGTCTTCTTATATACAAGGCCTTTTTTATACCGCTCGGTGAAGAATTTCTGTTCCCACCGGTAGTATTCCGGCGTACAGGt                                                                                                    >  1:607690/1‑251 (MQ=255)                                                                                                                                                    ttATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGcc                                                                                                >  1:269910/1‑249 (MQ=255)                                                                                                                                                      aTCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCt                                                                                            >  2:502783/1‑251 (MQ=255)                                                                                                                                                             tCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCt                                                                                     >  2:170727/1‑251 (MQ=255)                                                                                                                                                                 tGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTcc                                                                                   >  1:491239/1‑249 (MQ=255)                                                                                                                                                                       gcTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATCCCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATaa                                                                           <  1:431827/251‑1 (MQ=255)                                                                                                                                                                          cGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTAt                                                                                                                  >  2:150356/1‑209 (MQ=255)                                                                                                                                                                          cGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTAt                                                                                                                  <  1:150356/209‑1 (MQ=255)                                                                                                                                                                           gTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCaa                                                                       <  1:8626/251‑1 (MQ=255)                                                                                                                                                                                cGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGccc                                                                  >  2:408994/1‑251 (MQ=255)                                                                                                                                                                                                            cGGGATCTCTTTACGTTCCACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATAcg                                      >  1:114981/1‑251 (MQ=255)                                                                                                                                                                                                                        aCGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCt                                                                                                                                   <  1:537098/146‑1 (MQ=255)                                                                                                                                                                                                                        aCGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCt                                                                                                                                   >  2:537098/1‑146 (MQ=255)                                                                                                                                                                                                                         cgtcaaaCTTTGGTATCGCAGCGACAGCAACAGCCGTCGATAACTTTTTCGTTCGCCAGTACGGTCTGGTCGGTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTGTGAGCTGGTTTTTCATATACGCGATGTTGTCGTa                         <  2:114218/246‑1 (MQ=255)                                                                                                                                                                                                                           ttCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACg                       >  1:33847/1‑251 (MQ=255)                                                                                                                                                                                                                              aaCTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTcc                    >  1:484035/1‑251 (MQ=255)                                                                                                                                                                                                                                cTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTTAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCAc                  >  1:221292/1‑251 (MQ=255)                                                                                                                                                                                                                                  ttGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCAc                  <  1:408994/249‑1 (MQ=255)                                                                                                                                                                                                                                        tCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCgg          >  1:609037/1‑251 (MQ=255)                                                                                                                                                                                                                                               gcCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGgtg   >  2:245800/1‑251 (MQ=255)                                                                                                                                                                                                                                                ccAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGgtgt  <  2:491239/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     GTACAACCCATAAAGGTGTCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGATCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAGCAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGT  >  NC_000913/674014‑674503

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