Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,192	0	G	C	31.1%	127.4 / 23.7	95	intergenic (‑27/+8)	leuV/leuP	tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (18/44);  new base C (0/28);  total (18/77)
Fisher's exact test for biased strand distribution p-value = 5.58e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCC  >  NC_000913/4606043‑4606340                                                                                                                                                       |                                                                                                                                                      attcaactaataatcaacaaAGGAATTAATGATCTTTGTTGCCGGGTGTGGT‑ACTTGAACCCCGACGTCCGTAAGGACACTAACAACTGAAGCTAGCGCGTCTAGCTATTCCGCCACCTTCGCTCAGTCATCTTACTTTTTTTGATATCGc                                                                                                                                                      <  3:100350/137‑1 (MQ=255) aaTCAAGTACGAAACAAGAACGGAAGTAAACATATTTGGGGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                                                                                                                                                      <  8:846502/151‑1 (MQ=18) aaTCAAGGAACAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                                                                                                                                                      <  4:587625/151‑1 (MQ=255) aaTCAAGCAACAAACAAGAAAGGAAGGAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                                                                                                                                                      <  3:100343/151‑1 (MQ=255)  ataaatgaagaaAAAAGAATGCAAGTAAACATAATTGGTGCGAGGTGGGTT‑AGTTTAACCCACACGTCCGTAAGGACACAAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGcc                                                                                                                                                     <  1:1126495/145‑1 (MQ=18)  aTCACGGAACAAACAAGAACGGAAGTAAACATAATGGGGGCGCGGGGGGGC‑ACTGGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGcc                                                                                                                                                     <  5:1306050/151‑1 (MQ=21)  aTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGcc                                                                                                                                                     <  3:873617/151‑1 (MQ=255)    caagcaagAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCt                                                                                                                                                    <  7:137694/150‑1 (MQ=255)    cAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCt                                                                                                                                                    <  2:48364/150‑1 (MQ=255)    cAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc                                                                                                                                                   <  4:1553010/151‑1 (MQ=255)    cAAGGAAGAAACAAGAAAGGAAGTAAACATAATTGGTGCGAGGGGGGGG‑ACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc                                                                                                                                                   <  2:5634/151‑1 (MQ=255)    cAAGCAACAAACAAGAAAGGAAGTAAAGATATTGGGGGCGCGGGGGGGG‑CCTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc                                                                                                                                                   <  5:1257021/151‑1 (MQ=21)    aaatgacgaaACAACAACGGAAGTAACTATAATGGGGGGCGGGGGGGGC‑ACTTCAACCCCCACGGCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc                                                                                                                                                   <  3:794593/144‑1 (MQ=18)     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                                  <  4:136665/118‑1 (MQ=255)                                                                                                                                                      cGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                  >  3:135343/1‑118 (MQ=255)                                                                                                                                                      cGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGAGCCAGCAATATCAACGAGGTTTTTGCTGCGAGGGGCGGGACTTGAAcga  >  1:83564/1‑148 (MQ=21)                                                                                                                                                       |                                                                                                                                                      AATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCC  >  NC_000913/4606043‑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: ‑