Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,192	0	G	C	24.3%	124.4 / 10.1	76	intergenic (‑27/+8)	leuV/leuP	tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (13/40);  new base C (1/16);  total (15/61)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACC  >  NC_000913/4606042‑4606339                                                                                                                                                       |                                                                                                                                                      aaaTCAAGGAAGAAACAACAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCg                                                                                                                                                       <  7:411894/151‑1 (MQ=255)  aagctcggAAGAAACATCAAAGGAAGTAAAGATAATTGGGGGGAGGGGGGGCCCTGGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                                                                                                                                                      <  3:243461/145‑1 (MQ=21)  aaTCAAGGAAGAAACTAGAACGGAAGTACAGATAATTGGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGGAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                                                                                                                                                      <  1:1476276/151‑1 (MQ=255)   aTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGcc                                                                                                                                                     <  6:427044/151‑1 (MQ=255)     cAAGGAAGAAACAAGAAACGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCt                                                                                                                                                    <  4:884381/150‑1 (MQ=255)     aaaGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc                                                                                                                                                   <  4:450339/150‑1 (MQ=255)     aaaGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGGGCGAGGGGGGGGACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc                                                                                                                                                   <  7:935853/150‑1 (MQ=255)      aaGGAAGAACCAAGAACGGAAGTAAAGATAATTGGTGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                  <  7:1383613/151‑1 (MQ=255)      aaGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                  <  8:1015677/151‑1 (MQ=255)      aaGGAAGAAAAAAGAAACGAAGTAAACATAATTGGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                  <  8:143163/151‑1 (MQ=255)      aaGGAACAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                  <  8:3977/151‑1 (MQ=255)      aaGCAATAAACAACAACGTAAGAAAATAAAATTTGGTCTCGGGGGGGTACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                  <  4:822943/151‑1 (MQ=17)      aaGCAAGAAACAACAAAGGAATTAAAGATAATTGGGGCGCGGGGGGGCACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                                                                                                  <  3:1449733/151‑1 (MQ=21)       aGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGt                                                                                                                                                 <  8:1166516/151‑1 (MQ=255)        ggAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGGAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGtt                                                                                                                                                <  4:1460888/151‑1 (MQ=255)         gaacaaACAACAAAGGAACTAAAGATAATTGGTGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                                                                                                                                               <  6:1236828/151‑1 (MQ=255)          aagaaCCAAGAACGCAAGTAAATATAATTGGGGCGAGGGGGGGGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTg                                                                                                                                              <  2:198061/151‑1 (MQ=21)             aaaCAAGAAAGGAAGTAAACATAATTGCTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                                                                                                                                           <  1:715724/151‑1 (MQ=255)              aaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGc                                                                                                                                          <  8:142100/151‑1 (MQ=255)              aaCAAGAAAGCAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGc                                                                                                                                          <  8:1267408/151‑1 (MQ=255)                cAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGGAACGACACAAACACCTGAATCGAGACCCTATAGCGAATCCCGCAACTGCGCACCGGCATGATAAAATATTTGATCTCGCCGCGTGTGGTGCgg                                                                                                                                        >  6:1269568/1‑150 (MQ=9)                    aaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg                                                                                                                                    <  8:804886/151‑1 (MQ=255)                     aagcaacaaaaTATAAGTGGTCCGCGGGGGGGCACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggg                                                                                                                                   <  5:1328995/143‑1 (MQ=21)                     aaGGAAGTAAAGATACTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggg                                                                                                                                   <  7:581176/151‑1 (MQ=255)                     aaGCAAGTAAATATAATGGGTGCG‑CGGGGGGGCCTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                                                                                                                  <  4:900816/151‑1 (MQ=21)                      aGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                                                                                                                  <  8:1101896/151‑1 (MQ=255)                        gAAGTAAAGATAATTGCGGCGAGGGGGGGGCCTTGAACCCCCACGTCCGTAACGACACTAACACCTGAAGCTACCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggg                                                                                                                                <  7:92688/151‑1 (MQ=255)                         aattAAAGATAATTGGTTCTCGGGGGTGGACTTGAACCCCCACGTCCTTAATGACACTAACACCTGAAGCTAGCGCTTCTACCAATTCCGCCACCTTCTCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggg                                                                                                                                <  4:267101/147‑1 (MQ=25)                         aataaaaaaaaaTTGGGGCCAGGGGGGGGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGa                                                                                                                               <  2:132193/141‑1 (MQ=21)                            tAAAGATAATTGGGGCCAGGGGGGGGACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACt                                                                                                                             <  5:673756/150‑1 (MQ=255)                             acagacacttGGGGCGAGGGGGGGGCCTTCCACCCCCCCGGCCGGAACGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCCCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTg                                                                                                                           <  2:450770/143‑1 (MQ=21)                                      gggTCCGCGGCGGCGGACTTGTACCACCTCTTCCTTACGCACACTCACAGCTGAACCTACCTCCTTAACCTATTCCGCCACCCTAGGCGCTTCCGCTTTCTTTATTTGATTTATCCTCTTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCAc                                                                                                                 <  7:1275886/150‑1 (MQ=255)                                       ggggCGGGGGGGGGCACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTCGCGCGTCTACCAATTCCGCCACCTTCCCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCAc                                                                                                                 <  6:1223157/147‑1 (MQ=25)                                       gcggcccgcggggggCCTACAACCCCGACGTCCCTAACGACACTACCACCTGAAGCTCCCCCGTCTACCCAATCCGCCACCATCCCACACGCATCTTACATTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCAc                                                                                                                 <  1:71133/141‑1 (MQ=255)                                        gggcctgggcggggCCTAGCACCCCCACGGGCCTTACGAGACTAAATCCTGCAGCTCCCGCGTCAACCCATTCCGGCACCTTCCCACAAGCATCTAACATTTTTTGAACTCCCGTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGt                                                                                                               <  1:560228/145‑1 (MQ=255)                                        gggCGAGGGGGTGGCCTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTACCCCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGt                                                                                                               <  5:863355/149‑1 (MQ=25)                                                   gggccttccaCCCCCCCGGCCTGAACGACACTCCCCCCTGACGCTCGCCCGTCAACCCAATCCGCCACCTACGCACAGTCAACTAACATTTTTGGAAATCCCCTCCTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAaca                                                                                                    <  8:1281065/142‑1 (MQ=12)                                                   gggACTTGAACCCCCACGTCCGTAAGGACACTAAACCCTGAAGCTAGCTCGTCTAACAATTACGACACATTAGAAGAGTCATCTTACTTTTTTTGATATCTAATCTTTTGGTGCGAGGGGGGGGACTTG‑ACCCCCCACGTCCGTAAGAac                                                                                                     >  7:324488/1‑150 (MQ=21)                                                   gctatatcaccccccAGTTCCTTAAGCACATTAAAACTTTAAGTTATCTCTTATACAAATTCCGCCACTTTCCCACATTCATCTTACTTTTTTTTATATCCCCTTTTTTGTTGCGTGGGGGGGTACTTG‑AACCCCCACGTCCGTAAGAaca                                                                                                    <  4:565651/141‑1 (MQ=255)                                                    gcccgggaccccccACCGCCGGAACGACACTACCACCGGAAGCTAGCGCGTCTACCAAATCCGCCACCTTCCCCCCCTCAGCTTACTTTATTTGATATCGCCGGTTTTGGCGCGAGTGGGGGGACTTG‑AACCCCCACGTCCGTAAGAacac                                                                                                   <  6:323340/142‑1 (MQ=16)                                                      cctacaaCCCCCACGGCCGTAATGACACTAACACCTGAAGCTAGCCCGGCTACCAATTCCGCCACCTTCCCACAGGCATCTAACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTa                                                                                                 <  8:871137/146‑1 (MQ=37)                                                      ccgggAACCCCCACGGCCGCTAGTACACTCACACATGAAGCTCCCTATTATACCAATTCCGCCACCAAACCACACGCACCAAACTTTAATTGAACTCCCCTCCTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTa                                                                                                 <  8:88627/147‑1 (MQ=9)                                                         acacccccccCGGCCTTAAGGACACTCAAACATGAAGCTCCCCCTTCTACCAATTCCGCCACCTACCAACAGTCATCAAACTTTTTTTTATTTCGCCTGGTTTGGGGCGTGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAaca                                                                                              <  5:60122/147‑1 (MQ=16)                                                         acaaCCGCCACTGCCGAATGTACCCTAACACATGAAGCTCCCCCTTAAACCAATTCCGCCACCTACACACACGCATCAAACTTTATTTGAACTACCCACCTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAac                                                                                               <  8:1279082/148‑1 (MQ=9)                                                             cccccACGTCCTTATGGACACTAACACATGAAGATACCCCTTATACCAATTCCGCCACCATCCAACAGTCATCTAACTTTTTTTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTg                                                                                         <  2:210945/151‑1 (MQ=21)                                                               cccACGGCCTTAAGGACAGTCACACCTGACGCTCGCGCTTCAACCAATTCCGCCCCCTACAAACAGTAATCTTACTTTTTTTGATATCCCCTCTTTTGGGGCGCGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGaa                                                                                       <  8:725708/151‑1 (MQ=21)                                                                 cACGTGCGTTAGTAGACTAACTCCTGATGCTAGAGCTTATACCAATTCCGGGACCTTCCAACAGTCATGTAAATTTTTTTGTTATCTCTTCTTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGc                                                                                     <  8:899098/151‑1 (MQ=12)                                                                    ggcctttagtacacTAACACCTGATGATCACGCTTAAACCTATTCCGCCACCAACACACAGGCATCTAACATTATTTGAAATCCCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAg                                                                                  <  6:1006813/141‑1 (MQ=12)                                                                           atgACAGTAACACCTGCTGATACCAATTCTAGCAAATCCGGCACCTTACCACACGCATCTTACTTTATTTGATATCCCCTTTTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCt                                                                           <  5:96814/149‑1 (MQ=12)                                                                                           aaGCTAGCGCGTCTACAAATTACGCCAACTTCGAACAGTCATCTTAATTTTTTAGATATCGACTCGTTTGGTGCGAGGGCATGAACTAA‑AACCACCACGTCCGTAAGAACACTAACA‑CCTGCAGCTAGCGCGTCTAAGAGGTCCGCCAccg                                                           >  7:495498/1‑150 (MQ=255)                                                                                             gCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCATCGTTTGGTGCGAGGGCGGGGACTTG‑AACCCCCACTTCCGTAAGAACACAAACA‑CCTCAAGCTAGACCGTATACCAAGTCCGACCCCtac                                                         >  4:382835/1‑149 (MQ=25)                                                                                                   acttCTACCAATTCCGGCAGCTAACAACAGTCAACTAACTATTTTTGATATCCACTATATTGGGGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                                                   <  8:1454903/148‑1 (MQ=14)                                                                                                     gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGGAGCTAGCGCGTCTACCAATTACACCCCAATCGCCCacgg                                                 >  2:227745/1‑148 (MQ=25)                                                                                                     gTCTACCAATTCCGCCACCTTCGCACACTCAGATTACTTTTTTTGATATCGAAGCGTTTGGTGCGAGGGGGGGGACTTT‑CACCCCAACGTCCGTAATAAGACTAAAA‑CCTCAAGCTAGGGCCTCTACCAATTCCGCCCCCATCGCCGAGtc                                                 >  2:1077661/1‑150 (MQ=18)                                                                                                       cTACCAATTCCGCCACCTTCGCAGAGTCAGATTACTTTAATAGATATCCGAGCGGGTGGTGGGAGTGGGGGGCGTTGCAGCGCCCACG‑CCCTAAGCACTC‑AACACCCTTAAGCTAGTGCTACTACCAATTCCGCCACCTTCGccaacaaac                                               >  2:362144/1‑144 (MQ=255)                                                                                                             aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGGCTTTTTTTTTGATAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTAGCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                         <  4:468688/151‑1 (MQ=21)                                                                                                             aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGGCGCGGTTGGTGCGAGGGGGGGGCCTTG‑AACCCCCACGTCAGTAAGAACACTAACA‑CCTGAAGCTAGGGCGTCTAAAAATTCCGCCACCTTCGCACAgaagcaaaaa                                          >  5:1066835/1‑141 (MQ=18)                                                                                                             aTTCCGCCACCTTCGCACAGTCATCTTAATTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCACACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                         >  2:468117/1‑151 (MQ=255)                                                                                                                           actcaGGCATCAAACATTATTTGAAGTAGCCTGTTTTGGGGGGCGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                           <  2:1383137/148‑1 (MQ=17)                                                                                                                            aacaGTCATCTAACTTTTTTGTATATCCCCTCTTTGGGTGCGCGGGGGGGGACTTG‑ACCCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                          <  2:1074867/150‑1 (MQ=21)                                                                                                                             acatgcATCTAACTTTATTGGATCTACCGTTGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                         <  7:1185657/146‑1 (MQ=18)                                                                                                                             acagggaTCTTACTTTTTTTGATATCGCGTCTGTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                         <  8:414154/145‑1 (MQ=255)                                                                                                                              catgccacatACTTTTTTTGAAACACCCGCCGTTGGTGCGCGGGGGGGCACTTC‑ACCCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                        <  2:664164/142‑1 (MQ=17)                                                                                                                                  gaacaaaCTTTATGTGATATAGCCACCTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGCACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                    <  5:616926/145‑1 (MQ=18)                                                                                                                                   aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACAAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                    >  5:810363/1‑150 (MQ=255)                                                                                                                                    acatACATTTATTGATATACCCTAGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAg                  <  8:139801/148‑1 (MQ=18)                                                                                                                                     attACTTTTTTTTATCTCTACTGTTTTGGTGCGCGGGGGGGGACTTC‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgg                 <  6:739174/150‑1 (MQ=18)                                                                                                                                       tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACAAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                   <  1:811141/147‑1 (MQ=255)                                                                                                                                       tACTTTATTTGAATTCTAGTCTTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGTACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggg                <  3:74834/150‑1 (MQ=18)                                                                                                                                       aacttgatatGATTTCTCGTCGTTTGGTGCTAGGGGGGGGAGTTG‑ACCCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg               <  6:925474/142‑1 (MQ=18)                                                                                                                                        aCTTTTATTGATGTCGCGAGGATTGGTGCGAGGGGGGGGACTTG‑AACCGCCACGTCCGTAAGCACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggg              <  5:1434478/151‑1 (MQ=18)                                                                                                                                           ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc                     <  3:1197215/141‑1 (MQ=255)                                                                                                                                           ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc                     >  8:1191281/1‑141 (MQ=255)                                                                                                                                                 gATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGAGTGCCAGGAATATAAACGTGGTTATTGGTGCGAGGGGGGGGCCTTgg     >  2:615295/1‑150 (MQ=18)                                                                                                                                                   tatCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGACAGAAATATAAACGTGGTTTTTGGAGCGAggggcggtctttgg     >  2:792616/1‑142 (MQ=18)                                                                                                                                                    atCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGAGTGGCAGGAATATAAAAGAGGTTTTTCGTGGGAGGGGGGGGCCTTgcggg  >  8:1421490/1‑147 (MQ=17)                                                                                                                                                    atCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACAACGAGCAATATCAAAGATGTGATAGCTGCGAGCGGCGGGCCATGAAcc  >  3:980389/1‑151 (MQ=14)                                                                                                                                                       |                                                                                                                                                      AAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG‑AACCCCCACGTCCGTAAGAACACTAACA‑CCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACC  >  NC_000913/4606042‑4606339

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