Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,186,053 T→G 5.9% E136A (GAG→GCG)  aceK ← isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,186,0530TG 5.9% 246.6 / 5.3 81E136A (GAG→GCG) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (37/39);  new base G (2/3);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

GCATTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTGGATA  >  NC_002947/5185908‑5186196
                                                                                                                                                  |                                                                                                                                                
gcaTTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCAAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTcgc                                                                                                                                               >  6:270386/1‑149 (MQ=255)
gcaTTCGCCAACTTGTTCTTCCAGCCGCGCCAGCTCGCCTTCCAGGTCGCCATACGGCATGTCGAACGGGTAATAGGCGAACACCGCTCGAAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTAAGGCTTTCAGCGGCGCACTcgc                                                                                                                                               >  4:95359/1‑149 (MQ=255)
       cgcAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGAGCGCTCGCGGCcgcg                                                                                                                                       >  5:424758/1‑149 (MQ=255)
            ttgttCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATg                                                                                                                                   >  3:247950/1‑148 (MQ=255)
            ttgttCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATg                                                                                                                                   >  2:278566/1‑148 (MQ=255)
            ttgttCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCCGCGAGCGCGCCGCGCaaac                                                                                                                                  <  8:384596/149‑4 (MQ=255)
               ttcttcCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGaa                                                                                                                                <  7:182952/148‑1 (MQ=255)
                   tcCAGCCGCACCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCg                                                                                                                           >  1:265398/1‑149 (MQ=255)
                    ccAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCgg                                                                                                                          >  3:15456/1‑149 (MQ=255)
                    ccAGACGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGCTAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCgg                                                                                                                          >  5:292108/1‑149 (MQ=255)
                     cAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCggg                                                                                                                         >  8:286167/1‑149 (MQ=255)
                          gcgcCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTcgc                                                                                                                                               >  4:47731/1‑123 (MQ=255)
                              cAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGa                                                                                                                <  2:265398/149‑1 (MQ=255)
                               aGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGAt                                                                                                               <  1:278566/149‑1 (MQ=255)
                                 gTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGa                                                                                                             <  4:177132/149‑1 (MQ=255)
                                  tCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGaa                                                                                                            <  7:65741/149‑1 (MQ=255)
                                      cTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATg                                                                                                        <  7:87592/149‑1 (MQ=255)
                                         ccAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCa                                                                                                      <  4:336700/148‑1 (MQ=255)
                                          cAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGc                                                                                                    >  3:189806/1‑149 (MQ=255)
                                           aGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGcc                                                                                                   <  4:247950/149‑1 (MQ=255)
                                             gTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGt                                                                                                 >  3:258254/1‑149 (MQ=255)
                                                     aGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGAt                                                                                          >  2:141715/1‑148 (MQ=255)
                                                        gCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAgg                                                                                      >  8:52337/1‑149 (MQ=255)
                                                          aTGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCgg                                                                                                                          <  7:205240/111‑1 (MQ=255)
                                                          aTGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCgg                                                                                                                          >  8:205240/1‑111 (MQ=255)
                                                                    ggTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGa                                                                                                             <  3:319797/114‑1 (MQ=255)
                                                                    ggTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGa                                                                                                             >  4:319797/1‑114 (MQ=255)
                                                                    ggTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATg                                                                                                        <  7:187350/119‑1 (MQ=255)
                                                                    ggTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCCATAGGTACGGGTTTGCGCGGCGCGCTAGCGGCCCCGCATGGAACGCCCGCTATTCAGTATGAACATg                                                                                                        >  8:187350/1‑119 (MQ=255)
                                                                       aaTCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATc                                                                                         <  5:11274/131‑1 (MQ=255)
                                                                       aaTCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATc                                                                                         >  6:11274/1‑131 (MQ=255)
                                                                         tCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGGTCAGGTCATGGCTGAAGAGGCt                                                                     >  3:223908/1‑149 (MQ=255)
                                                                         tCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCt                                                                     <  4:258254/149‑1 (MQ=255)
                                                                              gAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaa                                                                 >  4:303358/1‑148 (MQ=255)
                                                                              gAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaa                                                                 >  3:81098/1‑148 (MQ=255)
                                                                                      cgcgcAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGGTCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAg                                                        <  4:223908/149‑1 (MQ=255)
                                                                                      cgcgcAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAg                                                        <  4:15456/149‑1 (MQ=255)
                                                                                      cgcgcAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAg                                                        <  4:141470/149‑1 (MQ=255)
                                                                                        cgcAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGtt                                                      <  7:286167/149‑1 (MQ=255)
                                                                                        cgcAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGtt                                                      <  4:351117/149‑1 (MQ=255)
                                                                                          cagcaAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGt                                                    <  7:194383/149‑1 (MQ=255)
                                                                                             caAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                                  >  8:275493/1‑148 (MQ=255)
                                                                                                gTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAgg                                              >  3:14284/1‑149 (MQ=255)
                                                                                                 ttCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGt                                             <  1:141715/149‑1 (MQ=255)
                                                                                                    cTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGtctc                                          >  5:188548/1‑149 (MQ=255)
                                                                                                     tGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCa                                                >  1:181678/1‑142 (MQ=255)
                                                                                                     tGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCa                                                <  2:181678/142‑1 (MQ=255)
                                                                                                       agaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATTAACATGCAGCCGTCGCTGATCAGGTCATTGCTGAAGAGTCTGCAAAACAGCGAGTTGTACCAGGTctct                                         <  8:247730/147‑1 (MQ=255)
                                                                                                       agaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTAAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTg                                        <  7:52337/148‑1 (MQ=255)
                                                                                                        gaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGc                                       <  3:95359/148‑1 (MQ=255)
                                                                                                          cagCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGACGGCTACTGTGGATTAACATCCAGCCGTCGCTGATCAGGTCAAGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGACTCTGCCa                                     >  1:115546/3‑148 (MQ=255)
                                                                                                              cgtgatgCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGTATGGAATGCCGGGTAGTGTGGATGAACATGCTGCCGTCGCTGATCAGGTCATGTCTGAAGAGTCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtcg                                <  2:115546/143‑1 (MQ=255)
                                                                                                               gTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtcgt                               >  5:111642/1‑149 (MQ=255)
                                                                                                                   aGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGGATGGAAGGCCGGGTAGTGTGGATGAAGATGCAGCCGTCGCTGAGCAGGGCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCgtcgtc                           <  7:360032/149‑1 (MQ=255)
                                                                                                                    gCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCgtcgtc                           >  6:324356/1‑148 (MQ=255)
                                                                                                                    gCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTcc                          <  4:349209/149‑1 (MQ=255)
                                                                                                                           ggTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgg                    <  4:189806/148‑1 (MQ=255)
                                                                                                                            gTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                                  <  6:208478/117‑1 (MQ=255)
                                                                                                                            gTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                                  >  5:208478/1‑117 (MQ=255)
                                                                                                                            gTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtcgtc                              <  4:118936/137‑1 (MQ=255)
                                                                                                                            gTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtcgtc                              >  3:118936/1‑137 (MQ=255)
                                                                                                                            gTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgggg                  >  6:184105/1‑149 (MQ=255)
                                                                                                                            gTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgggg                  <  6:424758/149‑1 (MQ=255)
                                                                                                                            gTACGGGGTTGCGCGTCGCGCTCGAGGCCGCGAATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgggg                  <  6:292108/149‑1 (MQ=255)
                                                                                                                             tACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCggggg                 <  6:188548/149‑1 (MQ=255)
                                                                                                                               cGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATgg                                                                                <  3:272567/84‑1 (MQ=255)
                                                                                                                               cGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATgg                                                                                >  4:272567/1‑84 (MQ=255)
                                                                                                                                   tttGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATc           <  5:324356/149‑1 (MQ=255)
                                                                                                                                   tctGCGCGGCGCGCTCGCGGCCGCGCATTGAAGGTCGTGTAGTGTGGATGAACATGCAGCTGTCGCTGATCAGGTCATGTCTTAAGAGGCTGCAAAACTGCGAGTTGTACCAGGTCTCTGCGAACTCGTCGTCCAGGCGGTGGTCGAt            <  1:56025/146‑1 (MQ=255)
                                                                                                                                   catGCGCGGCCGCGAGCGCGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                                  <  4:377925/108‑1 (MQ=255)
                                                                                                                                   catGCGCGGCCGCGAGCGCGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                                  >  3:377925/3‑110 (MQ=255)
                                                                                                                                   catGCGCGGCCGCGAGCGCGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgg                    <  3:77555/138‑1 (MQ=255)
                                                                                                                                   catGCGCGGCCGCGAGCGCGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCgg                    >  4:77555/3‑140 (MQ=255)
                                                                                                                                     tGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAg         >  6:9183/1‑149 (MQ=255)
                                                                                                                                      gcgcggcgcgCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAgg        >  3:394692/1‑149 (MQ=255)
                                                                                                                                       cgcggcgcgCTCTCGTCCGCGCATGGAAGGAAGCGTAGTGTCCATGAACATGCAGCCCACCATGATCATGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGCTCTCTGCCAACACGTCGTCCAGGCGGGGGTCTAt            >  2:56025/1‑144 (MQ=255)
                                                                                                                                       cgcggcgcgCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGt       <  3:303358/149‑1 (MQ=255)
                                                                                                                                       cgcggcgcgCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGt       <  7:275493/149‑1 (MQ=255)
                                                                                                                                              gcgcTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTGGata  >  5:146900/1‑148 (MQ=255)
                                                                                                                                                gcTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGt       >  5:159826/1‑141 (MQ=255)
                                                                                                                                                gcTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGt       <  6:159826/141‑1 (MQ=255)
                                                                                                                                                  tCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATc           <  8:144792/135‑1 (MQ=255)
                                                                                                                                                  tCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATc           >  7:144792/1‑135 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCATTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGATCAGGTGGATA  >  NC_002947/5185908‑5186196

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