Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351	A→G	L206L (CTA→CTG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351	0	A	G	98.2%	168.2 / ‑4.4	54	L206L (CTA→CTG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (28/25);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

ATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCC‑GGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTC  >  NC_000913/4104‑4592                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   aTAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTgggg                                                                                                                                                                                                                                                 <  2:225938/251‑1 (MQ=255)   aaGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGACGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGtt                                                                                                                                                                                                                                               >  1:353567/1‑251 (MQ=255)    aGCCAGTGACCATTCTTACCGCGACCTCCGGTGATTCCGGAGCTGCAGTGGCTCATGCTTTCTGCGGGTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGGGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTa                                                                                                                                                                                                                                              <  2:69892/251‑1 (MQ=255)    aGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTa                                                                                                                                                                                                                                              <  2:28287/251‑1 (MQ=255)           gACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCCCGAGGCAAAATCAGTCCACTGCAAGAAACACTGTTCTGTACATTGGGCGGCAATATCGAAAATGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGCAGAGCTGAACGTGGCGCTGGGGTTAAACTc                                                                                                                                                                                                                                         >  2:234080/1‑249 (MQ=255)             ccATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGGCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCt                                                                                                                                                                                                                                     <  1:338859/251‑1 (MQ=255)              catgcTGACCGCGACCTGCGGTGATACCGGAGCGGCAGTGGCGCATGTTTTCTACGGTTTACCTAATGTGAAAGGGATTATCGGCTATCCACGAGGCTAAATCAGTCCTCTGCTAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCTATCGCCGGCAATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTGTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTATACTCAGCTa                                                                                                                                                                                                                                    <  2:244474/247‑1 (MQ=255)              cATTCTGACCGCTACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCt                                                                                                                                                                                                                                     >  1:229784/1‑250 (MQ=255)                 tCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACt                                                                                                                                                                                                                                 >  2:314494/1‑251 (MQ=255)                    gACCGCGACTGCCGGTGATTCCGGCGCGGCAGTGGCTCATGTTGTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAACCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAATTGGCGCTGGGGTTAAACTCAGCTAACTCg                                                                                                                                                                                                                               <  1:190854/250‑1 (MQ=255)                       cgcgACCTCCGGTGGTACCGGTGCGGCAGTGGCTCAGGCTTTCTACGTGTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAAACAGGCCACTGCAATAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGGTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAACCTCAGCTATCTCGATTa                                                                                                                                                                                                                           <  2:16663/251‑1 (MQ=255)                            ccTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATc                                                                                                                                                                                                                      <  1:49830/251‑1 (MQ=255)                                    gATACCTGAGCGGCAGTGGCTCATGCTTTCTACGTTTTAACGAATGAGAAAGGGGTTATCCTCTATCCACGAGGCAAAATCAGTCCAATGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGCTTAACATCAGCCGttt                                                                                                                                                                                                              <  1:229863/251‑1 (MQ=255)                                       aCCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGc                                                                                                                                                                                                            <  1:190582/250‑1 (MQ=255)                                        ccGGAGCGGCAGTGGTTCATGCTTTCTACTGTTTACCGAATGTGAAAGGGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATAGTGAGACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTg                                                                                                                                                                                                          <  1:437941/251‑1 (MQ=255)                                              cggcatcggcgtaggCTTTCTACGGGTTTGCGAATGTGAAAGTGGTTATCTCCTAGCCACGAGGCAAAAGCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTTTTGCCATCGAAGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCTTTTGAGGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAg                                                                                                                                                                                                    <  2:261279/237‑1 (MQ=255)                                                 cAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGCGGAAAAATCAGTCCAATGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGAtt                                                                                                                                                                                                 >  1:311021/1‑251 (MQ=255)                                                          aTGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCtactac                                                                                                                                                                                        <  2:272945/251‑1 (MQ=255)                                                               ttCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGGAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTg                                                                                                                                                                                    >  1:263927/1‑250 (MQ=255)                                                                  tACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGc                                                                                                                                                                                >  1:271782/1‑251 (MQ=255)                                                                     ggTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCAACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGcagt                                                                                                                                                                             >  1:291551/1‑248 (MQ=255)                                                                      gTTTACCGAATGTGAAAGTGTTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                            <  1:372504/251‑1 (MQ=255)                                                                      gTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGACGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCCGATTTGCTACTACTTTGCAGCAGtt                                                                                                                                                                            >  2:82443/1‑251 (MQ=255)                                                                      gTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                            >  1:165690/1‑251 (MQ=255)                                                                         tACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTgcg                                                                                                                                                                         >  2:416985/1‑251 (MQ=255)                                                                         tACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAAGCCTGTCAGGCGCTGGTGAAGAAGGAGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTAAAACTCAGCGAACTCGATTAAAATCAGCCGTTTGCTGGCGCAGATTTGATACTACTTTGAAGCAGTTgcg                                                                                                                                                                         >  2:385582/1‑251 (MQ=255)                                                                             gAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGc                                                                                                                                                                     >  1:87499/1‑251 (MQ=255)                                                                                tgtgAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGGTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGc                                                                                                                                                                  >  1:25999/1‑251 (MQ=255)                                                                                               ccTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGgaagcgcgac                                                                                                                                                   >  2:315981/1‑243 (MQ=255)                                                                                                    aTCCACGAAGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAg                                                                                                                                              <  1:287994/251‑1 (MQ=255)                                                                                                        aCGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTgg                                                                                                                                          >  2:478783/1‑251 (MQ=255)                                                                                                              cAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtct                                                                                                                                    <  1:314494/251‑1 (MQ=255)                                                                                                                             gCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTGGCCATCGACGGCGATTTCGAGGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                     <  2:441470/251‑1 (MQ=255)                                                                                                                                              gTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGAc                                                                                                     <  1:378608/250‑1 (MQ=255)                                                                                                                                               tACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGAcgg                                                                                                   >  2:40017/1‑251 (MQ=255)                                                                                                                                                       gcgacAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGGCTGGCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGTCGTTTGCTGGCGCAGATTTGCTACTCCTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGAAACCTTCGGCGATTTGAcggcggtgtgt                                                                                           >  2:110777/5‑246 (MQ=255)                                                                                                                                                        cggcAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGctgc                                                                                          >  1:175353/1‑251 (MQ=255)                                                                                                                                                            aaTATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCTGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGc                                                                                      <  1:174879/251‑1 (MQ=255)                                                                                                                                                               atCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGaa                                                                                   <  1:315981/251‑1 (MQ=255)                                                                                                                                                                    aaCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCggggcgtctgggggagggac                                                                              >  2:424058/1‑234 (MQ=255)                                                                                                                                                                        ggtGCCAGCGACGGCGATTTCGATGCCTGGCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCGGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCgg                                                                          <  2:263927/249‑1 (MQ=255)                                                                                                                                                                         ttGCAATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTAGATGATGAAGAGCTGAAAGTGACGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCCGATTTGCTACTACTTTGAAGCCGTTGCGCAGCTGCCTCAGGAAGCGCGCAACCAGCTGGTAGTCTCGGTGCCAGGCAGACACTTCAGCGCTTTGCCGGCgggattatcgaccaagacaaccgcc                                                                         >  2:297593/1‑229 (MQ=255)                                                                                                                                                                                 gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCAccccggctgcccgg                                                                 >  1:394372/1‑238 (MQ=255)                                                                                                                                                                                  aCGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCGAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGCAGGCACTCGGTCTGCC‑gggg                                                               >  1:125055/1‑249 (MQ=255)                                                                                                                                                                                       gATTACGATGCCTGTCAGGCGCTGGTGAAGCAGGGGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACGACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACGCGGTCTGCC‑GGTGAAACg                                                          <  2:67633/251‑1 (MQ=255)                                                                                                                                                                                                      cAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGGTGGGGGAGTCACCCGGTCTGCC‑gggggaacggtttttttgttcgac                                           >  2:151445/1‑229 (MQ=255)                                                                                                                                                                                                        ggctcTGGTGAAGCAGGCGGTTGATGATGAAGAGCGGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTGGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCC‑GGTGAAACGTTTTATTGCTGCGACCa                                         <  2:415541/247‑1 (MQ=255)                                                                                                                                                                                                               gggaagaaGGCTTTTGATGATGAAGAGCTGAAAGTGGCGCGGGGGTTAAAATCAGCTAACTCTATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCC‑GGTGAAACGTTTTATTGCTGCGACCAACGTGaa                                  <  2:121494/249‑1 (MQ=255)                                                                                                                                                                                                                 gAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCC‑GGTGAAACGTTTTATTGCTGCGACCAACGTGAAc                                 <  1:205177/250‑1 (MQ=255)                                                                                                                                                                                                                           gatgattaataagagCTGAAAGTGGGGCTGGGTTTAAACTCAGCTAACTACATTAATACAAGCCTTTTGATGGAGCAGTTTTGCATCTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCACCCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCC‑GGTGAAACTTATTAGTGCTGCGTCCAACGTGAACGATACCGTGcc                      <  2:251499/241‑1 (MQ=255)                                                                                                                                                                                                                            ttGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCC‑GGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGcc                      <  1:393458/250‑1 (MQ=255)                                                                                                                                                                                                                                   tgaAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGGGGAGAGCCACGGTCCGCCCGGGGAAA‑GTTTTATTGCTGCGACCAAccggcacaaacacgcgccacgggtct              >  1:200525/1‑226 (MQ=255)                                                                                                                                                                                                                                           tGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGGCACTCGGGCTGCC‑GGTGAAACGTTTTATTGCTGCGGCCAACGAGAACGATACCGAGACACGGTTCCTGCacgac      >  1:386266/1‑251 (MQ=255)                                                                                                                                                                                                                                               aGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGGGAAGGAAACCGGGCTGGC‑GGCGAAACCTTTTATTGTTGCGACCAACGTGGACGATAACGCGCACCGGTTCCTGCACGGCGGTc  >  1:198318/1‑251 (MQ=255)                                                                                                                                                                                                                                               aGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGGAGGCACTCGGGCTGCC‑GGTGAAACGTTTTATTTCTGCGCCCCACGTGAACGAAACCGTGCCACGTTTTCTGCCCGACgggc  >  1:112693/1‑249 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   ATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCC‑GGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTC  >  NC_000913/4104‑4592

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