Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,492,131	A→C	F56C (TTT→TGT)	ppiA ←	peptidyl‑prolyl cis‑trans isomerase A (rotamase A)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,492,131	0	A	C	95.1%	140.5 / NA	41	F56C (TTT→TGT)	ppiA	peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (21/18);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGATTTCCTTAACAAACAGAGCAGAAAAAAAGCAAGTGAAGAGATTCTAAATAGCCTTCAAGAACAAG  >  NC_000913/3491905‑3492366                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             aaCGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGGGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑caca                                                                                                                                                                                                                                  <  1:735495/238‑1 (MQ=255)               aCTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGACGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCtt                                                                                                                                                                                                                     <  1:159708/237‑1 (MQ=255)                        tggCGCTGTCTTTGTCAGCGGTACGTGCCATCGGGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATc                                                                                                                                                                                                           <  1:35662/238‑1 (MQ=255)                             cTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCt                                                                                                                                                                                                      <  1:365790/238‑1 (MQ=255)                             cTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTTGTTTTTTTTTCTGTCTCCTCTC‑GGTGAAAACGCCGCCCTGAAACATAAAACCAGGAATG‑AGCCGGTGGAAGGAAGGGT‑TGTTTTAAAAAACACTGTTTCCATAATCGACAACGGTTTT‑TACAACAACTGGCGCGTTTTTTTTTTTCAGca                                                                                                                                                                                                      >  2:1129109/1‑237 (MQ=255)                                              aCGTGCCATCGCGATGTTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAg                                                                                                                                                                                    <  2:875767/239‑1 (MQ=255)                                                    cATCGGGAGGGTGCCACGCGTGTGGCGCAGGCCGTTATCGGCTTCATTGTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCGC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAAGCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAg                                                                                                                                                                               <  1:189704/238‑1 (MQ=255)                                                             ggTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATc                                                                                                                                                                                                           <  1:586953/201‑1 (MQ=255)                                                             ggTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGGTTGGTTTTTTTTTGTTCCTTTGCCC‑CGGGTAAAACCCCCCCCAAATCATAAAACCCCGAATGAACGCGGGGAAAGGAGAGGT‑GGTTTTAAAAACGCATGTTAAAAAAATCGAACAAGTTTTGACACAAacaagcg                                                                                                                                                                                                                          >  2:345202/1‑184 (MQ=255)                                                                 ccACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑caca                                                                                                                                                                                                                                  <  1:5108/174‑1 (MQ=255)                                                                      cgTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGCGGGGTTTGGTTTTTTTTCTTGCCTTTCCCC‑GGGGAAAAGGCCCCCCTCATACAAAATACCCGGAAAA‑ATGCGGGGAAAGGGAG‑GG‑TGTTTTAAAAAACACCGTTTTCATAAATGACACACGCTTG‑TACACACAACGCCGCCTTTTTGTTATCCCACACCCGTTGGCTGTTACACGCTGAGGGTGTGCACAATaaactcg                                                                                                                                                            >  2:221860/1‑233 (MQ=255)                                                                              gcAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTcccc                                                                                                                                                     <  1:15147/238‑1 (MQ=255)                                                                                       tATCGGCTTCATTTTTGATTGGC‑GGGTTTTGTTTTTTTTTCTTGATTTGCTCTGGGGAAACCCCCCCCCGCAATATTAAACCCGGGAA‑G‑ACGCGGGGGAAGGGGGGGTGTTTTTTAAAAACCCCTG‑TCACATAAATGACACAGGTTTG‑CACAGACACTGGAGCTTTTTGTTTTTTCAACGCCCCGTCGGTGTTATCTGCTGAGGGTGGTAACAAAAAGTGCGGGGGCCCTTTTGCttgc                                                                                                                                           >  2:1156631/1‑236 (MQ=255)                                                                                              ttCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATtgctg                                                                                                                                    <  1:705343/239‑1 (MQ=255)                                                                                                  tttttGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCggg                                                                                                                                <  1:1129109/239‑1 (MQ=255)                                                                                                      tGATTGGA‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGagaa                                                                                                                            <  2:1082377/239‑1 (MQ=255)                                                                                                                                 tgcATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGagaaagagcagaa                                                                                                                   >  1:652107/1‑222 (MQ=255)                                                                                                                                 tgcATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGagaaagagcagaa                                                                                                                   <  2:652107/222‑1 (MQ=255)                                                                                                                                          tC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATcgccg                                                                                          >  1:293375/1‑238 (MQ=255)                                                                                                                                             ggTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTAt                                                                                                                                                                                                            >  2:511562/1‑122 (MQ=255)                                                                                                                                             ggTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTAt                                                                                                                                                                                                            <  1:511562/122‑1 (MQ=255)                                                                                                                                              gTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAg                                                                                      >  2:919911/1‑239 (MQ=255)                                                                                                                                                       ccgccCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGAttg                                                                             >  2:418588/1‑238 (MQ=255)                                                                                                                                                         gccCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTTGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGGTTTGa                                                                           >  1:715731/1‑239 (MQ=255)                                                                                                                                                                       aGCCAGGAATG‑ACGCGGGGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGCATTCCt                                                             >  1:42520/1‑239 (MQ=255)                                                                                                                                                                         ccAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑caca                                                                                                                                                                                                                                  <  1:660228/72‑1 (MQ=255)                                                                                                                                                                         ccAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑caca                                                                                                                                                                                                                                  >  2:660228/1‑72 (MQ=255)                                                                                                                                                                         ccAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGCGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTAGATTTGCACATGAGATTTCCTTa                                                           >  2:271986/1‑239 (MQ=255)                                                                                                                                                                            ggAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGt                                                                                                                                                                                                 <  1:896224/102‑1 (MQ=255)                                                                                                                                                                            ggAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGt                                                                                                                                                                                                 >  2:896224/1‑102 (MQ=255)                                                                                                                                                                            ggAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGCTTTGAACATGAGATTTCCTTaaca                                                        >  1:882397/1‑239 (MQ=255)                                                                                                                                                                             gAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGATTTCCTTaacaa                                                       >  1:1157259/1‑239 (MQ=255)                                                                                                                                                                                   aCGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGATTTCCTTAACAAACaga                                                  >  2:212579/1‑239 (MQ=255)                                                                                                                                                                                    cgcgGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACt                                                                                                                                                                                                                            <  1:682082/68‑1 (MQ=255)                                                                                                                                                                                    cgcgGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACt                                                                                                                                                                                                                            >  2:682082/1‑68 (MQ=255)                                                                                                                                                                                        gTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGATTTCCTTAACAAACAGAGCAGa                                             >  1:139524/1‑239 (MQ=255)                                                                                                                                                                                                   tgt‑tgtTATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGATTTCCTTAACAAACAGAGCAGAAAAAAAGCAAg                                  >  1:154378/1‑239 (MQ=255)                                                                                                                                                                                                        gttATAAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGATTTCCTTAACAAACAGAGCAGAAAAAAAGCAAGTGaa                              >  1:997379/1‑239 (MQ=255)                                                                                                                                                                                                           ataAAAACCGCTGTTCACATAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGATTTCCTTAACAAACAGAGCAGaa                                            >  1:693706/1‑222 (MQ=255)                                                                                                                                                                                                                             aTAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCAt                                                                                               >  1:891538/1‑153 (MQ=255)                                                                                                                                                                                                                             aTAATCGACACAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCAt                                                                                               <  2:891538/153‑1 (MQ=255)                                                                                                                                                                                                                                    acacAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATcgccgc                                                                                         >  2:198731/1‑152 (MQ=255)                                                                                                                                                                                                                                    acacAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATcgccgc                                                                                         <  1:198731/152‑1 (MQ=255)                                                                                                                                                                                                                                    acacAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGATTTCCTTAACAAACAGAGCAGAAAAAAAGCAAGTGAAGAGATTCTAAATAGCCTTCAAGAACAAg  >  2:288692/1‑239 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             AACGTTGATAAAGAACTGGCTGGTGGCGCTGTCTTTGTCAGCGGTACGTGCCATCGCGATGGTGCCACGCGTGTTGCGCAGGCCGTTATCGGCTTCATTTTTGATTGGC‑GGGTTTGGTTTTTTCTGCTGCATCTGCTC‑GGTGAAACCGCCGCCCTGAATCATAAAGCCAGGAATG‑ACGCGGTGAAAGGTAGTGT‑TGTTATAAAAACCGCTGTTCACATAATCGACAAAGTTTTG‑CACAGACACTGGCGCTTTTTGTTTATCCAGCTCCAGTTCGATGTTACCAGCTGAGGTTGTCAACAATACGTGCGGGTCCCCTTTCGCTGCCATTGCTGCGGGAGAAAGAGCAGAAAGAGCGAAAACAGCAGCCATCGCCGCCAGGGTCGATTTGAACATGAGATTTCCTTAACAAACAGAGCAGAAAAAAAGCAAGTGAAGAGATTCTAAATAGCCTTCAAGAACAAG  >  NC_000913/3491905‑3492366

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