Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,434,882	0	T	G	41.9%	100.7 / 10.6	62	Y252S (TAC→TCC)	truA	tRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (7/29);  new base G (26/0);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 2.41e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434740‑2435031                                                                                                                                                 |                                                                                                                                                       atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAT‑ACCGGGCC‑GGGGAATCCAcc                                                                                                                                                >  1:327467/1‑151 (MQ=255) atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAA‑ACCGGGCA‑TGGGAAaccacc                                                                                                                                                >  2:817563/1‑151 (MQ=255)   atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAT‑ACCGGGCA‑GGGGAATCCAcccc                                                                                                                                              >  1:677961/1‑149 (MQ=255)   atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGAAGATCAT‑AACGGGCA‑GGGGAATCCACcgc                                                                                                                                              >  1:945237/1‑151 (MQ=255)    tttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAT‑ACCGGGCA‑GGGGAATCCAcccc                                                                                                                                              >  1:992563/1‑148 (MQ=255)      tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAT‑ACCGGTCA‑GGGGAATCCACcgc                                                                                                                                              >  1:989347/1‑148 (MQ=255)        ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCa                                                                                                                                         <  1:633278/151‑1 (MQ=255)           aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAT‑ACCGGGCA‑GGGGAATCCACCGCGAcc                                                                                                                                          >  1:844265/1‑147 (MQ=255)             gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑a                                                                                                                                    <  1:64896/150‑1 (MQ=255)             gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑Ac                                                                                                                                   <  1:677444/151‑1 (MQ=255)              ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTTGCAGGGTAA‑AACGGCCAGGGGTAATCCAcccggcccaata                                                                                                                                      >  2:142493/1‑139 (MQ=255)               cTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCAT‑ACCGGTCC‑GGGGAATACAcc                                                                                                                                                >  2:1143826/1‑137 (MQ=255)                 tGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGc                                                                                                                                <  1:350918/150‑1 (MQ=255)                  gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACCA‑AACGGGCA‑GGGGAATCCACCCCCACCCGAT‑AAAGcc                                                                                                                               >  1:837478/1‑150 (MQ=255)                        gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGAGCCCATTGGCGGTTTTGGAAGATCAT‑GCCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGc                                                                                                                        <  2:774241/151‑1 (MQ=255)                           aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGAAGAAACT‑ACCCGTCA‑GGGGAAACCACCGCCACCCGAG‑ACCACCCtttcccct                                                                                                                      >  2:703428/1‑144 (MQ=255)                             cTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCAT‑CCCGGTCC‑GGGGTATCCACCGCGACCCCAT‑AAAAGCCCTCCCCttttt                                                                                                                   >  1:1272035/1‑150 (MQ=255)                              tttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTGGAAGATTATAACCGGT‑A‑GGGGAATCACCCGCGGCCAGAT‑ACAGGCCTTTCGCcttt                                                                                                                    >  2:283451/1‑149 (MQ=255)                                 cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAAATCAT‑ACCCGGCA‑GGGGAATCCACCCCGACCAGAT‑ACAGCCCTTCCGCttttttcg                                                                                                                >  2:186627/1‑146 (MQ=255)                                     cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATTAT‑ACCCGGCA‑GGGGAATACACCCCGGCCAGAT‑ACAGCCCTTCCTCCTTTTCCGTtgc                                                                                                            >  1:741635/1‑150 (MQ=255)                                       gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCCT‑ACCGGTCA‑AGGGAATACACCCCCACCCGAT‑ACAGCCCTTCCGCCTTTTCCGTtgct                                                                                                           >  2:838346/1‑149 (MQ=255)                                         aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTgccg                                                                                                       <  1:1080607/151‑1 (MQ=255)                                               gTTGGCAACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATCAT‑AACCGTCC‑GGGGAAACCAACGCGCCCCGAT‑TAAGACCCTCCGCTTTTTCCGGTGCTGTCGCCAGc                                                                                                  >  2:791909/1‑150 (MQ=255)                                                  ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCAT‑AACGGGCA‑AGGGAATCCACCCCGGCCAGAT‑ACAACCCCTTCGCTTTTGCCGTTGCTGCCCCCCACGtt                                                                                               >  1:164039/1‑150 (MQ=255)                                                     tACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACCA‑ACCCGTCA‑GGGTAATCCCCCCCGACCAGATCACGCCCCTTCGCCTTTTTCC‑TTTTTGCCGCCAGCGTTCt                                                                                             >  2:101309/1‑149 (MQ=255)                                                        aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTAA‑ACCCGTGC‑TGGGAATACACCCCGGCCCCAT‑TCAGCCCCTCCGCCTTTTCCGCTGCTTGCCCCCGCGTTCCGttc                                                                                         >  2:330060/1‑148 (MQ=255)                                                        aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCt                                                                                        <  1:165421/151‑1 (MQ=255)                                                          ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAT‑ACCGGGCA‑AGGGAACCCCCCGCCCCCCGAT‑ACAGCGCTTCCGCTTTTGCCGTTGCCGCCGCCAGCGGTCTGTGCttt                                                                                      >  2:146002/1‑151 (MQ=255)                                                          ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAT‑ACCCGGCA‑GGGGAATACACCCCGACCAGAT‑AAAGCCCTTCCCCCTTTTTCGCTGGTGGCGCCAGCGTTCttttcctt                                                                                      >  2:1024050/1‑144 (MQ=255)                                                          ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                                      <  1:1114050/151‑1 (MQ=255)                                                            tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAT‑ACCCGTCA‑GGGGAATCCCCCGCGACCAGAT‑ACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCctt                                                                                      >  1:142799/1‑149 (MQ=255)                                                            tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCAT‑ACCCGTCA‑GGGTAATCCCCCGCGACCAGGA‑ACAACCCTTCCGCTTTTGCCGTTGCTGCCCCCAGCGGTCCGTCCtct                                                                                      >  1:686794/1‑147 (MQ=255)                                                             gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgcc                                                                                   <  1:48328/151‑1 (MQ=255)                                                               aaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAT‑ACCGGGCA‑GGGGAAACCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTTCCGGTTCTGCCGCCCGCGTTCTGTCCTTTgcc                                                                                   >  1:779898/1‑149 (MQ=255)                                                                   ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAT‑ACCGGTCA‑GGGGAATCCACCGCGACCCGAT‑AAAACCCCTCCCCTTTTTCCGTTTCTGCCGCCAGCGGTCTGTCCTTTgccgcc                                                                                >  1:717334/1‑148 (MQ=255)                                                                       gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAAct                                                                         <  1:1037645/151‑1 (MQ=255)                                                                            ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTAT‑AACCGGCA‑GGGGAATCCACCGCGACCAAAT‑AAAGCCCTTCCGCTTTTTCCGTTGCTGCCCCCAGCGGTCTTTCCTTTTCCCCCAGCAAc                                                                          >  1:309025/1‑145 (MQ=255)                                                                               ttCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCAT‑ACCGGTCA‑GGGGAATCCACCGCGACCAGAT‑ACAGCCCTTTCCCTTTTTCCGGTTCTGCCGCCCGCGGTCTGTCCTTTTCCGCCAGCAACTCTGCTAt                                                                  >  1:43226/1‑150 (MQ=255)                                                                                tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATc                                                                 <  1:1248386/150‑1 (MQ=255)                                                                                 ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGGCA‑GGGGAATCCACCGCGACCAGAT‑ACAGCCCTTCACCTTTTGACGTTGGTGCCCCCACCGGTCTTTGCTTTGCCGCCAGGCACACTGCTAtac                                                                >  2:505791/1‑148 (MQ=255)                                                                                      tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGctc                                                           <  1:93654/150‑1 (MQ=255)                                                                                                  aaTAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgt                                              <  2:946145/151‑1 (MQ=255)                                                                                                        gggCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTGTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGATATCCAGCTCTCCGGCTGGTTGTGTGCAc                                        <  2:842630/151‑1 (MQ=255)                                                                                                              aTTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACtt                                  <  1:1208676/151‑1 (MQ=255)                                                                                                                tGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                            <  2:821712/139‑1 (MQ=255)                                                                                                                tGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                            >  1:822252/1‑139 (MQ=255)                                                                                                                  gCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt                              <  1:789489/151‑1 (MQ=255)                                                                                                                   cGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                             <  2:99360/151‑1 (MQ=255)                                                                                                                   cGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                             <  2:142738/151‑1 (MQ=255)                                                                                                                    ggTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa                            <  1:626872/151‑1 (MQ=255)                                                                                                                       tttctaaGATCAT‑ACCGTTCA‑GGATAATCCAACGCGACCAGAT‑ACAGCCCTTCGGCTTTTGCCGTTGCTGACCCAAGCGTCCTGTCCTTTGCCGCCAGCAACTCTGATATCCTGCTATCCGGCTGGTTGTGGGAACCGGCTTCCATCAgg                          <  2:249412/145‑1 (MQ=255)                                                                                                                       tttGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCCCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc                         >  2:1000338/1‑151 (MQ=255)                                                                                                                                 cAT‑ACCGGTCA‑GGGTAATCCACCGCGCCCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGGCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCCCCGGCGGGTTGTGTGCACCGACTTCCATCAGGCTGCCGACaa                <  2:46511/150‑1 (MQ=255)                                                                                                                                  aT‑ACCTGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAata              <  2:329374/151‑1 (MQ=255)                                                                                                                                  aT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAata              <  1:838484/151‑1 (MQ=255)                                                                                                                                  aT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat               >  1:765326/1‑150 (MQ=255)                                                                                                                                     aCCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt            <  1:408901/151‑1 (MQ=255)                                                                                                                                       cGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTcc          <  2:59364/151‑1 (MQ=255)                                                                                                                                          tCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg        >  1:358800/1‑150 (MQ=255)                                                                                                                                           cA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAc      <  1:1144489/151‑1 (MQ=255)                                                                                                                                            a‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc     <  1:1060215/151‑1 (MQ=255)                                                                                                                                                gTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCata  <  1:876493/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                       ATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAT‑ACCGGTCA‑GGGTAATCCACCGCGACCAGAT‑ACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434740‑2435031

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