Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	718,669	0	C	T	42.6%	85.4 / 32.7	68	E598K (GAA→AAA)	speF	ornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base C (10/29);  new base T (29/0);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 4.87e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAT  >  NC_000913/718519‑718818                                                                                                                                                       |                                                                                                                                                       cATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTc                                                                                                                                                        <  2:871782/151‑1 (MQ=255)      ggCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGac                                                                                                                                                   <  1:906931/151‑1 (MQ=255)       gCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGAAACAAAAACTGCATTTTCTTac                                                                                                                                                   >  2:105753/1‑150 (MQ=255)       gCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGaca                                                                                                                                                  <  1:1258596/151‑1 (MQ=255)         aGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGGTACAAAACATGCTTTTTCtttcac                                                                                                                                                 >  1:259341/1‑145 (MQ=255)          gACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGAcaca                                                                                                                                                <  1:1008212/150‑1 (MQ=255)           aCGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAAc                                                                                                                                              <  2:858014/151‑1 (MQ=255)            cGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACt                                                                                                                                             <  2:1021497/151‑1 (MQ=255)             gAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGACTTTTCTGACCCAAcc                                                                                                                                             >  1:625380/1‑149 (MQ=255)                 cAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGc                                                                                                                                        <  2:1199500/151‑1 (MQ=255)                   gTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGTGAAACAAAACATATATTTTTTTACACAACTGAccc                                                                                                                                        >  1:953469/1‑147 (MQ=255)                     tCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAgg                                                                                                                                     <  1:1052620/150‑1 (MQ=255)                     tCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAgg                                                                                                                                     <  1:157209/150‑1 (MQ=255)                       cACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATTATGCATTTTCTTTCACAAATGAAGCAggg                                                                                                                                    >  2:855067/1‑149 (MQ=255)                        aCTTCACCGCGTAAATAGGGGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCGTACGGAACATCTCTTGTTGCAGTTGTATCACGTTGTGGCGGGCATACAAATCATGCATTGCCTGACACAACTGACGCAggg                                                                                                                                    <  1:983351/148‑1 (MQ=255)                         cTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGCTACAAAACATGCATTTCTTGAGACAACTAGCGCCGGGgat                                                                                                                                 >  2:907555/1‑150 (MQ=255)                         cTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAAAAAAAATGCATTTTCTTTCACAAATAACGCAaggg                                                                                                                                   >  2:682752/1‑147 (MQ=255)                             aCCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATTCATTTTCTTCCACAACAGACGCCgggggtaaa                                                                                                                              >  1:678255/1‑143 (MQ=255)                              ccGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTACAAATAATGCATTTTCTTGCACAACCGAACCAgggggtaaa                                                                                                                              >  1:706162/1‑142 (MQ=255)                               cgcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTACAAAAAATACACTTTCTGTCAAAACTGACGCCgggggtaaa                                                                                                                              >  2:1071440/1‑141 (MQ=255)                               cgcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGc                                                                                                                          <  2:74536/151‑1 (MQ=255)                                 cgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGt                                                                                                                        <  2:851511/151‑1 (MQ=255)                                   tAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACgg                                                                                                                           >  2:1214980/1‑146 (MQ=255)                                    aaaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCGAAAAAAAATGCATTTTCTTTCAAAAACGACGGCGGGGGTAAACGGCGTAgg                                                                                                                     >  1:45799/1‑150 (MQ=255)                                    aaaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATTCATTTTCTTACCCAACTTAAGCAGGGGATaa                                                                                                                               >  2:1236811/1‑141 (MQ=255)                                         ggCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCCTTTTCTGACACAAACTACGCAGGGGATAAACGGCGGAgggccc                                                                                                                 >  1:598818/1‑148 (MQ=255)                                              aGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTTGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATg                                                                                                           <  1:133919/151‑1 (MQ=255)                                              aGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATTCATTTTCTGACACAAATGAAGCAGGGGATAAAAGGCGGAGGGCTCtt                                                                                                               >  2:147789/1‑147 (MQ=255)                                                   gCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATTCATTTTCTTTCACAACATCCGCAGGGGGTAAACCGCGGAGCGCTCTTCTTGCTGtt                                                                                                      >  2:612440/1‑151 (MQ=255)                                                     ttcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCATTTTCTTTCCCAACAGACGCAGGGGATAACCCGCGGAGAGCTCTTCATGCTGttt                                                                                                     >  2:183566/1‑150 (MQ=255)                                                      tcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGTAGGGTATAACCGGCGTAGCGCTCTTCATGCTGttt                                                                                                     <  1:682038/149‑1 (MQ=255)                                                         tGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACGGCGGAGCGCTCTTCATGCTGttttt                                                                                                   >  2:966153/1‑146 (MQ=255)                                                          gCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACAAGCCTTTTCTTACACAAATTACGCCGGGGATAAACCGCGTAGCGCTCTTTCTGCTGttt                                                                                                     >  2:21463/1‑145 (MQ=255)                                                                 cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATATAAAAAATGCCTTTTCTTACACAAACAACCCAAGGGATAAACGGCGCTGCGGCCTTTATTCTGTTTGTTGATGGGAgg                                                                                         >  2:1040740/1‑150 (MQ=255)                                                                     cTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAAAAAAAAAGCAATTTCTTACACAAATCAAGCAGGGGGTAAAACGCGTAGCGCCCTTTTTGCTGTTTTTTGATGGAAgg                                                                                         >  1:1162152/1‑146 (MQ=255)                                                                          gcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCCAAAAATACTGCATTTTCTTACACAACTCACGCAGGGGATAACCGGCGTAGCGCCCTTCTTGCTGTTTTTTGATggggggaaggaccttt                                                                               >  1:832400/1‑137 (MQ=255)                                                                              gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTa                                                                           <  1:195579/151‑1 (MQ=255)                                                                                gAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATGCATTTTCTGACACAAATGACGCCGGGGATAAACCGCGGAGCGCCCTTCATGCTGTTTGTTGATGGGAGGGAGGACTTCTGCtta                                                                          >  1:1110443/1‑148 (MQ=255)                                                                                 aaGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATg                                                                        <  1:255472/151‑1 (MQ=255)                                                                                      cTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGGCAAAAAAAATTTCTTTTCTGTCCCAACCGACGGAGGGGGGAAACAGCGGAGAGGCCTTTCTGCTGGTTTTAGAAGGAAGGGAGGAAGTCTGCTAATgggg                                                                     >  1:740605/1‑147 (MQ=255)                                                                                        ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTc                                                                 <  2:575506/151‑1 (MQ=255)                                                                                          cttccggaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCg                                                                <  1:958572/143‑1 (MQ=255)                                                                                           tACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTGACACAACTGACGCAGGGGATAAACGGCGGAGCGCGCTTCATGCTGTTTGTTGATGGGAGGGAGGACTTCTTCTAAATGCGcggcgg                                                               >  2:92837/1‑150 (MQ=255)                                                                                            aCGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCTAAAAAATCATTCGTTTTCCCATAAAACTGTCGCAGGGGATAAACCGCGCAGCGCTTCTTATGCGGTTTTTAGAAGGGAAGGAGCACCTTTTCTAATTGCGCGTCgg                                                               >  2:873670/1‑149 (MQ=255)                                                                                               gAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCa                                                         <  1:361756/151‑1 (MQ=255)                                                                                                aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCaa                                                        <  2:1273762/151‑1 (MQ=255)                                                                                                 aCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAg                                                       <  2:130933/151‑1 (MQ=255)                                                                                                   aTCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGc                                                      <  1:166818/150‑1 (MQ=255)                                                                                                      tctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATAATGCATTTTCTTACACAAATGACGCAGGGGATAAACGGCGGAGCGCTCTTCCTGCTGTTTGTAGATGGAAGGGAGGACTTCTGCTAATTGGGCGGCGGGC‑TCa                                                         >  2:850796/1‑144 (MQ=255)                                                                                                       ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAAAGCATTTCCTGACACAAAAGACGCCGGGGGTAAACGGCGGAGCGGTCTTTATGCTGGTTTTAGATGGGAGGGAGGACTTCTGCTAATGGGGCGTCGGAC‑TCAAACAGgtt                                                 >  2:399775/1‑151 (MQ=255)                                                                                                            tGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATTCATTTTCTGATACAACTGACGCCGGGGATAAACCGCGGAGCGCTCTTCATGCTGTTTGTAGATGGGAGGGAGCACTTCTGCTAATGGCGCGGCGGGC‑TCAAGCAGgtttt                                               >  1:245167/1‑148 (MQ=255)                                                                                                             gCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGaa                                           <  1:231620/151‑1 (MQ=255)                                                                                                             gCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGaa                                           <  1:76872/151‑1 (MQ=255)                                                                                                              cAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGCATTTCCTGACACAAACGACGCAGGGGGAAAAAGGCGGAGCGCCCTTCATTCTGTTTGTAGATGGAAGGGAGGACTTCTTCTAATAGCGCGGCGGGG‑TCAAGAAGgtttt                                               >  1:1055293/1‑146 (MQ=255)                                                                                                              cAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGGATTTTCTGAGACAAACGACGCAGGGGGTAAACGGCGGAGCGCTCTTCATGCTGTTTGTAGATGGGAGGGAGGAATTTTGCTTATTGGGCGGCGGGC‑CCAAGCAGgttttt                                              >  1:483629/1‑146 (MQ=255)                                                                                                                gttgttTCACGTTGTGGCGGGCATACAAAACATGCATTTCCTTACACCAAATACGCCGGGGATAACCGGCGGAGCGCCCTTCCTGCTGTTTTTTGGTGGGAGGGAGCAACTCTGCTAATGGCGCGGCGGGCTTCAACAGGTTTTTCGag                                           >  1:611596/1‑148 (MQ=255)                                                                                                                         cGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACcag                                  <  1:1297403/148‑1 (MQ=255)                                                                                                                          gTTGTGGCGGGCATACAAATAATGCATTTTCTTCCACAACTGACGCAGGGGGTAAACGGCGGAGCGCCCTTCCTTCTTTTTGTTGATGGAAGGTAGCACTTCCTCTTATTGCGCGCCGGGC‑TCTAGCAGGTTTTTGTAACGTAcccaccgg                              >  1:454482/1‑145 (MQ=255)                                                                                                                             gtgGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGc                            <  2:215116/150‑1 (MQ=255)                                                                                                                                  gggCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca                        >  2:739516/1‑149 (MQ=255)                                                                                                                                   ggCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGt                     >  2:329362/1‑151 (MQ=255)                                                                                                                                          aaaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa               >  2:162709/1‑150 (MQ=255)                                                                                                                                           aaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAg             <  2:799485/151‑1 (MQ=255)                                                                                                                                              cATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGttt          >  2:599194/1‑151 (MQ=255)                                                                                                                                              cATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGttt          <  1:584449/151‑1 (MQ=255)                                                                                                                                              cATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGtt           <  1:254218/150‑1 (MQ=255)                                                                                                                                              cATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGtt           <  2:86074/150‑1 (MQ=255)                                                                                                                                                 gCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGc       >  1:255113/1‑151 (MQ=255)                                                                                                                                                  cATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGcc      >  1:862611/1‑151 (MQ=255)                                                                                                                                                       ccTGACACAACTGAGGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCatat  <  2:72672/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGAC‑TCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAT  >  NC_000913/718519‑718818

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