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

TTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGGGTTCGGGACT  >  NC_000913/2434636‑2435128                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             ttAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATACCGGGC‑AGGGGAATc                                                                                                                                                                                                                                                          >  1:385400/1‑251 (MQ=255)  tAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATc                                                                                                                                                                                                                                                          <  1:159188/250‑1 (MQ=255)     gAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑Acc                                                                                                                                                                                                                                                     <  2:77161/251‑1 (MQ=255)        cAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGTAGATAATTCCGGGC‑AGGGGAATAC‑Accccc                                                                                                                                                                                                                                                  >  2:373269/1‑248 (MQ=255)            aTGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGAGATTATACCGGTC‑AGGGGAATCC‑ACCGCGAccc                                                                                                                                                                                                                                              >  1:364682/1‑250 (MQ=255)                cATAAACCCAGGCGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATa                                                                                                                                                                                                                                          <  2:410944/251‑1 (MQ=255)                             gCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCt                                                                                                                                                                                                                             <  1:381808/251‑1 (MQ=255)                                         gACCAGCTCTGCCATATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGGTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTt‑gct                                                                                                                                                                                                                <  1:176870/251‑1 (MQ=255)                                         gACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTt‑gct                                                                                                                                                                                                                <  1:91420/251‑1 (MQ=255)                                            cAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGAGAATATACCGGTC‑AGGGGAAACC‑ACCGCGACCAGAAAAAGCCCCTCCGCTTTTTCCGGt‑gctggc                                                                                                                                                                                                             >  1:309645/1‑249 (MQ=255)                                            cAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑GTTTTGAAGATCATACCCGGC‑CGGGGAATCC‑CCCGCGGCCAGGTCCAGACCTTCCCCTTTTTCCGtttcctcc                                                                                                                                                                                                              >  2:366054/1‑245 (MQ=255)                                             aGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATACCAGTC‑AGGGGAAACC‑CCCCCCAACAAATACAACCCTTCCCCttttttcgttcgtgtcc                                                                                                                                                                                                             >  2:250405/1‑238 (MQ=255)                                                   gCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACTGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAg                                                                                                                                                                                                        <  1:294390/249‑1 (MQ=255)                                                              gTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCttt                                                                                                                                                                                           <  2:55825/251‑1 (MQ=255)                                                              gTCAATGGGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGA‑GCACATACCCGGC‑AGGGGAATCCCACC‑CCACCACAGAAAGCCCCTCCCCTTTTTCCGTT‑GTTgccgccccactgctggtcttt                                                                                                                                                                                           >  2:337591/1‑236 (MQ=255)                                                                  aTGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTgccg                                                                                                                                                                                       <  2:197284/251‑1 (MQ=255)                                                                      gCAGGATAAAATCGCTGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTTTAGATCATACCGGTCCAGGGGAAACC‑‑CCCCGGCCCGAAACAACCCTCCCCCCTTTTCCGTT‑GCTGCCTCCAGCGTTCTGTCCTTTgccgcccc                                                                                                                                                                                   >  2:273856/1‑249 (MQ=255)                                                                             aaaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAActc                                                                                                                                                                             <  1:373269/250‑1 (MQ=255)                                                                                       aGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGAACATAACGGTC‑AGGGGAATTC‑ACCGCGACCAAATACAACCCTTCCGCTTTTTCCCGT‑GCTGCCGCCCGCCGTCTGTTCTTTTTCGCCCGCAACTCTTCCATCCAGc                                                                                                                                                                  >  1:379474/1‑251 (MQ=255)                                                                                          aaataaatCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTG‑AGGGGAATCC‑ACCCCGACCCGACACAGCCCCTCCGCTTTTTCCGTT‑GGTGCCGCCCGCGGTCTGTCCTTTTCCGCCACCAAACCTGCTATCCAGcctt                                                                                                                                                               >  1:384190/1‑248 (MQ=255)                                                                                                atCAGGTCCATATATGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCt                                                                                                                                                        <  2:233729/251‑1 (MQ=255)                                                                                                    ggTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATACCGGGC‑AGGGGAAACC‑ACCGCGAACAGATACAGCCCCTCCGCTTTTGCCCGT‑GCTGCCGCCAGCGGTCTGTCCTTTTCCGCCAGCAAATCTGCCATCCAGCTC‑TCCGGccgggt                                                                                                                                                    >  1:284708/1‑246 (MQ=255)                                                                                                    ggTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGtt                                                                                                                                                    <  2:336912/251‑1 (MQ=255)                                                                                                          tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGA‑GATCATACCGGGC‑AGGGGAATCC‑ACCGCGGCCAGAAACAGCCCTTCCGCTTTTTCCGTT‑GCTGCCGCCAGCGGTCTTTCCTTTGCCGCCCGCAACTCTGCTATCCAGCCC‑CCCGGCGGGTTGT‑gggg                                                                                                                                             >  2:223171/1‑250 (MQ=255)                                                                                                                   aTGCATTGTACTCGTCATCTTTCATGCTGTAACTACGTTGGCTACATCCTGAAATGTTCCGAGCATTATGCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCGTTTGCCGTT‑GCTGCCGCCAGCGTTCTTTCCTTTGCCGCCAGCAACTCTGCTATCCTTCTC‑CCCGGCTGGTTGT‑GGGCACCTACtt                                                                                                                                     <  1:315297/250‑1 (MQ=255)                                                                                                                            aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGc                                                                                                                            <  1:4891/250‑1 (MQ=255)                                                                                                                             cTCGTCATCTTTCAGGCCGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATAACGGTC‑AGGGGAATCC‑ACCGCGCACAGATAAAGCCCTTCCCCTTTTGCCGTT‑GCTGCCCCCAGCGTTCTGTCCTTTGCCGCCAGCAACCCTGCTTTCCAGCTC‑TCCCGCTGGGTGT‑GGGCACCGACCTCCATCAAGCTg                                                                                                                          >  1:89595/1‑251 (MQ=255)                                                                                                                                                  aCTACGGTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGAC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCCGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCat‑a                                                                                                   <  1:6418/251‑1 (MQ=255)                                                                                                                                                       gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGGTAATACCGGTC‑AGGGGAATCC‑ACCCCGACCAAATACAGACCTTCCGCTTTTTCCGGT‑GCGTCCGCCAGCGGTCTGTCCTTTTCCCCCAACAACACTGCTATCCCGCTC‑TCCCGCCGGTTGTTGGGCACC‑ACCTTCATCCGGCTGCCGACAATAtttcctaacatttgaat                                                                                                >  2:303916/1‑236 (MQ=255)                                                                                                                                                           gCTACAGCCTGAAATTTTCCGAGCATTATGCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATCCCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑TGGCACCGACTTCCAGCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACa                                                                                          <  2:239922/251‑1 (MQ=255)                                                                                                                                                           gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACa                                                                                          <  1:141299/251‑1 (MQ=255)                                                                                                                                                               cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCAGACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGTTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATAATCCT‑GACCAT‑ATGATGTACAAAAg                                                                                      <  1:346832/251‑1 (MQ=255)                                                                                                                                                                  ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCAt                                                                                   <  2:385400/251‑1 (MQ=255)                                                                                                                                                                        aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGAACATACCGGGC‑TGGGGAATCC‑ACCGCCAACAGATACAGCCCTTTCGCCTTTTCCGTTTGCTGC‑GCCCGCGGTCTGTCCTTTTCCCCCCGCACATCCCCTTTCCAGCTC‑TCCGGCTCGGTGG‑GGGCGCCCACTTCCCTCAGGCCGGCGGCAATATTTCT‑TACCAT‑ATAATGTACCAAAGCAATCGCtt                                                                             >  2:174155/1‑251 (MQ=255)                                                                                                                                                                                   attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGGTCATAACGGGC‑ACGGGAATCCCAAC‑CGACCAGATACACCCCCTCCCCTTTTTCCGTT‑TCTGCCGCCCCCGTTCTTTCCCTTTCCGCCCGCAACTCCTCTATTCCCCTCCTCCG‑CTGGTTGT‑GGGCCCCCACCTCCCTCAGGCTGCCGCCCATATTTCT‑GACCATTAT‑ATTTACAAAAACATTCGCTTTTATATCTaac                                                                  >  2:419734/1‑249 (MQ=255)                                                                                                                                                                                        tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTAccacca                                                              <  1:349874/250‑1 (MQ=255)                                                                                                                                                                                           cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACcaccac                                                             <  1:146662/248‑1 (MQ=255)                                                                                                                                                                                           cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACAt                                                           <  1:430664/250‑1 (MQ=255)                                                                                                                                                                                                      cAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCCTTACGGTC‑CGGGGTATCC‑ACCGCGACCAGATACCACCCCTCCGCTTTTTCCGTT‑TGTGCCCCCCGCGGTCTTTTCTTTTTCGCCCGCAACTCTGCTATTCCGCTC‑TTCCGCTGGTTTT‑TGGCCCCGACTTCCCTCCGGCTGTCTACCATATTCCTTCACC‑T‑ATTATTTACAAAAGCATTTGCTTTTATATCtttccccccacaatgcacctgc                                                >  1:282665/1‑229 (MQ=255)                                                                                                                                                                                                             aGCGGGCCCATTGGCGG‑TTTTTGAAGATCATACCGGGC‑AGGGGAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTT‑GCTTCCGCCAGCGTTCTGTCCTTTTCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAAAATTCCT‑GACCCT‑ATGATGTACAAAAGCATTCGCTTTAATATCTAACACCACATAAGGACCGTGACGCGTGAc                                        >  1:83482/1‑251 (MQ=255)                                                                                                                                                                                                                 ggCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTa                                    <  1:218122/251‑1 (MQ=255)                                                                                                                                                                                                                        tGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGc                              <  1:222131/250‑1 (MQ=255)                                                                                                                                                                                                                        tGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGGAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGCCAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGc                              >  2:355695/1‑250 (MQ=255)                                                                                                                                                                                                                                       gATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAAGATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCAc               <  1:367606/251‑1 (MQ=255)                                                                                                                                                                                                                                           aTACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACgggg           >  1:391447/1‑251 (MQ=255)                                                                                                                                                                                                                                               cGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGGGTTCg       >  2:150309/1‑251 (MQ=255)                                                                                                                                                                                                                                                ggTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATGTGCATAACGTTGCGCCACGGGGTTCgg      >  2:362758/1‑251 (MQ=255)                                                                                                                                                                                                                                                     aGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGGGTTCGGGACt  >  1:228688/1‑251 (MQ=255)                                                                                                                                                                                                                                                       ggTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCCGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAat                                                                         <  2:298860/178‑1 (MQ=255)                                                                                                                                                                                                                                                       ggTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCCGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAat                                                                         >  1:298860/1‑178 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             TTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTC‑AGGGTAATCC‑ACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTT‑GCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTC‑TCCGGCTGGTTGT‑GGGCACCGACTTCCATCAGGCTGCCGACAATATTCCT‑GACCAT‑ATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTGCGCCACGGGGTTCGGGACT  >  NC_000913/2434636‑2435128

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