Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 86,345 G→C C239S (TGC→TCC)  ilvI → acetolactate synthase/acetohydroxybutanoate synthase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091386,3450GC91.7% 160.2 / ‑1.5 62C239S (TGC→TCC) ilvIacetolactate synthase/acetohydroxybutanoate synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base C (24/31);  minor base T (5/0);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 2.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATATTCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTAGCGGCAAAAAAACCGGTTGTCTACGTAGGCGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAACGGTGGAGGCGTTGAATCTGCCCGTTGTTTGCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCC  >  NC_000913/86118‑86585
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
aGATATTCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTgggggggg                                                                                                                                                                                                                            >  1:359632/1‑248 (MQ=255)
   tatTCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAACAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAACTTAAAGAACTGGTGGAGGGGGTTGATGTGCCCGTTGTTTTCTCTTTGAGGTGGCTGGGGGGGTTTc                                                                                                                                                                                                                         >  2:179063/1‑251 (MQ=255)
      tCTTAATCCGGCGAACAAATTACCCTATGGCTGGCCGGAGTCGGTTAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGTGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCGGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCg                                                                                                                                                                                                                       <  2:440339/250‑1 (MQ=255)
        ttAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCCTCAGCAGTTGGAAGAACTGGGGGCGGCGTTTAATCTGCCCGTTTTTTCCTCATTGGTTGGGGTGGGGGCGTTTTCGGc                                                                                                                                                                                                                     >  2:113398/1‑250 (MQ=255)
                          tACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACt                                                                                                                                                                                                   <  1:277638/250‑1 (MQ=255)
                            cccTATGTCTGGCCGGAGTCGGTCAGTTAGCGTTCTTACAAGCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCATCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTggg                                                                                                                                                                                                <  1:92836/251‑1 (MQ=255)
                                          ggAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAAAAATGGTGGTGGCGTTGAATCTGCCCGTTGTTTCCTCATTTATGGGGCTGGGGGGGTTTCCGGCAACGCATCGTCAGGCACTGGGGATTCTGGGAATGc                                                                                                                                                                                 >  1:283674/1‑251 (MQ=255)
                                                 gTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCAt                                                                                                                                                                                                                                              <  1:627425/184‑1 (MQ=255)
                                                 gTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGGGGGGGCGTTGGATCTGCCCGTTGTTTTCTCAt                                                                                                                                                                                                                                              >  2:627425/1‑184 (MQ=255)
                                                                 aCAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATg                                                                                                                                                          <  1:125299/251‑1 (MQ=255)
                                                                  caagcccACTGCTACCGGACATAAAGGGCAAATTAATCGTGCGCTGCAAACGCTGGTGGCGGCAAAAAAACCGGGTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGGGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGa                                                                                                                                                         <  2:270761/247‑1 (MQ=255)
                                                                     gccctctacgaCCGGAAATAAAGGGCAAATTAAGCCGGCGCTGCAAATGCTGGGGGCGGCAAAAAAACTGGGTGTCGTCGTAGGTGGTGGGGCAATCACGGCGGGCTGCCAGTAGCGTTTGAAAGAAAGGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGCTGGGGCTGGGGGCGTTTCCGGCAACGCATCGTTAGGCACGGGGCATGCTGGGAATGCACGTTACCTGCGAGGCCAATATGACGa                                                                                                                                                      <  2:467868/246‑1 (MQ=255)
                                                                        cACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGGGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGc                                                                                                                                                   <  2:178211/251‑1 (MQ=255)
                                                                         actactACCGGACCAAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTTGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCa                                                                                                                                                  <  2:334546/251‑1 (MQ=255)
                                                                                 cGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCCTCAGCAGTTGAAAGGAATGGGGGGGGCGTTGAATCTGGCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCCACGCATCGTCAGGCACTGGGGCTGCTGGGGATGAACGGTACCTACGATGCCAGAAAGACGAGGCATaaacggg                                                                                                                                          >  2:100663/1‑246 (MQ=255)
                                                                                        aaaGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCGGCCATCCGCACTTGAAAAACATGGTGGGGGGGTTGAATCTGCCCGTTGTTTCTTCATTGTTGGGGCTGGGGGGGGTTCCGGCAACGCATCGTCCGGCACTGGGCATGCTGGGAATGGCAGGGACCCACACACCCAATATTACGATGCATAAAGCggtgggtgt                                                                                                                                   >  2:245958/1‑244 (MQ=255)
                                                                                        aaaGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGATGCCCTCAGCAGTTGAAAGAAAAGGGGGAGGCGTTGAATCTGCCCGTTGTTTTCTTATTGATGGGGCTGGGGGGGTTTCCGGCAACGCATAGTTAGGACCTGGGGATGATGGGAATGCACGGGAACCACGAAGCCAATATGTTGATGGATAACGCGGATGTGAt                                                                                                                                   >  2:161878/1‑251 (MQ=255)
                                                                                               aaaTTAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGGGGAGGCGTTGAATCTGCCCGTTGGTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCCTGCTGGGGATGCACGGTACCTACGAAGCCAATAAGAAGCTGCATAACGCGGATGTGATTTTAGcc                                                                                                                            >  1:562425/1‑251 (MQ=255)
                                                                                                  ttAAGCGTGCCCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCAACAGCAGTTGAAAGACATGGGGGAGGCGTTGGATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGAACTGGGAATGCTGGGAATGGACGGTCACTACGAAGCCAAAATGACGATGCATAACGCGGGAGTGATTTTAGCCggc                                                                                                                         >  1:682484/1‑249 (MQ=255)
                                                                                                   tAAGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGATGAAAGAAATGGGGGAGGCGTTGAATCTGCCCGTTGTTTTCTCATTGATGGGGATGGGGGCGGTTCCGGCAAAGCATCGTCAGGCACTGGGCCTGCTGGGAATGGACGGTACCTACGAAGCCAATATGACGAAGCATAACGCGGATGTGGTTTTCggcggcg                                                                                                                        >  1:662764/1‑245 (MQ=255)
                                                                                                     aGCGTGCTCTGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCCGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCCTGCTGGGAATGCCCGGTACCTACGAAGCCAATATGGCGATGCATAACGCGGATGTGATTTTCGCCGTCggg                                                                                                                      >  1:242352/1‑251 (MQ=255)
                                                                                                            tctGCAAACGCTGGTGGCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTTCGGCAACGCATCGGCAGGCACTGGGCATGCTGGGAATGCACGGTAACTACGAAGCCAATATGACGATGCCTAACGCGGCTTTGATTTTCGCCGTCgggggtcgtt                                                                                                               >  1:163593/1‑245 (MQ=255)
                                                                                                               gccaaCGCTGGCGGCGTCAAAAAAATCGGTTGGTTAAGAAGGTGGTGGGGCAACCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGGTGAATCTGCCCGTGGTTGCCTCATTGATGGGGCTGGGGGCGTTTCCGGCACCGCATCGTCAGGCACGGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGAtgtg                                                                                                            <  2:19721/248‑4 (MQ=255)
                                                                                                                        ggtggCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGGACGCTTTGATGACCGa                                                                                                    >  1:428038/5‑250 (MQ=255)
                                                                                                                         gtggCGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAac                                                                                                  <  1:387582/248‑1 (MQ=255)
                                                                                                                            cgatccaaaaacCCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCAGCAGTAGTTGAAAGAAACGGTGGAGGCGGTTAACCTGCCCGTTGTTTGCTCAGTGATGGGTCTGTGGGCGTTTCCGGCAACTCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTCCCTACGAAGCCAATATGACGATGCATAACGCTGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAacgac                                                                                               <  1:312540/245‑1 (MQ=255)
                                                                                                                             cGGCAAAAAAACCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGGGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATATGATTTTCGCCGGCGGGGTTCGATTTGATGACCGAacgacg                                                                                              >  1:507298/1‑251 (MQ=255)
                                                                                                                                  aaaaaaCCGGTTGGCTACGGAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACaa                                                                                         <  2:428038/251‑1 (MQ=255)
                                                                                                                                      aaCCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGGTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTg                                                                                     <  1:678704/251‑1 (MQ=255)
                                                                                                                                      aaCCGGTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGAGGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTg                                                                                     <  1:449626/251‑1 (MQ=255)
                                                                                                                                          ggTTGTCTACGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCaa                                                                                 <  1:480726/251‑1 (MQ=255)
                                                                                                                                              gtcgaCGTAGGTGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACTCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGATCAATCTGGCAAAGt                                                                              <  2:662764/246‑1 (MQ=255)
                                                                                                                                              gTCTACGTAGGTGGTGGGGCAATCCCGGCGGGCCGCCATCAGCAGTTGAAAGAAAGGGTGTAGGCGTTGAACCTGCCCGTTGTTTCCTCATTGAGGTGGCTGGGGGCGTTTCCGGCAACGCATCGGCAGGCACTGTGCATGCTTGGAATGCTCGGGACCTCCGAAGCCAATATGACGATGCATAACGCGGATGTGATGTTCGCCGTCGGGGTACGATTTGAGGACCGAACGTCTATCAATCTGGCAAAGTa                                                                             <  2:99567/251‑1 (MQ=255)
                                                                                                                                                  aCGTAGGTGGTGGGGTAATCACGGCTGGCTGCCAGCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAAACGGCCCGGTGTTTCCTCAGTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGGCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGc                                                                         <  2:527141/251‑1 (MQ=255)
                                                                                                                                                    gTAGGTGGTGGGGCAATCACGGCGGGCTGACATCAGCAGTTGCAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCAATGGGCATGCTGGGAATGCACGGGACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGGTTTTCGCCGTCGGGGGACGATTTGATGACCGAACGACGAACAATCttgccacatactgccc                                                                       >  1:571943/1‑236 (MQ=255)
                                                                                                                                                       ggtggTGGGGCAATCACGGCGGGCCGCCATCAGCAGTTTAAAGAAATGGTGGAGGCGGTGGATCTGACCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTTCCGCCACGCATCGACCGGGACTGGGCATGGTGGGGATGTACGGTACCTACGAAGCCAATATGACGATGCCTCACGCGGATGTGATTTTCGCCGTCGGGGTACGAGTTGATGACCGAAAGGCGAACAATCCTGCAATGTCCTGCCCaa                                                                     >  2:685228/4‑250 (MQ=255)
                                                                                                                                                         tggTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTTCTCATTGGTGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGGTTTTCGCCGTCGGGGGAAGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCaaa                                                                    >  1:423809/2‑249 (MQ=255)
                                                                                                                                                             ggggCAATCCCTGCGGGCGGCCATCGTCAGTTGAAAGAAATGGTGGGGGCGGTGAAGCTGCCCGTTGTTTCCGCATTGATGGGGCTGGGGTCGTTTCCGGCAACGGATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCTCAAATGCCAc                                                              <  2:507298/251‑1 (MQ=255)
                                                                                                                                                                   aTCACGGCGGGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCt                                                        >  1:219215/1‑251 (MQ=255)
                                                                                                                                                                      ccggcggGCTGCCAGCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCCCGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATTCCACTGTTCTGCa                                                     <  1:366408/250‑1 (MQ=255)
                                                                                                                                                                       cggcggGCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCat                                                    <  1:27326/251‑1 (MQ=255)
                                                                                                                                                                          cggGCTGCCATCATCAGTAGAAAGAAATGGTGGATGCGTTGAATCTGTCCGTTGTTCCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTTCATATc                                                 <  2:571943/251‑1 (MQ=255)
                                                                                                                                                                           gggCTGCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCCTAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGACCGACGAACAATCTGGCAAAGGACTGCCCAAATGCCAATGTTCTGCATATCg                                                >  2:322908/1‑251 (MQ=255)
                                                                                                                                                                                gCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGCATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCCCGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGCTTTTCGCCGTCGGGGGAAGATTTGATGACCGAACGCCGAACAATCTGGCAAAGGACTGCCCAAATGCCCCTGTTCTGCATTACGATAtt                                           >  2:703091/1‑251 (MQ=255)
                                                                                                                                                                                gCCATCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATAtt                                           <  1:11433/251‑1 (MQ=255)
                                                                                                                                                                                  cATCAGCAGTTGAAAGAAATGGTGGAGGTGGTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGGTTCCGGCAACGCATCGTCAGGCACTGGTCATGCTGGGAATGCTCGGTACCTACGAAGCCAATATGACTATGCATAACGCGTATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGATCGACGAACATTCTTGCAAAGTACTGCCCAAATGCCTCTGTTCTGCATATCGCTAttta                                         <  2:151873/251‑3 (MQ=255)
                                                                                                                                                                                    tCAGCAGTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATATGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGACCTGCCCAAATGCCACTGTTCTGCATATAGATATTGATc                                       >  1:173768/1‑251 (MQ=255)
                                                                                                                                                                                      ggcagTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTGGTTTCCCCATTGAGGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGAGGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACTACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCt                                     <  2:359632/250‑1 (MQ=255)
                                                                                                                                                                                         agTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGTGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACt                                  <  1:108453/251‑1 (MQ=255)
                                                                                                                                                                                          gTTGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCATGCTCTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACtt                                 <  1:659712/251‑1 (MQ=255)
                                                                                                                                                                                           ttGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGAACGATTTGATGACCGAACGACGAACAATCTGGCAAATTACTGCCCAAATGCCACTGATCTGCCTATcattatcactcctacttc                                >  2:62270/1‑234 (MQ=255)
                                                                                                                                                                                           ttGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGAACGGTACCTACGAAGCCAATATGACGATGCATAAAGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGCTGACCGAACGACGAACAATCTGGCAAAGTACGGCCCAAATGCCAGTGATATGCATATCGAAATTGATCCTACtac                                >  2:214611/1‑249 (MQ=255)
                                                                                                                                                                                            tGAAAGAAATGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTc                                <  1:113398/250‑1 (MQ=255)
                                                                                                                                                                                                 gaaaTGGTGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGACCGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCTTATCGATATTGCTCCTACTTCCCTTTc                          >  2:668585/1‑251 (MQ=255)
                                                                                                                                                                                                       gTGGAGGCGTTGAATCTGCCCGGTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAAc                    <  2:171861/251‑1 (MQ=255)
                                                                                                                                                                                                        tGGAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAAcc                   >  1:166213/1‑251 (MQ=255)
                                                                                                                                                                                                          gAGGCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGCAACCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGAACGATTTGATGACCGAACGACGCACAATCTGGACAAGAAATGCACAAATGCCCCTGTTCTGCATATAGCTATTGATCCTACTTCCATTTCTAAAACCGt                 >  2:136140/1‑251 (MQ=255)
                                                                                                                                                                                                           aGGCGTTGAATCTGCCCGTTGTTTGCTCATTGATGGGGCTGGGGGCGTTTCCGCCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGTTGACCGACCGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTCCTAAACCCGTg                >  1:45423/1‑251 (MQ=255)
                                                                                                                                                                                                             gCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGAc              <  1:81148/251‑1 (MQ=255)
                                                                                                                                                                                                             gCGTTGAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGAc              >  1:468549/1‑251 (MQ=255)
                                                                                                                                                                                                                   gAATCTGCCCGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCgg         >  1:106992/1‑251 (MQ=255)
                                                                                                                                                                                                                          cccGTTGTTTCCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATccc  <  1:179063/251‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
AGATATTCTTAATCCGGCGAACAAATTACCCTATGTCTGGCCGGAGTCGGTCAGTATGCGTTCTTACAATCCCACTACTACCGGACATAAAGGGCAAATTAAGCGTGCTCTGCAAACGCTGGTAGCGGCAAAAAAACCGGTTGTCTACGTAGGCGGTGGGGCAATCACGGCGGGCTGCCATCAGCAGTTGAAAGAAACGGTGGAGGCGTTGAATCTGCCCGTTGTTTGCTCATTGATGGGGCTGGGGGCGTTTCCGGCAACGCATCGTCAGGCACTGGGCATGCTGGGAATGCACGGTACCTACGAAGCCAATATGACGATGCATAACGCGGATGTGATTTTCGCCGTCGGGGTACGATTTGATGACCGAACGACGAACAATCTGGCAAAGTACTGCCCAAATGCCACTGTTCTGCATATCGATATTGATCCTACTTCCATTTCTAAAACCGTGACTGCGGATATCCC  >  NC_000913/86118‑86585

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