Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,434,882	0	T	G	32.6%	75.8 / 10.3	50	Y252S (TAC→TCC)	truA	tRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (3/30);  new base G (16/0);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 2.89e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGG  >  NC_000913/2434645‑2435118                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   aaaTGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGAc                                                                                                                                                                                                                                       <  1:329436/251‑1 (MQ=255)                 gACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑ATTTTGAA‑GAGGATACAGG‑TGCGGG‑GAAACCCCCGCGACC‑‑AGATAGAGACCTTcc                                                                                                                                                                                                                     >  1:329485/1‑251 (MQ=255)                    gCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGA‑GGTCATAACCG‑TCAGGG‑GAATCCACCCCGACC‑‑acaaaaaacccctcccct                                                                                                                                                                                                                  >  1:216862/1‑234 (MQ=255)                       gTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCtttt                                                                                                                                                                                                               <  1:230998/251‑1 (MQ=255)                            ccGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGcc                                                                                                                                                                                                            <  1:379652/249‑1 (MQ=255)                                          gCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGc                                                                                                                                                                                             <  1:145012/250‑1 (MQ=255)                                                           gtgCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGG‑‑GAACATCCCGG‑GCCGGG‑GAATACACCACGACC‑‑CGATACAAACCCTCCCCTTTTGTCGTGTCCGCCCCCCGC‑G‑TCCTGTCCTTTgacgcc                                                                                                                                                                         >  2:489591/1‑251 (MQ=255)                                                           gtgCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTgccgcc                                                                                                                                                                         <  2:216862/251‑1 (MQ=255)                                                             gCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAg                                                                                                                                                                       <  1:11336/251‑1 (MQ=255)                                                                  aTAAAATCGATGAGGAAATAAATCAGGACCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAAGA‑CAAAACGG‑GACGGG‑GGAGCCACCGCGGCC‑‑AGATACAACCCCTCCGCCTTTGCCGGTGGTGCCGCCCGC‑G‑GTATGTCCttttcacgccgcaac                                                                                                                                                                   >  1:153605/1‑238 (MQ=255)                                                                     aaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGA‑GATCACACCCG‑TCCGGG‑GAATCCCCCCCCACC‑‑AGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCGCCCGC‑G‑TTCTGTCCTTTTCCCCCAGCCActctt                                                                                                                                                               >  1:140063/1‑250 (MQ=255)                                                                          gATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑c                                                                                                                                                         <  1:321617/251‑1 (MQ=255)                                                                              aGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGc                                                                                                                                                     <  2:540265/251‑1 (MQ=255)                                                                               ggAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGA‑GAACATACCCG‑TCAGGG‑GAATCCACCCCGACC‑‑CCATAAAGCCCTTCCGCCTTTGTCGTTGGTGCCCCCCCC‑G‑TTCTGTCCTTTTCCCCCCGCAACCCTGccttcccgcc                                                                                                                                                     >  1:79985/1‑242 (MQ=255)                                                                                   ataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAA‑GATCATAACGG‑TCAGGG‑GAATCCACCGCGACCCAAGA‑ACA‑CCCTTTCGCCTTTTCCTTTGTTGCCGCCAGC‑G‑TTCTGTCCCTTgccgcccgaaacccggctttaccgctctcc                                                                                                                                                 >  1:396936/1‑227 (MQ=255)                                                                                     aaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCgg                                                                                                                                              <  2:22753/251‑1 (MQ=255)                                                                                         cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGTA‑AAACATAACGG‑GCCGGG‑GAATCCAACCCGGCC‑‑AGATACAACCCTTCCCCTTTTTCCGGTGGGGCCCCCCAC‑G‑GTCTGTCCTTTGCCGCCAGCAACCCTGCTAT‑CCAGCCCTcccgccgg                                                                                                                                          >  1:20741/1‑245 (MQ=255)                                                                                              ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGA‑AATCATACCGGGCCAGGG‑‑GATAAACCGCCACC‑‑AGATACAGCCCCTACGCTTTTGTCGGTGCTGCAGCCCGC‑G‑GTCTGTCCTTTGCCGCCAACAACACTGCTAT‑CCCACTCCCCGGCTGGGtgtg                                                                                                                                     >  1:81553/1‑251 (MQ=255)                                                                                                  atTTGCCTAGGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAAGGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGc                                                                                                                                  <  2:20741/250‑1 (MQ=255)                                                                                                                tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCa                                                                                                                   <  1:181541/250‑1 (MQ=255)                                                                                                                tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCa                                                                                                                   <  1:37114/250‑1 (MQ=255)                                                                                                                 gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAgg                                                                                                                 <  1:102774/251‑1 (MQ=255)                                                                                                                       gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCg                                                                                                           <  1:217084/251‑1 (MQ=255)                                                                                                                        cattagttCAGGCTGTAACTCAGTTGGCTACAGACTGAAATGTTCCGAGGTGTATGTGCTAGAGCCCGCCAGAAATAGCTGGCACATTGGCGG‑TTGTGGAATCAT‑ATACTGG‑TCAGGGGTTA‑CAACCGCGACC‑‑ATTTACAGCCTTTCCGCTTTTGCCTTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCACCAACTCTGCTAT‑CCAGCACTCCGGCTGGTTGTGGGCACCG‑TCTTCCATCAGGCTGCCGa                                                                                                          <  2:223786/245‑1 (MQ=255)                                                                                                                          aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑GTTTGGAA‑GGACATAACGG‑GCCGGG‑GAATACACCCCGACC‑‑AGATACAGCCCCTCCGCGTTTGCCCCTGCTGCCCCCACC‑GATTC‑GTCCTTTGCCGCCCGCAACCCCGCTAT‑CCCACTCTCCGGGTGGGTGCGGGGACCGCACT‑CCCTCAGGCTGCCGAc                                                                                                         >  2:111305/1‑250 (MQ=255)                                                                                                                            cTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGA‑GATCATACCGG‑GGCGGG‑GAATTCACCCCGGCC‑‑AGATACAACCCTTCCCCTTTTTCACTTGCTGCCCCCCGCCG‑TCCT‑TCCTTTTCCGCCAGCCACTCTGCTAT‑TCTGCTCTCCCGCTGGTTGTGGGCGCCG‑GCTTccctcccgccgcccgccat                                                                                                      >  2:66389/1‑234 (MQ=255)                                                                                                                              ttCAGGCTGTAACTTCGTTTGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCTG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCATCTCTCCGGCTGGTTTTGCGCACCG‑ACTTCCATCAGGCTGCCGACaa                                                                                                       <  1:424668/248‑1 (MQ=255)                                                                                                                                cAGGCTTTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTGGACCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATAtt                                                                                                   <  2:193144/250‑1 (MQ=255)                                                                                                                                   gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAA‑GATCATACCGG‑GCCGGG‑GTATCCAACCCGACC‑‑CGAAACAGCCCCTCCCCTTTTTCCGTTGCTGCCCCCCGC‑G‑TTATGTCCTTTTCCGCCAGCAACTCTGCTAT‑CCCGCCCTCCGGCTGGGTGTGGGCACCC‑CCTTCCCTCAAGCTGCCGACCAAATTccct                                                                                               >  1:308907/1‑249 (MQ=255)                                                                                                                                              gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATAACGG‑GCCGGG‑GAAACCACCGCGACC‑‑AGATACAGCCCCTCCGCTTTTTTCGGTGCTGCCGCCCCC‑G‑GTCTTTCCCTTTTCGCCAGCAACTCTTCTTT‑TCAACTCTCCGGCGGGTTGTGGGGACCG‑ACTTCCATCAGGGTGCCGACAAAATTCCTGACCATatgatg                                                                                    >  1:251510/1‑251 (MQ=255)                                                                                                                                                gggtggactccggaaaggttcaggtcATTTTTCGCGTTAGCCCGCCATACATAGCGGGCCCATTGGCGT‑TTTTGGAA‑GATCATACCGG‑TCATGG‑TAATCCACCGCGACC‑‑AGAGACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCTATCAGGCTGCCGACAATATTCCTGACCATAGGATGt                                                                                   <  2:381728/225‑1 (MQ=255)                                                                                                                                                  gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGA‑GATCATACCCG‑GCCGGG‑GAATCCACCGCGGCC‑‑CGAAAAAGCCCTTCCCCCTTTTCCCTTTCTGCCGCCACC‑G‑GTCTGTCCTTTTCCGCCCGCAACTCTGCTAT‑TCCGCCCTCCCGCTGGGTGTGGGGACCG‑ACTTCCATCCAGCTGCCGACAAACTTCCTGACCATATGATGTAc                                                                                 >  1:492759/1‑250 (MQ=255)                                                                                                                                                      cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAg                                                                            <  1:273484/251‑1 (MQ=255)                                                                                                                                                            gAAATGTTCCGAGCATTAGTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCGGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTc                                                                       <  1:74767/250‑1 (MQ=255)                                                                                                                                                             aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGA‑GAACATACCCG‑TCAGGG‑GAAACCACCCCGACCC‑AG‑TACAGCCCCTCCCCTTTTTCCGTTGCCCCCCCCAGC‑G‑CCCTGTTCTTTGTCCCCCCCACCTCGGCTAT‑CCAGCTCTCCCCCTGGGTTTGGGCACCG‑CCTTCCATCCAGCCGCCACCAATANTTCTGACCACATGATGTTAAAAAACATTccc                                                                     >  2:247602/1‑249 (MQ=255)                                                                                                                                                                   tCCGAGCATTATTCTCGTTAGTCCGCCAGANATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTaa                                                                <  2:251510/250‑1 (MQ=255)                                                                                                                                                                   tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAA‑GATCATACCCG‑GCAGGG‑GAATCCACCGCGACC‑‑AGATAAAACCCCTCCCCCTTTTCCGTTGCTGCCGCCCGC‑G‑CTCTGTCCTTTTCCGCCAGCAACTCCGCTTT‑CCCGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCCGACCACATGTTGTACCAAAACAtttccttttat                                                               >  1:257498/1‑242 (MQ=255)                                                                                                                                                                       aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑GTTTTGAA‑GATCATACCGG‑TCAGGG‑TAAACCACCCCGACC‑‑CCAAACAACCCTTCCGCTTTTGCCGCTGCTGCCCCCAGC‑G‑TCCTGTCCCTTGCCGCCAGCCAAACTGCTATACCA‑CTCTCCCGCTGGGTGTGGGCCCCG‑ACCTCCCCCAGGAAGCCGCCAACATTCCTTAACACATGATGttaaaaaaacttcgctttataaact                                                           >  2:530955/1‑227 (MQ=255)                                                                                                                                                                                  cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACcaccac                                                   <  2:435312/248‑1 (MQ=255)                                                                                                                                                                                   gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCCCCGCGACC‑‑AGATACAGCCCTTCCGCGTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATaa                                               <  2:485569/251‑1 (MQ=255)                                                                                                                                                                                      aGTCCGCCAGAAAGAGCGAGCCCATGGGAGG‑TTTTGGAA‑GATCAGTCCGG‑TCTGGG‑TGATCCACCGCGACC‑‑AGATCCAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAgg                                             <  2:328515/250‑1 (MQ=255)                                                                                                                                                                                         ccGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑GCAGGG‑TAATCCACCGCGACC‑‑AAATAAAACCCTTCCCCCTTTTCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCCCCAGCAACTCTGCTAT‑CCAGCTCTCCCGCTGGGTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAAAATTCCTGAACATATAATGTACAAAAGCATTCGCTTTTATATCTACCACCCCCTAAGGAccc                                         >  1:118803/1‑250 (MQ=255)                                                                                                                                                                                         ccGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCAACCCGG‑TCAGGG‑GAATCCACCGCGACC‑‑AGATACAGCCCCTCCGCTTTTGTCGTTTCTGGCGCCAGC‑G‑GTCTGGCCTTTTTCCCCCGCCACTCTGCTTT‑CCAGCTCCCCGGCCGGTTGTGGGGACCGGACT‑CCCTCAGGCTGCCGACCATATTCCTGACCCTATGAAGTACAAAAGCAATCGCTTTAATATTTAacaccaactaaagaacg                                         >  1:125840/1‑240 (MQ=255)                                                                                                                                                                                               gAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGc                                                                                                                                                     <  2:100302/138‑1 (MQ=255)                                                                                                                                                                                                    aGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGa                               <  2:369274/250‑1 (MQ=255)                                                                                                                                                                                                        ggCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑TCTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGTTTTAATATCTACCACCGCATAAGGACCGTTACGCGTGACATTa                          <  2:81553/251‑1 (MQ=255)                                                                                                                                                                                                                     ggtgtGGAA‑GATCATACCGT‑TCAGGG‑TAATCCCCCGCGACC‑‑AGTTACAGCCCTTCCGCTGTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCATCTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGTGCACCG‑ACTTCCATCAGGCTGCCTACAATATTCCTTACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGTACCGTGACGCGTGACATTAATATGCATAACg              <  2:153605/247‑1 (MQ=255)                                                                                                                                                                                                                       tttGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGt             <  1:395836/250‑1 (MQ=255)                                                                                                                                                                                                                        ttGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGtt            <  2:492759/250‑1 (MQ=255)                                                                                                                                                                                                                                  cATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAACGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGAAGCTTGACATTAATATGCATAACGTTGCGCCACggg  >  2:402580/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   AAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAA‑GATCATACCGG‑TCAGGG‑TAATCCACCGCGACC‑‑AGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGC‑G‑TTCTGTCCTTTGCCGCCAGCAACTCTGCTAT‑CCAGCTCTCCGGCTGGTTGTGGGCACCG‑ACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGG  >  NC_000913/2434645‑2435118

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