Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,223,638 T→C intergenic (‑10/‑190) iclR ← / → metH DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,6380TC97.7% 148.7 / NA 44intergenic (‑10/‑190)iclR/metHDNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/28);  minor base A (1/0);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTG  >  NC_000913/4223398‑4223881
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCTGTAATTGCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTa                                                                                                                                                                                                                                               <  1:216591/251‑1 (MQ=255)
   aaCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCg                                                                                                                                                                                                                                            <  1:274170/251‑1 (MQ=255)
          cGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑GTTCGCGGG‑AATGGGTGCGACCATGACAGACTCCCTTTTTTTTATCTT‑GGa‑aa                                                                                                                                                                                                                                    >  1:82687/1‑250 (MQ=255)
                 aaCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGTTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AAtcattttc                                                                                                                                                                                                                            <  1:62533/251‑1 (MQ=255)
                 aaCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AAtcattttc                                                                                                                                                                                                                            <  2:164779/251‑1 (MQ=255)
                             gCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGt                                                                                                                                                                                                                 <  1:303700/250‑1 (MQ=255)
                             gCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGt                                                                                                                                                                                                                 <  1:225927/250‑1 (MQ=255)
                                 cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGAAACCCGGCTTGTTGCGCCAGTTCCGGGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGAGGCAACGGCGGGGTTTTTGCCGCG‑GGTCTCGGGGCATGGG‑GGGGCCATGAAAGTATCCTTTTTCTGTATCGT‑GGGGCATCATTTT‑ATTTTTTTTGttttgc                                                                                                                                                                                                            >  1:137482/1‑247 (MQ=255)
                                 cGTGGTGAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTTGAGTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGc                                                                                                                                                                                                             <  2:245538/250‑1 (MQ=255)
                                     gTTAGCAGGCGGCTTGACGTGGAAGTGGGTAACCCGTCTTGTCGCGCCAGTTCCGGGAGTGCCACACAGCCTTTGGAGTCGGCAATCCACTACAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCTCTGGTGCTGTTGCAACGGCGGGTTTTCTGCCGCT‑TGTCGTGGG‑AATGTGTGCTATCATGTCAGTCTCCTTTATCTGTAACGT‑GTA‑AAGCATTTTCAATTTTATTTTTAGCTAAt                                                                                                                                                                                                         <  2:136408/250‑1 (MQ=255)
                                        aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCa‑a                                                                                                                                                                                                    <  1:136706/251‑1 (MQ=255)
                                             gCGGTGGGTCGGGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGtt                                                                                                                                                                                               <  2:257936/251‑1 (MQ=255)
                                               ggTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACGCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACATTCTTCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTac                                                                                                                                                                                             <  2:34175/251‑1 (MQ=255)
                                                 tGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCG‑GTTCTCGGG‑GATGGGGGCGGCCCTGACAGGCTCCTTTTTTTTTATAGT‑GGG‑AATCATTTTTCTTTTTATTTTTAGGTAATGCA‑AAAATTactt                                                                                                                                                                                           >  1:268353/1‑250 (MQ=255)
                                                  gggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTactga                                                                                                                                                                                          <  1:77054/251‑1 (MQ=255)
                                                         ggAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCG‑GTTCGCGGG‑AATGGGTGGGACCATGACAGGCTCCTTTTTCTTTATCGT‑GGG‑AATCATTTTTTTTTTTTTTTTTTACTAATGCA‑ATAgatcgggaaagcaa                                                                                                                                                                                    >  1:230138/1‑237 (MQ=255)
                                                         ggAATTGGGTAACCCGGCGTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑At                                                                                                                                                                                                   <  2:230138/235‑1 (MQ=255)
                                                              tGGGTATCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCGCG‑GTTCGCGGG‑GATGGGGGCGGACATGACAGGCTCCTTTTTTTGTATAAT‑GGG‑AAACATTTTCATTTTTTTTGTTTGCTTATGCA‑ATAGGTACTGAA‑CTGATcccatt                                                                                                                                                                             >  1:93164/1‑247 (MQ=255)
                                                               gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGa                                                                                                                                                                            <  1:240793/251‑1 (MQ=255)
                                                                 gTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCGCGGTTTC‑CGGG‑GATGGGTGCGACCCTGACAGTCTCCTTTTTTTGTGTTGT‑TGA‑AATTATTTTTATTTTTTTTGTTAGCTTATTCC‑ATTGTTactgta‑cttgtccgctttgt                                                                                                                                                                          >  2:75592/1‑239 (MQ=255)
                                                                   aaCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTa                                                                                                                                                                        <  2:93164/251‑1 (MQ=255)
                                                                     cccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCG‑TTTCGCGGG‑GAAGGGTGGGGCCATGACAGGCTCCTTTTTCTTTATTGT‑GGG‑AATCATTTTCCTTTTTTTTTTTTGCTAATTCA‑ATAGTTTCTGAA‑CTTATccactgggatta                                                                                                                                                                       >  1:242237/1‑240 (MQ=255)
                                                                      ccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTGTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATg                                                                                                                                                                     <  2:213434/251‑1 (MQ=255)
                                                                         gCTTGTTGCGCCAGTGCCGTGAGTGCCACACGGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACTGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGT‑GGA‑AATCATGTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTGCTGAA‑CTTGTCCTATGAGGTAATGTTg                                                                                                                                                                  <  2:273295/251‑1 (MQ=255)
                                                                            tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCG‑GTTCGCGGG‑GATGGGTGGGAGCATGACAGTCTCCTTTTTTTGTATCGT‑GGG‑AATCATTTTTATTTTTTTTTTTTGCTTATGCC‑ATAGTTACTGAA‑ATTATCCGATGGGGTAATGTTgta                                                                                                                                                                >  1:104211/1‑248 (MQ=255)
                                                                                        tCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGAGGCAACGGCGGGTTTTCTTCCGCG‑TTTCGCGGG‑GATGGGTGCGACCCTGACAGCCTCCTTTTTCTGTTTCGTTGGA‑AA‑CATTTTTATTTTTTTTGTTTGTTAATTCA‑ATAGTTACTGAA‑ATGATCCTCAGAGTTAATGTTGATAAAATCTCAtttt                                                                                                                                                   >  1:6232/1‑248 (MQ=255)
                                                                                             gAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtg                                                                                                                                              <  2:3383/251‑1 (MQ=255)
                                                                                                   cacaCTGCCATTGGATTCGGCAATCCACTCCAGTACTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCG‑TTTCGCGGG‑GATGGGGGCGGCCATGACAGGCTCCTTTTTTTGTATAGT‑GGG‑AATCATTTTCATTTTTTTTGTTAGCTAATGCA‑AAAGTTACTGAA‑CTGATACGATGAGTTAATGTTTAACAAAACTCATTTTTCGTgggggtg                                                                                                                                        >  1:236234/1‑250 (MQ=255)
                                                                                                               ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCCCGGG‑GATGGGTACGACCATGACCGTCTCCTTTTTCTTTATCGT‑GGA‑AATCCTTTTCATTTTTATTGTTATCTAATTGA‑ATAGTTACTTGA‑CTTATCCGATGAGTTAATTTTGGACAAATCTCATGTTGCGTGGTGGTCGCTTTTTCcac                                                                                                                             >  1:308604/1‑250 (MQ=255)
                                                                                                                                  aGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTAt                                                                                                         <  2:308604/251‑1 (MQ=255)
                                                                                                                                     aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGt                                                                                                      <  1:75592/251‑1 (MQ=255)
                                                                                                                                               cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGTGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAtt                                                                                            <  1:188716/251‑1 (MQ=255)
                                                                                                                                                       aaaGACTGACCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTGTTGCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACGAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCTTTTATGCCAGTATTGTTTGTTGAATTTTTATTaa                                                                                    <  2:67507/251‑1 (MQ=255)
                                                                                                                                                         aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCCGG‑CTTCTCGGG‑AAAGGGGGAGACCATGACCGTCTCCATTTTCTGTTTCGG‑GGG‑AATCATTTTTCTTTTTATTGTTATCTAATGCACAT‑GTTACTGAAACT‑ATCCGCTGAGTTAATGTTTTACAAAACTCCTTTTGCGCGGTGGTCGCTTTTATCCCAGATGCGGTTTTGCTCGTGTTGTTTGTTTAAtttttttttta                                                                                   >  2:311470/1‑244 (MQ=255)
                                                                                                                                                                           cGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCg                                                                 <  2:234019/250‑1 (MQ=255)
                                                                                                                                                                                                          ccGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTCGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACt                                 <  1:39724/251‑1 (MQ=255)
                                                                                                                                                                                                                 ttCGCGGG‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                            <  1:99491/250‑1 (MQ=255)
                                                                                                                                                                                                                      ggg‑AATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGAAAAGTGGAACAACTGCGAGCGCAGt                      >  1:228504/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCggg                                                               >  1:156870/1‑195 (MQ=255)
                                                                                                                                                                                                                                      cATGACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCggg                                                               <  2:156870/195‑1 (MQ=255)
                                                                                                                                                                                                                                         gACAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCATCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg    <  1:111848/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aCAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc             >  1:8957/1‑241 (MQ=255)
                                                                                                                                                                                                                                          aCAGTCTCCTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc             <  2:8957/241‑1 (MQ=255)
                                                                                                                                                                                                                                           cAGTCTCATTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGATGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTg  >  2:89777/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCG‑TTTCGCGGG‑AATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGT‑GGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCA‑ATAGTTACTGAA‑CTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTG  >  NC_000913/4223398‑4223881

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