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.4% 128.5 / ‑4.9 39intergenic (‑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 (12/26);  total (13/26)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGG  >  NC_000913/4223396‑4223879
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
tCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCt                                                                                                                                                                                                                                             <  1:580336/250‑1 (MQ=255)
 cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTg                                                                                                                                                                                                                                            <  1:3410/250‑1 (MQ=255)
          gACGCACGAACCCCTGCTGTTGCATCGTGGTTCGCAGGCGTTGGGTCGTGGAATTGGGTAACCCGGTTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATGTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAAAGTC‑T‑CCTTTTTCTGTATCGTGGa                                                                                                                                                                                                                                   <  2:38978/250‑1 (MQ=255)
          gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGa                                                                                                                                                                                                                                   <  2:1766/250‑1 (MQ=255)
          gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGGGTTTCGCGGGGATGGGGGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTgggg                                                                                                                                                                                                                                  >  1:1766/1‑249 (MQ=255)
            cGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAAt                                                                                                                                                                                                                                <  1:315868/251‑1 (MQ=255)
            cGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAAt                                                                                                                                                                                                                                <  1:233623/251‑1 (MQ=255)
                cGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAAtcatt                                                                                                                                                                                                                            <  2:508331/251‑1 (MQ=255)
                cGAAACCCTGCGGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAAtcatt                                                                                                                                                                                                                            <  2:86050/251‑1 (MQ=255)
                    aCCCTGCTGTTGCTTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGGGTTTCCCGGGGATGGGGAGGACCATGGCAGTC‑TTCtttttttgttatgagggaattacttttt                                                                                                                                                                                                                         >  2:199284/1‑229 (MQ=255)
                         tctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCAGGTGCGGTGGCAACGGCGGGTGTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGGC‑T‑GCTGTTTCTGTATCGTGGAAATcatgttcatttt                                                                                                                                                                                                                    <  2:278812/249‑1 (MQ=255)
                          ctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCCGGAATTGGTGCGGACATGAAAGTCAT‑CTTTTTT‑TTTATCGTGGGAATTatttttatttttt                                                                                                                                                                                                                  >  1:334753/1‑250 (MQ=255)
                            gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtt                                                                                                                                                                                                                <  1:112130/251‑1 (MQ=255)
                             ttGCATCGTGGTTAGCAGGCGGTGGGTCGTGGCATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCCCGTTTTGCGGGGATGGGGGCGACCATGACAGAC‑T‑CCTTTTTTTGTATCGTGGAAAtcatttttcttttttttt                                                                                                                                                                                                               >  1:545669/1‑240 (MQ=255)
                                cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTa                                                                                                                                                                                                            <  1:499605/251‑1 (MQ=255)
                                 aTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAg                                                                                                                                                                                                           <  1:585427/251‑1 (MQ=255)
                                   cGTGGTTAGCAGGCGGTGGGTCGGGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCCCACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTAGCGTGGAAATCATGTTCATTTTTATTGTTAGCt                                                                                                                                                                                                         <  2:399893/251‑1 (MQ=255)
                                          ggcatgtggTGTGTCGTGGAATTGGGGAACCCGGCTTGTTGCGCCAGTGCCGGGAGTGCCACACGGCCATTGGATGCGGCAATGCTCTCCAGTAATTTCAGGCCACGCGTTAAAGACTTATCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCTGCGTTTCGCGTGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTTTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCAAAT‑GCaa                                                                                                                                                                                                 <  2:439819/244‑1 (MQ=255)
                                            cAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATa                                                                                                                                                                                               <  1:441270/251‑1 (MQ=255)
                                                                tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGAt                                                                                                                                                                            <  1:580201/250‑1 (MQ=255)
                                                                     aaCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTa                                                                                                                                                                      <  2:205611/251‑1 (MQ=255)
                                                                                  gcgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGCAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGAGTTTCGGGGGGGTGGGTGGGACCCTGGCAGCA‑T‑CCTTTTTCTGTCTAGTGGAAATCATTTTCTTTTTTTTTGTTTGCTAAT‑TCCATAGTTACTGTAATGATCCGAGTAGGTAATGTTGAACAAAt                                                                                                                                                         >  2:97132/1‑251 (MQ=255)
                                                                                        gTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGGGGGACCCTTACCATC‑T‑CCTTTTTCTTTTTCGTGGGAATCATTTTCTTTTTTTTTTTTTGCTAAT‑GCAATAGTTACTTAACTGCTCCGCTGAGTTAATGTTGTACAAATCTCATg                                                                                                                                                   >  2:181787/1‑251 (MQ=255)
                                                                                            cGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGAAGGGTGCGGCCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCCTTTTCATTTTTTTTGTTTGCTAAT‑GCACTAGTTACTGAATTTGTCCGATGGGTTTATGTTGAACCAATCTCATGTTgg                                                                                                                                               >  1:311166/1‑250 (MQ=255)
                                                                                                  tGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGAGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGTGGTTCGCGTGGATGGGGGAGGCCATGACAGTC‑T‑CCTTTTTTTGTTTTGTGGAAATTCTTTTTATTTTTTTTTTTAGCTAATTACAA‑AGTTAATTAACTGATCCGCAGAGTTAAAGttttacaaaacttactttgtgtggtg                                                                                                                                         >  1:573601/1‑227 (MQ=255)
                                                                                                                      cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTTCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGc                                                                                                                     <  2:309847/251‑1 (MQ=255)
                                                                                                                      cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATg                                                                                                                      <  1:403195/250‑1 (MQ=255)
                                                                                                                         cAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGGCCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGGAATCATTTTCATTTTTTTTGTTTGCTTAT‑TCAATAGGTACTGAAATGATCCGCTGTGTTTATGTTGTACAAAACTCCTGTTTCGTGGGGGGCGCttttttcacaaagtcgtt                                                                                                                  >  2:121437/1‑237 (MQ=255)
                                                                                                                               aCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggtg                                                                                                                                         <  1:314805/222‑1 (MQ=255)
                                                                                                                                            cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgt                                                                                                <  1:121437/250‑1 (MQ=255)
                                                                                                                                                ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAt                                                                                           <  1:120670/251‑1 (MQ=255)
                                                                                                                                                 cACGCGTTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTCTCGGGAAGTGGTGCGACCATGACAGTC‑T‑CCTGTTTCTGTATCGTGGAAATCAGTTTCATTGTTATTGTTAGCTAAT‑GCAATATTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTGGCGTGGTGGTCGTTTTTACCACAGATGCTTTTATGCCAGTATGGTTTGTTGAAtt                                                                                          <  2:498864/251‑1 (MQ=255)
                                                                                                                                                            gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgttg                                                                                              <  2:188429/236‑1 (MQ=255)
                                                                                                                                                            gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGCTTCCCCGGAATGGGTGGGACCAAGACAGTC‑T‑CCTTTTTCTTGATCGTGGGAATCATTTTTATTTTTTTTGTTTGATAAT‑GACATAGGTTCTGGACTGATCCGATGAGTTAATGTTTAACAAATCTCATGTTGCGTGGGGGTCGCTTTTACCACAGATGCGTTTATGGCAGTATGGTtttttggggtcgggagaggac                                                                               >  1:188429/1‑236 (MQ=255)
                                                                                                                                                                        ccAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGGGCGACCCTGACAGTT‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTTATTTTTATTGTTAGCTAAT‑GCAATAATTACTGAACTGATCCGCTGAGTTAATGTTGAACAAATCTCATGTTTCGTGGTGGTCGCTTTTTCCACAGATGCGTTTTTGCCAGTTTGGGTTGTTGAATTTTTTTTAAATCTGGGTTTAGCgg                                                                   >  1:120896/1‑250 (MQ=255)
                                                                                                                                                                                                                gTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgt                           >  1:581653/1‑251 (MQ=255)
                                                                                                                                                                                                                gTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgt                           <  2:545669/251‑1 (MQ=255)
                                                                                                                                                                                                                         gAATGGGTGCGACCATGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                  <  1:485039/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aTGACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCt    >  1:106457/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGAGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg  >  1:325762/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
TCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTC‑T‑CTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGG  >  NC_000913/4223396‑4223879

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: ‑