Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,380,292 2 bp→A coding (225‑226/735 nt) frdB ← fumarate reductase iron‑sulfur protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,380,2920G.93.9% 255.2 / 14.3 69coding (226/735 nt)frdBfumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  major base . (33/29);  minor base A (1/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01
*NC_0009134,380,2930TA92.5% 190.5 / 5.4 67K75N (AAA→AAT) frdBfumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base A (33/29);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATC‑TGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTT  >  NC_000913/4380047‑4380533
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
ggcgTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATC‑TGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAAtttgg                                                                                                                                                                                                                                              <  1:359818/251‑3 (MQ=255)
  ggtagcccgaaCCCCAGTTGGTGCCTCCGGAGAACGGGGGATCCTTCGCCATCGTG‑GCCGGGGTCTGGATGTTAGTACCCTGAGCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGGTTCCAGACTTTCGATGAAGTGGTTCATAGCGACCACCAGATCGCGTTCAATCTGGAAGTTAGCGAACGCTTCAACCGTCATTCCGTCGGTGTAATCACGCAGGAAGTTTTTCCATGCCAAtttggca                                                                                                                                                                                                                                            <  2:158123/242‑5 (MQ=255)
     aGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATC‑TGCGCCGGGGTCTGGATGATAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTTG‑CACgt                                                                                                                                                                                                                                        >  1:306050/1‑251 (MQ=255)
      gCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATC‑TGCGCCGGGGTCTAGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACgtt                                                                                                                                                                                                                                       >  1:83961/1‑251 (MQ=255)
        acaAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATC‑TGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACgttgt                                                                                                                                                                                                                                     <  1:274501/251‑1 (MQ=255)
        acaAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATC‑TGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACgttgt                                                                                                                                                                                                                                     >  1:201197/1‑251 (MQ=255)
                 aGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATC‑TGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACcatc                                                                                                                                                                                                                             >  1:237240/1‑250 (MQ=255)
                  gTTGATGCATCCGGAGAACTGGTGATTCTCCGCCATC‑TGCGCCGGGGTCTTGATGTTAGTACCCTGATCCTCGGTGCGGGAGTTGCCGATGATGTTCGGTTTGATCGCTTCCAGACTTTCGATGAAGTGTGTCTTATCGACCACCAGATCTCGTTCAATCGTTTAGTTAGCTAACGCTTTAACCTTCATACCGTCGGTGTAATCACGCAGGAAGTTTTTACATGCCA‑ATTTGGCACTTTGTTAACcatcat                                                                                                                                                                                                                           <  1:22398/251‑1 (MQ=255)
                    tGATGCAACCGGAGAACTGGTGATACTTCGCCATC‑TGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATg                                                                                                                                                                                                                          <  1:260920/250‑1 (MQ=255)
                                      ggTGATACTTCGCCATC‑TGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCg                                                                                                                                                                                                       >  1:407267/1‑251 (MQ=255)
                                                          cgcCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCg                                                                                                                                                                                    >  1:506859/1‑251 (MQ=255)
                                                            ccGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGt                                                                                                                                                                                  <  1:371929/251‑1 (MQ=255)
                                                                ggTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGACGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGCCCAGCGGTAGcg                                                                                                                                                                              >  1:448899/1‑250 (MQ=255)
                                                                       atgatAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACAATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCCCAGGTcc                                                                                                                                                                       >  2:384953/5‑251 (MQ=255)
                                                                         gTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCACCCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACCGCAGGACCAGCGGTAGCTCAGGTCCgg                                                                                                                                                                     >  2:275434/1‑251 (MQ=255)
                                                                         gTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGCAGCTCAGGTCCgg                                                                                                                                                                     >  1:426120/1‑251 (MQ=255)
                                                                          ttaggccccTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCgg                                                                                                                                                                     <  2:83961/244‑1 (MQ=255)
                                                                            aGTACCCTGATCCGTGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGc                                                                                                                                                                  <  1:176235/251‑1 (MQ=255)
                                                                            aGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTACGGTGc                                                                                                                                                                  >  1:202726/1‑251 (MQ=255)
                                                                                   tGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTg                                                                                                                                                           >  1:365909/1‑251 (MQ=255)
                                                                                     aTCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGGCCGGGGCCAGGTTGTc                                                                                                                                                         >  1:441048/1‑251 (MQ=255)
                                                                                        cgcgGTGCGGGAGTTGCTGGTGTTGTACTGTTTGATCGCGTCCAGACGTTCGAGGAAGTTGGGCATTTCGACCTCCAGATCGGGTTCAATCGGGAAGTGAGCGAACGCTTCAACCTTCATACCGTCGGTGTAATCGCGCGGGAAGGTTTTACATGCCA‑ATTTGGCGCGTTGTTAACCGTGTTGCCGCGGGATTCGGAGATCGTCATAGGGCAGGAGCAGCGGGAGCTCAGGTCGGGTGCCATGTTGTCttt                                                                                                                                                      <  2:541217/251‑1 (MQ=255)
                                                                                              gCGGGAGTTGCCGATGATGTTCGGGTGGATCGCTTCCAGCCTTTCGATGAAGGGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGTCAGGACCATCGGTAGCTCAGGTCCGGTGCCAGGGTGTCTTTTATGTa                                                                                                                                                <  2:224531/251‑1 (MQ=255)
                                                                                                   aGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCa                                                                                                                                           <  1:67185/251‑1 (MQ=255)
                                                                                                     ttGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCCGCGGTAGCGCAGGGCCGGGGCCAGGGTGTGTTTGATGTAGCCCAgc                                                                                                                                         >  2:109902/1‑251 (MQ=255)
                                                                                                         cGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCAt                                                                                                                                     <  1:168877/251‑1 (MQ=255)
                                                                                                          gatgatGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACTCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATc                                                                                                                                    >  1:364331/1‑251 (MQ=255)
                                                                                                             gatggtcGGTGTGATCGCTTCCAGACTTTCGATGAAGGGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGGTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCa                                                                                                                                  <  2:399498/244‑1 (MQ=255)
                                                                                                                         gATCGCTTCCAGACGTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGt                                                                                                                     <  1:297162/251‑1 (MQ=255)
                                                                                                                         ccctccacgacgctctgccgagctgTGGGTCATACCGACCCCCAGATCCCGTTCAACCTGGAAGTTAGCTAACGCTTCAACCCTCATACCGTCGGTTTAATCCCGCAGGAAGGTTTTACATGCCAGTTTTTGCACGTTGTTAACCAGCATGCCGCAGGAACCACAAATCGCCATTCGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAAGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTa                                                                                                                       <  2:283215/226‑1 (MQ=255)
                                                                                                                              cTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGt                                                                                                                >  1:365388/1‑251 (MQ=255)
                                                                                                                               ttCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGCTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCCGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTc                                                                                                               >  1:422332/1‑251 (MQ=255)
                                                                                                                                tCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTc                                                                                                               <  1:541043/250‑1 (MQ=255)
                                                                                                                                                 gTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCACAGCGCATCCAGTAATGAGGTaaggtcggaagagaacaagtctgaac                                                                                              >  1:283215/1‑226 (MQ=255)
                                                                                                                                                     gTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGATCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGGAGTTGCGTCATAAGGCCCTTAATAGAATgcg                                                                                         >  1:372993/1‑251 (MQ=255)
                                                                                                                                                       cATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACAGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCt                                                                                       <  1:362826/251‑1 (MQ=255)
                                                                                                                                                             gACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTAtgcgg                                                                                 <  2:24493/251‑1 (MQ=255)
                                                                                                                                                             gACCACCAGATCGCGTTCAATCGGGAAGTTAGCCAACGCTTCAACCTTCATACCGACGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGt                                                                                                                                                 <  1:139849/188‑1 (MQ=255)
                                                                                                                                                               cccccaGATCGCGTTCAATCGGGAAGGTAGTGAACGCTTCAACCTGCATACCGTCGGTGTAAGCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGGATCCAGTAATGAGGTAGTTGCGGCATAAGGCACTTCATAGAATGCGCTAtgcggtg                                                                               <  2:509917/248‑1 (MQ=255)
                                                                                                                                                               ccaccaGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGGTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCCTAGAATGCGCTAtgcggtg                                                                               >  1:442533/1‑251 (MQ=255)
                                                                                                                                                                caccaGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTAtgcggtgc                                                                              <  1:170458/251‑1 (MQ=255)
                                                                                                                                                                    aGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCCCAGGGCCGGTGCCAGGCTGTCTTTGATGTAGCCCAGCGCCTCCAGTACTGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATgcggtgcggt                                                                           >  1:499819/1‑250 (MQ=255)
                                                                                                                                                                      aTCGCGTTCAATCGGGAAGTTAGCGAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGTCACTTCATAGAATGCGCTATGCGGTGCGGTATc                                                                        <  1:463706/251‑1 (MQ=255)
                                                                                                                                                                            ttCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTc                                                                  >  1:201102/1‑251 (MQ=255)
                                                                                                                                                                              cAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCg                                                                >  1:501839/1‑251 (MQ=255)
                                                                                                                                                                                       aGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGACA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGAGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGCCTTCCGGGTTATAgc                                                       >  1:100528/1‑251 (MQ=255)
                                                                                                                                                                                        gTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAgcg                                                      <  1:366259/251‑1 (MQ=255)
                                                                                                                                                                                         ttAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGATGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGAACCAGCGGAAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGGAGTTGCGTCATAAGGACCTTCATAGAATGCGCTATGCGGTGCGGTATCGTCTTACGGGTtatatccc                                                     >  2:136375/1‑247 (MQ=255)
                                                                                                                                                                                              tAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGcaccac                                                >  1:391929/1‑251 (MQ=255)
                                                                                                                                                                                                 cccttaaaCCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTc                                             <  2:288387/245‑1 (MQ=255)
                                                                                                                                                                                                  gCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCa                                            <  1:437816/251‑1 (MQ=255)
                                                                                                                                                                                                      cAACCTTCATACCGTCGGAGTACTCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGATCAGGGCCGGTGACAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGCGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGAATTCCGGGTTATAGCGCACCACCTCAAttt                                        >  1:431649/1‑251 (MQ=255)
                                                                                                                                                                                                            tCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAg                                    <  1:163129/249‑1 (MQ=255)
                                                                                                                                                                                                                aCCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGtttt                               <  1:423598/250‑1 (MQ=255)
                                                                                                                                                                                                                  cGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTc                             >  1:285775/1‑250 (MQ=255)
                                                                                                                                                                                                                   gTCGGTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCAt                           <  1:396080/251‑1 (MQ=255)
                                                                                                                                                                                                                     cGGTTTAATCCCGCAGGAATGTTTTACATGCCA‑ATTTGGCACGTTGTTCACCATCATGGCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGACCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCTCATCCAGTAATGAGGTAGGTGCGTCATAAGACACTTCATAGAATGCGCTATGCGGTGCGGTTTCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGCTTTTCAtc                          <  2:499819/250‑1 (MQ=255)
                                                                                                                                                                                                                      ggTGTAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCATGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCTCTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCAtctc                        <  1:54366/251‑1 (MQ=255)
                                                                                                                                                                                                                        tgtAATCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACAATCATGCCGCAGGAACCAAAAATCGCCATACGGCAGGGCCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGGAAAGAGGTAGGTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTACGGGTTCTAGCGCACCACACCAATTTTCAGGTTTTTCAttcc                        >  2:14324/1‑246 (MQ=255)
                                                                                                                                                                                                                         ttAATCACGCAGGAAGGTTTTACATTCCA‑AGTTGGCACGGTGTTAACCATCAGGCCGCAGGAACCACCAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCAGCCAGTAATGAGGTAGTGGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGTCTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAg                      <  2:426120/249‑1 (MQ=255)
                                                                                                                                                                                                                          taataaCGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAAACATCATGCCGCAGGAAGCACAAATCGCAATACGGCAGGACCAGCGGTAGCTCAGGTCCGGAGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCCGTGATGAGGTAGTTGCGTCATAAGGCACTTAATAGAATGAGATATGCGGTGCGGGATCGACTGCCGGGTTATAGCGAACCACCTCAATTTTCAGGTTTTTCGTCTCAGcc                    >  2:15971/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaTCACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGcc                    >  1:150885/1‑250 (MQ=255)
                                                                                                                                                                                                                              cACGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTc                >  1:256060/1‑251 (MQ=255)
                                                                                                                                                                                                                               aCGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTGTTCATCTCAGCCATTCg               <  2:92850/251‑1 (MQ=255)
                                                                                                                                                                                                                                cGCAGGAAGGTTTTACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGc              >  1:244548/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gCAGGAAGGTTTTACATTCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGACCCACAAATCGCCATACGGCAGGACCAGCGGTAAGTCAGGTCCGGGGCCAGGTGGTCTTTGATGTAGCCCAGCGCATACAGTAATGAGGTAGTTGCGTCATACGGAACTTCATAGAATGCGCATTGCGGTGCGGTAACGGCTTTCGGGTTATCACACACCACCTCAATTTTCAGGTTTTTCCTCTCCGCCATTCGcc             >  2:430507/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cAGGAATGGTTTACATGCCAATTTTG‑CACGGTATTAACCAGCATGCCGCAGGAACCACAATTCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGTGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTGTCAGGGTTTTCATCTCATCCATTCGcct            <  2:443774/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCcttctt  <  2:36383/250‑1 (MQ=255)
                                                                                                                                                                                                                                             tACATGCCA‑ATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCcttctt  >  1:432082/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
GGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATC‑TGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTT  >  NC_000913/4380047‑4380533

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