Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	2,606,884	C→G	14.8%	R438G (CGG→GGG)	hyfF →	hydrogenase 4 subunit F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	2,606,884	0	C	G	14.8%	95.3 / 2.4	54	R438G (CGG→GGG)	hyfF	hydrogenase 4 subunit F
Reads supporting (aligned to +/- strand):  ref base C (17/26);  new base G (2/6);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACC  >  NC_007779/2606647‑2607093                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                        ggCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACg                                                                                                                                                                                                                        <  1:1222999/239‑1 (MQ=255)      tGTACTGCTCAAGTACGGCACGCGCGATCTCAACGGCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑gg                                                                                                                                                                                                                  <  1:1614831/239‑1 (MQ=255)           tGCTCAAGTCCGGCCCGCGCGGTCTCACCGTCTTTTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTgcg                                                                                                                                                                                                             <  1:1482218/239‑1 (MQ=255)           tGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGG‑ACGGGT‑gggcgc                                                                                                                                                                                                              >  1:377727/1‑234 (MQ=255)             tcaagggctggCACGGGCGATCACGTGGTCGTCGGTGGGCTTCTCAAAATCGTGCTGTTTCCCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAGTTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGcgcg                                                                                                                                                                                                           <  2:75296/230‑1 (MQ=255)             cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACCTTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCCGGTGCTGGC‑GGG‑CCCGGG‑GCGGAg‑ggtggccgg                                                                                                                                                                                                           >  2:1866635/1‑232 (MQ=255)             cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCCGGG‑TCGGGT‑gggttctgt                                                                                                                                                                                                           >  1:608939/1‑232 (MQ=255)             cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCCGGG‑ATGGagtggtgtgcg                                                                                                                                                                                                            >  1:567360/1‑229 (MQ=255)             cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCAAGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAAGCTGGTGCTGGC‑GGG‑GATGGT‑ACgggcgtgcgc                                                                                                                                                                                                               >  2:971291/1‑228 (MQ=255)             cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCCCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCCCCGGGCTGGCACGTAATCTCCTGCTGATTATCGTCCTGCTGTTTTTGCCGTTAACCCGGGGGCTTGG‑GGG‑GCTCGT‑GAGCGT‑GGCTCcgcg                                                                                                                                                                                                           >  2:2035694/1‑239 (MQ=255)                       gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGGCCCTGGA‑AGGGAT‑‑GCTGCGCGGGGGTTAAtt                                                                                                                                                                                                 >  1:791754/1‑238 (MQ=255)                            cgcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGGT‑TGGTGCGCGGGGGTTAAT‑gtggaa                                                                                                                                                                                           >  1:2010156/1‑234 (MQ=255)                                           gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGCACTTATGACCC‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCCGTTAACGCTGGCGCTTGC‑GGG‑GCCGGTGACGGTT‑GGTCCTGCGGGGGTT‑TT‑GGCGGCACcgccgcccggccg                                                                                                                                                                            >  2:1701396/1‑232 (MQ=255)                                                              ccATGCCCTTTCCCGCCCTTCTG‑TTGGCGGACGGTGCCGTGGCGCTGGCGGGCCTCCCGCCCTCAAACATCTCTCTTAGGAAATTTATGACCT‑TTACCGCCTGACTGGCACGTAACCACCTGCTGATTATCGTCCTGCTGTTATTGGTGTTAACGCTGGTGCTGGG‑GGG‑CCGGGT‑AGGGTG‑GGCTGGGGGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGc                                                                                                                                                      <  2:232495/238‑1 (MQ=255)                                                                 tGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAAAGCTGGTGCTGGCGGGG‑CCGGGT‑AGGGAG‑GGC‑GCGCGGGGGGTTAA‑GGGGACAACCCCGCACGGGTGTTACACAGGGGTG‑‑TCTCGGCTg                                                                                                                                                    >  2:1862605/1‑239 (MQ=255)                                                                          cccccgTGGTTTTTGGCGG‑GGGGTGGGTGGGGTCGGCAGGTGTGCCCCCCCTTCATATTTTTTTTTGCGATTTTTAGCCCATTTCCCGCCGGCCTGGCTGGAAACCCCCTGCTGTTTATCGTCTGGTTTTTTTTGTTGTTACCGTGGGGGTTGGC‑GGG‑CCTGGT‑A‑GGGT‑GGGCGGGGGGGTTTTAATGGGCAAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGaccac                                                                                                                                           <  2:1334357/236‑1 (MQ=255)                                                                                                                        cgccCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGGAACCGCCGCC‑GGCCGGTA‑AACGGGGGGA‑TCCCGGCTGGGTGGCCACCCCGCCCAAGGTGGTTTTCGGGGTCAT‑GATGCTGGCgggggggac                                                                                             >  2:520232/1‑230 (MQ=255)                                                                                                                          cctttCAACATTTTTCTTAGCTAATTTATGCCCA‑TTACCGCCGGACTGGCACGTAATCCCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGTT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCa                                                                                          <  1:1414636/236‑1 (MQ=255)                                                                                                                                     ttttttTAGCGATTTTATGACCT‑TTACCGCCGGACTGGCACGTAATCCCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCTGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAAc                                                                               <  1:1419759/239‑1 (MQ=255)                                                                                                                                     ttttCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGGGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGGGA‑TCTCGGCTGGGTGACCACCCCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGGGATGGGAACGCATATTCCCCaa                                                                                >  1:829822/1‑238 (MQ=255)                                                                                                                                       ttCTTAGCGATTTTATCCCCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTACCGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCt                                                                             <  1:341396/239‑1 (MQ=255)                                                                                                                                            aGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGGGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGGGA‑TCTCgg                                                                                                                                                       >  1:222260/1‑161 (MQ=255)                                                                                                                                                   ttAGGACCA‑TTACCGCCGGACTGGCACGTAATCCCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGGGCTGGC‑GGG‑CCTGGT‑ACGGAG‑GGCTGCGCGGGTGTTAAG‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAAc                                                                                             <  1:68114/211‑1 (MQ=255)                                                                                                                                                             ttACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGGCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACCACCTGTCCTCAGGATCCTggcgggg                                                        >  1:1524881/1‑238 (MQ=255)                                                                                                                                                                  gccgGATTGGCACGTAATCCCCTGCTGATTATCGTCCTGCTGTTATTGCGGTTAACGCTGGTGCGGGC‑GGG‑CCTGGT‑ACGGGT‑GGCGGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTc                                                   <  1:1623376/239‑1 (MQ=255)                                                                                                                                                                   ccgGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGGT‑GGCTGCGCGGGCGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCgggggggtc                                                                                                                                                            >  1:37374/1‑130 (MQ=255)                                                                                                                                                                     ggACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑GCTGGT‑ACGGAT‑GTCTGCGCGGGGGTTAAT‑GGGGAAACCGCCGCC‑GGGCGGTT‑AACGGGGGGA‑TTTCGGCTGGGTGACCACCCCCCCCATGGGGGTTCTTCCGGGCATGGTTGCTGG‑GGTGGGGAGCCAAAATTCACAACCTTTCATCCGGGACCTtgcggggggctcca                                                 >  2:262341/1‑230 (MQ=255)                                                                                                                                                                      gACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACcgccgc                                                                                                                                                                                   <  1:1174084/111‑1 (MQ=255)                                                                                                                                                                                tAATCTCCTGCTGTTTATCGTCTGGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCGGGT‑TGGGAT‑GGCGGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGgcgc                                                      <  1:1479450/222‑1 (MQ=255)                                                                                                                                                                                tAATCACCGGCTGATTATCGTCCTGCTGTTTTTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGTT‑ACGGAT‑GGCTGCGGGGGTGTAAAG‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                     <  1:1701396/239‑1 (MQ=255)                                                                                                                                                                                   tccccTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCat                                                                                         <  1:1210819/181‑1 (MQ=255)                                                                                                                                                                                     cccTGCTGATTATCGTCCAGCTGTTATTGCGGTTAGCGCTGGGGCTGGC‑GGG‑CCTGGT‑GCGGAG‑GGCTGCGGGGGTGTTAGT‑GGAAAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTg                                                                            <  1:820917/194‑1 (MQ=255)                                                                                                                                                                                      ccggctgATTATCGTCCTGCTGTTTTTGCGGTTAACGCTGGTGCGGGC‑GGG‑CCTGGT‑ACGGAG‑GGCTGCGCGGGTTTTAAG‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGATCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGa                               <  1:104780/236‑1 (MQ=255)                                                                                                                                                                                          ctgATTATCGTCCTGCTGTTATTGCTTTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGc                                                                                                                                                      <  1:253315/117‑1 (MQ=255)                                                                                                                                                                                              ttttCGTCTTGTGGTTAGTGCTGTGACCGCTGGTGCGGGG‑GGG‑CCTGGT‑CGGTAT‑GGCTGCGCGGGTTTTAAG‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGa                        <  1:1347881/235‑1 (MQ=255)                                                                                                                                                                                              ttATCGTCGTGCTGTTTTTGCTGTTACCGTGGGTGGTGGG‑GGG‑CCTGGT‑CGGGAG‑GGCGGCGCGGGTGTTAAG‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGa                        <  1:484095/238‑1 (MQ=255)                                                                                                                                                                                              ttATCGTCCTTGTGTTTTTGCGGTAACCGCGGGTGCGGGC‑GGC‑CCTGGT‑CCGGTT‑GGCTGCGGGGGTTTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGa                        <  1:1862605/238‑1 (MQ=255)                                                                                                                                                                                               tATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCCGGG‑GCCGAT‑TGCTGTGCGGGTGTTAAA‑GGGGAAACCGCCCCC‑GGGCGCTA‑AACGGGGGGATTTTCGGGCGG‑TGACCACCCCCCCCAATGCGGTTCTGGTGGTCTT‑GCTGCGGGCGCTGGGGACGCATATTTCACAACCCGTCATAAgaac                                                                  >  2:820917/1‑192 (MQ=255)                                                                                                                                                                                                  cGTCCTGCTGTTTTTGCTGTTAACGCGGGTGCTGGG‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGTTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgc                   <  1:1573817/239‑1 (MQ=255)                                                                                                                                                                                                    tCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGGAAACCCCGCC‑GGGCGTTA‑ACCGGGGGGG‑TCTTGGCTGGGTggtca                                                                                                                                            >  2:253315/1‑113 (MQ=255)                                                                                                                                                                                                    tCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCCGGT‑ACGGAT‑GGCTGCGCGGGAGTTAAT‑GGGGAAAACGCCCCA‑GGGCGTTA‑ACCGGGGTGA‑TCTCCGCTGGGTGACCCCCTCCCCCATGGGGATTCTGCTGGt                                                                                                                   >  1:1457385/1‑142 (MQ=255)                                                                                                                                                                                                     ccTGCTTTTTTGGCTGTTAACGCTGGTGCTGGC‑GGG‑CGGGGT‑‑CGGATGGGCGGGGGGGGTGTTAAG‑GGCAAACCCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgcctg                <  1:855036/239‑1 (MQ=255)                                                                                                                                                                                                       tggtttttttGCTGTTAACTCTGTTGCGGGC‑GGG‑CCTGGT‑CGGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca              <  1:1803457/231‑1 (MQ=255)                                                                                                                                                                                                       tgctgTTATTGCTTTTACCGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGTAT‑GGCTGCGCGGGTTTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca              <  1:955563/239‑1 (MQ=255)                                                                                                                                                                                                       tgctgTTATTGCTGTTAACGCTGGTGGTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCAAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca              <  1:1273531/239‑1 (MQ=255)                                                                                                                                                                                                       tgctgTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGGGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGGGCTTCCACTATAGTCCTCTCAGGGACGCACGATctg                    >  1:515449/1‑233 (MQ=255)                                                                                                                                                                                                       tgctgTTATTGCTGTTAACGCTGGTGCGGGC‑GGG‑CCTGGT‑ACGGAT‑GTCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca              <  1:262341/239‑1 (MQ=255)                                                                                                                                                                                                        gctgTTATTGTTGTTAACGCTGGGGGGGGC‑GGG‑CTTGGT‑CGGGAG‑GGCTGCGGGGGTTTAAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca              <  1:1466065/238‑1 (MQ=255)                                                                                                                                                                                                         ctgTTATTGCTGTTACCGCGGGTGCGGGC‑GGG‑CCTGGT‑ACGGAT‑GGCGGCGCGGGTTTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCat                                                                                         <  1:356307/162‑1 (MQ=255)                                                                                                                                                                                                          tgTTATTGCTTTTACCGCGGTTGCTGGC‑GGG‑CCTGGC‑ACGGAG‑GGCTGCGCGGGTGTTAAT‑GGCGAACCCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca              <  1:1541420/236‑1 (MQ=255)                                                                                                                                                                                                          tgTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTTTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCat‑gatg                                                                                                           <  1:141579/143‑1 (MQ=255)                                                                                                                                                                                                            tattggCTGTTACCGGTGGTGCGGGC‑GGC‑CCTGGT‑ACGGTT‑GGCGGCGTGGGTGTAAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGgcgc                                                      <  1:2179637/189‑1 (MQ=255)                                                                                                                                                                                                                              gTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGGGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGGTGACCACCTCGCCAATGGGGATTCTGCTGGTCAT‑GATGCTGGCGATGGGGACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACtata                                            >  1:859064/1‑186 (MQ=255)                                                                                                                                                                                                                              gTGCTGGC‑GGC‑CCTGTT‑CCGGAT‑GGCGGCGCGGTTTTTAAG‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCAcc  <  1:462420/228‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                        GGCAATGTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGG‑CGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCA‑TTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGC‑GGG‑CCTGGT‑ACGGAT‑GGCTGCGCGGGTGTTAAT‑GGCGAAACCGCCGCA‑GGCCGTTA‑ACCGGGGTGA‑TCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCAT‑GATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACC  >  NC_007779/2606647‑2607093

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