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.2% 124.4 / ‑6.5 36intergenic (‑10/‑190)iclR/metHDNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (8/27);  total (9/27)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGG  >  NC_000913/4223397‑4223879
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGTGTTCGACCATGACAGTCT‑CCTTTTTCTg                                                                                                                                                                                                                                           <  1:393371/250‑1 (MQ=255)
 gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGGGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCTCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGTCAATCCACTCCAGTAAGTTCAGGCCACGCGGTAAAGACTTAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGTCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTGCTGTa                                                                                                                                                                                                                                         <  2:465671/251‑1 (MQ=255)
             cACGAAACCCTGCTGTTGCATCGTTGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAAtca                                                                                                                                                                                                                             <  2:164564/251‑1 (MQ=255)
                           gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCC‑CGTTCGCGGG‑GATGGGGGGGGCCCAGACCGGCT‑CCTTTTTTTGTATCGGGGAAATCATTTTCAtttttttt                                                                                                                                                                                                               >  1:430864/1‑248 (MQ=255)
                             tGCATCGTGGTTAGCAGGCTGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTAGCGTGGAAATCATTTTCATTTTTATTg                                                                                                                                                                                                              <  1:77642/250‑1 (MQ=255)
                               cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTa                                                                                                                                                                                                           <  1:301075/251‑1 (MQ=255)
                                                       ggTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg                                                                                                                                                                                   <  2:18599/250‑1 (MQ=255)
                                                       cGTGGAATTGGGTAACCCGGCTTGGTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACGGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTGTCATTTTTATTGTTAGCTAAGGCATTAGGTactgaactg                                                                                                                                                                                   <  2:418272/251‑1 (MQ=255)
                                                                              gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGTCGCCGTTT‑GCGGG‑GATGGGGGGGGCCCAGACAGGCTTCTTTTTT‑TTTATTGTGGGAATCATTTTTTTTTTTTTTTTTTACTAATTTAATAATTACTGTAATGATCCGATGGGTTAATTTTGAAc                                                                                                                                                             >  2:21032/1‑250 (MQ=255)
                                                                              gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGC‑GGTTCGCGGG‑GATGGGGGCGGCCATGACAGTCT‑CCTTTTTTTGTATCGTGGAAAACATTTTTCTTTTTTTTGTTTGCTTATGCAATAGTTAATGAACTGATCCGGTGAGGTAATGTTggac                                                                                                                                                             >  1:94084/1‑247 (MQ=255)
                                                                                     cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAGCCGATGAGTTAATGTTGAACAAATCTCa                                                                                                                                                     <  1:202234/251‑1 (MQ=255)
                                                                                                    cacaCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTc                                                                                                                                      <  1:366100/251‑1 (MQ=255)
                                                                                                     acacTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGC‑GGTTCGCGGG‑GATGGGGGCGGCCATGACAGTCT‑CCTTTTTTTGTATCGGGGGAATCATTTTTATTTTTTTTTTTTGCTTATGCAATAGTTACTGAACTGATTCGATGAGTTAATGTTTAACAAATCTCATGTTTCGtggtgggcg                                                                                                                                     >  1:285146/1‑248 (MQ=255)
                                                                                                          gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGGTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCttt                                                                                                                                 <  2:399860/250‑1 (MQ=255)
                                                                                                               tGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACca                                                                                                                            <  2:73036/250‑1 (MQ=255)
                                                                                                                ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcaca                                                                                                                          <  1:61909/251‑1 (MQ=255)
                                                                                                                         aaTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGTTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGAGAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCGAATGCAATAGTTACTGAACTGATCCGATTAGTTAATGTTTAACAAATGGCATGTTGCGTTGTGGTCGCTTTTACCACAGATTCGttt                                                                                                                 <  1:431772/251‑1 (MQ=255)
                                                                                                                            ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATg                                                                                                              <  1:366427/251‑1 (MQ=255)
                                                                                                                                   aGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTAt                                                                                                       <  2:285146/251‑1 (MQ=255)
                                                                                                                                    gTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTGTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTGTACCACAGACGCGTTTATGCCAGTATg                                                                                                      <  2:306388/251‑1 (MQ=255)
                                                                                                                                      aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGt                                                                                                    <  2:140101/251‑1 (MQ=255)
                                                                                                                                         ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGTGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGttt                                                                                                  <  2:346189/250‑1 (MQ=255)
                                                                                                                                                cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGaa                                                                                            <  2:369864/249‑1 (MQ=255)
                                                                                                                                                cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCCC‑GTTTTGCGGG‑GATGGGGGCGACCCTGACAGTCT‑CCTTTTTTTGTATCGGGGGAATCATTTTCATTTTTTTTGTTAGCTTATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAAAAAATCTCATGTTGCGGGGTGGGCGCTTTTACCACAGAAGCGTTTTTGCCAGTATGGTTTGTTGaaa                                                                                           >  1:369864/1‑249 (MQ=255)
                                                                                                                                                  cgcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAtttt                                                                                        <  2:84355/251‑1 (MQ=255)
                                                                                                                                                                     gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAg                                                                      <  1:126495/250‑1 (MQ=255)
                                                                                                                                                                     gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAg                                                                      <  1:31854/250‑1 (MQ=255)
                                                                                                                                                                        cAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGTCCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATCGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCgtg                                                                  <  1:123983/251‑1 (MQ=255)
                                                                                                                                                                         aGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGC‑GGTTCGCGGGGAATGG‑TGGGACCATGACAGTCT‑CCTTTTTTTGTATCGTGGGAAACATTTTTATTTTTTTTGTTTTCTCATGCAAAAGTTAATGAACTGATCTGATGAGGTTATGTTGAAAAAAACTCATGTTGCGGGGTGGTCGCTTTTTCCACATATGCGTTTATGCCAGTATGGGTTGTTGTATTTTTATTAAATCTggggtgggggggg                                                                 >  1:49720/1‑241 (MQ=255)
                                                                                                                                                                             gCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGC‑GGTTCGCGGG‑AATGGGGGGGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGGAATCATTTTTCTTTTTTTTTTTAGCTAATGCCATAGTTACTGAAATGATCCCATGAGTTAATGTTGAACAAATATCATGTTGCGTGGGGGTCGCTTTTTCCACAAATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTTTTAAATCTGGGTTGAGGGTGTCggg                                                             >  1:389362/1‑251 (MQ=255)
                                                                                                                                                                                                     tttCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTgg                                      >  1:297397/1‑250 (MQ=255)
                                                                                                                                                                                                          gCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaacaac                                <  1:21032/251‑1 (MQ=255)
                                                                                                                                                                                                              c‑gTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcg                            <  1:437939/251‑1 (MQ=255)
                                                                                                                                                                                                              c‑gTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGAGAGCAGCAAAGTGGAACAACtgcg                            <  1:150635/251‑1 (MQ=255)
                                                                                                                                                                                                                    cgcgGG‑AATGGGTGCGACCATGACAGTCT‑CCTTTTTCTATATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTgcgc                       >  1:309472/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gACAGTCT‑CCTTTTTCTATATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg  <  2:309472/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
CGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGG‑AATGGGTGCGACCATGACAGTCT‑CTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGG  >  NC_000913/4223397‑4223879

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