Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,006	G→A	S163S (AGC→AGT)	satP ←	acetate/succinate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	10,006	0	G	A	97.5%	113.7 / ‑4.2	41	S163S (AGC→AGT)	satP	acetate/succinate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/22);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

ACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGC  >  NC_000913/9760‑10246                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  aCGTTGTGGTACACGGTATGTTTGCCAGCGTACCGTACTGCATTCAGTTGGTGGAAGGGCCATCCGGTGAAACGGCACCAGAAGTGGTTGCTTCATTCAGACCGAAGAGTGCAAGACGCGACGTGAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAACCAATTAGTGGGA‑TTCACCAATCGGCCGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACTCATCGCCTGATAGATTGCACTgg                                                                                                                                                                                                                                               <  1:470116/251‑1 (MQ=255)           ccactgtagttttGCCAGCGTACCGTACTGCATTCAGTTGCTGAAATGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCTAGCGGGGGTATCTCAAAACAATTAGTGGGA‑TTCACCAATCGTCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGa                                                                                                                                                                                                                                      <  1:38289/241‑1 (MQ=255)                                 ccgtgctgcagtccattgattgAAGGGCCATACGTTGAAACTGCACCATAAGTGGTTTCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAATCAATTAGTGGGA‑TTCACCACTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCCCCCCAGATCAGCCCAATCCAGCCGGCAAAg                                                                                                                                                                                                               <  2:156575/230‑1 (MQ=255)                                       tGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGACGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGCTCAGCCCAATCCATCCGGCAACGTGgatg                                                                                                                                                                                                         >  1:38041/1‑251 (MQ=255)                                                  gCTGGAAGGGCCATACGTTGAAACTGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCAGCGCCAGATAGATTGCACTGGCACCACAGAGCAGACCAATCCAGCCGGCAAAGTGTATGATTGCGGCGt                                                                                                                                                                                               <  1:422348/250‑1 (MQ=255)                                                             cATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGt                                                                                                                                                                                   >  1:313888/1‑251 (MQ=255)                                                              aTACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGtt                                                                                                                                                                                  <  1:191727/251‑1 (MQ=255)                                                                         cGGCACCAGAAGTGGTTGCTGCATGCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACATATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCca                                                                                                                                                                       <  2:291351/251‑1 (MQ=255)                                                                          ggCACCGGAAGTGGTTGCTGCATTCAGACCGCAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACAACTCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGcc                                                                                                                                                                        >  2:214286/1‑249 (MQ=255)                                                                             accagAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGTATAAAACAATCCCCCCGAGCGGGGGGATCTAAAAACAATTAGTGGGA‑TTCACAAATCGCCAGGACGTTGCGACCAAACTGCTCGTTCAGGACTTCACCCATCTCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGCCAAAGTGGATGATTGCGGCGTACCCGGCACTGTTACCGATCGCcagc                                                                                                                                                                     >  1:38951/5‑249 (MQ=255)                                                                                      ggTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAAc                                                                                                                                                            >  2:423885/1‑249 (MQ=255)                                                                                          gctgcaTTCAGCCCGAAGGGTTCAAGACTCGTCGTTAGCGAATAATAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAGTAGTGGGACGTC‑CCAATCGGCAGAACGGGGCTACCAAACTGCTCGTTCAGGACTTCGCCCATCGCCAGATAGATGTCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGGGGATGATTGCGGCGTTACCGGCAATGTTCCCGATCGCCAGCAGGGCAAACAGCAc                                                                                                                                                       <  2:259614/247‑1 (MQ=255)                                                                                                         aaGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGaa                                                                                                                                        >  1:39659/1‑250 (MQ=255)                                                                                                          agagTGCAAGACGGGACGGTAGCGAATAAAAAACCCCCCCCGAGCGGGGGGATCACAAAACCATTAGTTGGA‑TTCCCCAATCGGCAGAACGGTGCGACCAAACTGCTCTTTCAGTACTTCTCCCAGCGCCAGATAGAGTGCACTGGCCCCACAGCTCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAATaaaa                                                                                                                                      <  2:377805/251‑1 (MQ=255)                                                                                                            agTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGAAAACAGCACAGTCAGGCTAAAGAAAACg                                                                                                                                    >  2:359168/1‑251 (MQ=255)                                                                                                            agTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGAGGGATCTCAACACAATTAGTGGGT‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCCGACCAATCCAGCCGGCAAAGTGGATGATTACGGCGTTACCGGCCATGTTACCGATCGCCAGCCGTGCAAACAGCCCAGTCAGGCTAAAGAAAAc                                                                                                                                     >  1:129139/1‑250 (MQ=255)                                                                                                                 aaGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTg                                                                                                                               <  1:342872/251‑1 (MQ=255)                                                                                                                          cGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAAGTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAAcgcg                                                                                                                      <  1:158091/251‑1 (MQ=255)                                                                                                                            ttAGCGAGTAAAAAATTCCCCCCGAGCTGGGGGATCGCAAAACAATTTGTGGGC‑TTCACCAACCGGAAGAACGGTCCGACCAAACTGCACGTTCAGTACTCCACCCATCGCCAGATAGATTGCACTGGCACCACAGAGCAGACCAATCCAGCCGGCAAAGTGGATGAGTGCGGCGTTACCGGCAATGTTACCGAGCGCCAGCAGGTCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAAcgcgg                                                                                                                     <  2:342341/250‑2 (MQ=255)                                                                                                                                          aaTCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTg                                                                                                     <  1:121847/251‑1 (MQ=255)                                                                                                                                               ccccGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCgaa                                                                                                <  1:359168/251‑1 (MQ=255)                                                                                                                                                       ggggggATCTAAAAACAATGAGTGGGA‑GTCACTAAGCGGCAGAACGGTGCGATCAAACGGATAGGTCAGTACTGCACCCATCGCCAGATAGAATGCACGGGCTCCACAGATGAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTGACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAAGTGCAGAACGCGTGCGGCTTTCAG‑AGTGCCGAAGAACATaa                                                                                        <  2:38951/251‑1 (MQ=255)                                                                                                                                                        gggggATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATaa                                                                                        >  1:415084/1‑250 (MQ=255)                                                                                                                                                        gggggATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATaa                                                                                        >  1:198434/1‑250 (MQ=255)                                                                                                                                                                     aCAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAg                                                                                                                                                 <  2:404201/181‑1 (MQ=255)                                                                                                                                                                     aCAATTAGTGGGA‑TTCACCAATCGGAAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAg                                                                                                                                                 >  1:404201/1‑181 (MQ=255)                                                                                                                                                                       aaTTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGcc                                                                        <  2:362781/251‑1 (MQ=255)                                                                                                                                                                            gTGGGA‑TTCACCATTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTCCCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCGAAAGAAAACGAATTGCAGAACGCGAGCGCCTTTCAG‑CGTGCCGAAGAACATACACAGCGTAAATACGCCCcac                                                                    >  2:433393/1‑250 (MQ=255)                                                                                                                                                                               ggA‑TTCACCAAGCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAAGCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGAc                                                                <  2:313888/251‑1 (MQ=255)                                                                                                                                                                                         aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGt                                                          >  2:53009/1‑248 (MQ=255)                                                                                                                                                                                         aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGa                                                       <  1:53009/251‑1 (MQ=255)                                                                                                                                                                                            cGGCAGAACGGTGCGCCCAAACTGCTCGTTCAGTACTTCACCCATCTCCAGATAGATTGCCCTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAACCAGCACAGTCAGGCTAGAGAACACACATTGCAGCACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCCGGTAGaccc                                                    >  2:260812/1‑249 (MQ=255)                                                                                                                                                                                                     ggTGCGACAAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACt                                           <  1:186499/251‑1 (MQ=255)                                                                                                                                                                                                               aaCTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTgg                                 <  1:246031/251‑1 (MQ=255)                                                                                                                                                                                                                    cTCGTTCAGTACTTCACGCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGACAACGAATTGCAGAACGCGTGCGCCTTT‑AGTCGTGCCGAAGTACATAAACAGCGTAACTCCGCCACACAGACCCCGGTAGCCACCAAGGAAcggggcctttagcgcac                            >  2:94721/1‑235 (MQ=255)                                                                                                                                                                                                                         tCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGTGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGt                        <  2:391422/250‑1 (MQ=255)                                                                                                                                                                                                                         tCAGTACTTCACCCATCGCCAGATAGATTGCACTGGCACCACAGATAAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTc                       >  1:196256/1‑251 (MQ=255)                                                                                                                                                                                                                                          gCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCAt      >  1:23020/1‑251 (MQ=255)                                                                                                                                                                                                                                          gCCAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCAt      <  2:415084/251‑1 (MQ=255)                                                                                                                                                                                                                                            cAGATAGATTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATca    <  1:302250/251‑1 (MQ=255)                                                                                                                                                                                                                                              gatagatTGCACTGGCACCACAGATCAGACCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACAGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATcagc  >  1:280108/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  ACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGA‑TTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCGCTGGCACCGCAGATCAGCCCAATCCAGCCGGCAAAGTGGATGATTGCGGCGTTACCGGCAATGTTACCGATCGCCAGCAGGGCAAACAGCACGGTCAGGCTAAAGAAAACGAATTGCAGAACGCGTGCGCCTTTCAG‑CGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTCGGCATCAGC  >  NC_000913/9760‑10246

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