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.3% 114.5 / NA 37intergenic (‑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 (9/27);  minor base A (0/1);  total (9/28)
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.63e-01

CAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTG  >  NC_000913/4223392‑4223881
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
cAGTGCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGTTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCctttt                                                                                                                                                                                                                                                    <  2:344575/251‑1 (MQ=255)
           ccTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTg                                                                                                                                                                                                                                         <  2:400541/251‑1 (MQ=255)
                 gCACGAAACCCTGCTGTGGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCGTGTTGCGCCAGTTACGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAt‑c                                                                                                                                                                                                                                  <  2:310824/251‑1 (MQ=255)
                         accTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCATACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCCCCTTTTGCGGGGAGAGGTGCGGACACGACAGGCTCCTTTTTTTGTATCGTTGaacatatttttct                                                                                                                                                                                                                           >  2:230125/2‑240 (MQ=255)
                                gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAG‑CA‑TTTGCATTTTTAtt                                                                                                                                                                                                                  <  2:183028/251‑1 (MQ=255)
                                       cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGc                                                                                                                                                                                                            <  1:288851/250‑1 (MQ=255)
                                             tAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGAGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TGTTCATTTTTATTGGTAGCTAATGc                                                                                                                                                                                                      <  2:183760/250‑1 (MQ=255)
                                                         ggttggtggaTTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGTCATTGGATTCGGCAATCCACTCCAGTAATGTCGGGCCACGCTTTAAAGGCTGAACCTGTTCTGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGTCGCGTTTCGCGGGAAAGGGTGCGACTATGACAGTCGCATTTTTCTGTATCGTTGAAGT‑CA‑TTTTCGTTTTTATCGTTAGCTTTTGGAATAGTTaatgaa                                                                                                                                                                                         <  2:265119/242‑1 (MQ=255)
                                                         ggTCGTGGAAGTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaa                                                                                                                                                                                         <  1:6155/251‑1 (MQ=255)
                                                            cGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg                                                                                                                                                                                      <  1:282957/251‑1 (MQ=255)
                                                                    tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGAt                                                                                                                                                                              <  1:141060/250‑1 (MQ=255)
                                                                      ggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAAGAGTTACTGAACTG‑ATCCGATGAg                                                                                                                                                                           <  1:81851/251‑1 (MQ=255)
                                                                       gTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCCTGACAGTCTCCTTTTTTTGTATCGTGGGAAT‑CA‑TTTTTATTTTTTTTGTTAGCTAATGCAATAGTTACTTAAATG‑ATCCGATGAgg                                                                                                                                                                          >  1:113836/1‑250 (MQ=255)
                                                                           cccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAAt                                                                                                                                                                      <  2:359629/251‑1 (MQ=255)
                                                                                  tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACTCGTTAAAGACTGAACCTGGCCAGTCGCTGTTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGaa                                                                                                                                                                <  1:216413/250‑1 (MQ=255)
                                                                                          cAGTTCCGTGAGTGCCCCACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCa                                                                                                                                                       <  2:113836/251‑1 (MQ=255)
                                                                                          cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTAAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCTCGGGGATGGGGGCGGCCATCACAGTCCCCTTTTTCTTTTTCGTTGAAAA‑CC‑TTTTTCTTTTTTTTTTTAGCTACAGCAATTGTTAATGAACTGTATC‑GATGAGTTaaatttttaaaaatatcc                                                                                                                                                       >  2:331160/1‑235 (MQ=255)
                                                                                                               gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGTGCGACCATGTAAGTCTCCTTTTTTTGTATCGTGGAAAT‑CA‑TTTTTATTTTTTTTTTTTGATAATTCAATAGTTTCTGGACTG‑ATCCGATGTGGTAATTTTGAACAAATTTTCTGTTGCGTGGTGGGCGCtttt                                                                                                                                  >  1:325861/1‑251 (MQ=255)
                                                                                                                 cATTGTATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑AGCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAc                                                                                                                                <  1:133234/251‑1 (MQ=255)
                                                                                                                  aTTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc                                                                                                                               <  2:98917/251‑1 (MQ=255)
                                                                                                                              aaTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCCCGGTTCGCGGGGATGGGGGCGGCCCTGACAGTCTCCTTTTTTTGTATTGGGGGAAT‑CAATTTTC‑TTTTTTTTTTTTGCTAATGCAATAGTTACTGAAATG‑TTCCGGTGGGGTTATGTTTAACAAATCTCATGTTGCGTGGGGGTCGCTTTTTCCACAGATGCggtt                                                                                                                   >  2:168989/1‑248 (MQ=255)
                                                                                                                                          gtaTTTCAGGCGACGCGGTAAGTACTGAACCTGTCCAGGCGCAGGAGCGGTGGCAACGGCGGGTTGTCGTCCGCGGTTCGCTGTAATGGGTGCGACCATGACAGTCTCCTATTTCTGTATCTTGGAAAT‑CA‑TTTTCATTTGTAATGTTAGCAAATGCAATAGTTACTGATCTG‑ATCCGATGAGTTAATGTTGAACAAATCGCATGTTGCGTCGTGGTCGCTGTTACCACAGATGCGTTCATGCCTGTAttg                                                                                                       <  2:272022/249‑3 (MQ=255)
                                                                                                                                           aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGt                                                                                                      <  2:277931/251‑1 (MQ=255)
                                                                                                                                                cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTTGCGGGGATGGGTGGGGCCCTGACAGGCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTTTTGTTAGCCAATGCAATAGTTAATTAACTG‑ATCCGCTGAGTTAATGTTTAACAAATCTCATGTTTCGTGGGGGTCGCTTTTTCCACAGATGCGttttttcccggttgggttttt                                                                                                 >  2:332724/1‑233 (MQ=255)
                                                                                                                                                 aGGCCACGCGTTAAAGACTAAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgttg                                                                                                <  1:331160/251‑1 (MQ=255)
                                                                                                                                                    ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAt                                                                                             <  1:109533/251‑1 (MQ=255)
                                                                                                                                                    ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGCGCGGTGGCAACGGCGGGATTTTTGCCCCGTTTCCCGGGGATGGGGGGGGCCATGACAATCTCCTTTTTCTGTATCGTGGGAATACT‑TTTTC‑TTTTTTTTTTTAGCTAATGCCATTATTACTTAACTT‑ATACAGTGAGTTAATTTTTTACAAATCTCCTGTTGCGTGGTGGGCGCCTTTACCCCAGATGCCGTTATTCCCGTATGGGTTGTTGatt                                                                                             >  2:59271/1‑249 (MQ=255)
                                                                                                                                                             aaaGACTGAACCTGGCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCTTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttatta                                                                                     <  1:269341/250‑1 (MQ=255)
                                                                                                                                                              aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTaaa                                                                                   <  1:168989/251‑1 (MQ=255)
                                                                                                                                                               aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAAt                                                                                  <  1:313273/251‑1 (MQ=255)
                                                                                                                                                                   tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGATTCGCGGGAAAGGGGGTGACCATGACAGTCTCCTTTTTTTGTATCGGGGGGAT‑AC‑TTTTCATTTTTTTTGTTTGCTAATGCAATAGTTAATTAACTG‑ATTCGATTGGGTAATGTTTAACAAATTTCATTTTGCGTGGGGGTCGCTTTTACCACAGAGGAGTTTATTCCAGTATTGTTTATTGAATTTTTTTTTAATCTgg                                                                              >  1:371866/1‑251 (MQ=255)
                                                                                                                                                                                                            tCTGCCGCGTGTCGCGGGAATGGGTGCGACCATTACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGTTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAAGCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaa                                      <  2:347949/250‑1 (MQ=255)
                                                                                                                                                                                                               gCCGCGTTTCGCGGGAATGGGTTCGACCATGACAGTCTCCTTGTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTTCCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaacaac                                  <  1:78807/251‑1 (MQ=255)
                                                                                                                                                                                                                 cgcgTTTCGCGGGAATTGGTGTGACCATGCCAGTCTCCTTTTGCTGTTTCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCATATGTTGCGTGGTGGTCGCTTTTACCTCAGATGCGTTTATGCCAGTATGGTTTGTTTAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACATCtg                                <  2:75765/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGACCAACtgcgt                             >  2:344933/1‑251 (MQ=255)
                                                                                                                                                                                                                             gAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                    <  1:332724/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aCAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt   <  2:216659/251‑1 (MQ=255)
                                                                                                                                                                                                                                               cAGTCTCCTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTg  <  1:304438/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAT‑CA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTG‑ATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTG  >  NC_000913/4223392‑4223881

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