Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,134 G→A G139G (GGC→GGT)  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,1340GA97.7% 121.7 / NA 45G139G (GGC→GGT) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (25/17);  minor base C (1/0);  total (26/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.82e-01

TTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAAC  >  NC_000913/65886‑66352
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                           
ttACGCAGATAGTGTTTATTCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCg                                                                                                                                                                                                                          >  2:38959/1‑251 (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‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccc                                                                                                                  >  1:157992/1‑251 (MQ=255)
                                                                                                        ggCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGaa                                                                                                                                                                                    >  1:220875/1‑184 (MQ=255)
                                                                                                        ggCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGaa                                                                                                                                                                                    <  2:220875/184‑1 (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)
                                                                                                                                                                                                               aaGGTTTCTTCGATGACGTTTCCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGTTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCTTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCgc           <  1:324286/251‑2 (MQ=255)
                                                                                                                                                                                                                        aCGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAAc  >  1:385792/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                           
TTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCA‑GGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAAC  >  NC_000913/65886‑66352

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: ‑