Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 82,909 C→T L207L (TTG→TTA leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091382,9090CT91.7% 165.2 / ‑1.6 60L207L (TTG→TTAleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base T (30/25);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCAC  >  NC_000913/82664‑83154
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
ccAACAATGGCTTTGTTGTGCGGGACCGGAATATCACAAAACTGGTTAACGAACTGGCTGGTGCGCCATATCCCCTGGTGATAAATGGCGGGGTGGACTTTGAGAATATCCTTACGAACTGTGATCTCCATGATGACTTCTCCCAGGGAACAGTTTCCGGCCCGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCCCCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAAtc                                                                                                                                                                                                                                                     <  2:158316/251‑1 (MQ=255)
     aaTGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGGCGGTGAAAATAACCTTACGAACTTTTATCGCCATGATGACTTTTTTCAGGGGACAGATTCCGGCACACTCGCCGGGACCGGTTATTGGGGCTTTCACCCGGCGGGCACCGGCATGGACCGCCGCCCGTGGGTTGCCGACCCGCCAGGCGCAAtcgtcgg                                                                                                                                                                                                                                                >  1:108774/1‑250 (MQ=255)
        ggCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAATTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCCGTGAGTTGCCGACCCCCAGGCCTCAATCGTCGTggg                                                                                                                                                                                                                                            >  1:477979/1‑251 (MQ=255)
                      gggATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAAAATATCCTTACGAAATTTGGTCCACCTGATGACTTCTTCCAGGGCACAGTTTCCGGCACGCCCGCCGATCCCGTTTCTTGCGCCCTCCCCCTGGCGGGCACCGGCCTGTACCGACGCCCATGAGGTGCCGCCCGCCAGGCCTAAATCGTCGGTGGGttgggagagccc                                                                                                                                                                                                                                >  1:419830/1‑240 (MQ=255)
                                tatTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGGAGTTGAGAATATCCTTACGAACTTTTATCGCCATGATGCCTTCTTCCAGGGAACAGTTTCCGGCCCGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCCCCGCCCGTGAGTTTCCGACCGCCCGGCCTACATCGTCGTGGGTATGTACgggaattaatgctttt                                                                                                                                                                                                                     >  1:334020/1‑237 (MQ=255)
                                   tACAAATCTGGTTAACTAACTGGCTGGTGCGCCATATCTCCTGGGGATGAATGGCGGTGTGGACGTTGAGAATATCTTTACGACCTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATGGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGACTGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTc                                                                                                                                                                                                                  <  1:493918/250‑1 (MQ=255)
                                                aactaactGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGGTGAGAATATCCTTACGGACTTTGATCGCCATGGTGGCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTTCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCCCCCGTGAAGTGGCCACCGCCCGGCCTAAATCGTACTGGGTGTGTTCGGAAAAAATGGCttttttggtgtttgggacgg                                                                                                                                                                                                    >  1:562203/1‑234 (MQ=255)
                                                     actGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGGGTGGGCGTTGAGAATATCCTTTCGAACTTTGATCGCCCTTATGACTTCTTCCAGGGGACAGTTTTCGGCACCCTCGCCGATCCCGTTCATTGGGCCTTCCACCTGGGGGGCCCCGGCATGGACCGCCGCCAGTGGGTTGCCGACCGCCAGGGCCAAATCGTCGTGGGGAGATACGGGAGAAATGGCTTTGTCGAAGTTatgtacgggttc                                                                                                                                                                                                >  1:13109/1‑239 (MQ=255)
                                                        ggctggTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCAt                                                                                                                                                                                              <  1:264602/249‑1 (MQ=255)
                                                        ggcgtgTGCGCCATACCTCCTGGTTATTAATGGGGGGGTGGGCGTTGAGAATATCCTTACGATCTTTGACCGCCAAGATGACTTCTTCTAGGGACCAGTTTCCGGCACGCTCGCCGATCCCGTTCATGGCGCCTTCCCCCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCTATGTTAGGTACGCGTTCAt                                                                                                                                                                                              <  2:149375/244‑1 (MQ=255)
                                                          ctggTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGAAGTTGAGAATATCCTTACGAACTTTTATCGCCATGATGCCTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGGTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCACTGAGGTGCCGACCGCCAGGCCTAAATCCTCGTGGGTCTGTTCGGAAAAAATGGCTTTTTCGATGTTGGGGACGCGTTCATACa                                                                                                                                                                                           >  1:327290/1‑250 (MQ=255)
                                                                 gcCATATCTCCTGGTGATTAATGGCGGTGGGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATACCGTTCATTGCGCCTTCCACCTGGCGTGCCCCGGCATGTAACGCCGCCAGTGAGTTGCCGACCGCCCGGCCTAAATCGGCGTGGGTATGTCCGGAAAAAATGGCTTTGTCGATGGTAGGGACGCGTTCATACAGGCCGCtg                                                                                                                                                                                   >  1:592660/1‑251 (MQ=255)
                                                                   cATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCtgat                                                                                                                                                                                 <  1:484385/251‑1 (MQ=255)
                                                                    atatCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCtgatg                                                                                                                                                                                <  1:72922/251‑1 (MQ=255)
                                                                    atatCTCCTGGTGATTAATGGCGGTGTGGAAGTTGAGAATATCCTTACGGACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGAGCCTTCCACCTGGGGTGCACCGGCATGTAACGCCGCCCGGGAGGTGCCGACCGCCCGGCCTAAATTGTCGTGGGGATGTAAGGGAATAATGGGTTTGTTCATGTTAGGTACGCGGTCATACAgggcggtgggg                                                                                                                                                                                >  1:647212/1‑242 (MQ=255)
                                                                       tctcCTGGTGATTAATGGCGGTGTGGACGTTGAGAAAATCCTTACGCACTTTGATCGCCATGATGACTTTTTCCAGGGAACCGGTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCCGTGAGGTGGCGACCGCCAGGCCCAAATCGGCGGGGGTATGTACGGGAATAATGGCTTTTTCGGTGGTAGGTACGCGGTCATAAAGGCCGCtgatgttt                                                                                                                                                                             >  1:168196/1‑248 (MQ=255)
                                                                                  taaaTGGCGGGGTGGACGTTGAAAATTTCCTTCCGACCTTTGATCCCCATGATGTCTCCTTCCAGGGAACAGTTTCCGGCACGCGCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACt                                                                                                                                                                   <  2:506947/248‑1 (MQ=255)
                                                                                        gtgttgtgtaagttaagaATATCCATACGAATGTGTATCGCCATGATGACTTCTTCGAGGTAACTGTTTCCGGCACGCCCGCCGATCCCGTTCACTGCGCCTTCCACCTGGCGTGCACCGTCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACgg                                                                                                                                                            <  2:164090/236‑1 (MQ=255)
                                                                                          ggTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCg                                                                                                                                                                                                                 >  2:525151/1‑196 (MQ=255)
                                                                                          ggTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCg                                                                                                                                                                                                                 <  1:525151/196‑1 (MQ=255)
                                                                                          ggTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAATGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGc                                                                                                                                                           >  1:680482/1‑250 (MQ=255)
                                                                                                gACGTTGAGAATATCCTTACGAACTTTGACCGCCATGATGACTCCTTCCAGGGAACAGTTTCCGGCACGCCCGCCGATCCCGTTCTTTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGgtgt                                                                                                                                                    <  1:572654/251‑1 (MQ=255)
                                                                                                  cGTTGAGAATATCCTTACGAACTTTGATCGCCATGACGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGATTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTGAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAg                                                                                                                                                  <  1:283031/251‑1 (MQ=255)
                                                                                                        gaATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACg                                                                                                                                            <  1:147707/251‑1 (MQ=255)
                                                                                                            atCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGACCAGTTTCCGGCACACTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCCCCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGGTGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAAATCAAACGGCATGGTGTAGCCCACggagt                                                                                                                                        >  2:424831/1‑248 (MQ=255)
                                                                                                              ccTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTcc                                                                                                                                     <  1:233896/251‑1 (MQ=255)
                                                                                                                        tttGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTCTGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGGTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATCGTGTAGCCCACGCGGTCCGGCATGTTg                                                                                                                            >  2:327453/1‑250 (MQ=255)
                                                                                                                             aTCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGAtggg                                                                                                                       >  1:662928/1‑250 (MQ=255)
                                                                                                                                 ccATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGAtggtggtg                                                                                                                   >  1:41399/1‑251 (MQ=255)
                                                                                                                                  cATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGACTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGAtggtggtgg                                                                                                                  >  1:298919/1‑251 (MQ=255)
                                                                                                                                           ttcttcCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCAt                                                                                                         >  1:37263/1‑251 (MQ=255)
                                                                                                                                             cttcCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCCGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTACGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCAtt                                                                                                        >  1:441902/1‑250 (MQ=255)
                                                                                                                                                    ggAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTACATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGCTGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATTGTGTAGCCCACGGTGTCCGGAATGTTGCTGGTGGTGGCACCGGCATTAAccgccg                                                                                                >  2:349086/1‑251 (MQ=255)
                                                                                                                                                     gAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCCGTGAGTTGCCAACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAACTAATAGCTTTGTCGCTGTTAGGTACGCGTTCATACAGGCAGCTGATGATTCCGTCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGTTGGTGGTGGCACCGGCATTAATcgccgc                                                                                               >  1:420152/1‑251 (MQ=255)
                                                                                                                                                         ggttcccggccggccctccGATCCCGTCATTTGCGCCTTCCCCGGGGCGGGCCCCGGGATGTACCGCCGCCAGTGAGTTTCCGACCGCCATGCCTAAATCGTCGTGTGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCg                                                                                           <  2:669302/234‑1 (MQ=255)
                                                                                                                                                            ttCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGCATGTACGGAAATAATGGCTTTGTCGATGTTAGGGACGCGTTCATACAGGCCGCTGCTGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCCCCGGCATTAATCGCCGCTTCGAcc                                                                                        >  2:572364/1‑251 (MQ=255)
                                                                                                                                                                ggCACGCTCGCCGACCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTc                                                                                    <  1:675963/251‑1 (MQ=255)
                                                                                                                                                                  cACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCCTGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTAGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGAATGGTGTAGCCCACGGTGTCCGGCATGTTGATGGCGGTGGCACCGGCATTAATCGACGCTTCGGCCACTcgc                                                                                  >  2:604572/1‑251 (MQ=255)
                                                                                                                                                                     gCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGACATGTACCGCCGCCAGTGAGTTGCCGGCCGCCAGGACTAAATCGTCGTGGGTATGGACGGACATAATGGCCTTGTCGATGTTAGGTACGCGTTCCTACAGGCCACTGATGATTCCGGCGATCACAAACGGCTTGGAGTAGCCCACGGCGGTCCGGGTGGTGATGGTGGTGGCGCCGGCATTaaccaccgactcgaacaatcacccc                                                                               >  2:688845/1‑228 (MQ=255)
                                                                                                                                                                          ccGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATc                                                                          <  1:173255/251‑1 (MQ=255)
                                                                                                                                                                          ccGATCCCGTTCATTGCGCCTCCCTCCTGGCGTGCACCGGCATGTACCGCGGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACGCAAACGGCATGGTGTAGCCCACGGTGCCCGTAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATc                                                                          <  1:710583/251‑1 (MQ=255)
                                                                                                                                                                           cGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCCCGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTCATCGCCGCTTCGACCACTCGCGCCAGGTCg                                                                         >  2:21121/1‑251 (MQ=255)
                                                                                                                                                                            cgttcccttgcgtggcgcCTTCCCCCTGGCGTGCCCCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATTTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCg                                                                         <  1:651896/237‑1 (MQ=255)
                                                                                                                                                                                  gTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATgg                                                                  >  1:39286/1‑251 (MQ=255)
                                                                                                                                                                                  gTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATgg                                                                  >  1:239218/1‑251 (MQ=255)
                                                                                                                                                                                     cATTGCGCCTTCCACCGGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGgt                                                                <  1:278022/250‑1 (MQ=255)
                                                                                                                                                                                                         cGTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACCGTGTCCGGTATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGGTCGGCAATGGGTGTACGCCCGGCTACTTCGcc                                            >  2:174577/1‑250 (MQ=255)
                                                                                                                                                                                                          gTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCa                                                                              <  1:328279/216‑1 (MQ=255)
                                                                                                                                                                                                          gTGCACCGGCATGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGGGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCAACCGCGCCa                                                                              >  2:328279/1‑216 (MQ=255)
                                                                                                                                                                                                             cacccgcaggcaCCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTTGGTATGTACGGAAAGAATTGTTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGATCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGaa                                        <  2:359843/240‑1 (MQ=255)
                                                                                                                                                                                                                ggccttgtcccgccgacgtggAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACTCGTTCATACAGGCCGCTGATGATTCCGGCGAACGCAAACGGCATGGGGTAGCCCACGGTGTCCGGAATGGTGATGGTGGTGGCACCGGCATTAATCGCCGCTACGACCACACGCTACAGAGCGGCAATGGGCGTACGCCCGGCATCTTTGCAAGAAAAt                                     <  2:73885/231‑1 (MQ=255)
                                                                                                                                                                                                                 ggcctGTACCGCCGCCAGTGAGTTGCCGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGAACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCACGAAAAtt                                    >  1:283726/5‑251 (MQ=255)
                                                                                                                                                                                                                   caggTACCGCCGCCAGTGAGGTGCCGACCGCCAGGCCTAAATCGTCGTGGGTTAGTACGGAAATAATGGCTTTGTCGATGGTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGAGAACACAAACGGCATGGTGTAGCCCACGGTGGCAGGTATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCGTCGACAACTCGCGCCAGATCGGAACTTGGTGAAAGCCCAGCATATTCGCAAGACAAtacc                                  >  2:259906/4‑248 (MQ=255)
                                                                                                                                                                                                                           cgccgctattgagtggccGACCGTCAGGGCTAAATCTTCGGGGGTATGTACGGAAATAATGTCTTTGTCGATGTTAGTTTCGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGGAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGTCATCTTCGCAAGAAAATTCAAcatcat                           <  2:664650/235‑1 (MQ=255)
                                                                                                                                                                                                                                       gttgccGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTAc               <  1:21121/247‑1 (MQ=255)
                                                                                                                                                                                                                                       gttgccGACCGCCAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTAc               <  1:683944/247‑1 (MQ=255)
                                                                                                                                                                                                                                          tCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGAACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGCGAAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCACTGGGTGTACGCCCGGCATCTTCGCTAGAAAATTCacaataatcggtcaaattaaggg            >  2:361362/1‑229 (MQ=255)
                                                                                                                                                                                                                                                    cAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTc    >  1:96507/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    cAGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACACAAACGGCATGGTGTAGCCCACGGTGTCCGGCATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTc    >  2:274896/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     aGGCCTAAATCGTCGTGGGTATGTACGGAAATAATGGCTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGGAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCAc  <  1:42867/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
CCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAACTTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCATGTACCGCCGCCAGTGAGTTTCCGACCGCCAGGCCCAAATCGTCGTGGGTATGTACGGAGATAATGGCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCAC  >  NC_000913/82664‑83154

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