Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	159,258	G→T	21.7%	V285V (GTC→GTA)	gluQ ←	glutamyl‑Q tRNA(Asp) synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	159,258	0	G	T	21.7%	82.1 / 6.3	58	V285V (GTC→GTA)	gluQ	glutamyl‑Q tRNA(Asp) synthetase
Reads supporting (aligned to +/- strand):  ref base G (18/26);  new base T (11/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 9.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAGCGAGATTTGCCTTACTGTTGGT  >  NC_000913/159032‑159495                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    aTCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCGCCCAGTTTTT‑AAGGcggactga                                                                                                                                                                                                                                         >  1:173883/1‑239 (MQ=255)        tCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TCAGCCG‑CCCGGTTTTTACGGGTGAACGAAgagttt                                                                                                                                                                                                                                  >  2:129313/1‑234 (MQ=255)          aCCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTTCCGACTCGGGCACGGCGG‑TAAGCCG‑CCCGTTTTTGACGGCTTGCTGAAGGATTTGc                                                                                                                                                                                                                                >  1:374885/1‑239 (MQ=255)              gtgGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAACCG‑CCAGTTTTTTAAGGCTCACCTAAGGGTTTtctcc                                                                                                                                                                                                                             >  1:291216/1‑237 (MQ=255)                ggACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑GAAACCC‑CCCGGTTTTGTCGGCTGGCCAAAGGATTTgtcccggc                                                                                                                                                                                                                          >  2:505020/1‑232 (MQ=255)                            gcttttttaGCCTCTTTTTCCCTGCGTCGTAGGCCCTTGCGGAAAAAATTAGGCACTGGGGTAAAAATAGTACACCTC‑GGT‑AGTGTAAAACCTCATTCAGGACAAAAAAATAGCGGCCAATCATAGCTCCGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑Gc                                                                                                                                                                                                            <  2:137578/231‑1 (MQ=255)                                                    cgGTTTAGCCCCTTGCGGAAAAAATTAGCGACTCGGGTAAAAATAGTACCCCTC‑GGT‑AGTGTCAAACATCATTCAGGAAAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGcccc                                                                                                                                                                                   <  2:291216/239‑1 (MQ=255)                                                      gCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGCCGG‑TAAGCCG‑CCCGTTTTTTACGGCTCACTGAAGGGTTTTCTCGCCCCGTAAAT‑CCT‑GCCAAGTG‑GCTTCTTCCGCCTGCCCCAg                                                                                                                                                                                 >  1:507532/1‑239 (MQ=255)                                                       tttgggcccTTGGGGAAAAAATTGGCGACTCGGGTAAAAATAGTCCCCCCC‑GGT‑AGTTGCAAAAATGATTCAGGACAAAAAAATAGCGTGTAATCATAGCCCAGCATGACGCATTTGAGAAATTTGAATTTAAAATTGCCGACTCGGGCACGGAGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGa                                                                                                                                                                                <  2:700223/232‑1 (MQ=255)                                                         taaccccTTGCGGCAAAAATTAGCGAATCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGaaa                                                                                                                                                                              <  1:129313/234‑1 (MQ=255)                                                                 tGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAACCG‑CCAGGTTTTTACGGCTGACTTAAGGGTTTTGCCCAGGCCGAAAA‑CCC‑GCC‑AGGGGGGTTTTGCCTCCTGCCCCCGaaaattaaatg                                                                                                                                                                      >  2:181107/1‑231 (MQ=255)                                                                      aaaaatTTAGGCACCCGGCTAAAAATGGTACCCCAG‑GGTGGGTGTC‑AACCTCTTGCAGGAAAAAAAAATGGCGGCTAAGCAGAGCTCGGAATGACTCTTTTGAGAATGTGGAATTGACAATTGCAGACTCGGGCACGGTGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGAGGCTGAAAT‑CCT‑GCC‑AGTGGGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑cc                                                                                                                                                                   <  2:642018/235‑1 (MQ=255)                                                                              aGCCACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCGAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTAc                                                                                                                                                       <  1:732943/239‑1 (MQ=255)                                                                                          aaaaaTAGTACACCGC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGAt                                                                                                                                            <  1:181107/238‑1 (MQ=255)                                                                                           aaaaaaGTACCCCTCGGGT‑AG‑GTCAAAACGTAATCAGGACAAAAAAATAGCGGGTAAGCAAAGCTCAGGATTGAGCATTTGAGAATTGTGAATTTACAATTGCCGACTCGGGCACGGTGG‑TAAGCCG‑ACAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCg                                                                                                                                          <  1:457678/239‑1 (MQ=255)                                                                                              aTAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑GAAACCG‑CCAGTTTTTGTCGGCCGCCGCAAGGAATTTGTCGGCGCTGAAAA‑CCT‑TCC‑CATGGGTCTTTGCCTGTGGCCCCCGAAATAGAA‑TG‑ACGCCGATTATAACCGGGGGGGGagagc                                                                                                                                         >  2:722752/1‑233 (MQ=255)                                                                                                  tACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGGTTAACCCG‑CCCGTTTTTTA‑GGCTGACCTAAGGCTTTGCTCGGCGATGTAATCCCT‑GCC‑A‑TGGGCTTTCCCCTCCTGCCCCCCAAATTTAAATGCCCG‑CG‑TTGGAACCGGGCGTgggaccccattg                                                                                                                                    >  2:382724/1‑228 (MQ=255)                                                                                                  tACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCCGAAATTGAAGTG‑CCG‑CGATTTGTACCTGGCGTGGATCCCCttttg                                                                                                                                   >  1:882130/1‑237 (MQ=255)                                                                                                  tACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGGTTTTTACAGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCCAAGTG‑GCTTCTGCCTGCGGCCCCAGCAATTTGAGTG‑CCG‑CGATTAGGACCGGGGCTGGTACCCCTTTCg                                                                                                                                   >  1:390476/1‑239 (MQ=255)                                                                                                   accccTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCAGAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCgg                                                                                                                                  <  1:828460/236‑1 (MQ=255)                                                                                                       cTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGGTTTTTACGGCTGACTGAAGGGTTTTCTCGACGCTGGAAT‑CCT‑GCC‑AGTGTGCTTCTTCCTGCTGCCCCAGAAATTGAAAGG‑CCG‑CGCTTTATAACCGGCGGGAGTCGCCTTTCGGCaa                                                                                                                               >  1:345646/1‑238 (MQ=255)                                                                                                       cTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGGTTTTGACGGCTGACTGAAGGAGTTGTCCGACGCTGGAAA‑CCT‑TCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTTAAGTG‑CCG‑CGGTTAGTACCGGGCGTGGATCCCCTTTCGGcca                                                                                                                               >  1:401375/1‑236 (MQ=255)                                                                                                         c‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTTACCGCTGACTTAAGGATTTGCTCGCCGCTTAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCCGAAATTTAAGTG‑CCG‑CGATTACAACCGGGCGTGGACTGCCATTCGGCacagc                                                                                                                            >  2:190364/1‑235 (MQ=255)                                                                                                          tggT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGc                                                                                                                            <  1:473027/238‑1 (MQ=255)                                                                                                             t‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑CGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGGTCGCCTTTCGGCAACGCAg                                                                                                                          >  1:199665/1‑238 (MQ=255)                                                                                                             t‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTTGAGTG‑CCG‑CGATTATTACCGGGCGTGGATCCCCTTTCGGCCACGCAgg                                                                                                                         >  1:4679/1‑239 (MQ=255)                                                                                                                 tgtCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCCCGACGCTGAAAT‑CCT‑TCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGGG‑CCG‑CGGTTAGGACCGGGCCTGGATCCCCTTTTGGCAACGCAGgggc                                                                                                                      >  1:846874/1‑239 (MQ=255)                                                                                                                 tgtCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAACCG‑CCCATTTTTGACGGCTGATCGAAGGGTTTGCTCGCACCCGAAAC‑CCT‑GCC‑CGGGGGCTTTTCCCTGCTTCCCCCGAAAATAAAAGG‑CCG‑CGGTTAGTAACGGGCGGGGAGTCCCTTTCGGGaaaacaagggg                                                                                                                      >  2:235419/1‑230 (MQ=255)                                                                                                                           catTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAAGG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGGGCATGGTTCt                                                                                                              >  1:547653/1‑237 (MQ=255)                                                                                                                             ttCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCCACGCAGGCGCATGATTTTGCtt                                                                                                          >  1:708118/1‑239 (MQ=255)                                                                                                                                aGGAAAAAAAAAATGCGGATAATCACAGCCCAGCATGACTCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGa                                                                                                       <  2:547653/239‑1 (MQ=255)                                                                                                                                  gACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGaaa                                                                                                     >  1:167141/1‑239 (MQ=255)                                                                                                                                  gACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAAGG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCCACGCAGGGGCATGATTCTGCTTGGaaa                                                                                                     >  1:710151/1‑239 (MQ=255)                                                                                                                                       aaaaaTAGCGGCTAATCATAGCTCGGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑cc                                                                                                                                                                   <  1:197616/173‑1 (MQ=255)                                                                                                                                       aaaaaGCGCGTCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTa                                                                                                <  1:668638/239‑1 (MQ=255)                                                                                                                                        aaaaTAGCGGCTAAGAAAAGACCAGCATGCCGCATTTAAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAg                                                                                               <  1:427961/239‑1 (MQ=255)                                                                                                                                         aaaTGGCGGTTAAACATAGCCCAGCAGGCAGCATTTGATAATGTTTAATTTCCAATGGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACGGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAgc                                                                                              <  2:456173/239‑1 (MQ=255)                                                                                                                                               cGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTTACGGCTGACTGAAGGATTTGCTCGACGCTGGAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCGGCCCCCGAAATTGAAGTG‑CCG‑CGGTTATTACCGGCCGTGGGTCGCCTTTCGGCAACGACGCGGCAGGATTCTCGCTGGAAAGATTAACGCCttt                                                                                        >  2:98434/1‑238 (MQ=255)                                                                                                                                                                gCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTTTGGATTAAGGGCCAACg                                                                       >  1:168813/1‑239 (MQ=255)                                                                                                                                                                 gtggACGTATTTGAAAATTGTGAATTTACAATTGCCGGCTGGGGTACGGAGG‑TAAGCCG‑CCGGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTTAAAT‑CCT‑TCC‑AGTGTGCTACTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCgg                                                                      <  2:184561/236‑1 (MQ=255)                                                                                                                                                                 cATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAACGTTTAGCGCCTTTTGGGTTAAAGGCCAGCgg                                                                      >  1:399549/1‑239 (MQ=255)                                                                                                                                                                 cATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGGACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGATTATCGCCTTTTGGATTAAGCGCCAGCgg                                                                      >  1:810868/1‑239 (MQ=255)                                                                                                                                                                  atggcggttttGAGAATTTTGAATTTACAAATGCCGTCTCGGGCCCGGCGG‑TAAGCCG‑CCAGTTTTGGCCTGTTGATGGAAGGATTTGCTCGAAGCGGAAAT‑CCT‑GCC‑AGTGTGCGTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCgg                                                                      <  2:173883/230‑1 (MQ=255)                                                                                                                                                                      cGCATTTGAGAATGTTGAATTTCCAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGa                                                                  <  1:118418/238‑1 (MQ=255)                                                                                                                                                                           ttGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGGAACGCAGGCGCATGATTCTGCTTGGAAAGTTTTGCGCCTTGGGGATTAAGGGCCAGCGGCAGGTGAATg                                                            >  1:489655/1‑239 (MQ=255)                                                                                                                                                                               gaATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCCGGCGGGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGTTTAAGCGCCAGCGGCAGATGAATGTaa                                                         >  1:851717/1‑238 (MQ=255)                                                                                                                                                                                            aCAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑ACAGATTTTTACGGCTGGCTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCCGAAATTGAAGTG‑ACG‑CGATTAGTACCGGGCGTGGATCGCCTTTCCGCAACGCCGGAGCATGGTTCTGCCTGGAAAGTTTAGCGCCTTGTGGTATAAGCGCCACCGGCAGATACATGTAAACTGGGACTTT‑CCa                                          >  2:169156/1‑239 (MQ=255)                                                                                                                                                                                                 tGCCGACTCGGGCGCGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCaa                                     <  1:190364/239‑1 (MQ=255)                                                                                                                                                                                                              aCGGGGG‑TAAGCCG‑CGAGTTTTTGGAGGGTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAg                        <  2:199665/239‑1 (MQ=255)                                                                                                                                                                                                               cggggg‑TAAGCCC‑GCAGTTTTGGACGGGTGACTTAAGGATTTGGTCGAAGCTGAAAT‑CCTGGCG‑AGT‑TGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAg                        <  2:374885/238‑1 (MQ=255)                                                                                                                                                                                                                ggcgg‑TAAGCCG‑CCAGGTTTTTACGGCTGACTGAAGGGTTTGCTCGACCCTGAAAT‑CCT‑GCC‑AGGGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGGGGATCGCCTTTCGGCAAAGCAGGCGCATGATTCTGCTTGGACAGTTTAGCGCCTTGTGGGTTAAGAGCCAGCGGCCGATGCATGTAATCTGGCCCTTTCCCCGCCaaaaaacgggaaaacg                       >  2:33688/1‑228 (MQ=255)                                                                                                                                                                                                                         gccg‑ccATTTTTGGGCGGCTGTGTGAAGGTTTGGCTGCACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGGCTGGTGCCCCAGAAATGTAAGTG‑CCG‑CGATTAGGACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAGCGAGATTTg               <  2:345646/238‑1 (MQ=255)                                                                                                                                                                                                                          ccg‑ccAGTTTTTTACGTCTGACTGAAGTATTTGCCCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAGCGAGATTTGcc             <  2:168813/239‑1 (MQ=255)                                                                                                                                                                                                                            cccccgTTTTTTGCGGTTTACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAGCGAGATTTGCCt            <  2:791552/234‑1 (MQ=255)                                                                                                                                                                                                                            ccccAGTTTTTGGCGGCTGACTGAAGGATTTTCCCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAGCGAGATTTGCCt            <  2:882130/237‑1 (MQ=255)                                                                                                                                                                                                                             cccATTTTTTGCCGGTTGACTGAAGGTTTTGCTCGACGTGGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAGCGAGATTTGCCtt           <  2:167141/238‑1 (MQ=255)                                                                                                                                                                                                                                ggTTTTTGCCGGTTTACTGAAGGATTTGCTGGCAGGTGAAAT‑CCT‑GCA‑AGTGGCCTGCTGGCTGCTTACCCAGAAAGTGAACTG‑CCG‑CGATTAGTACCTGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGGGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAGCGAGATTTGCCTTACt        <  2:219389/238‑1 (MQ=255)                                                                                                                                                                                                                                      tGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAG‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAGCGAGATTTGCCTTACTGTTGGt  <  1:315131/239‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ATCACCGTCACCTGTGGACGGGCGACTGCCTGTTCAGCCTCGCTTTCCTCGCGGCTTAGCACCTTGCGGCAAAAATTAGCGACTCGGGTAAAAATAGTACACCTC‑GGT‑AGTGTCAAACATCATTCAGGACAAAAAAATAGCGGCTAATCATAGCTCAGCATGACGCATTTGAGAATGTTGAATTTACAATTGCCGACTCGGGCACGGCGG‑TAAGCCG‑CCAGTTTTTGACGGCTGACTGAAGGATTTGCTCGACGCTGAAAT‑CCT‑GCC‑AGTGTGCTTCTGCCTGCTGCCCCAGAAATTGAAGTG‑CCG‑CGATTAGTACCGGGCGTGGATCGCCTTTCGGCAACGCAGGCGCATGATTCTGCTTGGAAAGTTTAGCGCCTTGTGGATTAAGCGCCAGCGGCAGATGAATGTAATCTGGCACTTT‑CCAGCCAAAAAGCTGGTACAGCGAGATTTGCCTTACTGTTGGT  >  NC_000913/159032‑159495

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