Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	718,659	0	T	A	44.6%	106.1 / 36.3	74	D601V (GAT→GTT)	speF	ornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base T (5/36);  new base A (33/0);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 4.47e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAG  >  NC_000913/718509‑718807                                                                                                                                                       |                                                                                                                                                       cGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAt                                                                                                                                                        <  1:781890/151‑1 (MQ=255)    ccTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCAt                                                                                                                                                     <  1:401309/151‑1 (MQ=255)          gCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCAttt                                                                                                                                               <  1:1083807/151‑1 (MQ=255)          gCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCAtt                                                                                                                                                <  1:591443/150‑1 (MQ=255)            aTCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTcc                                                                                                                                             <  2:593700/151‑1 (MQ=255)                   aGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGCTACAAAACCTGCGTTTTCTCaccc                                                                                                                                       >  1:58125/1‑148 (MQ=255)                   aGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTACAAAACATGCATTTTCTGac                                                                                                                                         >  2:885965/1‑148 (MQ=255)                    gACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACaa                                                                                                                                     <  1:1098355/151‑1 (MQ=255)                       gAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAACAAACATGCAttttttttcacaa                                                                                                                                     >  2:65752/1‑138 (MQ=255)                        aaCCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGa                                                                                                                                 <  1:307877/151‑1 (MQ=255)                          ccAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACg                                                                                                                               <  2:323188/151‑1 (MQ=255)                               tCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTGCACAAATTACGC‑Cggg                                                                                                                         >  1:906698/1‑151 (MQ=255)                                    ttCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGAATAAAAAACATTCATTTCCTTACACAACTGAAGA‑AGGGgata                                                                                                                     >  1:335259/1‑150 (MQ=255)                                     tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTGACACAACTGACGC‑CGGGGATaa                                                                                                                    >  2:517450/1‑150 (MQ=255)                                         cgcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAAAATGAATTTTCTTACACAACCGACGC‑AGGGGATaaa                                                                                                                   >  2:1145243/1‑146 (MQ=255)                                         cgcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCgtc                                                                                                               <  2:339244/151‑3 (MQ=255)                                         cgcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGc                                                                                                               <  1:731617/151‑1 (MQ=255)                                           cgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTTCACAACCGACGC‑AGGGGATAACCCGc                                                                                                               >  2:189914/1‑149 (MQ=255)                                              aaaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTGACACAACAGACGC‑CGGGGATAACCGGCGGAg                                                                                                           >  1:857439/1‑150 (MQ=255)                                               aaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGCATTTTCTTACACAACTGACGC‑AGGGGATAACCCGCGTAgcg                                                                                                         >  1:315931/1‑151 (MQ=255)                                                  aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGct                                                                                                       <  1:147092/150‑1 (MQ=255)                                                   ggCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAAAATGCATTTTTTTACACAACACAAGC‑AGGGGATAACCGGCGTAgcgcg                                                                                                       >  2:712794/1‑148 (MQ=255)                                                    gCGTAGTTGGCTTCTTGCGTATTCATTCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTCTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCtt                                                                                                    <  2:677982/151‑1 (MQ=255)                                                     cGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTc                                                                                                   <  2:272511/151‑1 (MQ=255)                                                      gTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCa                                                                                                  <  2:652296/151‑1 (MQ=255)                                                         gTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAAAAAACATTCATTTTCTTTCACAACCTACGC‑AGGGGGTAAACAGCGGAGGGCCCTTcct                                                                                                 >  2:1077026/1‑147 (MQ=255)                                                          ttGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTTCTTTCACAACCGACGC‑AGGGGATAAACGGCGGAGCGCCCTTCCTg                                                                                                >  2:522455/1‑149 (MQ=255)                                                            ggCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAACTGACGC‑AGGGGATAAATTGCGTAGCGCTCTTCATGCTGt                                                                                            >  1:530046/1‑151 (MQ=255)                                                              cttcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTTACACAAATGACGC‑AGGGGATAAACGGCGGAGCGCTCTTCATTCTGttt                                                                                          >  1:964303/1‑151 (MQ=255)                                                                   tGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGCATTTCCTTACACAAATGACGC‑AGGGGATAACCGGCGGAGCGCCCTTCTTGCTGttttttg                                                                                      >  1:959937/1‑146 (MQ=255)                                                                       gATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATGCATTTCCTGACCCAAAAGACGC‑AGGGGATAACCGGCGGAGCGCTCTTCATGCTGTTTTTAg                                                                                      >  1:945373/1‑146 (MQ=255)                                                                           cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAgg                                                                              <  2:1174475/150‑1 (MQ=255)                                                                            aTGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAACATTCCTTTTCTTACACAACTGCCGC‑AGGGGGTAAACCGCGCGGCGCCCTTCATGCTTTTTTTTGATGGGAgg                                                                              >  2:309226/1‑149 (MQ=255)                                                                            aTGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATACATTTTCTTTCACAAATGACGC‑AGGGGATAAACGGCGGAGCGCTCTTCATGCTGTTTTTTGATgg                                                                                  >  1:675183/1‑145 (MQ=255)                                                                                 gACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑t                                                                      <  1:913670/151‑1 (MQ=255)                                                                                    gcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACAAGCATTTTCTGACACAACTCACGC‑AGGGGGTAAACGGGGGGGGGCCCTTCATGCTGTTTTTTGATgg                                                                                  >  2:204016/1‑137 (MQ=255)                                                                                     cgTGGGAAGTGCTCCTTACAGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACAATCATTTCCTGACACCAAAGACGC‑AAGGGATAAACGGCGGAGCGCTCTTTATTCTGTTTGTAGATGGAAGGGAGCAC‑TTCt                                                                   >  1:544785/1‑150 (MQ=255)                                                                                       tGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCAAACCAAACATGGCTTTTCTTAAACAAATGACCA‑CGGGGATAAACAGCGGCGGGGCCCTCATGCTGTTTTTAGGTGGGAGGGAGGAcat                                                                      >  2:498849/1‑144 (MQ=255)                                                                                        gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCt                                                                <  2:366628/150‑1 (MQ=255)                                                                                        gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTa                                                               <  2:81552/151‑1 (MQ=255)                                                                                        gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTa                                                               <  2:335254/151‑1 (MQ=255)                                                                                           aaGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATg                                                            <  2:717670/151‑1 (MQ=255)                                                                                             gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGgc                                                          <  2:885425/151‑1 (MQ=255)                                                                                             gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGgc                                                          <  1:972186/151‑1 (MQ=255)                                                                                               gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGgcgc                                                        <  1:1115315/151‑1 (MQ=255)                                                                                                cTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATTCATTTTCTTACACAACTGACGC‑AGGGGATAAACCGCGGAGGGCCCTTCCTGCTGTTTGTTGATGGGAGGGAGCAC‑TTCTTCTTATGGcgcg                                                       >  1:356256/1‑151 (MQ=255)                                                                                                 tCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATTCAAAACATGCATTTCTTTTCTCAAACGACGC‑AGGGGGTAAACCGCGGAGGGCTCTTCCTTTTGTTTTTTGATGGGAGGGAGGAC‑TTCTTCTAATgggg                                                         >  2:874680/1‑146 (MQ=255)                                                                                                   cTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCACACAAAACATGCATTTTTTGTCACAACAGAAGC‑AGGGGATAACACGGGTAGCGCTCTTCTTGCTGTTTGTATATGGGAGGTAGGAC‑TTCttcctattgcg                                                         >  2:595405/1‑136 (MQ=255)                                                                                                    ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGAAACAAAAAAAGGATTTTCTGTCACAAATGACGCAAGGGAAAAACCGG‑GGAGCGGCCTTTATTCTTTTTTTTGATTGAAGGGAGGACATTCTgtaaagggggggcgg                                                    >  1:247117/1‑137 (MQ=255)                                                                                                      aCGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCAAAAAAACATTCCTTTTCTGTCACAACTGACGC‑AGGGTGTAAACGGCGGAGCGGTCTTCCTGCTGTTTGTTGATGGGAGGTGGAAC‑TTCTGCTAATGGCGCGGCGGa                                                  >  2:944539/1‑150 (MQ=255)                                                                                                         gAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCa                                              <  1:65751/151‑1 (MQ=255)                                                                                                          aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTGGCACAACTCACGC‑CGGGGGTAAACGGGGGAGGGCGCTTTCTTCTGTTTTTTGATGGGAGGGAACAC‑TTCTGCTAATGGCGcggcgg                                                   >  2:1090847/1‑145 (MQ=255)                                                                                                                 ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCa                                          <  1:217982/147‑1 (MQ=255)                                                                                                                 ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAg                                         <  2:693430/148‑1 (MQ=255)                                                                                                                   ttttGCAGTTGTTTGACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGtt                                       <  2:226902/148‑1 (MQ=255)                                                                                                                    tttGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGttttt                                    <  1:992877/150‑1 (MQ=255)                                                                                                                    tttGCAGTTGTTTCACGTTGTGGCGGGCACAAAAAAAATGTATTTCCTTACACAACCGACGC‑AGGGGATAAACCGCGTAACCCCCTTCATTCTTTTTGTAGGTGGGAGGAAGAAA‑ATCTTCTCATGGCGCGGCCGGCTCAAACAGattttt                                   >  1:191801/1‑150 (MQ=255)                                                                                                                     ttGCAGTTGTTTCACGTTGTGGCGGGCGTACAAAACATGCATTTTCTGACACAAACGACGC‑CGGGGGTAACCGGCGGAGGGGCCTTCTTTCTGTTTTTTGATGGGAGGGAGAAC‑CTTTTCTTATGGGGCGGCGGGCTCCAGCAGgtttt                                    >  2:118290/1‑149 (MQ=255)                                                                                                                     ttGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATGCATTTTCTTACACAAAATACGC‑CGGGGATAACCCGCGTAGCGCTCTTCTTGCTGTTTTTAGATGGAAGGGAGCAC‑TTCTGCTAATGGCGCGTCGGGCTCAAGAAGtttttt                                   >  1:278383/1‑149 (MQ=255)                                                                                                                      tGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCg                                  <  1:885987/150‑1 (MQ=255)                                                                                                                      tGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATTCATTTTCTGACACCACCCACGC‑AGGGGATAAACGGCGTAGCGCTCTTCATGCTGTTTGTTGATGGAAGGGAGCAC‑TTCTGTTTATGGCGGGGCGGGCTCCAACAGgttttt                                   >  2:1024067/1‑148 (MQ=255)                                                                                                                       gCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGaa                                <  2:426568/151‑1 (MQ=255)                                                                                                                        cAGTTGTTTCACGTTGTGGCGGGGGTAAAAAACATGCCTTTTCTTTCCCCACAGACGC‑GGGGGATAAACCGCGGAGCGCGCTTCATTTTGGTTTTTGGTGGAAGGGAGCAC‑TTCTGCTAATGGCGCGTCGGGCTCCAACAGttttt                                    >  2:873150/1‑146 (MQ=255)                                                                                                                        cAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAg                               <  2:427059/151‑1 (MQ=255)                                                                                                                            tgttTCACGTTGTGGCGGGCCTACAAAATAAGCATTTTCTCACACAAAACACGG‑AGGGGATAAAACGCGTGCCCCTCTTTTTGCTTTTTTTTGATTGTAGGAAGCAC‑CTCTTCTAAAGTGGCGTCCGGCTCAAGAAAtttttttatgcgt                            >  2:1071403/1‑142 (MQ=255)                                                                                                                                tCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACca                        <  1:1064510/150‑1 (MQ=255)                                                                                                                                 cACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACcag                       <  1:158177/150‑1 (MQ=255)                                                                                                                                      tgtgGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAggg                  <  2:802136/150‑1 (MQ=255)                                                                                                                                           cGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGTGcaa               <  2:1064033/148‑1 (MQ=255)                                                                                                                                           cGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca             >  2:188994/1‑150 (MQ=255)                                                                                                                                             ggCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGt          <  2:600361/151‑1 (MQ=255)                                                                                                                                                aTACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGc       <  2:246157/151‑1 (MQ=255)                                                                                                                                                aTACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGc       >  1:1052537/1‑151 (MQ=255)                                                                                                                                                 tACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCt      >  2:819447/1‑151 (MQ=255)                                                                                                                                                      aTCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAg  >  1:1013379/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGC‑AGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCAC‑TTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAG  >  NC_000913/718509‑718807

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