Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,703,805	G→A	G168G (GGC→GGT)	ydgJ ←	putative oxidoreductase YdgJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,703,805	0	G	A	98.1%	165.1 / NA	52	G168G (GGC→GGT)	ydgJ	putative oxidoreductase YdgJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/30);  minor base C (0/1);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTT  >  NC_000913/1703558‑1704043                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               cATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGaccg                                                                                                                                                                                                                                             >  1:273354/1‑251 (MQ=255)                                        aCCGTGCTGGCGCAGCGACGTCCAGCATTGTACCGTGTAAAATGACTTGCCGCAGTGGATAGGAAAAGATGGCGTGGAAACAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCTTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATTCCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCGGCGGACGGAAg                                                                                                                                                                                                      <  2:332365/251‑1 (MQ=255)                                         ccgggCTGACTCAGCGGCTGCCAGCATGGTTCCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGCCTGTGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGc                                                                                                                                                                                                     <  1:215155/247‑1 (MQ=255)                                         ccGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCgcgccgatcggcagagcgtc                                                                                                                                                                                                     >  2:338606/1‑235 (MQ=255)                                         ccGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCgcgc                                                                                                                                                                                                                     <  1:338606/235‑1 (MQ=255)                                                gACTCAGCGGCTGCCAGCTTGGTCCCGTGTAAAATGACTCGCCGCTGTGTATATGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCTGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCCGGACCGCCCTGTTCACGCCAACGATCTCGCACCTGCGGACGGAAGCgggca                                                                                                                                                                                                <  2:300134/249‑4 (MQ=255)                                                          gggcagtgaagtaaccgggcaaaaTGTCGCGCCGCTGTGGGTAGGACAAGATGGCGTGGACATCACCGGCCGACTGCGCTCCGGGCCGTAACGGTGCCAAATCTCCCGTCATGCTGACCAGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATgaggcgca                                                                                                                                                                                    <  2:72411/231‑6 (MQ=255)                                                             ccAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAg                                                                                                                                                                                 <  1:281370/251‑1 (MQ=255)                                                              cAGCTTGGTACCGTGTAAAAGGACTCGCCGCTGGGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCCCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGt                                                                                                                                                                                <  2:334855/251‑1 (MQ=255)                                                                   tGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGcg                                                                                                                                                                                                              <  2:99438/216‑1 (MQ=255)                                                                   tGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGcg                                                                                                                                                                                                              >  1:99438/1‑216 (MQ=255)                                                                                 atgtttctcctcgttgGATAGGAACAGATGGCGTGGAAATAATCGGGCGACGGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTAATGCTGACCGGTAAACCAAAAAGCGGAATGGTCTGATCAAGAAGATGTGGTGCTAAATCGTACCTGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAATTTCACCCAGCAc                                                                                                                                                             <  2:156564/237‑1 (MQ=255)                                                                                       cgccgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCtt                                                                                                                                                        <  1:246824/250‑1 (MQ=255)                                                                                           gcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCcgcg                                                                                                                                                   >  1:183819/1‑251 (MQ=255)                                                                                           gcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCcgcg                                                                                                                                                   <  1:245645/251‑1 (MQ=255)                                                                                             tgtgGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCAACCAGCAAGCCTTCcgct                                                                                                                                                   >  1:322751/1‑248 (MQ=255)                                                                                             ggtgGATAGGACAAGATGGCGTGGAAATAATCGGTCGACGGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTGCACCCAGCACGCCTTCCGCGAg                                                                                                                                                 <  2:157997/250‑1 (MQ=255)                                                                                                ggATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTaa                                                                                                                                              <  1:79478/251‑1 (MQ=255)                                                                                                 gATAGGACAAGAGGGAGGGGAAAAAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTTACCGGTAAACCAAAAAGCGTAAGGGCCTGATCAAGAAGATGTGGGGCTAAATCGTACCCTATACCGCTGCCTGGACCTCCCTGTTCACGCCAACGAGCTCTCACCTGCTGACGGAAGCGGGCTAAATGAGACGCAAAGTCAGCAACTTCACCCattacgtctttcgcgattaac                                                                                                                                             <  2:49186/251‑21 (MQ=255)                                                                                                      gACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCttt                                                                                                                                        >  1:283450/1‑251 (MQ=255)                                                                                                      gACAAGATGGCGTGGAAAGAATCGGTTGATGGCGCTCCGGGCCGTAAATGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTTGTGCTAAAGCGTACCAGATACCGCTGCCTGGCCCGCCCTGTGCACGCCAACGATCGCGCACCTGCGGACGGGCGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCttt                                                                                                                                        <  1:232495/251‑1 (MQ=255)                                                                                                                cGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCa                                                                                                                                >  1:299639/1‑249 (MQ=255)                                                                                                                   ggAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGaa                                                                                                                            <  1:271686/250‑1 (MQ=255)                                                                                                                    gAAATAATCGGTCGCCTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGATGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAAt                                                                                                                          <  1:127205/251‑1 (MQ=255)                                                                                                                        taatCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCt                                                                                                                      >  1:226719/1‑251 (MQ=255)                                                                                                                                cGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGc                                                                                                              <  1:53816/251‑1 (MQ=255)                                                                                                                                   cTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCCCCCAGCACGCCTTCCGCGAGTAAAcc                                                                                                                                           >  2:1320/1‑219 (MQ=255)                                                                                                                                   cTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAAcc                                                                                                                                           <  1:1320/219‑1 (MQ=255)                                                                                                                                                      ggtgCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGAGGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAATTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGGCAAGAAAGCGCTATCCCAGCGACGGTTATGGAATacagacag                                                                                        <  2:88743/250‑1 (MQ=255)                                                                                                                                                                  aCCGTCATGATGACCGGTAAACCAAAAAGCGTAATGGCATGATCAAGAAGAGGTGGTGCTAAACCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGTCGGTTATGGAATACAGACAGCACACGCCCCAg                                                                            <  1:319712/251‑1 (MQ=255)                                                                                                                                                                   ccGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGGGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCACCTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACTGTTATGGAATACAGACAGCACACGCCCCAgg                                                                           <  1:356945/251‑1 (MQ=255)                                                                                                                                                                             gACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCGGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTAACCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCa                                                                 <  1:270884/251‑1 (MQ=255)                                                                                                                                                                              aCCGGTAATCCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGACTCCTCTGCCTGGACCGCCCGTTTCACGCCAACGATCGCGCACCTTCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTTTCCCGGCACGCCTTCCGCGAGTAATCCTTTTACCGTCAAGAAAGCGCTATCCCAGCGTCGGTTTTGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAg                                                                <  1:402497/251‑1 (MQ=255)                                                                                                                                                                                cGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAgcg                                                              >  1:271905/1‑251 (MQ=255)                                                                                                                                                                                 ggtaaaaaaaaaaGCGTAATGGCCTGATCACGACGATGTGGAGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCAAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGGGACTCAAAGAAAGCAACTTAACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGCAATCGCTCTCCCAGCGACGGTTATGGAATACAGCCAGCACCCCCCCCCGGCTTTTTTCCAgcgc                                                             >  2:373151/9‑251 (MQ=255)                                                                                                                                                                                 ggTAAACCAAAAAGAGTAATGGCCTGATCAAAAAGATGTGGTGCGTAATCGTACCAGATTCCTCTGCCTGGACCCCCCTGCTCACGCCACCGATCGCGCACCTTCGCACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCCTCAAGAAATCGCTATCCCTGCGACGGTTATTGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAgcgc                                                             <  2:238589/251‑1 (MQ=255)                                                                                                                                                                                        cAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGATACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGc                                                      >  1:265092/1‑251 (MQ=255)                                                                                                                                                                                             aGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATCCCGCTGCCTGGACCGCACTGTTCACGCCAACGATCGCGCACCTGCGGACGGCAGCGGTCAAAATGAGCCTCAAAGTAAGCAACTTCCCCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAACTCGCTATCCCAGCGACGGTTATGGAATACAAACAGCACACTCCCCAGGCTTTTTGCCAGCGCAGCCAtctcccg                                                 >  2:220396/1‑248 (MQ=255)                                                                                                                                                                                                 aaaTGGCCTGAGCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCt                                             <  2:299639/250‑1 (MQ=255)                                                                                                                                                                                                          gATCAAGAAGATGTGGTGCTAAATCTTCCCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCGGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCATTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCCGACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACa                                      <  1:360599/249‑1 (MQ=255)                                                                                                                                                                                                                gaagaTGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTCTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACg                              >  1:211229/1‑251 (MQ=255)                                                                                                                                                                                                                 aagaTGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCACAGTAAGCAACTTCACCCAGCACGCCTTCCGCGCGTAAACCTTTTAGCGTCAAGAACTCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTtgtgcaggagtcacgg                             >  2:18122/1‑239 (MQ=255)                                                                                                                                                                                                                 aagaTGTGGTGCTAAATCGTACCAGATACCGATGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAAATTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCACGAAATAGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTATCACgg                             >  2:415891/1‑251 (MQ=255)                                                                                                                                                                                                                    aTGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTa                           >  1:360538/1‑250 (MQ=255)                                                                                                                                                                                                                        ggtgCTAAATCGTACCAGATCCCGCTGCCTGGACCGCCCTTTTCACGCCAACGATCGCGCCCCTGCGGACGGAAGCGGTCAAAATGATACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAgg                       <  2:265092/250‑1 (MQ=255)                                                                                                                                                                                                                                  cGTACCAGATACCGCTTCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGTTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGTTTTATCAAc             <  2:183819/250‑1 (MQ=255)                                                                                                                                                                                                                                   gTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATacagac                                                                                          >  1:312195/1‑172 (MQ=255)                                                                                                                                                                                                                                   gTACCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGGCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATacagac                                                                                          <  2:312195/172‑1 (MQ=255)                                                                                                                                                                                                                                     aCCAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGCGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGCCAGTGTCACGGTACAGGGTTTATCACCGAcc         >  1:290660/1‑251 (MQ=255)                                                                                                                                                                                                                                      ccAGATACCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGAcc         >  1:104739/1‑250 (MQ=255)                                                                                                                                                                                                                                            aCCGCTGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGtt  >  1:223051/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTa                                                                                                      >  2:103591/1‑146 (MQ=255)                                                                                                                                                                                                                                                 tGCCTGGACCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTa                                                                                                      <  1:103591/146‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               CATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTT  >  NC_000913/1703558‑1704043

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