Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,197	0	G	T	37.5%	105.3 / 21.4	92	intergenic (‑32/+3)	leuV/leuP	tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (13/37);  new base T (0/30);  total (14/78)
Fisher's exact test for biased strand distribution p-value = 1.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAC  >  NC_000913/4606047‑4606338                                                                                                                                                        |                                                                                                                                              aaGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCg                                                                                                                                               <  8:866743/151‑1 (MQ=255)  acgaacaaACAAGAAAGCACGAAAAGATAATTGGTGCCAGGGGGGGGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGt                                                                                                                                              <  4:449022/149‑1 (MQ=21)    gaacaaAAAACAACGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                            <  5:336823/151‑1 (MQ=255)    cacgaaACAAGAACGGAATTAAAGATAATTGGTGCCAGGGGGGGGCCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                            <  7:1467925/148‑1 (MQ=255)     aagacACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTg                                                                                                                                           <  5:585157/151‑1 (MQ=255)     aagaaACAACAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTg                                                                                                                                           <  8:1396856/151‑1 (MQ=255)     aacaaaCACGAAAGGAAGTAAACATAATGGGGGCGCGGGGGGGGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTCTTGATATCGCCTCGTTTg                                                                                                                                           <  7:1202346/148‑1 (MQ=21)     aacaaaCAACAAAGTAAGAAAACATAATTGGTGCGAGGGGGGGGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTg                                                                                                                                           <  2:859119/148‑1 (MQ=21)      acaaaaaaTAACGGAAGTAAATATATTTGGTGCGCGGGGGGCTAGTCCACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTgg                                                                                                                                          <  8:1566498/149‑1 (MQ=18)       gaaACATGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                                                                                                         <  6:103941/151‑1 (MQ=255)       caaaCAAGAACCGAAGTAAAGATAATTGGGGCGCGGGGGGGCACTTCACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                                                                                                         <  4:826216/150‑1 (MQ=21)       caaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                                                                                                         <  6:98951/150‑1 (MQ=255)        aaaCAATATAGGATGTTAAGTTTTTTGGTGCGAGGGGGGTGAGTTGTACGCCGAGGTGCGTAAGGAGAGTATGTGGTGTAGGTAGCGCGTCTAGCAATTGCGCCAGCTTCGCTGAGTCA‑TGTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                        <  1:1031400/151‑1 (MQ=255)        aaaCAATAACGTAATTAAAGTTTTTTGTGGCGTGGGGGGTGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                        <  3:788669/151‑1 (MQ=18)        aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                        <  1:60799/151‑1 (MQ=255)        aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTACGCCACCTGCGCAGAGTCA‑TCATACTTTTTTTCATATCTCCTCGTTTGGTg                                                                                                                                        >  4:1160957/1‑151 (MQ=21)           ccacaaaGTAATTAAATATACTGCGTGCGTGGGGGGTCACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                                                                                     <  2:1362735/147‑1 (MQ=18)             atactggAAGTAAAAATAATTGGGGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgg                                                                                                                                   <  4:479632/146‑1 (MQ=255)              caacggAAGAAACGATAATTGGGGCGAGGGGGGCGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                                                                                                  <  4:1197708/147‑1 (MQ=255)              aaaaGGAAGTAAAGATAATTGGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                                                                                                  <  5:739255/150‑1 (MQ=255)               aacggAAGTAAAGATAATGGGGGCGCGGGGGGGCACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                 <  7:719960/148‑1 (MQ=25)               aacggAAGAAAAGATAATTGGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                 <  5:979253/148‑1 (MQ=255)               aaaGGACTTAAAGATAATTGGTGGGAGGGGGTGGACTTGAACCCCCACGTCAGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGCTGCGAgggg                                                                                                                                 <  8:648975/151‑1 (MQ=21)               aaaGGAATTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                                                                                                  <  1:980275/150‑1 (MQ=255)               aaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                 <  8:691185/151‑1 (MQ=255)               aaaGGAAGTAAAGATAATTGGGGCGAGCGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGATGCTAGGGCGGCTACCAATTCCCCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                 <  5:1007578/151‑1 (MQ=25)                aaGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggg                                                                                                                                <  8:959394/151‑1 (MQ=255)                       tAAACATAATTGGTGCGCGGGGGGGCCCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGGCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACt                                                                                                                          <  7:1260663/150‑1 (MQ=255)                        aatcttatttGGTGCGTGGGGGGTGACTTTCAGCCCCACGTCCGTAAGGTCACTAACACATGAAGCTTTCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTg                                                                                                                        <  4:162178/143‑1 (MQ=21)                             aatttgcggccaGGGGGGGTCCTTCCACCCCCACGTCCCTAAGGACACTAACACCTGCAGCTAGCGCGTCTACCAATTCCGCCACCTACGCACACTCA‑TCTTACTATTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAAccc                                                                                                                   <  8:1568932/140‑1 (MQ=21)                                tttttgctacgggggccACTTCAACCACCACTTAAGTAATCACACTTACCCCTCAAGCTCGCGCTTCTAACAATTCCGGCATCAAAAAACATTCA‑TCTTACTTTATTTGATATACCCTACTTTGGTGCG‑CGGGGGGGACTTGAACCCCCa                                                                                                                <  4:959364/140‑1 (MQ=255)                                     gcctggcgtctgccTTCAACCGCCACGTCCCTATCTTCACTACCACCTGAAGCTACCGAGGCAACAAATTCCGCGACCTTCGCACAGTGA‑TCTAACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                                                                                           <  2:146886/138‑1 (MQ=19)                                        cggggggggCCTGTACCCCCCACGTCCGATTTTTCACTACAACATGAAGCTCCCGCTGATAGGAATTCCGCCACCTTCCCACAGGCA‑ACTTACTTTTTTTGATTTCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                        <  8:1346728/150‑1 (MQ=12)                                         ggggggggACTTCAACCCCCACGTCCTTAAGTTCACTAACACCTGAAGCTAGCCCGTCAACCAATTCCGCCACCTTCCCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTaa                                                                                                       <  5:1547236/151‑1 (MQ=25)                                         gggcggttccTTCAACCCCCACGTCCGTAAGTACACTCACACCTGAAGCTCCTCCGGCAACCAATTCCGCCACCTACCCACAGTCA‑CGTTACTTTTTTTGATCTCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTaa                                                                                                       <  8:624795/142‑1 (MQ=21)                                          gggttggACTTTAACCCCCACGTCCTTAAGTAAACTAAAACCTGAAGCTAGAGCGGAAACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                      <  5:201991/146‑1 (MQ=25)                                          ggggtgtccTTCAACCCCCACGGCCTTTTTTACAGTAAATCCTGATGCTAGCGCGGATACCAATTCAGCCACCTTCCCACAGGCA‑CCTAACATTTTTTGATCTCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                      <  8:342304/143‑1 (MQ=16)                                          gcgcgggcctcccaCCCCCACGGCCCCATCGCCACTCCCACCTGATGCTCCCCAGTCAACCCAATCCGCCACCAACCCACACGCA‑CCTAACTATTTTTGAACTCCCCACCTTTGGTGCCAGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                      <  1:284816/138‑1 (MQ=2)                                              tttccTTCAACCCGCACGTCCTTATTTACACTAACTCATGAAGATACCCATTAAACCAATTCCGCCACCAAACAACATTCA‑TCTAACATTTATTGATCTCCCCTCCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                                                                                                   <  8:1084292/146‑1 (MQ=12)                                              gttactacaaCCCCCACGTGCGTATATACACTAAAACATGAAGCTTGCACGGCTACCAATTCCGGCACCTTCCCACATGCA‑TGTTACTTTTTTTGATATCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                                                                                                  <  1:304299/143‑1 (MQ=21)                                               tgACTTCAACCCCCACGTCCGTTAGGACACTAACACCTGATGCTCCCGCTTCTACCAATTCCGCCAGCTTACCACAGGCA‑TCTAACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                                                                                                  <  5:741201/149‑1 (MQ=25)                                                  atacaaCCCCCACGTCCTTAATTACACTAAAACCTGAAGATAATCATTAAATCAAATCCGCCACCTTCTAACAGTCA‑TCTAACTATTTTTGATATATCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                              <  2:1323116/147‑1 (MQ=16)                                                    agAACCCCCCCGGCCCTAACTACACTAACACCTTACTCTCTCTCTTCAAACAATTCAGCCACCTTATAACATTCA‑CCAAACAATAATTTAACTCCCCTCCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                                                                                            <  1:286923/150‑1 (MQ=2)                                                             aCGGCCGTAAGGACACTAACACCTGAAGCTAGCCCGTATACCAATTCCGCCACCTTCCAACAGTCA‑TCTAAATTTTTTTGATATCCCCTCTTTTGGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                                                                   <  7:1137471/151‑1 (MQ=25)                                                                    tatttacacTAACACATGAAGATAAACCGTATATCAAATCCGGCACCAAAAAACAGTCA‑TCTAACATTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcgc                                                                              <  5:775458/144‑1 (MQ=12)                                                                      aggacagtccaacatGATGCTACCGAGTATACAAATTCCGCCACCTTACAACAGGCA‑TCATACTTTTTTTGATTTCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                                                           <  6:311287/136‑1 (MQ=17)                                                                      aGGACACTAAAACATGCAGCTCCCGCGTATACCAATTCCGCCACCTTCCCACAGTCA‑TCTTACTTTTTTTGATATCCCGTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                                                          <  7:1074905/151‑1 (MQ=21)                                                                       gtacagTAACACATGAAGATATCGCGTATACCAATTCCGCCACCTTAGCACAGGCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                                                          <  5:720019/148‑1 (MQ=25)                                                                       ggacattcccacCTGATCCTAGCGCTTATTCGAATTCCGCCACGAACCACCAGGCA‑TCAAACTTTTTTTGATGTCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                                                          <  6:189508/141‑1 (MQ=9)                                                                          cacTAACCCCTCACGCTAGCGGGTCTACCATTTCCGCCCCGTGCGACCCGTAA‑TCTAACTTTTTTTGAAATCCCGTGGTTTGGTGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                      <  7:1164728/151‑1 (MQ=17)                                                                          cacTAACACCTGAAGATCGATCTTCTACGAATTCCGGCACCATCCAACATGCA‑ACTAACTTTTTTTGAAATAGCCTCCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                      <  8:1343553/151‑1 (MQ=17)                                                                              aatacCTGAAGCTCGCGATTATACAAATTCCGCCACGAAAAAACATTCA‑ATAAACTTTTTTTTATTTCCCCTCTTTTGGTGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                                  <  7:823413/148‑1 (MQ=9)                                                                              aaCACCTGAAGCTCCCGCGGATACCAATTCCGCCACGTAACCACAGTCA‑TCTAACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                                  <  5:819876/151‑1 (MQ=21)                                                                               aaacCTGCAGCTACCGCGTAATACAATTCCGCCACCTTCCCCAACGCA‑CAAAAACTTTTTTGATATCCCCCTCTTTGGTGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                 <  8:470755/149‑1 (MQ=2)                                                                                aacattatgcttacgctTCTACAAATTCCGCCACTTTCCAACA‑TCATTGTAACTTTTTTTGATTTCCCCACTTTTGGTGGGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                                <  7:96842/138‑1 (MQ=9)                                                                                         cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCGGAAGCTAGCGCGTCAGCCAAAACCGCccccaaac                                                       >  5:1539911/1‑147 (MQ=25)                                                                                         atgcacattatACAAAATCCGGCAGCTACACACAGGCA‑TCTAACTTTTTTTGAATTCCCGTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                       <  8:376360/141‑1 (MQ=14)                                                                                             ctcgtaaacaaaTTCCGCCACCTTCCCACAGTCA‑TCTAACTTTTTTTTATATCCCCTGTTTTGGGGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                        <  8:1572606/136‑1 (MQ=17)                                                                                             cgcggatACCAATGCCGCCACCTACCCACAGTCA‑TCTAACTTTTTTTTATATCCCGTCTTTTGGGGCGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                                                   <  7:1138158/145‑1 (MQ=17)                                                                                                gtatagctaTTCCGCCTCCTAAGCACAGGCA‑ACAAACTTTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCACCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                                                  <  7:573919/141‑1 (MQ=17)                                                                                                 tatACCAATGCCGCCACCAACCAACAGCCA‑TGAAAATTTTTTTGATTTCCCCTGTTTGGGGGGGGGGGGGGGGACTTCACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                <  7:851514/148‑1 (MQ=9)                                                                                                 tCTACCTATTCGTCCATATTTTCACATTCA‑TCTAACTTTTTTTTATTTATAGTTTTTTGGTGGGTGGGGGGGGATTTTAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                <  4:604061/150‑1 (MQ=2)                                                                                                   tACCAATTCCGCCACCTTCGCACAGTCA‑TATTACTTATTTTGATTTCGGAACCTTTGGTGCGAGGGGGGGGCCTTGAACCCGCACGTCAGTAAGTAAAAATCCACATACCGGTAGAGCCTCTACCAATTCCGACACCTACGctaccttcc                                              >  5:58255/1‑142 (MQ=255)                                                                                                    tccAATGCCGCCACGTTACAACATGCA‑TCTTACATTTTTGTATTTAGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                              <  6:541522/148‑1 (MQ=17)                                                                                                    acaaaTTCCGCCACCTTACCACAGGCA‑TCTAACTTTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                            <  8:702268/148‑1 (MQ=18)                                                                                                     cctaTTCCGCCACCTAATAACATGCA‑TCTTACTTTTTTTGATTTCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc                                           <  1:1564933/148‑1 (MQ=18)                                                                                                     ccAATTCCGCCACCTTCGCACAGTCA‑GCTTACTTTTAATGATATAGCTGGGTTTCGGGCCCGGGGGGGCACGAGAACCACCACCGCCGTAAGAACACTAACACCTGAAGCTAGCCCGTCTACCCATTCCGCCCCCTTCGCGCAgggggcgg                                           >  7:65095/1‑144 (MQ=255)                                                                                                      cAATTCCGCCACCTTCGCACAGTCA‑GCTTACTTTTATTGATATCGCCTCGTTTGGGGCGAGGCGGGGGCCTTGAACCACCACGTCCGTAAGAACACTACCACCTGCCGCTAGCGCCTCTACCGAATCCGCCCCCTTCGcaaagccagaca                                           >  8:1274103/1‑139 (MQ=255)                                                                                                      aaaTTCCGCCACCAAACAACATGCA‑ACTAACATTTATTGATTTACCCTCCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCa                                          <  7:1256961/150‑1 (MQ=14)                                                                                                        aTTCCGCCACCAACCCACATTCA‑ACAAACTTTTATTGAATTCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                        <  7:501701/151‑1 (MQ=17)                                                                                                              gcacaTTCCCCCAGTCA‑TGTAACTTTTTTTGATATACCTTCTTTGGGGCCTCGGGGGGGGACTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                  <  4:711260/150‑1 (MQ=17)                                                                                                              ccACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACGGCGACCAATATaaa                                   >  6:641201/1‑150 (MQ=255)                                                                                                                   atatctcaGTGA‑TCTTACTTTTTTTGATTTCCCCTTATTTGGTGGGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTgg                              <  4:562620/144‑1 (MQ=18)                                                                                                                   aaaacacaCTCA‑TCTAACAATTTTTGAACTCCCCTCCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                             <  5:28957/147‑1 (MQ=18)                                                                                                                    acccacaGGCA‑TCTTACTTTTTTGGATATCGCGTTGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                            <  8:514223/148‑1 (MQ=255)                                                                                                                     aacacaGGCA‑TCAAACATTATTGGAAATCCCGTCTATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttt                           <  5:1201706/149‑1 (MQ=17)                                                                                                                      gCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGGGAGACAAATATAAAAGTGTtttt                          >  6:1140102/1‑151 (MQ=18)                                                                                                                        cacatgcatcaaaCTTTATTTGAAATCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                         <  7:1045852/139‑1 (MQ=18)                                                                                                                        atcattcttCTTACTTTTTTTGATTTATCGTTTTTTGGTGGGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                          <  4:1152533/142‑1 (MQ=18)                                                                                                                         caGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGTGAGACCAATATCAACGAGGTTTTTgc                       >  5:215728/1‑150 (MQ=255)                                                                                                                         acatgcaccaaaCTTTTTTGGATGTCACCTGTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                        <  7:300118/140‑1 (MQ=18)                                                                                                                          aattgatgtaaCTTTTTTTGATTTATTCTTTTTTTGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTTAGAAGACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                       <  3:1362275/141‑1 (MQ=18)                                                                                                                          aattcaacaaacattttttGATATATACTTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                       <  4:443809/138‑1 (MQ=17)                                                                                                                            tCA‑TGTTACTATTATTGATTTCTCTTAATTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                     <  5:1221528/150‑1 (MQ=18)                                                                                                                             cA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGACGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGGGGCAGCAATATCAACggcgttgtaggagcg                   >  5:96262/1‑137 (MQ=18)                                                                                                                              gaacaaaCTTTATTGGATTTACCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc                    <  3:31431/144‑1 (MQ=18)                                                                                                                               aacaaacattttttGAATTACCCTTGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                  <  8:1377254/143‑1 (MQ=18)                                                                                                                                 cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCAGAGGGTGAGCTATTTAAAAGACGTTATTGGTGCGAgg                >  8:218034/1‑151 (MQ=17)                                                                                                                                   aagttttttgGATTTACCGTCATTTGGGGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTATCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                  <  8:378660/144‑1 (MQ=18)                                                                                                                                      tttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGACACACGAATAACAACGAGGTGGTAGCTGCGAggcgggg           >  4:746553/1‑151 (MQ=14)                                                                                                                                        tttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCCTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGC‑‑CCGTCGCCCAGTGCGTGCAATATGAACGTGGTTTGTGGGCCGGGGGGGGCGagg       >  3:1331477/1‑149 (MQ=18)                                                                                                                                            tGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAAGTGGTTTTTGGTGCGAGGGGGGGGACTTg     >  4:93763/1‑151 (MQ=255)                                                                                                                                               tatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGAGAGCGACCAATATCAAAGTGGTTTTTGGTGCGAGGCGGGGGCCTTGCAc  >  7:32546/1‑151 (MQ=18)                                                                                                                                                        |                                                                                                                                              AAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCA‑TCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAC  >  NC_000913/4606047‑4606338

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: ‑