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→AATfrdBfumarate 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

GGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTT  >  NC_000913/4380047‑4380533
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
ggcgTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAAtttgg                                                                                                                                                                                                                                              <  1:359818/251‑3 (MQ=255)
  ggtagcccgaaCCCCAGTTGGTGCCTCCGGAGAACGGGGGATCCTTCGCCATCGTGGCCGGGGTCTGGATGTTAGTACCCTGAGCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGGTTCCAGACTTTCGATGAAGTGGTTCATAGCGACCACCAGATCGCGTTCAATCTGGAAGTTAGCGAACGCTTCAACCGTCATTCCGTCGGTGTAATCACGCAGGAAGTTTTTCCATGCCAAtttggca                                                                                                                                                                                                                                            <  2:158123/242‑5 (MQ=255)
     aGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGATAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTTGCACgt                                                                                                                                                                                                                                        >  1:306050/1‑251 (MQ=255)
      gCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTAGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACgtt                                                                                                                                                                                                                                       >  1:83961/1‑251 (MQ=255)
        acaAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACgttgt                                                                                                                                                                                                                                     <  1:274501/251‑1 (MQ=255)
        acaAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACgttgt                                                                                                                                                                                                                                     >  1:201197/1‑251 (MQ=255)
                 aGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACcatc                                                                                                                                                                                                                             >  1:237240/1‑250 (MQ=255)
                  gTTGATGCATCCGGAGAACTGGTGATTCTCCGCCATCTGCGCCGGGGTCTTGATGTTAGTACCCTGATCCTCGGTGCGGGAGTTGCCGATGATGTTCGGTTTGATCGCTTCCAGACTTTCGATGAAGTGTGTCTTATCGACCACCAGATCTCGTTCAATCGTTTAGTTAGCTAACGCTTTAACCTTCATACCGTCGGTGTAATCACGCAGGAAGTTTTTACATGCCAATTTGGCACTTTGTTAACcatcat                                                                                                                                                                                                                           <  1:22398/251‑1 (MQ=255)
                    tGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATg                                                                                                                                                                                                                          <  1:260920/250‑1 (MQ=255)
                                      ggTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCg                                                                                                                                                                                                       >  1:407267/1‑251 (MQ=255)
                                                          cgcCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCg                                                                                                                                                                                    >  1:506859/1‑251 (MQ=255)
                                                            ccGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGt                                                                                                                                                                                  <  1:371929/251‑1 (MQ=255)
                                                                ggTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGACGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGCCCAGCGGTAGcg                                                                                                                                                                              >  1:448899/1‑250 (MQ=255)
                                                                       atgatAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACAATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCCCAGGTcc                                                                                                                                                                       >  2:384953/5‑251 (MQ=255)
                                                                         gTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCACCCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACCGCAGGACCAGCGGTAGCTCAGGTCCgg                                                                                                                                                                     >  2:275434/1‑251 (MQ=255)
                                                                         gTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGCAGCTCAGGTCCgg                                                                                                                                                                     >  1:426120/1‑251 (MQ=255)
                                                                          ttaggccccTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCgg                                                                                                                                                                     <  2:83961/244‑1 (MQ=255)
                                                                            aGTACCCTGATCCGTGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGc                                                                                                                                                                  <  1:176235/251‑1 (MQ=255)
                                                                            aGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTACGGTGc                                                                                                                                                                  >  1:202726/1‑251 (MQ=255)
                                                                                   tGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTg                                                                                                                                                           >  1:365909/1‑251 (MQ=255)
                                                                                     aTCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGGCCGGGGCCAGGTTGTc                                                                                                                                                         >  1:441048/1‑251 (MQ=255)
                                                                                        cgcgGTGCGGGAGTTGCTGGTGTTGTACTGTTTGATCGCGTCCAGACGTTCGAGGAAGTTGGGCATTTCGACCTCCAGATCGGGTTCAATCGGGAAGTGAGCGAACGCTTCAACCTTCATACCGTCGGTGTAATCGCGCGGGAAGGTTTTACATGCCAATTTGGCGCGTTGTTAACCGTGTTGCCGCGGGATTCGGAGATCGTCATAGGGCAGGAGCAGCGGGAGCTCAGGTCGGGTGCCATGTTGTCttt                                                                                                                                                      <  2:541217/251‑1 (MQ=255)
                                                                                              gCGGGAGTTGCCGATGATGTTCGGGTGGATCGCTTCCAGCCTTTCGATGAAGGGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGTCAGGACCATCGGTAGCTCAGGTCCGGTGCCAGGGTGTCTTTTATGTa                                                                                                                                                <  2:224531/251‑1 (MQ=255)
                                                                                                   aGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCa                                                                                                                                           <  1:67185/251‑1 (MQ=255)
                                                                                                     ttGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCCGCGGTAGCGCAGGGCCGGGGCCAGGGTGTGTTTGATGTAGCCCAgc                                                                                                                                         >  2:109902/1‑251 (MQ=255)
                                                                                                         cGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCAt                                                                                                                                     <  1:168877/251‑1 (MQ=255)
                                                                                                          gatgatGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACTCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATc                                                                                                                                    >  1:364331/1‑251 (MQ=255)
                                                                                                             gatggtcGGTGTGATCGCTTCCAGACTTTCGATGAAGGGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGGTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCa                                                                                                                                  <  2:399498/244‑1 (MQ=255)
                                                                                                                         gATCGCTTCCAGACGTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGt                                                                                                                     <  1:297162/251‑1 (MQ=255)
                                                                                                                         ccctccacgacgctctgccgagctgTGGGTCATACCGACCCCCAGATCCCGTTCAACCTGGAAGTTAGCTAACGCTTCAACCCTCATACCGTCGGTTTAATCCCGCAGGAAGGTTTTACATGCCAGTTTTTGCACGTTGTTAACCAGCATGCCGCAGGAACCACAAATCGCCATTCGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAAGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTa                                                                                                                       <  2:283215/226‑1 (MQ=255)
                                                                                                                              cTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGt                                                                                                                >  1:365388/1‑251 (MQ=255)
                                                                                                                               ttCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGCTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCCGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTc                                                                                                               >  1:422332/1‑251 (MQ=255)
                                                                                                                                tCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTc                                                                                                               <  1:541043/250‑1 (MQ=255)
                                                                                                                                                 gTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCACAGCGCATCCAGTAATGAGGTaaggtcggaagagaacaagtctgaac                                                                                              >  1:283215/1‑226 (MQ=255)
                                                                                                                                                     gTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGATCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGGAGTTGCGTCATAAGGCCCTTAATAGAATgcg                                                                                         >  1:372993/1‑251 (MQ=255)
                                                                                                                                                       cATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACAGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCt                                                                                       <  1:362826/251‑1 (MQ=255)
                                                                                                                                                             gACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTAtgcgg                                                                                 <  2:24493/251‑1 (MQ=255)
                                                                                                                                                             gACCACCAGATCGCGTTCAATCGGGAAGTTAGCCAACGCTTCAACCTTCATACCGACGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGt                                                                                                                                                 <  1:139849/188‑1 (MQ=255)
                                                                                                                                                               cccccaGATCGCGTTCAATCGGGAAGGTAGTGAACGCTTCAACCTGCATACCGTCGGTGTAAGCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGGATCCAGTAATGAGGTAGTTGCGGCATAAGGCACTTCATAGAATGCGCTAtgcggtg                                                                               <  2:509917/248‑1 (MQ=255)
                                                                                                                                                               ccaccaGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGGTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCCTAGAATGCGCTAtgcggtg                                                                               >  1:442533/1‑251 (MQ=255)
                                                                                                                                                                caccaGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTAtgcggtgc                                                                              <  1:170458/251‑1 (MQ=255)
                                                                                                                                                                    aGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCCCAGGGCCGGTGCCAGGCTGTCTTTGATGTAGCCCAGCGCCTCCAGTACTGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATgcggtgcggt                                                                           >  1:499819/1‑250 (MQ=255)
                                                                                                                                                                      aTCGCGTTCAATCGGGAAGTTAGCGAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGTCACTTCATAGAATGCGCTATGCGGTGCGGTATc                                                                        <  1:463706/251‑1 (MQ=255)
                                                                                                                                                                            ttCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTc                                                                  >  1:201102/1‑251 (MQ=255)
                                                                                                                                                                              cAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCg                                                                >  1:501839/1‑251 (MQ=255)
                                                                                                                                                                                       aGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGACAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGAGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGCCTTCCGGGTTATAgc                                                       >  1:100528/1‑251 (MQ=255)
                                                                                                                                                                                        gTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAgcg                                                      <  1:366259/251‑1 (MQ=255)
                                                                                                                                                                                         ttAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGATGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGAACCAGCGGAAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGGAGTTGCGTCATAAGGACCTTCATAGAATGCGCTATGCGGTGCGGTATCGTCTTACGGGTtatatccc                                                     >  2:136375/1‑247 (MQ=255)
                                                                                                                                                                                              tAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGcaccac                                                >  1:391929/1‑251 (MQ=255)
                                                                                                                                                                                                 cccttaaaCCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTc                                             <  2:288387/245‑1 (MQ=255)
                                                                                                                                                                                                  gCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCa                                            <  1:437816/251‑1 (MQ=255)
                                                                                                                                                                                                      cAACCTTCATACCGTCGGAGTACTCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGATCAGGGCCGGTGACAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGCGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGAATTCCGGGTTATAGCGCACCACCTCAAttt                                        >  1:431649/1‑251 (MQ=255)
                                                                                                                                                                                                            tCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAg                                    <  1:163129/249‑1 (MQ=255)
                                                                                                                                                                                                                aCCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGtttt                               <  1:423598/250‑1 (MQ=255)
                                                                                                                                                                                                                  cGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTc                             >  1:285775/1‑250 (MQ=255)
                                                                                                                                                                                                                   gTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCAt                           <  1:396080/251‑1 (MQ=255)
                                                                                                                                                                                                                     cGGTTTAATCCCGCAGGAATGTTTTACATGCCAATTTGGCACGTTGTTCACCATCATGGCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGACCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCTCATCCAGTAATGAGGTAGGTGCGTCATAAGACACTTCATAGAATGCGCTATGCGGTGCGGTTTCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGCTTTTCAtc                          <  2:499819/250‑1 (MQ=255)
                                                                                                                                                                                                                      ggTGTAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCATGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCTCTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCAtctc                        <  1:54366/251‑1 (MQ=255)
                                                                                                                                                                                                                        tgtAATCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACAATCATGCCGCAGGAACCAAAAATCGCCATACGGCAGGGCCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGGAAAGAGGTAGGTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTACGGGTTCTAGCGCACCACACCAATTTTCAGGTTTTTCAttcc                        >  2:14324/1‑246 (MQ=255)
                                                                                                                                                                                                                         ttAATCACGCAGGAAGGTTTTACATTCCAAGTTGGCACGGTGTTAACCATCAGGCCGCAGGAACCACCAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCAGCCAGTAATGAGGTAGTGGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGTCTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAg                      <  2:426120/249‑1 (MQ=255)
                                                                                                                                                                                                                          taataaCGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAAACATCATGCCGCAGGAAGCACAAATCGCAATACGGCAGGACCAGCGGTAGCTCAGGTCCGGAGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCCGTGATGAGGTAGTTGCGTCATAAGGCACTTAATAGAATGAGATATGCGGTGCGGGATCGACTGCCGGGTTATAGCGAACCACCTCAATTTTCAGGTTTTTCGTCTCAGcc                    >  2:15971/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaTCACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGcc                    >  1:150885/1‑250 (MQ=255)
                                                                                                                                                                                                                              cACGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTc                >  1:256060/1‑251 (MQ=255)
                                                                                                                                                                                                                               aCGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTGTTCATCTCAGCCATTCg               <  2:92850/251‑1 (MQ=255)
                                                                                                                                                                                                                                cGCAGGAAGGTTTTACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGc              >  1:244548/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gCAGGAAGGTTTTACATTCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGACCCACAAATCGCCATACGGCAGGACCAGCGGTAAGTCAGGTCCGGGGCCAGGTGGTCTTTGATGTAGCCCAGCGCATACAGTAATGAGGTAGTTGCGTCATACGGAACTTCATAGAATGCGCATTGCGGTGCGGTAACGGCTTTCGGGTTATCACACACCACCTCAATTTTCAGGTTTTTCCTCTCCGCCATTCGcc             >  2:430507/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cAGGAATGGTTTACATGCCAATTTTGCACGGTATTAACCAGCATGCCGCAGGAACCACAATTCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGTGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTGTCAGGGTTTTCATCTCATCCATTCGcct            <  2:443774/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCcttctt  <  2:36383/250‑1 (MQ=255)
                                                                                                                                                                                                                                             tACATGCCAATTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCcttctt  >  1:432082/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
GGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGATGTACGGTTTGATCGCTTCCAGACTTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTT  >  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: