Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,092 A→G F153F (TTT→TTC 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,0920AG98.2% 178.7 / ‑4.2 57F153F (TTT→TTCaraDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (32/24);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

ATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCT  >  NC_000913/65846‑66336
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ataCAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGt                                                                                                                                                                                                                                                  <  1:356135/251‑1 (MQ=255)
    aGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCt                                                                                                                                                                                                                                              <  1:12742/251‑1 (MQ=255)
              gCCCGTAATATGCCTTCGCGCCATGCTTACGCAGCTAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCCCGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGt                                                                                                                                                                                                                                    >  2:438626/1‑251 (MQ=255)
                  gTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTAcc                                                                                                                                                                                                                                >  1:310626/1‑251 (MQ=255)
                     atatGCCTTCGCGCCATGCTTACGCAGATAGAGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATAGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGt                                                                                                                                                                                                                             >  1:412876/1‑251 (MQ=255)
                                 gcaaTGCTTACGCAGATAGTGTTGATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGAGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCata                                                                                                                                                                                                                 <  1:150033/248‑1 (MQ=255)
                                 gcCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCata                                                                                                                                                                                                                 <  1:64720/251‑1 (MQ=255)
                                     tGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTc                                                                                                                                                                                                              >  1:199514/1‑250 (MQ=255)
                                      gCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTca                                                                                                                                                                                                             <  1:5153/250‑2 (MQ=255)
                                      gCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGAATATCCGTAAACTGCGGCGCTAACTGACGGCAGAATACCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTGCGGCATTTTTGCCCATGCAAATGGGCCGTTAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTca                                                                                                                                                                                                             <  2:358278/250‑2 (MQ=255)
                                       cTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTcacc                                                                                                                                                                                                           >  1:235404/1‑248 (MQ=255)
                                        ttACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCg                                                                                                                                                                                                          >  2:38959/1‑251 (MQ=255)
                                          aCGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCAGAACCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGAGACCTCTTCCAGCACGATGGCGTTATGCACCGAATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTAGAAGGTTTCTACGATGACGTTACCGGTTTCCCGCTCATATTcacagtt                                                                                                                                                                                                        >  2:41979/1‑245 (MQ=255)
                                              aGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGAtt                                                                                                                                                                                                    <  1:428399/251‑1 (MQ=255)
                                                   gtgtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCCTCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGc                                                                                                                                                                                               >  2:124492/1‑251 (MQ=255)
                                                    tgtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGTTACCCTGTTTTTCGCAGGTTTCTACGATTGCGTTACCGGTTTCCCACTCCTATTCACCTTTGATTTCTGcc                                                                                                                                                                                              >  2:412342/1‑250 (MQ=255)
                                                    tgtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCAGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCg                                                                                                                                                                                              >  1:129173/1‑251 (MQ=255)
                                                         attcagcagCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGTCAGAATTTCACCATATAAGCGACCTCTTCCACCACGATGGCGTTATGCACCGCATCTTCGACCTTTTTGCCCATGCAAATGGGCCGCGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGCCGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGt                                                                                                                                                                                         >  1:263057/4‑251 (MQ=255)
                                                           tcagcagCGTTTGCTGAATATCCGGTAACTGCGGCGCTAACTGACGGCAGACTATCCACATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACAGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTc                                                                                                                                                                                        >  1:297911/2‑250 (MQ=255)
                                                                       gctgcATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTg                                                                                                                                                                            <  1:38959/250‑1 (MQ=255)
                                                                               tCCGGTAACTGCGGCGCTAACTGACGGCAGCATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCCCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCGGGGaa                                                                                                                                                                    >  2:325874/1‑249 (MQ=255)
                                                                                              gcTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTCGAAAt                                                                                                                                                     >  1:251870/1‑249 (MQ=255)
                                                                                                 aactgcgagccaacgatCCCCATATAAGCGACCTCTTCCAGCACGACGGAGTTATGCACCGCATCCTCGGCATTGTTGCCCATGCAAATGGGCCGTGAGAAGGGTCCAGGACGCCTGGCAATTGCGCTGCCTCGATACCCTGTTATTCGAAGGTATCTACGATGACGCGATCGGGTTCCCGCGCATATTCACCGTTGATTTCTACGTCTGTCATTTTGCGGAGCAAGGGAATGGTACAGTAGAAATGGTCg                                                                                                                                                >  2:106129/16‑251 (MQ=255)
                                                                                                 aaCTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGt                                                                                                                                                  <  1:47674/249‑1 (MQ=255)
                                                                                                 aaCTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTc                                                                                                                                                 <  1:445002/250‑1 (MQ=255)
                                                                                                     gACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGc                                                                                                                                              <  1:139186/250‑1 (MQ=255)
                                                                                                         gCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAACTAGTCGGCgtgggg                                                                                                                                        >  2:194296/1‑250 (MQ=255)
                                                                                                           aGAAAATCACCATATAAGCGACCTCTGCCAGAACGATGGCGTGATGCACCGCAGCTTCTGCATTTTTGCCCATGCAAATGTGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCTCTCATATTCACCGTTGATTGCTGCGTCGGTCATTTTGCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCttgggtgt                                                                                                                                      <  2:93031/251‑2 (MQ=255)
                                                                                                             aaTATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGgtggtg                                                                                                                                    >  1:226186/1‑251 (MQ=255)
                                                                                                              atatCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGc                                                                                                                                   >  1:228097/1‑251 (MQ=255)
                                                                                                                           aGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCGTCGGCATTTTTGCCCATGCATATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTGTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGGCGGCGTGGGTGGTGCCGGTTGCTGGAAt                                                                                                                      <  2:328883/251‑1 (MQ=255)
                                                                                                                             cgtccTCTTCCAGCACGGTGGCGTTTTGCACCGCACCTTCGGCATTTTTGCCCAGGCAAATGGGCCGTGAGAATGGACCAGTACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGAGGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTTCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCCGCGTGGGTGGTGCCGTTTGCTGGAATCg                                                                                                                    <  2:341356/248‑1 (MQ=255)
                                                                                                                                ccTCTTCCAGCACGATGGCGTTATTCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACt                                                                                                                 <  2:415002/251‑1 (MQ=255)
                                                                                                                                       cAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGc                                                                                                          >  1:170752/1‑251 (MQ=255)
                                                                                                                                               tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCa                                                                                                                                                                                                                   <  1:41307/139‑1 (MQ=255)
                                                                                                                                               tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCa                                                                                                                                                                                                                   >  2:41307/1‑139 (MQ=255)
                                                                                                                                               tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccc                                                                                                  >  1:157992/1‑251 (MQ=255)
                                                                                                                                                ggCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGaa                                                                                                                                                                    <  2:220875/184‑1 (MQ=255)
                                                                                                                                                ggCGTTATGCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGaa                                                                                                                                                                    >  1:220875/1‑184 (MQ=255)
                                                                                                                                                        gCACCGCATCTTCGGCATTTTTGCCCATGCAAATGGGCCGGGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCGAGAtagtgg                                                                                         >  2:413474/1‑251 (MQ=255)
                                                                                                                                                           ccGCATCTTCGGCATTTTTGCCCCTGCAAATGGGCCGTGCGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGCAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTTCCGTAGGAATAGTCTGCGGGGGTGGAGGCGGTTGATGGAATCGACTGGCCCGCCTGCGCCCAgtgagtggcgg                                                                                      >  2:54415/1‑241 (MQ=255)
                                                                                                                                                                 ccttcgtcattttGGCCCATGCAAATGGGCCGTGAGAATTGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTTCCGTTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGgcgc                                                                                 <  1:272575/242‑1 (MQ=255)
                                                                                                                                                                   tCGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATGGTCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCCCTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGTAATCGACTGGCCCGCCTGCGCCCAGATAGtggcgtggcg                                                                                  <  1:364440/247‑1 (MQ=255)
                                                                                                                                                                    cGGCATTTTTGCCCATGCAAATGGGCCGTGAGAATCGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAg                                                                              <  1:206653/250‑1 (MQ=255)
                                                                                                                                                                        aTTTTTGCCCATGCAAATGGGCCGTGAGAATGGACCCGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtggagggagggggt                                                                         >  1:443435/1‑240 (MQ=255)
                                                                                                                                                                               gcccATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCTTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAAGGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTGGCTGGAATCGTCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGca                                                                    <  2:334298/249‑1 (MQ=255)
                                                                                                                                                                                    tGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCAGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTgtgggtgcgcaatgc                                                              >  2:155069/1‑237 (MQ=255)
                                                                                                                                                                                      aaaatgggCCGGGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATTCCATGTTTTGCGACGGTTTCTACGATGACGTTACCGGTTTCCCACTCAAATTCTCCTGTGATTTCTGCGTCGGTCATTTCGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGTTTGCTGGAACCGACGGGTCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATgccg                                                            <  1:71902/246‑1 (MQ=255)
                                                                                                                                                                                              ccGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGGGTGTGTGCACAATGCCGCCCATgg                                                     >  1:260598/1‑250 (MQ=255)
                                                                                                                                                                                                    gAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATGCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCCCAATGCCGCCAATGGACGGGa                                               >  1:418595/1‑250 (MQ=255)
                                                                                                                                                                                                        ggACCAGGACGCCGGGCATTTGCGCTTCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTGCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGcc                                          <  1:335307/251‑1 (MQ=255)
                                                                                                                                                                                                                   ccGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAAAGCCGCCAATGGACGGGAATGCCTGATAGagc                                 >  1:193923/1‑249 (MQ=255)
                                                                                                                                                                                                                                cTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGGTGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGCCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGCATCGACTGGCCCGCCTGCGCCCAGATGGTGGCGAGGCGGGAGTGTGTGTGCACAATGCCGCCAATGTACGGGCATTCCTGATAGAGCAGTCGGTGAGTTggg                  >  1:259780/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aTCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACTTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgtat              <  1:128284/251‑3 (MQ=255)
                                                                                                                                                                                                                                    aTCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGTTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGGGCACAATGCCCCCAATGGACGGGAATGCCTGATAGAGGCGGCGGTGAGTTggggtat              >  2:198964/1‑246 (MQ=255)
                                                                                                                                                                                                                                       gATACCCTTTTTGTCGAAGTTTTCTACGATGCCTGTACCGTTTTCACACTCAGCTGCACCGTTGATTGCTGCGTCGGTTATTGTGCGGGTGCAGGGAATGGTGCCGTAGAAAGAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGGGGCGTGGCGCGAGTGTGTGGGCACAATGCCGTCAATGGACTGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCgt           <  2:36664/251‑2 (MQ=255)
                                                                                                                                                                                                                                                tttttCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGCGGGCt  >  1:429037/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCT  >  NC_000913/65846‑66336

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