Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	653,049	T→G	21.5%	V272G (GTG→GGG)	citA →	sensory histidine kinase in two‑component regulatory system with CitB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	653,049	0	T	G	21.5%	120.4 / 3.0	66	V272G (GTG→GGG)	citA	sensory histidine kinase in two‑component regulatory system with CitB
Reads supporting (aligned to +/- strand):  ref base T (32/17);  new base G (10/6);  total (43/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 = 1.00e+00

TCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAA  >  NC_000913/652813‑653286                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  tCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAgggg                                                                                                                                                                                                                                                 >  1:712772/1‑239 (MQ=255) tCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCCGGGGCTGAGTTCACCCAGACGGCAATGGTTGGGAAAACCAATTGTGGAAgggg                                                                                                                                                                                                                                                 >  2:655455/1‑239 (MQ=255)        ttGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGc                                                                                                                                                                                                                                          >  1:417537/1‑238 (MQ=255)        ttGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACACGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCGGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCCGAGTTCCCCCGGCCGGCAATGGTTGGGTAAACCCATTGTTGAAAGGGTAA‑GGc                                                                                                                                                                                                                                          >  2:218625/1‑238 (MQ=255)                        tGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCACAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCCATGGTTGGGTAAACCCCTTGTTG‑AGGGGTGACGGCCCGCCGCTTTTTTTaac                                                                                                                                                                                                                         >  2:44401/1‑237 (MQ=255)                        tGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCGGGGGGTGAGTTCCCCCGGACGACAATGGTTGGGAAAAACCATTGTTGAAGTGGTTG‑GGCCCGCCGCTTCTTTTAcg                                                                                                                                                                                                                         >  2:222583/1‑238 (MQ=255)                            cgcgCATACCCGTCGGAAGTTGAGGGGCAGGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGaa                                                                                                                                                                                                                      <  1:284784/238‑1 (MQ=255)                             gcgcATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAGGGGGTCA‑GGCCCGCCGATTTCTTTAGCggaaa                                                                                                                                                                                                                    >  2:530208/1‑235 (MQ=255)                             gcgcATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACa                                                                                                                                                                                                                    >  1:547047/1‑239 (MQ=255)                              cgcATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGCTTAGTTGGGTTTATGAAGGGCTGATTGCGGGGGATCCGCACGTTTACATTACAGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGCAGCAAATGGATGGGGAAACCCATTTGTAAAGAGGTCA‑GGCCCCCGGCTTTTTTTACCGaaaaa                                                                                                                                                                                                                   >  2:349099/1‑236 (MQ=255)                                               ctgatgGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGGTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCCAGAAAGATGCTGCGGCTGGGTTCCCCCGGACGGCACTGGTAGGGTAAAACCATTGTTGAAGGGGTAA‑GGCCCGCCGATTTCTTTACCGAAAAAATAGATGaaaaagataa                                                                                                                                                                                                  >  2:1001341/2‑234 (MQ=255)                                               atgatgGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGAGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCGGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGACCGGATTGATGAAAAACGTCa                                                                                                                                                                                                  >  1:616327/1‑239 (MQ=255)                                                    gggCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGCTTTCTTTACCGAACAGATTGTTGAAAAACGCCAGGCtg                                                                                                                                                                                             >  1:955750/1‑239 (MQ=255)                                                             gcaaaaGCAAACCGCACGCGTGGGCCGTCAGCAAGGGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATGGGGGTGGTCCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCg                                                                                                                                                                                     <  1:900595/235‑1 (MQ=255)                                                                     aaaCCGCACGGGGGGCCCGTCAGCAAGGGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCGGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATATGGTGGCGAACTTTa                                                                                                                                                                              <  1:218625/237‑1 (MQ=255)                                                                      aaTCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACg                                                                                                                                                                           >  1:387058/1‑239 (MQ=255)                                                                        tgccgcgggTGGTCGGTCAGGAAGAGGGGCTGTTTGGTTCGTTTTAAGGAGGGCGGGTTGGGGGAGATCCGCTGGGTTAGTTTACCGCCGT‑TAGTCGTAGCGCAAGAAAGTTGGTGGGGTTGGGTCCCCCGGGAGGGCAATGGTTGGGTAAACCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGt                                                                                                                                                                         <  2:176181/231‑1 (MQ=255)                                                                        tcccACGCGTGGTCGGTCAGCAAGAGGCGCTGTTTAGTCGGGTTTATAAAGGTCTGATTGCGGTGGATCCGCATGGTTACATGTCCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAACCCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGGCAGGATGTGGTGGCGAACTTTAACGGt                                                                                                                                                                         <  1:503516/236‑1 (MQ=255)                                                                             cgcgTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTAGTGAAAAACGTCAGGATGTGGGGGGGAACTTTAACGGTCTGAg                                                                                                                                                                    >  1:514066/1‑239 (MQ=255)                                                                              gcgTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTTCCGACAAAATTGATAAAAAACGTCAGGAGGGGGTGGCGAACTTAAAGGGTcgaagg                                                                                                                                                                   >  2:503516/1‑234 (MQ=255)                                                                                        tCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGACCAGATTGATGAAAAACGTCAGGTGTTGGTGGCGAACTTTAACGGACTAAGCGTTATTGC‑Ca                                                                                                                                                        >  1:417681/1‑239 (MQ=255)                                                                                              agagGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCCGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTTTTGC‑CAACCggg                                                                                                                                                  >  1:824523/1‑239 (MQ=255)                                                                                                     cTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTAt                                                                                                                                           >  1:492816/1‑239 (MQ=255)                                                                                                      tGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGGTGAAAAACGTCAGGATGTGGTGGCCAACTTTAACGGCCTGAGCGTTATTTC‑CAACCGGGAAgacaa                                                                                                                                           >  2:83733/1‑234 (MQ=255)                                                                                                          tAGTTCGGTTTAAGAAGGTCTGTTTGGGGTGGACCCGAATGGTTCAATTACCGCCAT‑CAATCGAAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGtt                                                                                                                                      <  1:530208/239‑1 (MQ=255)                                                                                                               cGGTTTATGAAGGGTTGAGTGGGGTGGTCCCGCATGGTTACATTACCGCCAT‑CAAGCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCATTGTGTGGGTAACCCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGt                                                                                                                                 <  1:218396/239‑1 (MQ=255)                                                                                                                  tttATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGtgat                                                                                                                              >  1:606379/1‑239 (MQ=255)                                                                                                                  tttATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAATTCCCCCGGACGGCAATGGTGGGGAAAAACCATTGTTGAGGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGACAAACGTCGGGCAGTGGTGGCGAACTTTAAAGGAGTAAGCGGTATTGC‑CAACCGGGAGGCAACTAGTTTAgatgat                                                                                                                              >  2:29588/1‑239 (MQ=255)                                                                                                                  tttATGAAGGGCTGATTGCGGTGGATCCGCATGGTTAAATTACCGCCAT‑CAATCGTACCGCAAGAAAAATGCTGGGGCTTAGTTCCCCAGGACGGCAATGGTTGGGTAAACCCATAGTTGCAGGGGCCA‑GGCCCGCCAATTCTTTTACCGAACAGCTTGTTGCAAAACGTACGGATGTGGTGTCGAACTTTAAGTGTCTGCACATTATTAC‑CAACCCGGACGCAATTCGGTCGGGtggt                                                                                                                              >  2:507655/1‑237 (MQ=255)                                                                                                                       gAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGGTAAGGGGAT‑GAttt                                                                                                                        >  2:11572/1‑239 (MQ=255)                                                                                                                           ggCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAATGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGTTTGGGTAAACCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCCGGTGTTGGGGGCGAAATTTAACGGTATGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGGT‑GGttggtgg                                                                                                                    >  2:656093/1‑234 (MQ=255)                                                                                                                            gCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAGTCGTAACGCAAGAAAGGTGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGGAAACCCATTTTTGAAGGGGTCA‑GGCCCGCCAATTTCTTTACCGAACAGATGTATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGCTGGGGCGTTATTGC‑CAACCGTGAAGCCTTTCGTGCAGGTGAT‑GATTTgtggg                                                                                                                   >  2:843488/1‑235 (MQ=255)                                                                                                                            gCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGGGCGTTATTGC‑CAACCGGGAAGCTATTCGTTTAGGGGAT‑GATTTGCTgg                                                                                                                   >  2:716386/1‑239 (MQ=255)                                                                                                                                aTTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGGGAT‑GATTTGCTggggg                                                                                                                >  1:96983/1‑238 (MQ=255)                                                                                                                                 tggtgggggATCCGGTTTGTTTCATTTCCCCCATGCAAGTGTAACCTAGGAATGGAGCTGGGGCCGGTTTCCCCCGGACGGAAATGGTTGGGTAAACCGATTGTGGAGGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGc                                                                                                               <  2:616327/232‑1 (MQ=255)                                                                                                                                          aaccGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCAGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGCGGTGGTGGCGCCCTTTAACGGTCTGAGCGTTATTGC‑CCACCGGGGAGCACTTCGTTACGGTGAT‑GATTTTCTGGGGGCCATTAtaa                                                                                                       >  2:765198/3‑235 (MQ=255)                                                                                                                                                   gTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGACCAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGGCTGAGCGTTATTGC‑AAACCGGGAAGCTATTAGGTCAGGTGAG‑GTTTTGTTGGGGGCCATTATCAGC‑tttgtaggt                                                                                           >  2:691453/1‑233 (MQ=255)                                                                                                                                                       cATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAg                                                                                       >  1:832949/1‑239 (MQ=255)                                                                                                                                                          tACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGCT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACg                                                                                    >  1:484024/1‑239 (MQ=255)                                                                                                                                                             cGCCAT‑CACCCGTAACGCAAGACAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGTGTGGGTAAACCCATTGTTGCAGGGGTCC‑GGCCCCCCGATTTCTTTACCGAACAGATTGTTGAAAAACACCGGGATG‑GGTGGTCGACTTTTAAGGACCGGGGGTTTATGCACAACCGGGAAGATATTC‑ATCCGGGGATGGATTTGC‑GGGGGGCACTGTCAGCAGTTGGTAGAAAAGAGGaaa                                                                                 >  2:1029196/1‑239 (MQ=255)                                                                                                                                                             cGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGaaa                                                                                 <  1:44401/239‑1 (MQ=255)                                                                                                                                                              gCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGGAGCTATTCGTTCAGGGGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTGGTAAAGCCGAAAt                                                                                >  1:549944/1‑239 (MQ=255)                                                                                                                                                                 aT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCAATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGGGAT‑GTTTTGCTGTGGGCCATTATCGTG‑TTTGCTATTAAAAACGAAAttt                                                                              >  2:227084/1‑238 (MQ=255)                                                                                                                                                                  t‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGCCGAAATTTc                                                                             >  1:835742/1‑238 (MQ=255)                                                                                                                                                                        cGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGCTGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCt                                                                       >  2:595992/1‑239 (MQ=255)                                                                                                                                                                             cgaaagaaagaTGTTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAACCCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATg                                                                  <  1:929341/236‑1 (MQ=255)                                                                                                                                                                                  gaaagaTGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGGGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCa                                                              >  1:344203/1‑238 (MQ=255)                                                                                                                                                                                        tgttGGGGCTTGTTTCCCCCGGGCGGAATTGGTTGGGTAAACCCGTTGTAGAAGGGGTAA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGAc                                                        <  2:96983/235‑1 (MQ=255)                                                                                                                                                                                               gtttgttccccccGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATa                                                 <  1:655455/230‑1 (MQ=255)                                                                                                                                                                                                             gtcgtcattttttGGGTAACCCCTTTTTTGAGGTGGTCA‑GGCCCGTCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGGGGTTGCTAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTgag                                  <  2:514066/228‑1 (MQ=255)                                                                                                                                                                                                             gACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTAt                                                                                                         <  1:554037/169‑1 (MQ=255)                                                                                                                                                                                                             gACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGATGGGGGCCATTAt                                                                                                         >  2:554037/1‑169 (MQ=255)                                                                                                                                                                                                              aCGGCATTGGTTGGGTAACCCCATTTTTAATGTGGCCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGGCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGTAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTgag                                  <  1:545622/238‑1 (MQ=255)                                                                                                                                                                                                              aCGGCAATGGTTGGGTAAACCCATTTTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTgag                                  <  1:691453/238‑1 (MQ=255)                                                                                                                                                                                                                tgtattGGTTGAGTAAACCTTTTTTTGAA‑TGGGCACGCCCCTCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGagag                                <  2:835742/233‑1 (MQ=255)                                                                                                                                                                                                                gtaaagggTTGGGTAAACCCATTGTTGAA‑TGGGCAGGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGagag                                <  1:1001341/232‑1 (MQ=255)                                                                                                                                                                                                                                 cccATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGc‑c                                                                                                                                                         <  1:528258/102‑1 (MQ=255)                                                                                                                                                                                                                                 cccATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGc‑c                                                                                                                                                         >  2:528258/1‑102 (MQ=255)                                                                                                                                                                                                                                  ccATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGacac             >  1:108207/1‑239 (MQ=255)                                                                                                                                                                                                                                   ggttgttgAATTGGCCA‑GCCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCGCCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACg            <  1:519164/237‑1 (MQ=255)                                                                                                                                                                                                                                     ttttagAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAg          <  1:716386/234‑1 (MQ=255)                                                                                                                                                                                                                                      tgttgAAGTGTTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGCGCGGTATTGC‑CAACCGGGAAGCTATTCGTTCAGGGGCT‑GATTTGCTAGGGTCCATTATCCGC‑TTTCGTAGTAAAAGCGAAATTTCCACCATCAATGCGTAACTTACGCACATAAAACAAAAAGATAAGACCCTTCGTACTTTGCGccaataca         >  2:945858/1‑234 (MQ=255)                                                                                                                                                                                                                                      tgttgAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGCGAGCCTTCTTCCATTGCGACACGAGc         >  1:767929/1‑239 (MQ=255)                                                                                                                                                                                                                                      tgttgAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACCCGAGc         >  1:293664/1‑239 (MQ=255)                                                                                                                                                                                                                                      tgttgAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGACCCTTCGTCCATTTGGACACGAg          >  1:398278/1‑238 (MQ=255)                                                                                                                                                                                                                                        ttgAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCAt       <  1:11572/239‑1 (MQ=255)                                                                                                                                                                                                                                        tgaaaGTGTTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCAt       <  1:227084/236‑1 (MQ=255)                                                                                                                                                                                                                                             gTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGCCGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAACAAAAACGTTGAGAGCCTTCGTACATAGCGACCCGAGccacacaa  >  2:267598/1‑232 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  TCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAA  >  NC_000913/652813‑653286

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