Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,024 Δ1 bp coding (527/696 nt) araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,0240C.98.2% 245.8 / ‑3.7 54coding (527/696 nt)araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (34/19);  total (35/19)
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.99e-01

GCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCAC  >  NC_000913/65781‑66271
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
gctgtaattCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCaaat                                                                                                                                                                                                                                                 >  1:64194/9‑250 (MQ=255)
 ctgtaattCCATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCaaatg                                                                                                                                                                                                                                                >  2:323573/8‑249 (MQ=255)
            tCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGACCGTAATATGCCTTCGCGCCATGATTACGCAGATAGTGTTTATTCAGAAGCGATTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTgaga                                                                                                                                                                                                                                      >  1:24618/1‑247 (MQ=255)
                   aCCAGGCTTGATTATAGCCTGGTGTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCa                                                                                                                                                                                                                              <  1:16897/251‑1 (MQ=255)
                      aGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATTCAAATGGGCCGTGAGAATGGACCAgg                                                                                                                                                                                                                            >  2:265252/1‑249 (MQ=255)
                           tgattataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGc                                                                                                                                                                                                                        <  2:69513/245‑1 (MQ=255)
                             attataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCggg                                                                                                                                                                                                                    >  2:260021/3‑251 (MQ=255)
                                ataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGCAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCAt                                                                                                                                                                                                                 >  1:347545/1‑251 (MQ=255)
                                    ccTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTgc                                                                                                                                                                                                             <  2:64194/251‑1 (MQ=255)
                                               ttGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGAt                                                                                                                                                                                                  >  1:86959/1‑251 (MQ=255)
                                                           ggTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTc                                                                                                                                                                                      <  1:34048/251‑1 (MQ=255)
                                                                  taCAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGtt                                                                                                                                                                               >  2:57695/1‑251 (MQ=255)
                                                                          ttACTGCCCGTAATATGCCTTCGCGCCCTGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGAt                                                                                                                                                                       >  2:224861/1‑251 (MQ=255)
                                                                           tACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGAt                                                                                                                                                                       <  1:305233/250‑1 (MQ=255)
                                                                            aCTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGa                                                                                                                                                                     <  1:303328/251‑1 (MQ=255)
                                                                             cTGCCCGTAATCTGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGAc                                                                                                                                                                    >  1:142972/1‑251 (MQ=255)
                                                                             cTGCCCGTAATATGCCTTCGCGCCATGCTTCCGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACAGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGAc                                                                                                                                                                    >  2:287080/1‑251 (MQ=255)
                                                                                cccGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGtt                                                                                                                                                                 <  1:117054/251‑1 (MQ=255)
                                                                                cccGTAATACGCCTTCGCGCCATGCTTACGCAGATAGCGCTTATTCAGCAGCGTTTGCTGCATATACGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATCAGCGAACACTTCCAGCACGATGGCGCTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGAAGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGtt                                                                                                                                                                 >  1:241349/1‑251 (MQ=255)
                                                                                                  gcCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCt                                                                                                                                                                                             <  2:168696/205‑1 (MQ=255)
                                                                                                  gcCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCt                                                                                                                                                                                             >  1:168696/1‑205 (MQ=255)
                                                                                                       gCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTcac                                                                                                                                          >  1:219303/1‑249 (MQ=255)
                                                                                                       gCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATACGGTAACTGCGGCGCTACCTGAAGGCAGAATATCCACATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGCGAATGGACCAGGACGCCGGGAATTTGCGCTGCATCGATACCCTGTTTTACGAAGGTTTCTAAGATGACGTTACAGGTTTCCCACTCATATTcac                                                                                                                                          >  1:164567/1‑249 (MQ=255)
                                                                                                               aGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGAtt                                                                                                                                  >  1:232347/1‑251 (MQ=255)
                                                                                                                   aGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCt                                                                                                                               <  1:285842/250‑1 (MQ=255)
                                                                                                                                   cGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATa                                                                                                                                                                                                 <  2:49862/168‑1 (MQ=255)
                                                                                                                                   cGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATa                                                                                                                                                                                                 >  1:49862/1‑168 (MQ=255)
                                                                                                                                       tgctgcATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGt                                                                                                           <  1:95276/250‑1 (MQ=255)
                                                                                                                                         ctgcATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGc                                                                                                         >  1:123996/1‑250 (MQ=255)
                                                                                                                                                      aaCTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCg                                                                                           >  1:356052/1‑251 (MQ=255)
                                                                                                                                                           cggcgCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGTGGCCGTAGaa                                                                                      <  1:112971/251‑1 (MQ=255)
                                                                                                                                                           cggcgCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGa                                                                                       >  2:362134/1‑250 (MQ=255)
                                                                                                                                                            tgcgcTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGTATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGaa                                                                                      >  2:253222/2‑250 (MQ=255)
                                                                                                                                                               gcTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGAt                                                                                                                                   <  1:50592/202‑1 (MQ=255)
                                                                                                                                                               gcTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGAt                                                                                                                                   >  2:50592/1‑202 (MQ=255)
                                                                                                                                                                  aaCTGACGGCACAATATCCCCATATAAGCGACATCTTCCAGAACGATGCCGTTATGCACCGCATCTTCGGCATTTTTCCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTGGAAATAGTCg                                                                               >  2:209706/1‑251 (MQ=255)
                                                                                                                                                                      gACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTTCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCg                                                                            <  2:186925/250‑1 (MQ=255)
                                                                                                                                                                       aCGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCCTTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCTACTCATATTCACCGTTGATTTCGGCGTCGGTCATTGTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtg                                                                          <  2:177576/251‑1 (MQ=255)
                                                                                                                                                                         ggCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtggg                                                                        >  1:7504/1‑251 (MQ=255)
                                                                                                                                                                            aGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCgtgggtgg                                                                     >  2:169395/1‑251 (MQ=255)
                                                                                                                                                                             gAATATCCGCATATAAGCGGCCTATTCACGCAAGATGGCGTTATGCACCGCATCTTAGGCCTTTTTGCCCATGCAAATGGGCCGTGAGCATGGACTAGGACGACGGGCATTTGCGCTGCCTCGATACCCTGTTTTTCGAAGGTTTCTACGATACCGTTACCGGTTTCCAACTCATCTTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGCATGGTGCCGTAGAAATCGTCGGCgtgggaggg                                                                    >  2:321440/1‑247 (MQ=255)
                                                                                                                                                                               atatCCCCATATAAGCTACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCAGCATTTTTGCCGATGCAAATGGGCCGGGAGAATGGACCAGGACGCCGGGCATTTACGCTGCCTAGATACACTGTTTTTGGAAGGTTGCTACGATGACGTTACCGGATTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGTGATGGTTCCGTCGAAATCGTCGGCGTGGGGGGTGc                                                                  >  1:159285/1‑251 (MQ=255)
                                                                                                                                                                                   ccccATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGt                                                              >  1:208378/1‑251 (MQ=255)
                                                                                                                                                                                     ccATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCACATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTg                                                            >  2:9213/1‑251 (MQ=255)
                                                                                                                                                                                           aaGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGa                                                       >  1:113131/1‑250 (MQ=255)
                                                                                                                                                                                                        cAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGc                                         >  1:292438/1‑251 (MQ=255)
                                                                                                                                                                                                        cAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGCGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGCATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGc                                         >  1:235897/1‑251 (MQ=255)
                                                                                                                                                                                                                  gcggtaggCACCGCGTCTTCGGCATGTTTGCCCATGCAAATGGGCCGGGAGAAGGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAg                               <  2:274976/244‑1 (MQ=255)
                                                                                                                                                                                                                          cACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggc                       <  1:182558/251‑1 (MQ=255)
                                                                                                                                                                                                                           aCCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGAGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTAATTTTGCGGGTGCAGGGACTGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcg                      >  2:18870/1‑251 (MQ=255)
                                                                                                                                                                                                                            ccGCAGCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTGGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgt                     <  1:253222/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGtggcgtggcggc               >  2:356007/1‑249 (MQ=255)
                                                                                                                                                                                                                                      ggCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTg           >  1:124452/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTGCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGCGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCGGGAATCGACTTGCCCGCCTGCGCCCAGATAGTGGCTTGGCGCGAGTgt          <  2:24618/251‑2 (MQ=255)
                                                                                                                                                                                                                                                gcccATGCAAATCGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGCTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGtggcgtggcgggggtgtgggtgcac  >  2:25530/2‑236 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
GCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCAC  >  NC_000913/65781‑66271

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