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→TTC) araDL‑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

ATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCT  >  NC_000913/65846‑66336
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ataCAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGt                                                                                                                                                                                                                                                  <  1:356135/251‑1 (MQ=255)
    aGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCt                                                                                                                                                                                                                                              <  1:12742/251‑1 (MQ=255)
              gCCCGTAATATGCCTTCGCGCCATGCTTACGCAGCTAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCCCGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGt                                                                                                                                                                                                                                    >  2:438626/1‑251 (MQ=255)
                  gTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTAcc                                                                                                                                                                                                                                >  1:310626/1‑251 (MQ=255)
                     atatGCCTTCGCGCCATGCTTACGCAGATAGAGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATAGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGt                                                                                                                                                                                                                             >  1:412876/1‑251 (MQ=255)
                                 gcaaTGCTTACGCAGATAGTGTTGATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGAGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCata                                                                                                                                                                                                                 <  1:150033/248‑1 (MQ=255)
                                 gcCATGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCata                                                                                                                                                                                                                 <  1:64720/251‑1 (MQ=255)
                                     tGCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTc                                                                                                                                                                                                              >  1:199514/1‑250 (MQ=255)
                                      gCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTca                                                                                                                                                                                                             <  1:5153/250‑2 (MQ=255)
                                      gCTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGAATATCCGTAAACTGCGGCGCTAACTGACGGCAGAATACCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTGCGGCATTTTTGCCC‑ATGCAAATGGGCCGTTAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTca                                                                                                                                                                                                             <  2:358278/250‑2 (MQ=255)
                                       cTTACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTcacc                                                                                                                                                                                                           >  1:235404/1‑248 (MQ=255)
                                        ttACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCg                                                                                                                                                                                                          >  2:38959/1‑251 (MQ=255)
                                          aCGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCAGAACCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGAGACCTCTTCCAGCACGATGGCGTTATGCACCGAATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTAGAAGGTTTCTACGATGACGTTACCGGTTTCCCGCTCATATTcacagtt                                                                                                                                                                                                        >  2:41979/1‑245 (MQ=255)
                                              aGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGAtt                                                                                                                                                                                                    <  1:428399/251‑1 (MQ=255)
                                                   gtgtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCCTCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGc                                                                                                                                                                                               >  2:124492/1‑251 (MQ=255)
                                                    tgtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGTTACCCTGTTTTTCGCAGGTTTCTACGATTGCGTTACCGGTTTCCCACTCCTATTCACCTTTGATTTCTGcc                                                                                                                                                                                              >  2:412342/1‑250 (MQ=255)
                                                    tgtTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCAGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCg                                                                                                                                                                                              >  1:129173/1‑251 (MQ=255)
                                                         attcagcagCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGTCAGAATTTCACCATATAAGCGACCTCTTCCACCACGATGGCGTTATGCACCGCATCTTCGACCTTTTTGCCC‑ATGCAAATGGGCCGCGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGCCGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGt                                                                                                                                                                                         >  1:263057/4‑251 (MQ=255)
                                                           tcagcagCGTTTGCTGAATATCCGGTAACTGCGGCGCTAACTGACGGCAGACTATCCACATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACAGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTc                                                                                                                                                                                        >  1:297911/2‑250 (MQ=255)
                                                                       gctgcATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTg                                                                                                                                                                            <  1:38959/250‑1 (MQ=255)
                                                                               tCCGGTAACTGCGGCGCTAACTGACGGCAGCATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCCCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCG‑GGGaa                                                                                                                                                                    >  2:325874/1‑249 (MQ=255)
                                                                                              gcTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTCGAAAt                                                                                                                                                     >  1:251870/1‑249 (MQ=255)
                                                                                                 aactgcgagccaacgatCCCCATATAAGCGACCTCTTCCAGCACGACGGAGTTATGCACCGCATCCTCGGCATTGTTGCCC‑ATGCAAATGGGCCGTGAGAAGGGTCCAGGACGCCTGGCAATTGCGCTGCCTCGATACCCTGTTATTCGAAGGTATCTACGATGACGCGATCGGGTTCCCGCGCATATTCACCGTTGATTTCTACGTCTGTCATTTTGCGG‑AGCAAGGGAATGGTACAGTAGAAATGGTCg                                                                                                                                                >  2:106129/16‑251 (MQ=255)
                                                                                                 aaCTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGt                                                                                                                                                  <  1:47674/249‑1 (MQ=255)
                                                                                                 aaCTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTc                                                                                                                                                 <  1:445002/250‑1 (MQ=255)
                                                                                                     gACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGc                                                                                                                                              <  1:139186/250‑1 (MQ=255)
                                                                                                         gCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAACTAGTCGGCgtgggg                                                                                                                                        >  2:194296/1‑250 (MQ=255)
                                                                                                           aGAAAATCACCATATAAGCGACCTCTGCCAGAACGATGGCGTGATGCACCGCAGCTTCTGCATTTTTGCCC‑ATGCAAATGTGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCTCTCATATTCACCGTTGATTGCTGCGTCGGTCATTTTGCGGGGGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCttgggtgt                                                                                                                                      <  2:93031/251‑2 (MQ=255)
                                                                                                             aaTATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGgtggtg                                                                                                                                    >  1:226186/1‑251 (MQ=255)
                                                                                                              atatCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGc                                                                                                                                   >  1:228097/1‑251 (MQ=255)
                                                                                                                           aGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCGTCGGCATTTTTGCCC‑ATGCATATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTGTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGGCGGCGTGGGTGGTGCCGGTTGCTGGAAt                                                                                                                      <  2:328883/251‑1 (MQ=255)
                                                                                                                             cgtccTCTTCCAGCACGGTGGCGTTTTGCACCGCACCTTCGGCATTTTTGCCC‑AGGCAAATGGGCCGTGAGAATGGACCAGTACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGAGGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTTCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCCGCGTGGGTGGTGCCGTTTGCTGGAATCg                                                                                                                    <  2:341356/248‑1 (MQ=255)
                                                                                                                                ccTCTTCCAGCACGATGGCGTTATTCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACt                                                                                                                 <  2:415002/251‑1 (MQ=255)
                                                                                                                                       cAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGc                                                                                                          >  1:170752/1‑251 (MQ=255)
                                                                                                                                               tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCa                                                                                                                                                                                                                   <  1:41307/139‑1 (MQ=255)
                                                                                                                                               tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCa                                                                                                                                                                                                                   >  2:41307/1‑139 (MQ=255)
                                                                                                                                               tGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccc                                                                                                  >  1:157992/1‑251 (MQ=255)
                                                                                                                                                ggCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGaa                                                                                                                                                                    <  2:220875/184‑1 (MQ=255)
                                                                                                                                                ggCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGaa                                                                                                                                                                    >  1:220875/1‑184 (MQ=255)
                                                                                                                                                        gCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGGGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCGAGAtagtgg                                                                                         >  2:413474/1‑251 (MQ=255)
                                                                                                                                                           ccGCATCTTCGGCATTTTTGCCCC‑TGCAAATGGGCCGTGCGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGCAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTTCCGTAGGAATAGTCTGCGGGGGTGGAGGCGGTTGATGGAATCGACTGGCCCGCCTGCGCCCAgtgagtggcgg                                                                                      >  2:54415/1‑241 (MQ=255)
                                                                                                                                                                 ccttcgtcattttGGCCCATGCAAATGGGCCGTGAGAATTGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTTCCGTTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGgcgc                                                                                 <  1:272575/242‑1 (MQ=255)
                                                                                                                                                                   tCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGTCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCCCTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGTAATCGACTGGCCCGCCTGCGCCCAGATAGtggcgtggcg                                                                                  <  1:364440/247‑1 (MQ=255)
                                                                                                                                                                    cGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATCGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAg                                                                              <  1:206653/250‑1 (MQ=255)
                                                                                                                                                                        aTTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCCGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtggagggagggggt                                                                         >  1:443435/1‑240 (MQ=255)
                                                                                                                                                                               gcccATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCTTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAAGGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTGGCTGGAATCGTCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGca                                                                    <  2:334298/249‑1 (MQ=255)
                                                                                                                                                                                    tGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCAGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTgtgggtgcgcaatgc                                                              >  2:155069/1‑237 (MQ=255)
                                                                                                                                                                                      aaaatgggCCGGGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATTCCATGTTTTGCGACGGTTTCTACGATGACGTTACCGGTTTCCCACTCAAATTCTCCTGTGATTTCTGCGTCGGTCATTTCGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGTTTGCTGGAACCGACGGGTCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATgccg                                                            <  1:71902/246‑1 (MQ=255)
                                                                                                                                                                                              ccGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGGGTGTGTGCACAATGCCGCCCATgg                                                     >  1:260598/1‑250 (MQ=255)
                                                                                                                                                                                                    gAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATGCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCCCAATGCCGCCAATGGACGGGa                                               >  1:418595/1‑250 (MQ=255)
                                                                                                                                                                                                        ggACCAGGACGCCGGGCATTTGCGCTTCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTGCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGcc                                          <  1:335307/251‑1 (MQ=255)
                                                                                                                                                                                                                   ccGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAAAGCCGCCAATGGACGGGAATGCCTGATAGagc                                 >  1:193923/1‑249 (MQ=255)
                                                                                                                                                                                                                                cTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGGTGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGCCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGCATCGACTGGCCCGCCTGCGCCCAGATGGTGGCGAGGCGGGAGTGTGTGTGCACAATGCCGCCAATGTACGGGCATTCCTGATAGAGCAGTCGGTGAGTTggg                  >  1:259780/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aTCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACTTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgtat              <  1:128284/251‑3 (MQ=255)
                                                                                                                                                                                                                                    aTCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGTTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGGGCACAATGCCCCCAATGGACGGGAATGCCTGATAGAGGCGGCGGTGAGTTggggtat              >  2:198964/1‑246 (MQ=255)
                                                                                                                                                                                                                                       gATACCCTTTTTGTCGAAGTTTTCTACGATGCCTGTACCGTTTTCACACTCAGCTGCACCGTTGATTGCTGCGTCGGTTATTGTGCGGGTGCA‑GGGAATGGTGCCGTAGAAAGAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGGGGCGTGGCGCGAGTGTGTGGGCACAATGCCGTCAATGGACTGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCgt           <  2:36664/251‑2 (MQ=255)
                                                                                                                                                                                                                                                tttttCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGCGGGCt  >  1:429037/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCT  >  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: ‑