Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	252,387	G→A	13.3%	K29K (AAG→AAA)	yafO →	mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	252,387	0	G	A	13.3%	169.0 / 3.3	70	K29K (AAG→AAA)	yafO	mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (18/42);  new base A (8/2);  total (26/44)
Fisher's exact test for biased strand distribution p-value = 4.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

ACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAA  >  NC_000913/252156‑252621                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                aCATGCTTTGGGAAAAAGAAGAAAAAAAGCCCAAAAAGGCGCTGTTCCTGCCAAGTGCTGCAAGATTTGAGGAAATTACAATGCGCGCTGAAAAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGAc                                                                                                                                                                                                                                          <  2:239670/238‑1 (MQ=255)        tgctGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGTTTTTATTTCCTTTAAGC‑GTGCGGGGGtttt                                                                                                                                                                                                                                  >  1:236845/1‑239 (MQ=255)         gttggaCAAGGGGGGAAAAACCCCATAAAGGCGCGCTTCCGTTCAAGGTGTGCAAGATTAGAGTAAATTACACGCCGCGCTGAACAATATCGTAATGATATGACTGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTg                                                                                                                                                                                                                                 <  2:1229707/234‑1 (MQ=255)          ctGAATAAGGGGAAAAAAAGGCCATAAAGGCCCGCTTCCGTCCAAGTTCTGCAAGATTAGAGGAAATTACCCGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGc                                                                                                                                                                                                                                <  2:247653/239‑1 (MQ=255)           tGAACAAGGGGAGAAAAAGGCCATAAAAGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTGCACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGcc                                                                                                                                                                                                                               <  2:410846/239‑1 (MQ=255)           tGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCTGCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGc                                                                                                                                                                                                                                <  1:799812/238‑1 (MQ=255)            gAACAAGGGGAGAAAAAGCCCAAAAAGGCGCGCTTCCGTCCAAGTGCTTCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGcc                                                                                                                                                                                                                               <  1:666822/238‑1 (MQ=255)            gAACAAGAGGAAAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGcc                                                                                                                                                                                                                               <  1:1022153/238‑1 (MQ=255)                          aaaGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCCATAAAC‑GTGACGGGGTTTTTCCCACTATATTTG‑GTc                                                                                                                                                                                                                >  1:143322/1‑238 (MQ=255)                            aGCCCATAAAGGCGGGTCTCCGTCCAAGTTCTGAAAGATTAGAGGAAATTACACTGCGCGCTGAACAATATCTTAATGATAGGCAGGATGATGATGTCAATGACTTTAAGGAAGAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTcgcg                                                                                                                                                                                                             <  2:1194678/239‑1 (MQ=255)                               ccATAAAGGGGCGGTTTCGGTCAAGTTATTCAAAATTAGAGGAAAGTACACGCCGGGGTGAAAAATAACTTAATGAAATGAAGGATGAGGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑Tg                                                                                                                                                                                                         <  2:1149664/239‑1 (MQ=255)                               ccATAAAGGCTCCCTTCCGTCCAAGGGCTGTAAGATTTGAGGAAATTACCCGCCGCGTTGAACAATATCTGAATGTTATGACGGATGATGATTTCAATTACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATGCGTCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑Tg                                                                                                                                                                                                         <  2:793516/239‑1 (MQ=255)                               ccATAAAGGCGCGCTTCCGTCAAAGTGCTGCAATATTAGAGGAAATTCAACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑Tg                                                                                                                                                                                                         <  2:1033236/239‑1 (MQ=255)                                 aTAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCa                                                                                                                                                                                                       <  2:218148/239‑1 (MQ=255)                                          cgcTTCCGGCCAAGTGTTGCAATATGAGAGGAAATTAACCGCCGCGCTGAAAAATATCTTAAGTAGATGACGGATGATGATTTCAATGTCTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATGCGCCTGCAACTTACAGCAGAGGAACGTTATGCGTTAACGGCGGATTTTATTTCCTATAAAC‑GTTAAGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTacgac                                                                                                                                                                                              <  2:814984/239‑1 (MQ=255)                                           gcTTCCGTCCAATGGCTGCAAAATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTacgac                                                                                                                                                                                              <  1:329126/238‑1 (MQ=255)                                           gcTTCCGTCCAAGTGCCGCAAAATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGAAGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTacgacg                                                                                                                                                                                             <  1:744360/239‑1 (MQ=255)                                            cTTCCGTCCATGTGCTGCAAAATTAGAGGAAATTACCCGCCGCGCTGAACAATATCTTAATGGTATGACGGATGATGATTTCAATGCCTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGTATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTT‑GTCGCGA‑TGCACTCTAcgacga                                                                                                                                                                                            <  2:842600/239‑1 (MQ=255)                                            cTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGGCGCGCTTAACAATATCTTAATGATATGACGGATGATGATTTCAATGGCTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTacgacg                                                                                                                                                                                             <  1:334604/238‑1 (MQ=255)                                                             cAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTacg                                                                                                                                                                                                >  1:1023644/1‑218 (MQ=255)                                                                          aaTTACACGCCGCGCTGAACCATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGa                                                                                                                                                             <  1:127013/239‑1 (MQ=255)                                                                          aaTTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGGATTTATTTTTCCAAAGA‑GTGAAGGGGTTTTTCCAGATATGTTTT‑TACGGTAGTGC‑CTCTTCAGCGACCCCTTTTCCTGGCCCTTAT‑ACAAAATTGa                                                                                                                                                             >  2:837845/1‑239 (MQ=255)                                                                           aTTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTTCTCTAAGC‑GTGACGGGGTTTTTACACGAATATTTT‑GTCTCGA‑TGCACACTCCAGCAATCCCTTTACCTGGCCCTTAA‑TTAAAttttga                                                                                                                                                            >  2:1060076/1‑236 (MQ=255)                                                                               cacGCCTGGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAg                                                                                                                                                       <  1:1060076/239‑1 (MQ=255)                                                                               cacGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAg                                                                                                                                                       <  1:301421/239‑1 (MQ=255)                                                                               cacGCCGCGCTGAACAATATCTTAATGATATGACGGATGAGGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAg                                                                                                                                                       <  1:184269/239‑1 (MQ=255)                                                                                accccgccatGAACAATATCTTAATGAAATGACGGAAGATGATTTCAATGACTTTAAGGAATAAGGCTGCGGGTATTCAAAACAAAACTTATTGTCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGt                                                                                                                                                      <  2:444591/230‑1 (MQ=255)                                                                                      cgatGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCTTGCAACTTACAGCAGAGGAAATTGATGCGTTAACGGCGGATTTTATTTCCTATTAAA‑GTTTAGCGGTTTTGTCCGTTTTAATTGTGTC‑ATT‑TGCTCTCTACTAACACGCCTTTACCTTGCCATTAA‑TAAAAATTTA‑GCAAGtgagtca                                                                                                                                                >  2:1144725/4‑233 (MQ=255)                                                                                                 atCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTgg                                                                                                                                     <  1:1212048/239‑1 (MQ=255)                                                                                                  tCTTAAGGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGGATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTTATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGc                                                                                                                                    <  1:1299412/239‑1 (MQ=255)                                                                                                     tAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGGGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTTGCCCTTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTTGCaaa                                                                                                                                 >  1:388869/1‑239 (MQ=255)                                                                                                            atGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGGGTTTTGCCAGATATATTTT‑GTCGCGA‑TGCACTCT‑‑‑ACGACTCCTTTACCCTGCCATTAA‑TCACATTTTA‑GCGAGTTGCTCCTATTCCACTGGCCAATGAGaataa                                                                                                                        >  1:388879/1‑238 (MQ=255)                                                                                                            atGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGGGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTTT‑GCGCGTTGCTCcta                                                                                                                                              >  1:904536/1‑219 (MQ=255)                                                                                                                cGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAAC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGGGA‑TGCACTCTACGACGACTCCTTTCTCTGGCCATTAA‑TTAAATTTGA‑GCGAGTTGCTCATATTCATTTTGCAAATGAGAAAAATc                                                                                                                      >  1:612023/1‑239 (MQ=255)                                                                                                                     gatgatTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTTCTATAAGC‑GTGACGGTGGTTTTCCAGATATATTTT‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TTAAATTTTA‑GCGAGTTGCTCCTATTCATTTTGCCAATGAGAATAATCCCtt                                                                                                                  >  1:207398/1‑238 (MQ=255)                                                                                                                             cAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTTTTTCCTATAAAC‑GGGACGGTGTTTTGCCAGCTATATTTG‑GTCGCGA‑TGACCTCTACGAGGACTCCTTTACCCGGCCATTAATTAAAATTTGA‑‑CGAGTTGGTCCTATTAATTCGCCCAATTAGAATAACACATTTCCCCccca                                                                                                         >  2:1015559/1‑239 (MQ=255)                                                                                                                                    tttAAGGAATAAGGATGGGGGTTTTAAAAACAAAATTTTTTCGCCTGACACTTACAGCAGGGGAACTTGATGCTTTAACGGCGGATTTTATTTCCTATAAGC‑TTGACGGTTTTTTGCCAGATTTATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTgcgc                                                                                                  <  2:612023/239‑1 (MQ=255)                                                                                                                                            aTAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAAC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTTGCAAATGAGAATAATCCATTTCCCCACCAATTGGGCCAATTca                                                                                           >  1:1243994/1‑238 (MQ=255)                                                                                                                                              aaGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTTCTATAAAC‑GTGACGGGGTTTTTCCAGAGATAATTG‑GTCGGGA‑TGCACTCTACGAGCACGCCTTTACTTGGGCATTAA‑TCCAAATTTA‑CAAGGTTTGTCTTATTCACTTTGCAAATGAGAAAAAACCCATTCGCCCACAATTTCGCCCAtttacaa                                                                                        >  2:1288715/1‑234 (MQ=255)                                                                                                                                                 gATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGa                                                                                      <  1:640760/238‑1 (MQ=255)                                                                                                                                                  aagCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGa                                                                                      <  1:22359/235‑1 (MQ=255)                                                                                                                                                     cGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTT‑GTCGCGA‑TGCCCTCTACGACGACTCCTTTACCTTGCCCTTAA‑TCAAATTTTATGCG‑GTTGCTCCTATTCATTTGGCAAATTGGAATAATCCATTTCCCCCCCCGTTTGCCCAATTTCACAGAACCaa                                                                                 >  1:625453/1‑239 (MQ=255)                                                                                                                                                     cGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTT‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCCTTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATTAGAATAATCCATTTCCGCCCCAGTTTCGCCAATTCAGCAGAACGaa                                                                                 >  1:911593/1‑239 (MQ=255)                                                                                                                                                     cGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCCCAGTTGCGCCAATTCACCACAACCaa                                                                                 >  1:1297649/1‑239 (MQ=255)                                                                                                                                                     cGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAAC‑GTGACGGTGTTTTTGCAGATATATTTT‑GTCGCGG‑TTCACTCTAAGACGAATCCTTTACCTGGCCCTTAA‑TCAAATTTTA‑GCGAGTTGCTCATATTCATCTTGCACATGAGAATAATCCATTTCCGCcaca                                                                                                         >  1:1068893/1‑215 (MQ=255)                                                                                                                                                           ttCAAAACAAAACTTATTCGCATGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGA‑GTGAAGGGGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCAATCTACGACGACTCCTTTACCTGGCCCTTAA‑TCAAATTTGA‑GTGAGGTGATCATTTTCATCTGGCAAATGAGAGTAATCCCTTTTCGCCACGTTTGTGCCAATTCACCAGAACGAATAACGa                                                                           >  1:326037/1‑239 (MQ=255)                                                                                                                                                             caaaacaaaaCTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAg                                                                         <  1:1288715/239‑1 (MQ=255)                                                                                                                                                             aaaaaaaaaacTTATTCGCCTTCAACTTAAAGCGGAGGAACTTTATGCTGTAACGGCGGATTTTATTTCCTATAGGC‑GTGAGGGTGTTTTGCCAGATATTTTTG‑GTCGCGA‑TGCACTTTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAg                                                                         <  2:754249/238‑1 (MQ=255)                                                                                                                                                                 acaaaacTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGGGAATAATCCATTTCCGCCACAGTTGCGCCAATTTCGCAGAACGAATGACGAAgcg                                                                       >  1:754249/1‑237 (MQ=255)                                                                                                                                                                  caaaccGTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCCCTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCAt                                                                    <  1:1094417/239‑1 (MQ=255)                                                                                                                                                                  caaaacTTATTCGCCTGCAACTTACAGCAGGGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTTTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCAt                                                                    <  1:487816/239‑1 (MQ=255)                                                                                                                                                                  caaaacTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCAt                                                                    <  1:1263774/239‑1 (MQ=255)                                                                                                                                                                  aaaaacTATTTCTCCTTCAATTTAGAGCAGGGGACCTTGATCGGTTAACGGCGGATTTTATTTCCTAAAAGC‑GGGACGGTGTTTTGCCAGATAGATTTG‑GTCGGGT‑TGCACACTACGACGACTCCTTTACCTGGGCATTAA‑TCAAATTTGA‑GCGGGTTGCTCATATTCATCTGGCAAAGGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCAt                                                                    <  2:236845/238‑1 (MQ=255)                                                                                                                                                                    aaacTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCAttt                                                                  <  1:1015559/239‑1 (MQ=255)                                                                                                                                                                    aaacTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCAttt                                                                  <  1:443962/239‑1 (MQ=255)                                                                                                                                                                     aacTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTTTTTCCTATAAGC‑GTGAAGGCGTTTTGCCAGATATATTTG‑GTCGCGA‑TTCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑ACGAGTTGCTCCAATTTATCTTGCAAATTAGAATAATCCATTTCCGCCACAGTTGCGCCACTTCAGCAGAACGAATGACGAAGCGCAtttt                                                                 >  1:339741/1‑238 (MQ=255)                                                                                                                                                                     aacTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTT‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTTAGCAGAACGAATGACGAAGCGCCTTTg                                                                 >  1:902023/1‑239 (MQ=255)                                                                                                                                                                     aacTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGGTGCTCATATTCATCTTGCAAATGAGGATAATCCATTTCCGCCACAGTTTCGCCAATTTACCAGAACGCATTACGAAACGCATTTg                                                                 >  1:1133789/1‑239 (MQ=255)                                                                                                                                                                        ttttgcccttGCAACTTTCAGCGGGGGAACTTGATGGTTTAACGGCGGTTTTTTTTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGGGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGta                                                              <  2:911593/230‑1 (MQ=255)                                                                                                                                                                                              gtagGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGAGGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTtgatg                                        <  1:325844/237‑1 (MQ=255)                                                                                                                                                                                                agGAACTTGATGCGTTAACGGCGGATTTTATTTTCTCTAAAG‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTAACTTGCCATTAA‑TTCAATTTTG‑GGGAGGTGCTCAAATTCATCTGGCAAATGAGAATAATACATTTCCCCCCCAGTTTGGCCCATTTCACCGAACGAAAGGCGAAACCCATTTTGGATATTGTTAGGGGGCGTTtgagtaa                                      >  2:720385/1‑235 (MQ=255)                                                                                                                                                                                                agGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGGGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTTA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGTTGag                                      >  1:830754/1‑239 (MQ=255)                                                                                                                                                                                                    actttgtggcttAACGGCGGGTTTTTTTTTCTATAAGC‑TGTACTGTTGTTTGCCAGATATATTTG‑GCCGCGA‑TGCCCTCTTCGAAGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                  <  2:562574/229‑1 (MQ=255)                                                                                                                                                                                                    acgtGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCGGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaa                                   <  1:1001951/235‑1 (MQ=255)                                                                                                                                                                                                     tttgtttcGTTAACGGCGTATTTTTTTTCC‑ATAAACGGTGACGGGGTTTTGCCAGATATATTTG‑GTCCGGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGGTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                  <  2:902023/231‑1 (MQ=255)                                                                                                                                                                                                     cTTGATTCTTTAACGGCGGATTTTAATTACTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGATGACTCCTTTCCCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                  <  1:46491/238‑1 (MQ=255)                                                                                                                                                                                                         aTGGGTTAACGGCGTATTTTATTGCCTATAAGC‑GTGACGGTGTTTAGCCAGATATATTGG‑GCCGCGA‑TGCCCTCTACGAGGACACCTTTACCTGGCCATGAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGAAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATgg                             <  1:828816/239‑1 (MQ=255)                                                                                                                                                                                                                     ggttttttA‑TTCCTATAAGC‑GTGTCGGGTTTTTTCCAGATAGAGTTG‑GTCGCAA‑GGTGCGTTACGACGTCTCCGTTACCTGGCCAGTAA‑TCAAATGTGA‑GCGATTTGCGCATATTCATCTGGCAAATGAGAATAAGCCAGTGCCGCCACAGTTGCGACAGTTCAGCAGAACGAATGGCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                 <  2:143322/234‑1 (MQ=255)                                                                                                                                                                                                                                   taaaaGG‑GTGACGGTTTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGa   <  1:148754/236‑1 (MQ=255)                                                                                                                                                                                                                                    ataAGC‑GAGACGGTGTTTTGACAGATATATTTG‑GTCGCGT‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTCGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCACTTCAGCCGAACGAATGCCGAAGCGCATTTGGTATAGTGTCAGGGGGCGTTTGTTGGGCCAGCCGGGTGGCTCCTTGCACTTCTTAAAACTGac  >  2:57832/1‑238 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                ACATGCTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGC‑GTGACGGTGTTTTGCCAGATATATTTG‑GTCGCGA‑TGCACTCTACGACGACTCCTTTACCTGGCCATTAA‑TCAAATTTGA‑GCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAA  >  NC_000913/252156‑252621

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