Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,197	0	G	T	36.1%	91.8 / 15.0	80	intergenic (‑32/+3)	leuV/leuP	tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (17/29);  new base T (1/25);  total (18/62)
Fisher's exact test for biased strand distribution p-value = 1.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCC  >  NC_000913/4606047‑4606340                                                                                                                                                        |                                                                                                                                                 aaGGAAGAAACAAGAACGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                  <  4:1119171/151‑1 (MQ=255) aaGGAAGAAACAAGAAAGGAAGTAAATATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                  <  8:1237789/151‑1 (MQ=255)   ggAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGCGGGGGGGC‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGtt                                                                                                                                                <  6:344236/151‑1 (MQ=255)   cgacgaaACAACAAAGGACGTAAAGATAATTGGTGCGAGGGGGGGG‑CCTGGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGtt                                                                                                                                                <  1:799696/147‑1 (MQ=21)    caagaaACAACAAAGCAATTAAAGAAAATGGGGGCGAGGGGGGGG‑ACTCCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                               <  5:679315/150‑1 (MQ=21)    caacaaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                               <  8:1110911/147‑1 (MQ=255)    caacaaaCAACAACGCAAGTAAAGATATTGGGGGCGCGGGGGGGC‑ACTTCACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                               <  5:671896/147‑1 (MQ=18)      agaaACAAGAAAGGAAGTAATGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTgg                                                                                                                                             <  8:1303554/151‑1 (MQ=255)      agaaAAAAGAAAGGAAGTAAACATAATTGGTGCGCGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTgg                                                                                                                                             <  3:1294310/151‑1 (MQ=255)       caaaCACGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                                                                                                            <  5:212817/150‑1 (MQ=255)        atatcattaaGGATGTAATGATATTTGGTTGTTGGGGTTGT‑ATTTCAACCCCATCTTTCTTAAGGACACTAACTACTGCAGCTAGCGCGTCTACCAATTCCGCCAGCTTCGCACATTCATCTTACTTTTTTTGATATCGCGTCGTTTGCTg                                                                                                                                           <  8:532763/143‑1 (MQ=12)        aatcAAGAATGCAAGTAAATTTATTTGGTGCTAGGGGGGGG‑ACTTTAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                           <  1:891045/148‑1 (MQ=18)        aaaCAAGAAAGGAATTAAAGATAATTGGGGCGCGGGGGGGGCACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                                                                                                            <  3:1063262/151‑1 (MQ=255)        aaaCAAGAAAGGAAGTAAAGATAATTGGGGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                           <  8:440555/151‑1 (MQ=255)           caacaaCGGAAGAAAAGAAAATTGGGGCGCGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                                                                                        <  7:1055223/151‑1 (MQ=21)            aacaacggAAGTAAAGATAATGGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑g                                                                                                                                      <  3:129751/145‑1 (MQ=255)               aaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑gggg                                                                                                                                   <  2:1331536/151‑1 (MQ=255)                 aGGAAGTAAACATAATTGGGGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑gggggg                                                                                                                                 <  3:1448249/151‑1 (MQ=255)                  gcaaGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑ggggggg                                                                                                                                <  2:1545610/149‑1 (MQ=255)                  gcaaGAAAAGATATTTGGGGCGCGGGGGGGC‑ACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑ggggggg                                                                                                                                <  5:1166848/149‑1 (MQ=21)                    aaGTAAATATAATTGGTGCTTTTGGTGTT‑CGTTGAAGCCCCACGTCCTTAATTACACTATCACCTGAAGCTATGGCGTCTACCAATTCCGCCAGCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTTGTTTGGTGCGA‑gggggggg                                                                                                                               <  7:1187567/150‑1 (MQ=255)                     aGTAATGATAATTGGTGGGCGGGGGGGG‑ACTGGACCCCCCACGTCCGTAAGGACACTAACACCTGCAGCTAGCGCGTCTACCAATTCCGCCACCTTCCCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGAc                                                                                                                             <  2:738138/151‑1 (MQ=25)                           aaaacttGGGGCGTTGGGGGTG‑ACTACAACCCCCACGTGCGAAAGTACACTAACACCTGAAGCTCCCCCTTAAACCAATTCCGCCACCTAACCACACTCATCAAACTTTTTATGATCTCCCCACGTTTGGTGCGA‑GGGGGGGGACTTGAAc                                                                                                                       <  5:1543190/146‑1 (MQ=12)                                   gggCGTGGGGGTGG‑ACTGGAACCCCCACGTGCTTATCTACACTACCACATGAAGATATCCATTATTCCTATTCCGGCAGCTTACCACAGTCAACTTACTTTTTTTGATCTAGCCTCCTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGt                                                                                                               <  2:940478/149‑1 (MQ=255)                                     gctatggggtgg‑ACTTGAACCCCCACTTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTcc                                                                                                             <  2:961131/146‑1 (MQ=25)                                      cGAGGGGGGGT‑ACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTACCGCGTCTACCAATTCCGCCACCTACGCACAGTCATCTAACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCg                                                                                                            <  6:282764/151‑1 (MQ=255)                                            gggttgccTTCAACCCCCACGGCCTTAATGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTACCCACAGTCATCTTACTTTTTTTGATATCCCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGa                                                                                                       <  5:533142/144‑1 (MQ=25)                                             gggtccctccaccccccACGTCCGTAAGTACACTAACACCTGATGCTACCGCTTATACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGAAATCCCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGaa                                                                                                      <  1:1162983/139‑1 (MQ=25)                                                  ccTTCAACCCCCACTTCCTTAAGTACACTAACACCTGAAGCTACCGCGTATAGCAAATCCGCCACCTTCCCACAGTCATCTTACTTTTTTTGATATCCCGTCCTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                                  <  6:1138731/149‑1 (MQ=25)                                                       aatccccACGTCCTTTAGTAGAGTAACTCCTGATGATTGCTAGTATTCTTATTCCGCGACCTTCTAACATTCATCTAACTTTTTTTGATATCCCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                                                            <  1:324859/148‑1 (MQ=12)                                                        ccccccACGTCCGTAATGACAGTAACACCTGAACCTACCACTTATACAAATTCCGCCACCTACAAACAGGCATCTAACTTTTTTTGATTTCCCCTCTTTTGGGGCGC‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                                                           <  5:1133489/150‑1 (MQ=21)                                                        acccaaacctcctaaaGCACACAATCTCCTCAACCTCGCGCGTCTACCTATTCCGGCATCTAATAACATTCATCATACTTTTTTTGATTTATCGTATTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                                                           <  8:587101/137‑1 (MQ=9)                                                         cccccACGGCCTTAAGGACACTAAAACCTGAAGCTAGCTCGTCTACCAATTCCGCCACCTTCTCACAGTCATCTTACTTTTTTTGATCTCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                                                          <  3:991583/151‑1 (MQ=37)                                                           ccgactgccGTAAGGACACTAACACATGAAGCTAGCACGTATACAAATTCCGCCACCAACCCACAGGCATCTTACTTTTTTTGATTTCCCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                                                                        <  3:497614/144‑1 (MQ=21)                                                              actgcctttattacagTAAAACCTGAAGATATCTCGTATTCCAATTCCGCCACCTTCCAACACTCAACTTACATTTTATGATCTACACTGGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                                                                     <  3:1424067/139‑1 (MQ=12)                                                               cGTCCTTAAGTACAGTAACACCTGATGCTAGAGCGGCTACCAATTCCGGCACCATCCCACATGCATCATACTTTTTTTGATATCCCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                    <  2:1592841/151‑1 (MQ=21)                                                                 gccctAAGTACAATAACACCTGCTCCTAACACGTAAACAAACTCCGCCACCTAAAAACATTCATCTAACATTTTTTGAAATACACTCCTTTGGTGCGC‑GGGGGGGGACTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                                                                  <  6:1316042/147‑1 (MQ=9)                                                                  ccGTAAGGACACTAACCCCTGACCCTCCCGCTGCTACCAATTCCGCCACCTTCGCACAGTCATATAACTTTTTTGTATTACCCCTCGTTTGGGGCGT‑GGGGGGGGACTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                 <  5:90156/151‑1 (MQ=21)                                                                     tAAGTACAGTAACTCCTGATGCTATTGCTTCTACCTATTCCGCTACTTTCTCACATTCATGTTACTTTTTTTGATATCGCCTCTTTTGGTGGGA‑GGGGGGGTACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                                                               <  4:1590323/150‑1 (MQ=18)                                                                     tAAGGACACTAACACCTGAAGATAGCGCGTCTACCAATTCCGCCACCTTCGCACAGGCATCTTACTTTTTTTGATATCGCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcgc                                                                                <  2:454691/149‑1 (MQ=255)                                                                      aactacaatacatccTGAAGATATCGCGTATTCCAATTCCGCCAGGTACCCACAGTCATCTTACTTTTTTTGATTTAGCCTCATTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                                                             <  2:770553/138‑1 (MQ=18)                                                                        gcacacTAACACCTGAAGCTAGAGCGTCAACCAAATCCCCCCCCTACCAACATCCATCTAACTTTTTTTGAAATCCCCTCTTTTGGGGCGC‑GGGGGGGCACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                                                            <  5:806561/148‑1 (MQ=18)                                                                            attAACACATGAAGCTCGCTCTGATTCCTATTACGCTACCTTCTCTCAGTGATCATACATTTTATGATATACCCATTTTTGGTGCGA‑GGGGGGGGACTTGAACCCGCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                       <  7:830977/149‑1 (MQ=9)                                                                             atacaacatGAAGCTAACACGTCAACCCAATCAGCCACCCAACCACAAGCATCTAACATCTATGCATATACCCACATTTGGTCCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                                                      <  5:985728/143‑1 (MQ=2)                                                                               aaCACCTGAAGATACCCCTTATACCAATTCCCCCACAAAACCACACTAATATTACTTTTTTTTATATCCCCTTTTTTGGTGCGT‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                                    <  2:893454/151‑1 (MQ=14)                                                                                 aacCTGAAGCTCGCCCGTCTACCCATTCCGCCACCTACCAACAGTCATCTTACTTTTTTGGATCTCCCCTCGATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                                  <  6:578497/150‑1 (MQ=18)                                                                                  acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTAATGATATCGAAGCGTTTAGTGCGA‑GGCGAGGGGCTTGCACCACCACGTCCGTAAGAACACAAACACCTAAAGCTAGAGCCTCTAACAATTACg                                                                  >  4:535384/1‑150 (MQ=21)                                                                                      taatctttaccttCTACCAATTCCGGCAGCTTTTCAGATTTATCTTACTTTTTTTTATTTATAGTTTTTTGGTGCGT‑GGGGGGGGACTTTAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                              <  7:737119/138‑1 (MQ=2)                                                                                              cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTATTGATATCGAATCGTATGGTGCGA‑GGCGGGGGACTACAACCACCACGTCCGAAAACACACAAACACCTGAAGCTAGAGCGTCTAACAATTAATCAACCTTCGcccc                                                     >  8:1253990/1‑150 (MQ=21)                                                                                                agtattctaaTTCCGCCAGCAAAACACAGGCATCTAACTTTTATTGATATCCCCTCATTTGGGGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                   <  2:602534/143‑1 (MQ=17)                                                                                                  tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGACTCGTTTGGTGCGA‑GGGGTGGGACTTGAAACACCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCACAAccgaagcgcacacc                                                 >  3:947999/1‑139 (MQ=21)                                                                                                   cTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGGGA‑GGGGTGGGACTTGCACCTCCACGTCCGTAAGAACACTAACACCTGAAGCTACAGCGTCTACCAATTCCGCCACCTACGcacacggcc                                                >  5:96170/1‑146 (MQ=21)                                                                                                     aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTAATTGATATCGACGCTATTGGTGCGA‑GGCGGGTGCATTTCACCACCACGTCAGTAATAAAACTAACACCTAAAGCTACATCGTGTAACAATTACGCAACCtaagctaacgagc                                                >  7:1568041/1‑137 (MQ=255)                                                                                                       cAATTCCGCCACCTTCGCACAGTCAGCTTACTTTTTTTGATATCGACGCGGTTGGTGCGA‑GGGGGGGGCATTGAACCACAACGTCCGTAATAACACTTACACCTGCCGCTAGCGGGTCTACCAATTACGCCACCGTCGcaaacgagcagg                                             >  1:713107/1‑139 (MQ=255)                                                                                                        aaTTCCGCCACGATCTCACATTGATCTAACATGAATGGATTTCCCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCa                                            <  4:788197/150‑1 (MQ=14)                                                                                                        aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAggtcgatct                                            >  5:1167420/1‑142 (MQ=255)                                                                                                         aTTCCACCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGACGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCaccggcacctat                                          >  2:1622746/1‑140 (MQ=21)                                                                                                              gCCACCTTCGCACAGTCATCATACTTTTTTAGATATCGCCTCGTTTGGTGCGA‑GGCGCGGGACTTGAAACCCCACGTCCGTAAGAACAATAACACATGAAGCTAGAGCGTCTAACAATTACGCCACCGTCGcacacgtccataaaaaaaa                                      >  4:19357/1‑135 (MQ=255)                                                                                                              cctacctaaattcattCATCTAACTTTTTTTGAATTACCTTTTTTTGGTGCGT‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa                                      <  1:1146029/135‑1 (MQ=17)                                                                                                               ggattaaaaaacaTTCATCAAACATTTATTTATTTCCCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                    <  1:1074666/141‑1 (MQ=17)                                                                                                               ggatcttataacaATCATCTAACTTTATTGGATTTACACTTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                    <  4:1230300/141‑1 (MQ=14)                                                                                                                cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAACGGCGATAAATATaaa                                     >  2:289106/1‑150 (MQ=18)                                                                                                                  gattattaaacgcATCTAACTTTTTTTGATTTACCCTATTTTGGTGCGT‑GGGGGGGGACTTTAACCCCCACTTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACg                                   <  2:1173375/137‑1 (MQ=17)                                                                                                                     acaaacaTTCAACATACATTTTTTGATCTCCCCTTCTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                              <  6:996471/147‑1 (MQ=18)                                                                                                                      aactcattcaagaaaCTTTAATTGATGTATCCTCTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttt                             <  2:695679/137‑1 (MQ=17)                                                                                                                       acacaATCATCAAACTTTATTTGATTTCCCGATTATTGGTGCGA‑GGGGGGGGACTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                            <  5:425606/150‑1 (MQ=17)                                                                                                                         acaGGCATCTAACTTTATAGGATATACCGTTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                          <  5:997075/151‑1 (MQ=17)                                                                                                                          catttattaaaCTTTTTTTGATTTTACTTTTTTTGGTGCGA‑GGGGGGGGAGTTTAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                          <  1:570700/140‑1 (MQ=17)                                                                                                                           agggatgtaaCTTTTTTTGATCTCGAGTATTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                        <  6:447555/142‑1 (MQ=18)                                                                                                                             ttaatttACTTTTTTTGATTTATACTTTTTTGGTGCGT‑GGGGGGGTACTTTAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                       <  1:1018698/145‑1 (MQ=18)                                                                                                                             ggatgaaaattttttgGATATACACTGCTTTGGTGCGA‑GGGGGCGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc                      <  6:82381/142‑1 (MQ=18)                                                                                                                                agtaacattttttGATTTACCCATGTTTGGGGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAg                   <  3:555769/144‑1 (MQ=18)                                                                                                                                   aacattttttGATGTCCACATGATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg                <  3:647870/147‑1 (MQ=18)                                                                                                                                       ttttttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACGCAGACCAATATAAACGACGTTATACGACCGAggggggtt            >  5:448159/1‑149 (MQ=17)                                                                                                                                           ttGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGAGGGAGGAATATAAAAGAGTTTTTTGCTCCGAgggggggtcctt        >  3:500265/1‑146 (MQ=17)                                                                                                                                             gATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGGCAGCTCTATAAAAGTTTTTTTAGGAGCGAGGCGGGGTAGTTgg      >  1:526061/1‑150 (MQ=17)                                                                                                                                               tatCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGACTCCCAGATATATCAACGAGTTTATTGGTGCGAgggggtgtcgtgggcc    >  3:1021168/1‑140 (MQ=17)                                                                                                                                               tatCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACAGCCAGCAATATCAACGTGGTTTTTCGACCGAGGGGGCGGAATTGaa     >  4:796438/1‑150 (MQ=18)                                                                                                                                                atCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACCCCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGACCCACGAATATAAAAGTATTTATTCCAGCGAGGGGGGGGACATgcgcc   >  6:1239546/1‑147 (MQ=14)                                                                                                                                                  cGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                  >  3:135343/1‑118 (MQ=255)                                                                                                                                                  cGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                  <  4:136665/118‑1 (MQ=255)                                                                                                                                                  cGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGAGCCAGCAATATCAACGAGGTTTTTGCTGCGAGGGGCGGGACTTGAAcga  >  1:83564/1‑148 (MQ=21)                                                                                                                                                        |                                                                                                                                                 AAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCC  >  NC_000913/4606047‑4606340

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