Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	718,690	0	T	G	42.2%	104.3 / 35.8	72	T591P (ACC→CCC)	speF	ornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base T (12/29);  new base G (30/0);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 1.20e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGC  >  NC_000913/718546‑718835                                                                                                                                                  |                                                                                                                                                  cACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAAcc                                                                                                                                             <  2:752836/151‑1 (MQ=255)    cgcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTGTAACCGGc                                                                                                                                          <  2:772811/151‑1 (MQ=255)         aaaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgc                                                                                                                                     <  1:853547/151‑1 (MQ=255)           aTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGTGATACAAATC‑ATGCATTTTCTTGCACAAATGACGCAGGGGATAAACGGCg                                                                                                                                         >  1:440879/1‑145 (MQ=255)              ggCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGctc                                                                                                                                 <  2:746117/150‑1 (MQ=255)                  tAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTTCTTACACAACTGACGCAGGGGATAAACAGCGGAGCGctc                                                                                                                                 >  1:88934/1‑146 (MQ=255)                  tAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAA‑ATTCATTTCCTTACACAACTGACGCAGGGGGTAACCGGCGGAGGGCCCTTCAt                                                                                                                            >  2:325520/1‑151 (MQ=255)                   aGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATg                                                                                                                           <  2:580341/151‑1 (MQ=255)                   aGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAC‑ATGCATTTTTTTCCACAAACGAAGCAGGGGGTAAACGGCGGAgcgcccctcctt                                                                                                                           >  2:970753/1‑143 (MQ=255)                      tGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATT‑ATGCATTTTCTTCCACAACTGACGCAGGGGATAAACCGCGTAGCGCTCtttcttct                                                                                                                         >  1:787278/1‑144 (MQ=255)                      tGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑AAGCATTTTCTGTCACAAATGACGCAGGGGATAAACAGGGTAgcg                                                                                                                                    >  1:587818/1‑139 (MQ=255)                             ttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAC‑ATTCATTTTCTTACACAAATGACGCAGGGGATAAACCGCGGAGCGCCCTTTCTGCTGgtttt                                                                                                                   >  1:429454/1‑147 (MQ=255)                                   aTTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGAAAAAAAAC‑ATGCATTTTTTGTAACAACAGACGCAGGGGATAAACGGCGTAGAGCTTCTCATCCTGTTTGTTGa                                                                                                                >  1:942786/1‑146 (MQ=255)                                    ttCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAATT‑ATTCATTTTCTTACACCAACGAACCCGGGGATAAACCGCGGAGCGCTCTTCATTCTGtttttt                                                                                                                  >  1:143753/1‑140 (MQ=255)                                      cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGt                                                                                                        <  1:205354/151‑1 (MQ=255)                                      cATGCTGACGCGTGGGAAGTGCTCCCTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAACTC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAgg                                                                                                         <  1:362922/150‑1 (MQ=255)                                        tGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGCAAAAAAAC‑AAGCATTTCCTGGCACAACAGAAGCCGGGGGTAAACGGGGGTGGGCTCTTTTTGCTGTTTGTtggtggga                                                                                                           >  2:492350/1‑144 (MQ=255)                                         gCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGc                                                                                                     <  2:787425/151‑1 (MQ=255)                                            gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAA‑ATTCATTTTCTTACACAAACGACGCAGGGGATAAACGGCGTGGCGCCCTTCCTGCTGTTTTTGGATGGGAGGGAGCa                                                                                                    >  2:593514/1‑149 (MQ=255)                                              cgcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTc                                                                                                <  1:243623/151‑1 (MQ=255)                                                cgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTg                                                                                              <  2:373832/151‑1 (MQ=255)                                                 gTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGGTGCAAAAA‑AAGAATTTTCCTACACAAATTACGCAGGGGGTAAACCGCGGAGCGCCCTTCCTGCTGTTTTTTGATgggaaggaagactttt                                                                                               >  1:178391/1‑135 (MQ=255)                                                   gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATTCAAAAA‑ATGCATTTCCTGACACAAACGACGCAGGGGAAAAACAGCGGAGGGCTCTTCTTGCTGTTTGTATGTGGGAGGGAGCAATTCTGCt                                                                                            >  2:690797/1‑150 (MQ=255)                                                   gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTa                                                                                           <  2:291788/151‑1 (MQ=255)                                                    ggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGCTAAAAAAA‑ATGTATTTTCTGACAAAAATCACGCAGGGGATAAACCGCGGGGGGCCCTTTTTTCTGTTTTTAAATTGAAGGGAGCACTTTTGCt                                                                                            >  2:507311/1‑149 (MQ=255)                                                        gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgc                                                                                      <  2:352768/151‑1 (MQ=255)                                                        gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCACACAAAAC‑ATGCATTTTCTGCCCCAACTGACGCAGGGGATAAACGGCGGCTGGCTCTTCATGCTGTTTTTTTGTGGGAGGTAGTGCTTCTTCTAAtag                                                                                       >  2:574806/1‑148 (MQ=255)                                                          gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgcgc                                                                                    <  1:737408/151‑1 (MQ=255)                                                          gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgcg                                                                                     <  2:937255/150‑1 (MQ=255)                                                             ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGcgcg                                                                                   <  1:625504/149‑1 (MQ=255)                                                             ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGcgcg                                                                                   <  1:965082/149‑1 (MQ=255)                                                              cTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAC‑ATTCATTTTCTGACACAAATGACGCAGGGGATAACCGGCGGAGCGCTCTTCCTGCTGTTTGTAGATGGGAGGGAGCACTTCTGCTAAttgggc                                                                                    >  1:524327/1‑142 (MQ=255)                                                               ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAAA‑ATGCGTTTTTTGTCACAAATTACCCAGGGGATAAACGGGGGTGCGCCCTTCTTTCTGTTTTTTGATGGGAAGGAGGACCTCTGCTAAttg                                                                                       >  2:673729/1‑141 (MQ=255)                                                               ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCAACAAAAC‑CTGCTTTTTCTTACACAACATACGCCGGGGCAAAACGGCGTAGCGCCCTTCTTCCTGTTTGTTGATGGGGGGAAGCACCTTTTCCAATgg                                                                                       >  2:870381/1‑143 (MQ=255)                                                               ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAC‑ATGCATTTTCTTCCACCACTGACGCCGGGGGTAAACAGCGTAGCGCTCTTTATGCTGTTTTTAGATGGGAGGGAGAAACTCTTCTCATGgggccgcgg                                                                               >  1:514268/1‑146 (MQ=255)                                                                   ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCAAACAAATC‑ATTGCTTTTCTGGCACAAATGACGCAGGGGAAAAACCGCGTAGCGCTCTTCCTCCTGTTTTTTGAAGGAAGGGAGAACCTTTGCTAATTGCGCGTCgg                                                                               >  1:699485/1‑147 (MQ=255)                                                                     aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGCACAAAAA‑ATGCCTTTTTTGTCACCACTGACGCAGGGGATAAAACGGGTAGCGCCCTTCATGCTGTTTGTAGATGGGAGGGAGAAATTCTCCTAATGGCGGGGCGGGCTc                                                                           >  2:729520/1‑149 (MQ=255)                                                                     aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCaa                                                                         <  2:877925/151‑1 (MQ=255)                                                                       cATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGAGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGGATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCaa                                                                         >  2:299445/1‑149 (MQ=255)                                                                         tctcTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAATA‑ATGTCTTTTCTGGCACAAATGACGCCGGGGGTAACCCGGGGGGCGCTCTTCTTTCTGTTTTTTGATGGAAGGGGGGAATTCTTCTTATGGCGcggcgg                                                                               >  2:612692/1‑141 (MQ=255)                                                                            ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAA‑ATGCATTTTCTTACACCACCGACGCCGGGGGTAAACCGGGTAGGGCCCTTCATGCTGTTTTTAGATGGAAGGGAGCAAATCTGCGAATGGCGCGGCGGGCCCAAGccggtt                                                                  >  2:456704/1‑146 (MQ=255)                                                                                   cAGTTGTTTCACGTTGTGGCGGGGATAAAAAAC‑ATGCCTTTTCTGGCACAACATACGCAGGGGATGACCCGCGGAGCGCCCTTTCTTCTGTTTTTTGATGGGAGGGAGGACCTCTTCTAATTGCGCGGCGGGCTCAAACAGatttt                                                                >  2:180671/1‑146 (MQ=255)                                                                                     gttgttTCACGTTGTGGCGGGCCTAAAAAACAATTCCTTTCTTGA‑ACAACTGACGCGGGGGATAAACGGCGCAGCGCTCTTCCTTCTGTTTGTTGATGGGAGGTAGCAATTCTGCTAATGGCGCGTCGGGCTCAACTAGgtttttcatgc                                                          >  2:52065/1‑144 (MQ=255)                                                                                     gttgttTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTGTGCATGCTGTTGGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGCAga                                                          <  2:269874/150‑2 (MQ=255)                                                                                       tgttTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGt                                                        <  1:257918/150‑1 (MQ=255)                                                                                         tttCACGTTGTGGCGGGCATACAAAAC‑ATGTATTTCCTGACACAACAAACGCCGGGGGTAAACCGCGTAGCGCCCTTCTTGTTGTTTTTTGATGGAAGGGAGGACATTTTCTAATGGTGCGGGGGGCTCCAACAGTTTTTTTAAGCGTa                                                       >  2:124933/1‑149 (MQ=255)                                                                                           tCACGTTGTGGCGGGCATACAAAAC‑ATGGATTTTCTGGCACAACTGAAGCCGGGGATAACCAGCGGAAGGCCCTTTCTTCTGTTTTTAGATTGGAGGGAGCACATCTTCTTATGGCGCGCCGGGCTCCAACAGgttttttgaggggaccc                                                    >  2:509980/1‑138 (MQ=255)                                                                                            cACGTTGTGGCGGGCATACAAAAC‑ATGCATTTTCTGTCACAACTCACGCAGGGGATAAACGGCGGAGCGCTCTTTATGCTGTTTTTAGATGGGAGGGAGCACTTCTGCTAATGGGGCGGCGGGCCCAAGCAGTTTTTTGAAGGGGACcagc                                                  >  2:864207/1‑151 (MQ=255)                                                                                             aCGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACcagca                                                 <  1:340187/151‑1 (MQ=255)                                                                                                 tgtgGCGTGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGc                                             >  1:18360/1‑151 (MQ=255)                                                                                                  gtgGCGGGCATACAAAAA‑ATGCATTTCCTGACACAAATGACGCAGGGGATAAACGGCGGAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGGCGGGCTCAAGCAGGTTTTTGAAGCGTACCAACAggg                                              >  2:471855/1‑149 (MQ=255)                                                                                                   tgGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGca                                            >  1:926459/1‑150 (MQ=255)                                                                                                     gCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaac                                          >  2:502810/1‑150 (MQ=255)                                                                                                           aTACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGc                                   >  1:140568/1‑151 (MQ=255)                                                                                                             aCAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg                                 >  2:555565/1‑151 (MQ=255)                                                                                                               aaaTC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa                                >  2:441853/1‑150 (MQ=255)                                                                                                                aaTC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAg                              <  1:885516/151‑1 (MQ=255)                                                                                                                   c‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGttt                           <  1:676852/151‑1 (MQ=255)                                                                                                                        cATTTCCTGACCCAACTGACGCAGGGCATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGAcc                       >  2:298496/1‑151 (MQ=255)                                                                                                                           ttCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCat                     <  1:10997/150‑1 (MQ=255)                                                                                                                            tCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCatat                   >  1:241174/1‑151 (MQ=255)                                                                                                                             ccTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCatat                   <  2:144523/150‑1 (MQ=255)                                                                                                                             ccTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCatat                   <  1:630411/150‑1 (MQ=255)                                                                                                                                gACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtctt                <  2:851200/150‑1 (MQ=255)                                                                                                                                 acacAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtcttct              >  2:212654/1‑151 (MQ=255)                                                                                                                                  cacaACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTg             <  2:517120/151‑1 (MQ=255)                                                                                                                                         gACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGt      >  2:522123/1‑151 (MQ=255)                                                                                                                                         gACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGt      <  2:627338/151‑1 (MQ=255)                                                                                                                                          aCGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGt      >  1:281924/1‑150 (MQ=255)                                                                                                                                            gCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGt      >  2:203365/1‑148 (MQ=255)                                                                                                                                            gCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTcag   <  1:883665/151‑1 (MQ=255)                                                                                                                                              aGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGCTTGGCCATATCTTCTGCCGGAGTcagc  >  2:386607/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATC‑ATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGC  >  NC_000913/718546‑718835

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