Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,336	A→G	A53A (GCT→GCC)	satP ←	acetate/succinate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	10,336	0	A	G	98.1%	153.7 / ‑4.7	53	A53A (GCT→GCC)	satP	acetate/succinate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (21/31);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

AGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAGCGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCAGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAAT  >  NC_000913/10094‑10574                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               aGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAAt                                                                                                                                                                                                                                        >  1:264176/1‑251 (MQ=255)          aGGCTAAAGAAAACGAAATGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCTGCAGACCGGCAAAAATTTGCGCGAt                                                                                                                                                                                                                               <  2:565609/251‑1 (MQ=255)                aaGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATgccgcc                                                                                                                                                                                                                         >  1:232568/1‑251 (MQ=255)                aaGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGGACTCCAGCAGACCGGCAAAAATTTGCGCGATgccgcc                                                                                                                                                                                                                         <  1:288299/251‑1 (MQ=255)                   aaaaCGAATTGCAGAACGCGTGCGCCGTTCAGCGTGCCAAAGAACATAAACCGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTTAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTTCAGCAGACCGGCAAAAATTTGCGCGATGCTGCCGTa                                                                                                                                                                                                                      <  2:653212/251‑1 (MQ=255)                    aaaCGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGACTGCGGTTAAACCGAAAGTGTTGCCTTTTTTGGACTCAAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTa                                                                                                                                                                                                                      >  1:202593/1‑250 (MQ=255)                         aaTTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAAt                                                                                                                                                                                                                <  1:38884/251‑1 (MQ=255)                               aGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGCTGCCGCCGTAGAAAATGCCCAt                                                                                                                                                                                                          >  1:43913/1‑251 (MQ=255)                                 aaCGCGTGCGCCTTTCAGCGTGCCGAAGAACATATACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAGCGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCATTTTTTGTACTCAAGCAGACCAGCAAAACTTTGCGCGATGCCGCCGTAGAAAATGCCCCTgg                                                                                                                                                                                                        >  2:307868/1‑251 (MQ=255)                                 aaCGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATgg                                                                                                                                                                                                        >  1:40821/1‑251 (MQ=255)                                     cgTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCa                                                                                                                                                                                                      >  2:132254/1‑249 (MQ=255)                                      gTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCCTCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGAGCATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCa                                                                                                                                                                                                      >  1:650331/1‑248 (MQ=255)                                        gcgcCATTCAGCGTGCCGAAGAACATAAACAGCGTAAATACTCCCCACAGACCCAGGTATACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCAAGTTTCGGTATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGaat                                                                                                                                                                                                 <  2:191910/251‑1 (MQ=255)                                            cTTTCAGCGTGCAGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAACAATTTGCGCGATGCCGCCGTATAAAATGCCCATGGCAAGaataat                                                                                                                                                                                              >  1:125194/1‑250 (MQ=255)                                               tCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAATTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATAcc                                                                                                                                                                                           <  1:637237/250‑1 (MQ=255)                                                  gCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCAAAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGGCAGCTAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTcc                                                                                                                                                                                       <  1:622368/251‑1 (MQ=255)                                                   cGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCa                                                                                                                                                                                      <  1:132254/251‑1 (MQ=255)                                                              acaaaaacgacGTAAATCCGCCCCACAGCCCCAGGTAGACACCAAGGAACTTCGCATTTGGCGCATCTGTCAGACCCAGTTTCGGCAACAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAAGGCGGTAAAACCGAAAGTGTTGCCTTTGTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGAGGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATaa                                                                                                                                                                           <  2:396684/241‑1 (MQ=255)                                                                        cGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCAGCGGTCAGACCCAGTTGCGGCAACAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTAAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTtgt                                                                                                                                                                 <  2:264176/251‑1 (MQ=255)                                                                                  ccccACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCCGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGGCCAGAGCGAAATAACCCCCGTTGTGCAGGTTa                                                                                                                                                         >  2:509807/1‑249 (MQ=255)                                                                                                  aGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGAATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATg                                                                                                                                       <  2:207852/251‑1 (MQ=255)                                                                                                        cAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGCCCAGAGCGAAATAACCCACGTTGTGCCGGTTAAGCAGAATGGTGGTCATGCCGaa                                                                                                                                  >  2:119463/1‑250 (MQ=255)                                                                                                                           cgcATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTg                                                                                                              <  1:638544/251‑1 (MQ=255)                                                                                                                                  gTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCACCGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATTCCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGAtt                                                                                                       <  1:623128/251‑1 (MQ=255)                                                                                                                                   ccAGACCCAGTTTCGGCAACAGCATAATGGAAACCACCGTAAGCCAGAAAGACCCGTAAGAGGTAAATGCGGTTAAACCGAAAGTGTGGCTTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATATCCCACGTTGTGCAGGTTAAGCAGAATGGTGTTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTa                                                                                                      <  2:447557/250‑1 (MQ=255)                                                                                                                                       aCCCAGTTTCGGCAACAGCAAAATCGCAACCAACGTCAGCCACAAAGAACCGGAAGAGGAGAATGGGGTTAAACCGAAAGTGTTGCCTTTGTAGTACGCCAGCAGACCGACAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGGCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGTCCAGCGGTGCCGGATTAGCCa                                                                                                  <  2:110624/251‑1 (MQ=255)                                                                                                                                               ccggtaccagcagAATCGCAACCAACGACAGCCATAAAGACCCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGGTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAgtg                                                                                          <  2:43913/244‑1 (MQ=255)                                                                                                                                                  gCATCAGTAGAATCGCAACCAGCCTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGCGTTGCCTTTTTTGTACCCCAGCAGACCGGCAAAAATTTGCGCGATGCAGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCGCGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGGGTTg                                                                                       <  2:48425/251‑1 (MQ=255)                                                                                                                                                             aTCGCAACCAACGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAATTGTTTCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTcc                                                                            <  2:179995/251‑1 (MQ=255)                                                                                                                                                                    ccaacGTCATCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATACCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGTTGGTCATGcc                                                                                                                                     <  2:307713/183‑1 (MQ=255)                                                                                                                                                                    ccaacGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGcc                                                                                                                                     >  1:307713/5‑187 (MQ=255)                                                                                                                                                                     caacGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGAAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCAAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATATTTCTTCAAAAat                                                                    >  1:344785/4‑251 (MQ=255)                                                                                                                                                                       acGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCTTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAatca                                                                  <  1:587091/250‑1 (MQ=255)                                                                                                                                                                        cGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAatcat                                                                 <  1:4021/251‑1 (MQ=255)                                                                                                                                                                               cAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGGGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGAGGTCATGCCGAAGCCCATAAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCga                                                           >  1:679206/1‑250 (MQ=255)                                                                                                                                                                                agaaagaaCCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTGGCGCGATGCCGCCGTAGAAAATGCCCAAGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATTGTGGTCATGCCGAATGCCATCAGGCCCAGCGGTGCTGGAGTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAAGCATCATCgaat                                                         <  1:30489/251‑1 (MQ=255)                                                                                                                                                                                      aaCCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCgaat                                                         >  1:559585/1‑245 (MQ=255)                                                                                                                                                                                      aaCCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATGTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTATCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCgaat                                                         <  2:559585/245‑1 (MQ=255)                                                                                                                                                                                                ggTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCAt                                                                                                                                                                                                          >  2:573171/1‑90 (MQ=255)                                                                                                                                                                                                ggTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCAt                                                                                                                                                                                                          <  1:573171/90‑1 (MQ=255)                                                                                                                                                                                                  tGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGAGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCAACAGGCCCAGCGGTGCCGGATTAGCCAAATTAGTGTTTCCCATAATTCCCCAAAAATCATCCTCGAATAAATGGTGAAAAAAtttt                                       >  2:575848/1‑250 (MQ=255)                                                                                                                                                                                                  ggAATGCGGTTAAACCGTAAGTGGTGTCTTGGTTGGACCTCAGGAGACCGGAAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATTCCGTCCAGAGCGAAATAACCCACGGTGTGCAGGTTAAGCAGAATGGTTGTCATGCCGAAGCCCATCATGCCAAGCGGTGCCTGAGTAGCCAACTTAGTGTTGCCCATAATTCCGCAAAAATCAGCATCGAATGAATGGTGAAATAATTTc                                       <  2:125194/250‑1 (MQ=255)                                                                                                                                                                                                       gCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGa                                  <  1:429671/251‑1 (MQ=255)                                                                                                                                                                                                                ccGAAAGTGTTGCCTTTTTTGTACTCCAGAAGAACGGCAAAAAGGTGCGCGAAGCCGACGTAGAAAATGCCCATGGCAAGAATAATACCGTCAAGAGCGAAATAACCCACGTTGAGAAGGTTAAGCAGCATGGAGGTCATGCCGAAGCCCATCAGGCCCAGCAGTGCCGGATTCGCAAACTTAGTGTTGCCAATAATTCCTCAAAAAACATCATCGAATGAATGGGGAACAAATTTCACTGGATaccaga                          >  2:344835/1‑245 (MQ=255)                                                                                                                                                                                                                     aGTGTTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGtt                    >  1:255410/1‑251 (MQ=255)                                                                                                                                                                                                                       tgtTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGtttt                  >  1:529585/1‑251 (MQ=255)                                                                                                                                                                                                                        gtTGCCTTTTTTGTACTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTc                 <  1:129725/251‑1 (MQ=255)                                                                                                                                                                                                                           gCCTTTTTTGTACGCAAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAgg              <  2:255410/251‑1 (MQ=255)                                                                                                                                                                                                                                      aCTCCAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAaataat                                          >  1:436292/1‑212 (MQ=255)                                                                                                                                                                                                                                      aCTCCAGCAGACCGGCAAAAACTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAaataat                                          <  2:436292/212‑1 (MQ=255)                                                                                                                                                                                                                                         ccAGCAGACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCataat  <  1:300041/249‑1 (MQ=255)                                                                                                                                                                                                                                           agcagACCGGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAAtgg                                                                                                                                               <  1:182040/106‑1 (MQ=255)                                                                                                                                                                                                                                           agcagACCGGCAAAAATTTGCGCGATGCCGCAGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACGAAAGTTGTGCAGGATAAGCAGaaggg                                                                                                                                               >  2:182040/1‑103 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               AGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGCAGAATCGCAACCAGCGTCAGCCAGAAAGAACCGTAAGAGGTGAATGCGGTTAAACCGAAAGTGTTGCCTTTTTTGTACTCCAGCAGACCAGCAAAAATTTGCGCGATGCCGCCGTAGAAAATGCCCATGGCAAGAATAATACCGTCCAGAGCGAAATAACCCACGTTGTGCAGGTTAAGCAGAATGGTGGTCATGCCGAAGCCCATCAGGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAAT  >  NC_000913/10094‑10574

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