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,6380TC98.2% 192.6 / NA 56intergenic (‑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 (17/38);  minor base A (0/1);  total (17/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCT  >  NC_000913/4223397‑4223883
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
cGCCAACCTGACGCACGAATCCCTGGTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATGGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTCCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTg                                                                                                                                                                                                                                                <  2:339573/250‑1 (MQ=255)
cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTg                                                                                                                                                                                                                                                <  2:281836/250‑1 (MQ=255)
    aaCCTGACGCACGAAACCCTGCTGGTGCAGCGTGGTGAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCg                                                                                                                                                                                                                                           <  2:332060/251‑1 (MQ=255)
     aCCTGACGCACGAAACCCTGCTGTGGCAGCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGt                                                                                                                                                                                                                                          <  1:417954/251‑1 (MQ=255)
          ccGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGaaa                                                                                                                                                                                                                                     <  1:462675/250‑1 (MQ=255)
                gAAACCCTGCTGGTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAAtcattt                                                                                                                                                                                                                               <  1:739455/251‑1 (MQ=255)
                 aaccccTGCTGTTGCAACGTGGTTAGCAGGCGTTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAAtcatttt                                                                                                                                                                                                                              <  2:555343/248‑1 (MQ=255)
                      ctgctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTTGC‑CGGGATGGGTGCGGCCCTGA‑CAGGATCCTTTTTCTGTATCCTGGGAAtccatttcattt                                                                                                                                                                                                                         >  1:140313/1‑251 (MQ=255)
                       ggctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATcattttcatttt                                                                                                                                                                                                                        <  2:429913/250‑1 (MQ=255)
                        gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATcattttcatttt                                                                                                                                                                                                                        <  1:69475/250‑1 (MQ=255)
                        gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCAttttt                                                                                                                                                                                                                       <  1:796598/251‑1 (MQ=255)
                              gcttCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGtt                                                                                                                                                                                                                 <  2:436320/248‑1 (MQ=255)
                                aTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTGCGCCGTTCCC‑GGGAAAGGGTGGGGCCATGATCA‑TCTCCTTTTTTTTTATTGTGGTAATCATTTTTATTTTTTTTGtttg                                                                                                                                                                                                               >  1:577456/1‑249 (MQ=255)
                                                gtggggTCGTGTAATTGGGTAACCCGGCTTGTTGCGCCATTTCCGTGAGTGCCACACTGCCATTGGATTCTGCAATCCACTCCATTAATTGCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTac                                                                                                                                                                                               <  2:109600/248‑1 (MQ=255)
                                                            aaTTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGC‑GGGGATGGGGGCGACCAAGA‑CAGTCTCCTTTTTCTTTTTTGTGGTAATCATTTTTTTTTTTTTTTTTAACTAATTAACTAGTTacttacactattcc                                                                                                                                                                                   >  1:557466/1‑240 (MQ=255)
                                                                   aaaCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGtt                                                                                                                                                                           <  1:253172/250‑1 (MQ=255)
                                                                                  cgcCAGTGCCGTGAGTGCCCCAGTGCCATTGGAGTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCGAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTGCTGTATCGTGGAAAACATGTTCAATTTGATTGTTAGCTAATGCAATAGCTGCTGAACTGATCC‑GATGAGTTAATGTTGAACAAAtc                                                                                                                                                            <  2:309461/251‑1 (MQ=255)
                                                                                            gTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCg                                                                                                                                                  <  1:471304/251‑1 (MQ=255)
                                                                                                gTGCCCCACTGCCATTGGATTCGGCAATCCATCCCAGTAATTTCAGTGCACTTGTTAAAAACTGAACCTGTCCAGTCGCTGTTGCGGTGGCAACTGTGGGTTTTCGGTCGCTTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAACTAATACAATAGTTATTGAGCTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTGGCgtggt                                                                                                                                              <  2:38166/251‑1 (MQ=255)
                                                                                                  gCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGTTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCAGTTTCATTTTTATTGTTAGCGAATGCAATAGTTACAGAACTGATCC‑GATGAGTTaa                                                                                                                                                                         <  2:37424/222‑1 (MQ=255)
                                                                                                       acTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCGCGCTTCTC‑GGGGAGGGGGGCGACCCTGA‑CAGGCTCCTTTTTTTTTATCGTGGGGATCATTTTTTTTTTTTTTTTTTACTAATGCAATAGTTTCTGAACTGTTCC‑GATGAGTTAATTTTTAACaaaactaattattttgggtgggggct                                                                                                                                       >  2:155316/1‑228 (MQ=255)
                                                                                                          gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCTTTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTa                                                                                                                                                                                                <  2:551229/192‑1 (MQ=255)
                                                                                                            cATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCCGGGATTGGG‑GCGGCCATGT‑CAATTTCCCTTTTCTGTTTCGGGGGGAACATTTTTATTTTTTTTGTTCACTTATGCAATAGTTACTGTACTGATCCTGAT‑AGGTAATGTTGAACAAAtctccttttgctagggtgttgcctttat                                                                                                                                  >  1:117278/1‑227 (MQ=255)
                                                                                                                     cGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGc                                                                                                                         <  2:765525/251‑1 (MQ=255)
                                                                                                                                      aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATgg                                                                                                         <  2:140313/250‑1 (MQ=255)
                                                                                                                                           cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgt                                                                                                    <  1:254008/250‑1 (MQ=255)
                                                                                                                                             ggCCTCGCGTTAAACACTGTACCGCGCCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GTGAATGTGTGCGACCATGA‑CAGTCTCCTGTTTCTGTATCGTGGAAAGCATTTTCATTTGTATTGTTAGCTAATGCAATAGTGACTGAACTGATCC‑TATGAGTGAATTTTGAACAAATCTCATGTTGCTTGGTGGTCGCTTGTACCACAGAGGCGGTTATGCCAGTATGGTttgttta                                                                                                 <  2:376692/251‑3 (MQ=255)
                                                                                                                                               ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGC‑GGGGATGGGGGCGCCCCTGA‑AAGACTCCTTTTTTTTTATCGTGGGAAACATTTTCATTTTTATTTTTAGCTAATGCAATAGTTACTGAACTGATCC‑GGTGAGTTAATGTTGAACAAATCTCATGTTTCGGGGGGGGCGCTTTTACCACAGAAGCGGTTATGCCCGTGTGGttttttttat                                                                                               >  1:443877/1‑244 (MQ=255)
                                                                                                                                                   gcgTTAAAGACTGAATCTGTCCAGTCGCTGGTGCGGTTGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAttttt                                                                                           <  2:219159/251‑1 (MQ=255)
                                                                                                                                                    cgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTTATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGTATTTtta                                                                                          <  1:268093/251‑1 (MQ=255)
                                                                                                                                                             cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCGCGTTTTGC‑GGGGATGGGTGCGACCCTGA‑CAGTCTCCTTTTTTTTTATCGTGGGAAACATTTTCATTTTTTTTGTTTGCTAATGCAATATTTACTGAACTTCTCC‑GCTGAGGTTATTTTGGACAAATCTCATGTTTCGTGGGGGTCGCTTTTTCCACAGATGGGTTTATGCCAGGATGGTTTTTTTAGTTTTTTTTAAATCtt                                                                                 >  1:241294/1‑250 (MQ=255)
                                                                                                                                                                aaCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGt                                                                              <  1:696230/251‑1 (MQ=255)
                                                                                                                                                                     gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAAGCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGa                                                                           <  1:638132/249‑1 (MQ=255)
                                                                                                                                                                     gTCCAGGCGCTGGTGCGGTGGGAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGGCCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGGTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGa                                                                           <  1:743360/249‑1 (MQ=255)
                                                                                                                                                                             gCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTTCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCggg                                                                 <  2:683821/251‑1 (MQ=255)
                                                                                                                                                                             gCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCGCGTTTTCC‑GGGGATGGGGGCGGCCATGA‑CAGTCTCCTTTTTTTGTATCGTGGGAATTATTTTTATTTTTTTTGTTTGCTAATGCAATAGTTACTTAACTTATCC‑GATGGGTTAATGTTGAACAAATCTCATGTTGCGTGGGGGTCGCTTTTTCCACAAATGCGTTTATGCCCGTATGGGTTGTTTAATTTTTATTTAATCTGGGTTGAGGGTGTCggg                                                                 >  1:243893/1‑251 (MQ=255)
                                                                                                                                                                                      gTGGCAACGGCGTGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGAACATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGGTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAgtgt                                                        <  2:662729/251‑1 (MQ=255)
                                                                                                                                                                                             cggcggGTTTTCTGACGCGTTTCGC‑GGGAAGGGGTGCGACCATGA‑CAGTCTCCTTTTTCAGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGagcagc                                                 <  2:289095/251‑1 (MQ=255)
                                                                                                                                                                                              gtcggGTTTTCTGCCGCGTTTAGC‑GGGTATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATGTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCACGTGTGAgcagca                                                <  2:562390/249‑1 (MQ=255)
                                                                                                                                                                                                  ggTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                            <  1:800003/251‑1 (MQ=255)
                                                                                                                                                                                                          gCCGCGTGTCGC‑GGGAATGGGTGCGATCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaacaac                                    <  2:443877/251‑1 (MQ=255)
                                                                                                                                                                                                           ccGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACt                                   >  2:321812/1‑251 (MQ=255)
                                                                                                                                                                                                                     c‑gGGAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTgcgc                           >  1:639166/1‑249 (MQ=255)
                                                                                                                                                                                                                      ggggAATTGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGGTTATGCCAGTAGGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAg                         <  2:241294/250‑1 (MQ=255)
                                                                                                                                                                                                                        ggAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGtt                                                                                                       <  2:498776/171‑1 (MQ=255)
                                                                                                                                                                                                                        ggAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGtt                                                                                                       >  1:498776/1‑171 (MQ=255)
                                                                                                                                                                                                                        ggAATGGGGGCGAACATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGACAGGGTGAGCAGCAAAGGGGAACAACTGCGTgcgcgggt                       >  2:19971/1‑247 (MQ=255)
                                                                                                                                                                                                                         gAATGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGtt                       >  2:21862/1‑250 (MQ=255)
                                                                                                                                                                                                                          aaTGGGTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                      >  2:371107/1‑250 (MQ=255)
                                                                                                                                                                                                                               gTGCGACCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAGAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa                <  2:817702/251‑1 (MQ=255)
                                                                                                                                                                                                                                agCGACCATGA‑CAGTCTCATTATTCTGTAACGTGGAAAACATTGTCAATGTTATTGATAGCTAATGAAATAGGTACTGAACTGAGCC‑GAGGAGTTAATGTTGAACAAATCTCAGGTTGCGTGGTGGTAGAGTTTACCACAGAGGCGTTTATGCCAGTATGGTTTGTTGAATATTTATTAAATCTTGGTTGAGCGTGTGGGGATCAAGAGTGAGAAGCAAAGAGGAACAACTGCGTGCGCAGTTAAATGAAc               <  2:390957/250‑1 (MQ=255)
                                                                                                                                                                                                                                    aCCATGA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAt           <  1:266333/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gA‑CAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg      >  1:543543/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc                                                                                                                                 >  2:247506/1‑126 (MQ=255)
                                                                                                                                                                                                                                            cAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc                                                                                                                                 <  1:247506/126‑1 (MQ=255)
                                                                                                                                                                                                                                              gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGTTCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCt  >  2:758446/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
CGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGGGTGCGACCATGA‑CAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCC‑GATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCT  >  NC_000913/4223397‑4223883

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