Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,776,263	0	T	G	42.2%	90.3 / 16.4	65	V226G (GTG→GGG)	ydiF	putative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base T (8/29);  new base G (27/0);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 4.46e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTTG  >  NC_000913/1776114‑1776412                                                                                                                                                       |                                                                                                                                                        cAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTg                                                                                                                                                        <  1:845306/151‑1 (MQ=255)   aaGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCa                                                                                                                                                      <  2:1126064/151‑1 (MQ=255)     gCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGa                                                                                                                                                    <  1:825227/151‑1 (MQ=255)       gATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGCATTTTGATGATGT‑GGGGGCAGa                                                                                                                                                    >  1:209057/1‑149 (MQ=255)        aTTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTCTTTTGTTGATGC‑AGGGGCAGgaa                                                                                                                                                  >  2:361398/1‑150 (MQ=255)          tGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAAt                                                                                                                                                <  1:1225423/150‑1 (MQ=255)            cgcCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGATGC‑CGGGGCCGaaaaatg                                                                                                                                              >  2:1194781/1‑147 (MQ=255)               cAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTaa                                                                                                                                          <  2:66021/151‑1 (MQ=255)               cAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTaa                                                                                                                                          <  1:212227/151‑1 (MQ=255)                aGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTaa                                                                                                                                          <  1:1034122/150‑1 (MQ=255)                       gCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGcc                                                                                                                                   <  2:252368/150‑1 (MQ=255)                       gCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGcc                                                                                                                                   <  1:549834/150‑1 (MQ=255)                        ccTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCAc                                                                                                                                 <  1:369845/151‑1 (MQ=255)                         cTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCAc                                                                                                                                 <  1:1106653/150‑1 (MQ=255)                         cTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTTTGGTGTTGC‑AGGGGCGGAAAATGGGTTAGaaa                                                                                                                                      >  1:962817/1‑145 (MQ=255)                          ttCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTTTGATGATTC‑AGGGGGGGAAAAAGGTTAAGaaaaacaccc                                                                                                                               >  1:932965/1‑144 (MQ=255)                           tCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGATGC‑CGGGGGAGAAAAAGGGTAAGAAAACCACGc                                                                                                                               >  1:602291/1‑150 (MQ=255)                           tCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACg                                                                                                                                <  1:926296/149‑1 (MQ=255)                               tCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGGGGAGAAAAAGGTTAAAAAAGCCACGCTg                                                                                                                             >  1:1103340/1‑148 (MQ=255)                                cgcgCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTTGGGTGGTGT‑AGGTGGAGAAAAAGGGTAAGAAAAACACGCTGCCTc                                                                                                                         >  2:1014174/1‑151 (MQ=255)                                      accaccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGATGC‑AGGGGCCGGAAAAGGGTAAAAAAAACACCCTGCCTCCTcaa                                                                                                                    >  2:111519/1‑150 (MQ=255)                                       ccaccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGC‑AGGGGCAGAAAAAGGGTTAGAAAACCACGCCGCATTCCaaa                                                                                                                    >  1:232884/1‑149 (MQ=255)                                        caccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGTTGATGC‑AGGGGCAGGAAATGGGTAAGAAAACCACGCCGCAACCTAAAACt                                                                                                                 >  1:904968/1‑151 (MQ=255)                                          ccTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGt                                                                                                               <  2:790992/151‑1 (MQ=255)                                                      aaGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCggg                                                                                                 <  2:1281449/151‑1 (MQ=255)                                                       aGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATACCGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGt                                                                                                <  2:276547/151‑1 (MQ=255)                                                           tACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGTTGATTG‑AGGGGGAGAAAAAGGGTAAGAAAAACACGCGGCCTCCTAAATATGTTCGT‑ATTTCCGGGtttt                                                                                             >  2:953653/1‑149 (MQ=255)                                                           tACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGGGCACAATAACGGCGGCGTTGTGAATTTGC‑ACGGGGAGAAAATGGTTTAAAAAAACACCCTCCCTCCTTAATCTGGCCTT‑G‑ATCCgggg                                                                                                >  1:1216605/1‑146 (MQ=255)                                                                cACGTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATTC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCtggtg                                                                                       <  2:1211069/151‑1 (MQ=255)                                                                 aCTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTTTTTTGATGATGC‑AGGGGCAGGAAATGGGTAAAAAAACCACGCCGCATCCTAAAACTTTCCGT‑A‑TTTCGGGGTATTtggggg                                                                                      >  2:48818/1‑151 (MQ=255)                                                                  cTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGGGGTGCACAATAACGGCGGTATTGTGATGATTC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TGCCGGGTTATCtggtga                                                                                      <  2:230439/150‑2 (MQ=255)                                                                     ttGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGata                                                                                   <  2:506713/150‑1 (MQ=255)                                                                     ttGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGata                                                                                   <  1:974248/150‑1 (MQ=255)                                                                      tGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGGTGATGC‑AGGGGTAGAAAAAGGGTAAGAAAACCACCCTGCCGTCTAAAACTTTTCGG‑G‑TTTCGGGGTATATGGTGGGTAtt                                                                                 >  2:995424/1‑151 (MQ=255)                                                                       gAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGGGCAGAAAATGGGTAAGAAAACCACGCTGCATCCTAAATCTTTCCGT‑A‑TTTCGGGGTATCTGGTGGGTAtt                                                                                 >  1:375011/1‑150 (MQ=255)                                                                         aGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGC‑AGGGGGAGGAAAAGGTTAAGAAAGCCACCCCGGCTCCTAAATCTTTTCGC‑A‑TTCCGGGGTATCTGGGGGGTAtt                                                                                 >  1:1218936/1‑148 (MQ=255)                                                                             gAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGC‑AGGGGCAGAAAAAGGTTTAGAAAACCACGCTGCATCCTAAAACTGTCCGT‑A‑TTCCGGGGTATCTGGTGGGTATTGtggggg                                                                          >  1:36219/1‑151 (MQ=255)                                                                                    tGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATATTGAGGATGCGAGGTG‑AGAAAATTGCTTAAAAAACCACGCTTAATCCCAAAACCGTCCGG‑A‑ATTCGGGTTATTTGGTGGCTATTTTGGTGGGCGa                                                                      >  1:869517/1‑148 (MQ=255)                                                                                      tATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTTTTGTGGTGAAGG‑AGGGGCCGAAAATTGGTAAAAAAACCACGCCGCATACTAAATCTGTCCGT‑A‑TTCCGGGGTTTTTTGGGGGTATTGTGGGGGGCGGTcccg                                                                 >  2:589895/1‑149 (MQ=255)                                                                                         ctcGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGC‑AGGGGGCGAAAAAGGGTAAGAAAACCACGCTGCATCCTTAAACTTTCCGT‑A‑TTCCGGGGTATCTGGGGGGTATTGTGGTGGTCGATcccggt                                                               >  1:9683/1‑146 (MQ=255)                                                                                          tcGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGTTGA‑AGGGGGAGAAAAAGGTTAAGAAAAGCCCGCCGCATCCTAAAACTGTTCGT‑A‑TTCCGGGGTTTTTGGTGGGATTTTTGGTGGGCGAGCCGGGTc                                                              >  1:936338/1‑150 (MQ=255)                                                                                              cGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGGATTGTGATGATGG‑AGGGGCAGAAAAAGGGTAAGAAAGCCACGCTGCATCCTAAAATTTTTCGG‑A‑TTTCGGGGTATCTTGTGGGTATTGTGGGGGTCGATCCGGATCAAAcc                                                         >  1:706981/1‑151 (MQ=255)                                                                                               gCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTTATGGTGC‑AGGGGCAGAAAAAGGGTAAGAAAGCCACGCTGCATCCCAAATCTGTTCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGCTCAAAccc                                                        >  1:551114/1‑151 (MQ=255)                                                                                                 aTTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTTTGATGATGC‑AGGGGCAGAAAATGGGTAAGAAAACCACGCTGCCTCCTAAATCTTTCCCT‑A‑TTCCGGGGTATCTGGGGGGTATTGTTGTGGTCGATCCCGATCAAAccc                                                        >  1:215469/1‑149 (MQ=255)                                                                                                   tGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGGGCAGAAAAAGGTTAAGAAAGCCCCGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGGTATCTGGTGGGTATTTTGGTGGTCGATCCGGATCAAAACCCAc                                                     >  1:792699/1‑150 (MQ=255)                                                                                                      ttATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTGATGG‑AGGGGCAGAAAATGGGTAAAAAAACCACGCCGCATCCTAAATCTTTCCGT‑A‑TTCCGGGGTTTCTTGTGGGTATTTTGGGGGGCGATCCGGGTCAAAACCaa                                                      >  1:960917/1‑146 (MQ=255)                                                                                                        aTTGCCCAGGCGGTGCACAATAACGGCGGTCTTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTAt                                                <  2:34591/150‑1 (MQ=255)                                                                                                             ccAGGCGGTGCACAATAACGGCGGTATTTTGATGTTGG‑CGGGGGAGGAAAAGGGTAAAAAAACCCCGCCGACTCCTTAATATGTCCCTAA‑TTTTGGGTT‑TTTGGGGGGGATTGTGGTGGTCGGGCCCGATCAAACCCCAccattttggggg                                          >  2:752256/1‑140 (MQ=255)                                                                                                               aGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGc                                        <  1:175915/151‑1 (MQ=255)                                                                                                                        acaATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTaa                               <  2:1113782/151‑1 (MQ=255)                                                                                                                           aTAACGGCGGTATTGTGATGATTC‑AGGGGCAGAAAAAGGGGTAGAAAACCACGCCGCATCCTAAATCTGTCCGT‑A‑TTCCCGGTTATTTGGGGGGGATTGTGGTGGTGGATCCGGGTCCCAACCAACTGTATGGGGGGGGGCCGGGTAAcc                             >  2:1173100/1‑150 (MQ=255)                                                                                                                             aaCGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGc                           <  1:882749/150‑1 (MQ=255)                                                                                                                                  cggTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTAtt                     >  2:1283312/1‑151 (MQ=255)                                                                                                                                          tgatgaTGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTg              >  2:271609/1‑150 (MQ=255)                                                                                                                                            atgatgC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACt           <  2:1024468/151‑1 (MQ=255)                                                                                                                                            atgatgC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACt           <  2:722347/151‑1 (MQ=255)                                                                                                                                            atgatgC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACt           <  2:602266/151‑1 (MQ=255)                                                                                                                                             tgatgC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGTTGACtt          >  2:943717/1‑151 (MQ=255)                                                                                                                                             tgatgC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACtt          >  1:1103372/1‑151 (MQ=255)                                                                                                                                             tgatgC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACt           <  1:582494/150‑1 (MQ=255)                                                                                                                                              gatgC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACtt          >  2:158346/1‑150 (MQ=255)                                                                                                                                                 gC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCAcc      >  2:486108/1‑151 (MQ=255)                                                                                                                                                    aGGCGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGCGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCCTTATTTCTGGGGACTTCACCCt    >  2:54202/1‑151 (MQ=255)                                                                                                                                                     ggTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCtt   <  1:72017/151‑1 (MQ=255)                                                                                                                                                      gTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTtg  <  1:551961/151‑1 (MQ=255)                                                                                                                                                       tGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTtg  >  1:1039128/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                        CAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGC‑AGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGT‑A‑TTCCGGGTTATCTGGTGGATATTGTGGTGGTCGATCCGGATCAAACCCAACTGTATGGCGGTGCACCGGTTAACCGCTTTATTTCTGGTGACTTCACCCTTG  >  NC_000913/1776114‑1776412

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