Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,277,448	G→A	17.2%	E38K (GAA→AAA)	yhaV →	toxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,277,448	0	G	A	17.2%	144.4 / 10.8	73	E38K (GAA→AAA)	yhaV	toxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/11);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 5.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACA  >  NC_000913/3277215‑3277682                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               tttCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAAt                                                                                                                                                                                                                                           >  1:240794/1‑239 (MQ=255)              atatCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑GTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAC‑CCTACACCGCTTTAGTTGCCGAAGGCGAGCCTTTAAAGGGAAAAGATCCTGAAAATAAACAGaaaaaa                                                                                                                                                                                                                               >  2:516769/1‑238 (MQ=255)                 tCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGGGAAAAgcc                                                                                                                                                                                                                            >  1:429919/1‑238 (MQ=255)                 tCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAATCGAAACATTAAAGGGAAACAATCCTAAAAATAAT‑AAAGaaaaaacg                                                                                                                                                                                                                           >  2:514417/1‑234 (MQ=255)                   cAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAACGGG‑GTAATGGTTGGGCG‑CTATATG‑CCCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCGTGAAAATTATCAGAGAAAAgctgac                                                                                                                                                                                                                         >  1:799868/1‑235 (MQ=255)                           cccGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGCGCCATTAAAGGGAAAAGAACCTGAAAATTATCAGAGAAAACCCG‑CACaaaagta                                                                                                                                                                                                                >  2:423071/1‑235 (MQ=255)                           cccGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAggt                                                                                                                                                                                                                 >  1:685468/1‑237 (MQ=255)                           cccGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACACTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGtt                                                                                                                                                                                                                 >  1:733799/1‑239 (MQ=255)                           cccGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGCCATGGACGTCAACATTGATGATGTGATTGGCGACGACGACTAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCACTGTTTTCAGGAAA‑CCGACGACGCTTTACTTGCCGAAGTCAAGACGTTAAAGGGAAAGAATCATGAAAATTATCAAAGAAAAGCCGCCGCAAAttt                                                                                                                                                                                                                 >  2:58098/1‑239 (MQ=255)                               cAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CTTACGAAGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAACAGAACCTGAAAATTACCAGAGAAAAGCCGCCCAAAAGAAATTg                                                                                                                                                                                                             >  2:258093/1‑239 (MQ=255)                                aaaaacCTCGTCCATTAAACATTCAACAAGAAAAGAAACTTGCCGCTGGCAGGGACGTCAACATTGATGATGAGTTTGGCGCCGACGAATAATGAATTTTCCAAAAAGGT‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTgg                                                                                                                                                                                                            <  1:514417/239‑1 (MQ=255)                                aaaaaaCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CTTACGACGCTTTAGTTGCCGAAGTCGAGCACTTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCCCCACAAAGAAATTgg                                                                                                                                                                                                            >  2:171502/1‑239 (MQ=255)                                  aaatCTCCTCTATTAAACATGCAAAAAGGAAAGAACATTGTCGCTGGAATGGACGTAAATATGGATGATGAGGTTGGATACGACGAATAATGGATTTTCCACAAAGGT‑TTAAGGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGACACATTAAAGGGAAAAGATCCGGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCg                                                                                                                                                                                                          <  2:505708/239‑1 (MQ=255)                                  aaagcgtgTCCTTTCAACATTAAACAAGGAAAGAAGTGTGTCGGTGGCAGGGACGAACACATATATGATGAGATTGGGGGCGACTAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTTTTTTAAGGAAA‑CCTACGACGTTTTAGTTGCCGATGTCGAGACATTAAAGGAAAAAGATCCGGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCg                                                                                                                                                                                                          <  1:558725/232‑1 (MQ=255)                                    aaCTCGTCCATTAAACATTCAACAAGGAAAAAACCTTGGCGCTGGCATGGACGTCAACAGTGATGATGAGATTGGCGACGAAGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGt                                                                                                                                                                                                        <  1:1742/239‑1 (MQ=255)                                    aaCTCGTCCATTAAACATTAAACAAGGAAAGAAACTTGTCGCTGGCAGGGACGTCAAAATTGAGGAGGAGATTGGCGACGACAAATAATGGATTTTCAAAAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATAAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCTGt                                                                                                                                                                                                        <  1:171502/239‑1 (MQ=255)                                        cGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCCCAAAGTTATTGGCGGTA‑GTc                                                                                                                                                                                                   >  1:27557/1‑239 (MQ=255)                                          ctctttaaacattaaaaaGGGAAAAAAACGTTCCCGTGGAGTGGACTTCAAATTTGAGGGTGGGTTTGGGGACGAGAAATAAGGGATTTTCAAAAAAGGG‑TTAATGGTTGGGCG‑TTATAGG‑CTCATCCCTGTTTAGAGGAAA‑CCTACGACGCTTTAGTCGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCa                                                                                                                                                                                                 <  2:27557/232‑1 (MQ=255)                                                     tccagaaaGGAAAAAAGTGTGCCGCTGGCATGGACGTCAACATTGATGATGAAATGGGCGACGACGAATAATGGATTTTCAAAAAAGGG‑TTAAGGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTg                                                                                                                                                                                     <  1:423071/233‑1 (MQ=255)                                                     tacaacaaGGAAAGAAACTTGTCGCTGGCATGGCCTCAAACATTGATGATGAGATTGGCAACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATAGC‑CTCATCCCTGTTTTCAGGAAC‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTg                                                                                                                                                                                     <  1:725782/237‑1 (MQ=255)                                                      tcaacaAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑AGCCA‑TAAAATGATTga                                                                                                                                                                                    >  1:777067/1‑239 (MQ=255)                                                      tcaacaAGGAAAGAAACTTGTCGCTGGCATGGACGTAAACATTGATGATGGATTTGGCGACGAAGAATAATGGATTTCCCAAAAAGGG‑TTAAGGGTTGGGCG‑TTAATTG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTga                                                                                                                                                                                    <  1:482008/239‑1 (MQ=255)                                                       caaaagggAAAAAATCTTGTCGTGTGCATTGCAGGCAAAATGGAGGAGGAGATTGGCGAAGACGAATAATGGTTTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGAGGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTgag                                                                                                                                                                                   <  1:1004408/233‑1 (MQ=255)                                                         acaaGGAAAGAAACTGGTCGCTGGTGTGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTCCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGAAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTgagga                                                                                                                                                                                 <  1:454325/239‑1 (MQ=255)                                                           aaGGAAAAAAACTTGTCGGGGGCATGGACGCAAAAATTGAGGAGGAGGTGGGGCAAGAGAAATAATGGATTTTCCAAAAAGGG‑TTAATGGTTGGGTG‑CTTTATG‑CTCTCCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCATGAAAATTATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGc                                                                                                                                                                               <  2:240794/239‑1 (MQ=255)                                                             ggAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCC‑TAAAATGATTTAGGGGcct                                                                                                                                                                             >  1:67265/1‑237 (MQ=255)                                                             ggAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTACA‑TAAAATGATTTAGGGGCat                                                                                                                                                                             >  1:752415/1‑239 (MQ=255)                                                                     atttttccttGGCATTGAGGTCAAATTTGATGTTGGGTTGGGCGCCAAAAAATAAGGGTTTTTCCACAAAGGG‑TAAAGGGTTGGGCC‑TTATATG‑CTCATCCCTGTTTGCAGAAAA‑CCTACGACGCTTTAGTTGCCAAAGTCGAGACATTAAAGGAAAAAGATCCAGAAAATTATAAGAAAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGt                                                                                                                                                                     <  2:545023/230‑1 (MQ=255)                                                                        tGTCGCTGGCAAGGACTAAAACATTGATGTTGAGATTGGCGCCAGCAAATAATGGATTTTCCACAAGGGG‑TTAATGGTTGGGCG‑TTATATG‑CTCCTCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCaa                                                                                                                                                                  <  1:115085/239‑1 (MQ=255)                                                                                         taaaCATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCAAAAGGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑Taa                                                                                                                                                                                             <  1:829879/193‑1 (MQ=255)                                                                                          cAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCCGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GGCCA‑TAAAGAGAATTAGGGGCCTAACACGGGCAAACCATCATCACCgccct                                                                                                                                                 >  1:1006119/1‑234 (MQ=255)                                                                                          cAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAAGGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCAAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCAt                                                                                                                                                 <  1:258093/238‑1 (MQ=255)                                                                                                  atgatgaGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAg                                                                                                                                                                                                                   >  1:524401/1‑166 (MQ=255)                                                                                                      tgggaTTGGCGACGACGAATAATGGATTTCCCAAAAGGGG‑TTAAGGGTTGGGCG‑CTATATG‑CCCATCCCTGTTTCCAGGAAA‑CCTACGACGTTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGACCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCa                                                                                                                                    <  1:83978/235‑1 (MQ=255)                                                                                                       gaggtgggcggcGCCAAATAATGGATTGTCCACAAGGGT‑TAAAGGGTGGGGGG‑TTAAAGGACTC‑TCCTGGTTTTCGTGAAT‑CTTCCGACTTTTTATTGGCGTAAGTCGAGACATTAAAGGGAAATGATCCAGAAAATTATCAGGAAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCaa                                                                                                                                   <  2:752415/232‑1 (MQ=255)                                                                                                            tGGCGAAGAAGAATAAGTGATTCTCCAAAAAGTG‑TAAATGTTGGGGGT‑TTATATT‑CTCCTCCAGTTTTTCAGAAAT‑CCTACGCCGGTTGAGGTGCCAAAGCCGAGAATTTAAAGGGAAAAGTTCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGtt                                                                                                                             <  1:58098/239‑1 (MQ=255)                                                                                                                       aataatGTATTTACCAAAAAGGG‑TTAAGGGTTGGGCG‑TTATATG‑CCCACCCCTGTTTTCAGAAAC‑CCGACGACGTTTTAGTTGCCGAAGTCGAAACATTAAAGGGAAAAGATCCTAAAAATTACCGGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGa                                                                                                                  <  1:449470/239‑1 (MQ=255)                                                                                                                       aataatGGATTTTCCAAAAAGGG‑TTAATGGTTGGGCG‑TTATATG‑CTCATCCTTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGt                                                                                                                                                                     <  1:740106/189‑1 (MQ=255)                                                                                                                        atatgggATTTCCCACAAGGGG‑TTAAGGGTTGGGCG‑CTATATG‑CGCATCCCTGTTTCAAGGAAA‑CCTACGACGTTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGa                                                                                                                  <  1:962341/233‑1 (MQ=255)                                                                                                                        ataatGGATTTCCCACAAAGGG‑TTAATGGTTGGGCC‑GTATGGG‑CTGCTCCCTGTTTTCAGGAAA‑CCTACGACGTTTTAGTTGCCGAAGTCAAGAAATAAAAGGGAAAAGACCCTAAAAATGATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGa                                                                                                                  <  1:263335/238‑1 (MQ=255)                                                                                                                             tgtttttTCCCAAAGGGTTTAATGG‑TGGGCGTCATTATG‑CACTCCCCTGTTCTCAGAAAA‑CCTACCACGGTTTAGTGCCCGAAGACGAGACATAAAAGGGAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAAt                                                                                                             <  2:667045/236‑1 (MQ=255)                                                                                                                                        aaaGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCC‑ATAAGAGATTTAGGAGACAATACCGGGCAATCCATACTCACCCGCATTCCGTGCT‑GGACAGTAGGTAAACTCTGGGAaaaaaaaaaaagtgacc                                                                                                 >  2:378232/1‑230 (MQ=255)                                                                                                                                        aaaGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGCATGTCCA‑‑AAAAAGATTTAGGGGGAAATCACGGTCAATACATCATTAACGGCATTCCGTCCTGGGCAAGT‑GGTAGGGACTGGGCaaaaaaaaaaacgagga                                                                                                 >  2:979477/1‑230 (MQ=255)                                                                                                                                            gg‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTgag                                                                                                                                                                                   <  1:389250/154‑1 (MQ=255)                                                                                                                                               ttAAGGGTTGGGGG‑CTATATG‑CTAATCCCTGTTTTAAGGAAA‑CCTTCGACGCTTTAGTTGCCGAAGTCGAGAAATTAAAGGAAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTa                                                                                          <  1:516769/239‑1 (MQ=255)                                                                                                                                                tAATGGTTGGGAG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTGTATTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTAGCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTa                                                                                          <  1:163958/238‑1 (MQ=255)                                                                                                                                                tAAGGGTGGGGCG‑CTAATTG‑CCCATCCCTGTTTCAAGGAAC‑CCTACGACGCTTTAGTTGCCGAAGTCGAGAATTTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTa                                                                                          <  1:949228/238‑1 (MQ=255)                                                                                                                                                aaagggTTGGGCG‑CTATAGT‑CTCATCTCTGTTTGCAGGAAA‑CCTACGACGCTTTAGGTGCCGAAGTCGAGACTGAAAAGGGAAAAGATCCTGAAAATTTGCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTa                                                                                          <  1:385797/234‑1 (MQ=255)                                                                                                                                                       tGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTG‑GTCCACAAAAGTG‑TTGAGGGGACAATACCGGGACATCCACCCTCCCCCGCATTCCCTAAG‑TGGAAGACGTTAGGCTCTTGGAAAAAAAAAAAC‑AGGACACGGGAaaaaaatgt                                                                                  >  2:167854/1‑234 (MQ=255)                                                                                                                                                       tGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAATGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTTCGTCAT‑GGCAAGTCGTTTGGCTCTTGGGAAAATCAAGAC‑TGGTCACGGGGCAAAATTgg                                                                                  >  1:803006/1‑239 (MQ=255)                                                                                                                                                           gggtttgatG‑CCCACCCCTGTTTT‑AGGAAACCCTACGACGTTTTAGTTGCCGAAGACGAAAAATTAAAGGAAAAAGATCCTGAAAATTAAAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTg                                                                                <  1:351388/231‑1 (MQ=255)                                                                                                                                                           gccttttttG‑CTCATTCTTTTTTTCAGAAAA‑CCTGCGATGTTTTAGTTGCAGAAAGCGAGACATTAAAGGGAAAAGATCCGGAAAATTATCAGAGAAAGGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTg                                                                                <  2:429919/230‑1 (MQ=255)                                                                                                                                                           ccgttataTT‑CTCATCCCTGTTTACAGGAAA‑CCTACGACGTTTTAGTTGCCGAAGTCGGGATATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAGGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTg                                                                                <  1:198836/234‑1 (MQ=255)                                                                                                                                                                 tatG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGTTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCAcc                                                                                                                                                      <  1:391027/164‑1 (MQ=255)                                                                                                                                                                          tCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATc                                                                                                                                                                <  2:688930/146‑1 (MQ=255)                                                                                                                                                                          tCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAATGATTGAGGAGCATCTCCCGGTCAATc                                                                                                                                                                >  1:688930/1‑146 (MQ=255)                                                                                                                                                                          tCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAAGGATTGAGGGGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAAAA‑TGGTCACGGGGAAAAATTGGGGCTGGTGGTTATCg                                                                   >  1:370340/1‑237 (MQ=255)                                                                                                                                                                            ccTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGGGGCTATCACGGTCAATCCATCATCACCGGCCTTCCGTCAT‑GGCAAGTCGTTAGGGCCTGGGAAAAATAAAAACGGGGCCAGGGACAAAATTTGG‑GCTGGCGGTTTTCGtc                                                                 >  1:548174/1‑237 (MQ=255)                                                                                                                                                                              tGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAAGGGTTTTGGAGCATCACCCGGGTAATCCCTCCTCAACGGCATTCCGTCCT‑GGGAAGTCGTTAAGCGTTGGGGGAAATAAAGAC‑AGGACACCGGGAAAATATTGTGTTGTTCggattagtccct                                                              >  2:565562/1‑227 (MQ=255)                                                                                                                                                                                     aGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTAt                                                                                                                                                                                                                                        >  2:682489/1‑65 (MQ=255)                                                                                                                                                                                     aGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTAt                                                                                                                                                                                                                                        <  1:682489/65‑1 (MQ=255)                                                                                                                                                                                     aGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GGCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGGTCTGGGGAAAATAAAAAC‑TTGGCACGGGGAAAAATTTGTGCTTGGCGTTTTCGGTCCTTCTTTCGt                                                      >  1:140652/1‑239 (MQ=255)                                                                                                                                                                                                cgccgCTTGAGTTGCCGAAGTCGAGACATAAAAGGGAAAAGCTCCTAAAAATTACCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGaaa                                            <  1:461046/236‑1 (MQ=255)                                                                                                                                                                                                      tttAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAAGGGTTGTGGAGGATATAACGGTCAATCCATCCTCAACGGCATTCCGTACA‑GGCCGGTCGTTAGGCTCCTGGGAAAATAAAGAC‑TGGTGCCGGGGAAAATTTTGTGCGGGGTGGTATGGTCTCTTCTTTCTTTAAAAAGaaaaaaaaa                                      >  2:11759/1‑235 (MQ=255)                                                                                                                                                                                                      tttAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GGCAA‑TAAAATGATTTTGGGGGAAATAACGGTCAATCCACCATCACCGGCATTCCGTCAT‑GGCAAGTACTTAGGCTCCGGGGAAAAAAAAAAA‑TGGTCACGGGGCAAAATTGTGGCTGGTCGTTTTATGCTCTTTTTTTGTTTTAAGGaaaaaaaaa                                      >  2:696866/1‑235 (MQ=255)                                                                                                                                                                                                            taataagagacagaCATAAAAGGAAAAGGATCAAAAAAATTAACAAAAAAAAGCCGACACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATAGCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtca                                <  2:170519/228‑1 (MQ=255)                                                                                                                                                                                                            tGCCGAAGTAGAGAAATTAAAGGGAAAAGATCCTAAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtca                                <  1:754814/239‑1 (MQ=255)                                                                                                                                                                                                             cccGAAGGCGAGAAATTAAAGGGAAAAGTCCCGAAAAATTATGAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtcat                               <  1:979477/238‑1 (MQ=255)                                                                                                                                                                                                                       ctaCAATAAGGGAAAAAGAGCTTGAAAATTAACAGAGAAAAGCCGGCCCAAAGTTATTGGCGGTA‑GCCCA‑AAAAGTTATTTAGGAACATATCACGGTCAAGCCATCATCACCGGCAGTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGgga                     <  2:491648/237‑1 (MQ=255)                                                                                                                                                                                                                        aaCATTAAAGGGAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGgga                     <  1:11759/237‑1 (MQ=255)                                                                                                                                                                                                                                aGGGAAAAGATCCCTAAAATTAACAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAAc            <  1:106103/239‑1 (MQ=255)                                                                                                                                                                                                                                    aaaaGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGAGATTGAGGAGCATATCACGGTCAATCCATCATCCCCCGCATTCCGGCCT‑GGCAAGTCGTTAGGCTCTGGGGAAAATAAAGAC‑TTGTCACCGGGAAAATTTGGTGCTGGTCGTTATCGCCTCTTCTTTCGTTTTAAGGAAAAAAAGAAAGtcatccatttggg                      >  1:170519/1‑217 (MQ=255)                                                                                                                                                                                                                                          tCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAaca  <  1:696866/239‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               TTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGG‑TTAATGGTTGGGCG‑CTATATG‑CTCATCCCTGTTTTCAGGAAA‑CCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTA‑GTCCA‑TAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCAT‑GGCAAGTCGTTAGGCTCTGGGAAAAATAAAGAC‑TGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACA  >  NC_000913/3277215‑3277682

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