Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,489,174	0	G	T	46.8%	87.4 / 29.3	62	G315C (GGT→TGT)	aldA	aldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base G (6/27);  new base T (29/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 4.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGATATTATT  >  NC_000913/1489030‑1489324                                                                                                                                                 |                                                                                                                                                       gAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTaaac                                                                                                                                                  >  2:204690/1‑149 (MQ=255)      ggCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCg                                                                                                                                               <  1:7624/149‑1 (MQ=255)        cAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAACCCGc                                                                                                                                              >  2:862201/1‑148 (MQ=255)        cAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAAcccc                                                                                                                                                >  2:831118/1‑146 (MQ=255)           tCAAAGCCATCGTTGATTAACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGGTCAATTTTGTAAcccccgcg                                                                                                                                            >  1:205782/1‑143 (MQ=255)             aaaGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGGTTAACTTGGGAAACCCCCCGaa                                                                                                                                          >  1:141429/1‑147 (MQ=255)               aGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGAACGCaa                                                                                                                                     >  2:355858/1‑150 (MQ=255)                 ccATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTTACCCCGCCGaaag                                                                                                                                        >  1:604594/1‑143 (MQ=255)                  cATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGAACGCAACGAc                                                                                                                                 >  2:969412/1‑151 (MQ=255)                   aTCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCCATTTTGTAAACCCGCTGAACGCAACGAc                                                                                                                                 >  2:10781/1‑150 (MQ=255)                   aTCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACa                                                                                                                                <  2:597543/151‑1 (MQ=255)                               cgcgTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGc                                                                                                                    <  2:820947/151‑1 (MQ=255)                                gcgTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGAACGCCACGACATTTCGATgggggc                                                                                                                   >  1:368863/1‑149 (MQ=255)                                     aTTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAAATTTGTAAACCCGCCGAACGCAACGAAAATTCGATGGGGCCGTTg                                                                                                               >  1:1011025/1‑150 (MQ=255)                                       tAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGAt                                                                                                             <  1:581566/150‑1 (MQ=255)                                              gggCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAAcgccg                                                                                                     <  2:1078787/151‑1 (MQ=255)                                              gggCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAAcgcc                                                                                                      <  1:45108/150‑1 (MQ=255)                                                      gtgtAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTg                                                                                             <  2:658072/151‑1 (MQ=255)                                                          aaCTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGaaa                                                                                         <  2:717718/151‑1 (MQ=255)                                                           aCTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGAACGCAACGACAATTGGATGGGGGCGTTGGTTa                                                                                                           >  2:974179/1‑132 (MQ=255)                                                              gtgCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAACTTTGTAACCCCCCCTGACGCAACAAGATTTCGATTGGGCCGATGATTAAAGCCGCCGCGCTGGGAAggg                                                                                      >  2:859659/1‑150 (MQ=255)                                                                gCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAAATTTGGAAACCCGCTGAACGCAACGACATTGCGATGGGGGCGTTGATTAACGCCCCGGCGCTGGAAAggg                                                                                      >  1:49772/1‑148 (MQ=255)                                                                    aaCGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAACCCCGCTGAACGCAACGACATTGCGGTGGGGCCGTTTTTTTACGCCGCCGCGCTGGGAAGGGGCggg                                                                                 >  1:440653/1‑147 (MQ=255)                                                                      cGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAACCCCCCTGAACGCAACGACATTTCGGTGGGGCCGGTGGTTAACCCCCCGGCGCCGGGAAGGGGCggg                                                                                 >  1:1013424/1‑145 (MQ=255)                                                                         gtTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCaaaaa                                                                           <  2:795183/150‑1 (MQ=255)                                                                           ttATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCCCAACGCAACGACATTGCGCTGGGGGCGGTGATTAAAGCCCCGGCGCTGGGACGGGGCGGGCaaaa                                                                            >  1:619738/1‑147 (MQ=255)                                                                            tATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTg                                                                        <  2:925481/150‑1 (MQ=255)                                                                                 aCAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTTATTAACGCCGCGGCGCTGGGAAGGGTCGAGCAAAAAAGGgcgcg                                                                   >  2:142322/1‑150 (MQ=255)                                                                                  cAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAACCCCGCTGAACGCCACCACCTTTTCATGTGGCCCTTTATTTAACCCCCGGCGCTGGGAAGGGTGGAGGCAAAAAAggcgggcg                                                                 >  2:994362/1‑151 (MQ=255)                                                                                   aGAAAGGCATTTATGATTAGTTCGTCAATCTGCTAGGTGAAGCGATGCAGGCGGTTCAATTCGGTAACCCCGCTGAACGCATCGCCACCGCGATGGGCCCGTTGATTAACGCGGAGTCTCTGGAAAGGGTCGAGCAAAATGTGgcgctcgc                                                                <  1:177819/151‑1 (MQ=255)                                                                                       aGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGGACGCAACGACATTGCGGTGGGGCCGTTGATTAACACCGCGGCGCTGGGAAGGGTCGAGCGAAAAATGGCGCGCGCAg                                                              >  2:919156/1‑149 (MQ=255)                                                                                         gCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCCCGGCGCTGGGAAGGGTCGAGCAAAAAATGGCGCGCGCCGt                                                             >  2:655615/1‑148 (MQ=255)                                                                                           aTTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCGCTGAACGCAACGACATTTCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGGAAGGGGGGAGCAAAAAGTGGCGCGGGCAgga                                                            >  2:143221/1‑145 (MQ=255)                                                                                             ttATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAACCCCCCTTAACGCAACGACATTGCGGTTGGGGCGTTGATTAACGCCGCGGCGCTGGAAAAGGTCGGGCCAAAAGTGgggcgggccggagga                                                         >  1:879038/1‑141 (MQ=255)                                                                                             ttATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTTACCCCGCTGAACGCAACGGCAATGGGATGGGGCCGTTGATTAACCCCCCGGCGCTGGGAAAGGTCGAGGAAAAAATGGCGCGCGCAGTagaa                                                         >  2:771562/1‑148 (MQ=255)                                                                                                 gATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAAATTTGTAAACCCGCTGAACCCAACGGCATTGCGATGGGGCCGGTGATTAACCCCCCCGCGCTGGAAAGGGGCGAGCAAAAAGTGGGGCGCGCGATAGAAGAAgggg                                                  >  2:223051/1‑151 (MQ=255)                                                                                                  aTCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGTAAACCCGCCGAACGCGACGACATTGCGCTGGGGCCGTTGGTTAACGCCGCGGGGCTGGAAAAGGTCGGGCAAAAAGTGGCGCGCGGCGGAGAAGAAgggg                                                  >  2:175136/1‑150 (MQ=255)                                                                                                   tCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGAACGCCACCACCTTGCGATGGGGCCGTTGATTAACGCCGCCGCGCTGGGAAGGGTCGAGCCAAAAATGGCGCGCGCCGTAGAAGAAGGGGCg                                                >  2:811118/1‑151 (MQ=255)                                                                                                     aGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGATCAAATTGGTAACCCCCCGGAAAGCCACGGGAATGAGATGGGGCCGGTGGTAAAAGCCGCGGCGCTGGGAAGGGGCGAGCGAAAAATGGCGCGGGCAGTcggaga                                                       >  2:1001013/1‑142 (MQ=255)                                                                                                      gTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTTAATTTTGGAACCCCGCTGAAAGCAACGACATTTCGATGGGGCCGTTGGTTAACGCCGCGGCGCCGGAAAGGGGCGAGCAAAAAGGGGCGCGCGCAGGAGAAGAAgggggggg                                              >  1:674200/1‑146 (MQ=255)                                                                                                      gTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCCATTTTGGAACCCCGCCGAACGCAACGACATTGCGAAGGGGCCGTTGGTTAACGCCCCGGCGCTGGGAAGGGGCGAGCAAAAAATGGCGCGGGCCGGAGAAGAAGGGGCggg                                              >  1:148122/1‑148 (MQ=255)                                                                                                         cGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGt                                           <  1:763955/150‑1 (MQ=255)                                                                                                                cGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTc                                    <  1:970764/150‑1 (MQ=255)                                                                                                                cGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTc                                    <  2:779129/150‑1 (MQ=255)                                                                                                                cGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCg                                   <  2:1018320/151‑1 (MQ=255)                                                                                                                 ggCTGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAACCCCCCTGAACGCAACGACATTGCGGTGGGGCCGTTGATTAACGCCGCGGCGCTGGGAAGGGGCGGGCAAAAAGGGGCGCGCGCCGTAGAAGAAGGGGGGGGGGGGGCGGTCgg                                  >  2:145985/1‑151 (MQ=255)                                                                                                                    tGGGTGAAGCGATGCAGGCGGTTCAATTTTGTAAACCCGCTGAACGCAACGACATTGCGATGGGGCCGGTGATTAAAGCCGCGGCGCTGGAAAGGGGCGAGCAAAAAATGGCGCGCGCAGTAGAAGGAGGGGGGAGAGGGGGGGTCGgggg                               >  2:214272/1‑151 (MQ=255)                                                                                                                    tGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGGCGTTGATTAACGCCGCGGCGCTGGGAAGGGTCGAGCAAAAAAAGGCGCGCGCAGTAGAAGAAGGGGGGAGAGGGGGGGTCGgggg                               >  1:76396/1‑151 (MQ=255)                                                                                                                     gggTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTg                                <  1:832360/149‑1 (MQ=255)                                                                                                                        tGAAGCGATGCAGGCGGTTCTATTTTGTAACCCCGCTGAACGCAAACACATTGCGGTGGGGGCGTTGATTAACCCCGCCGCGCTGGGAAGGGTCGAGCAAAAAATGGCGCGCGGCGTAGGAGAAagggggggggggggg                                       >  1:799755/1‑128 (MQ=255)                                                                                                                          aaGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGc                         <  2:873783/151‑1 (MQ=255)                                                                                                                              gATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTa                     <  1:870835/151‑1 (MQ=255)                                                                                                                                  cAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAgg                 <  1:259796/151‑1 (MQ=255)                                                                                                                                  cAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAgg                 <  2:911972/151‑1 (MQ=255)                                                                                                                                   aGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAggg                <  2:260206/151‑1 (MQ=255)                                                                                                                                     gCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGa              <  2:423477/151‑1 (MQ=255)                                                                                                                                     gCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGa              <  1:804218/151‑1 (MQ=255)                                                                                                                                       ggTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGATCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGaaa            <  2:1074699/151‑1 (MQ=255)                                                                                                                                           cAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGat        <  2:161551/151‑1 (MQ=255)                                                                                                                                           cAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGat        <  1:969347/151‑1 (MQ=255)                                                                                                                                           cAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAACAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGAGGCAAAGCGGGAGGGGGGAAAGGa         >  1:130772/1‑150 (MQ=255)                                                                                                                                            aaTTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGata       <  1:21029/151‑1 (MQ=255)                                                                                                                                                 ggTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGAtattatt  <  1:332122/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                       GAACTGGCAGTCAAAGCCATCGTTGATTCACGCGTCATTAATAGTGGGCAAGTGTGTAACTGTGCAGAACGTGTTTATGTACAGAAAGGCATTTATGATCAGTTCGTCAATCGGCTGGGTGAAGCGATGCAGGCGGTTCAATTTGGTAACCCCGCTGAACGCAACGACATTGCGATGGGGCCGTTGATTAACGCCGCGGCGCTGGAAAGGGTCGAGCAAAAAGTGGCGCGCGCAGTAGAAGAAGGGGCGAGAGTGGCGTTCGGTGGCAAAGCGGTAGAGGGGAAAGGATATTATT  >  NC_000913/1489030‑1489324

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