Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,651,601 C→T 5.5% R290H (CGT→CAT)  fba ← fructose‑bisphosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,651,6010CT 5.5% 275.3 / 3.5 90R290H (CGT→CAT) fbafructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (42/43);  new base T (3/2);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

CCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCCGCCGTACTGGTTGATGATCGCCAGC  >  NC_002947/5651454‑5651747
                                                                                                                                                   |                                                                                                                                                  
ccTTCCAGAGAGATCGGCTTGATATTCGAGGCATTGCCCGCGGTGCCGAAGGATTCGTAACTCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCAAAGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACg                                                                                                                                                   >  2:185403/1‑149 (MQ=255)
    ccAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGgc                                                                                                                                                 <  7:197235/147‑1 (MQ=255)
     cAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGgcgc                                                                                                                                               <  4:144657/148‑1 (MQ=255)
      agagagATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCAt                                                                                                                                             >  8:12432/1‑149 (MQ=255)
       gagagaTCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATg                                                                                                                                            >  6:292590/1‑149 (MQ=255)
        agagaTCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATgg                                                                                                                                           >  6:284581/1‑149 (MQ=255)
           gaTCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCAc                                                                                                                                        <  5:255544/149‑1 (MQ=255)
             tCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCg                                                                                                                                      >  5:176195/1‑149 (MQ=255)
                gCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATgg                                                                                                                                           <  4:5688/141‑1 (MQ=255)
                gCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATgg                                                                                                                                           >  3:5688/1‑141 (MQ=255)
                    gATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAgg                                                                                                                               >  5:315345/1‑149 (MQ=255)
                      tCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGc                                                                                                                              >  5:280664/1‑148 (MQ=255)
                       cTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGc                                                                                                                              >  3:268689/1‑147 (MQ=255)
                       cTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGc                                                                                                                              >  3:132593/1‑147 (MQ=255)
                        ttCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCACCATCAGACGGCGCATGGCACCGGTAGAGGc                                                                                                                              >  5:299580/1‑146 (MQ=255)
                         tCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGa                                                                                                                          >  1:191285/1‑149 (MQ=255)
                            aGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGc                                                                                                                       <  4:12027/149‑1 (MQ=255)
                              gCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAg                                                                                                                     >  1:126155/1‑149 (MQ=255)
                               cATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAgg                                                                                                                    >  3:304204/1‑149 (MQ=255)
                                    cccGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGt                                                                                                               <  2:191285/149‑1 (MQ=255)
                                     ccGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAg                                                                                                                     <  7:225351/142‑1 (MQ=255)
                                     ccGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAg                                                                                                                     >  8:225351/1‑142 (MQ=255)
                                     ccGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTa                                                                                                              <  3:76184/149‑1 (MQ=255)
                                      cgcgGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTAt                                                                                                             <  4:304204/149‑1 (MQ=255)
                                          gTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGAt                                                                                                         <  5:292590/149‑1 (MQ=255)
                                             ccGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGtt                                                                                                      <  7:12432/149‑1 (MQ=255)
                                             ccGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGtt                                                                                                      >  6:24895/1‑149 (MQ=255)
                                               gAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGAt                                                                                                         <  2:215480/144‑1 (MQ=255)
                                               gAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGAt                                                                                                         >  1:215480/1‑144 (MQ=255)
                                                   gCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCtt                                                                                                >  2:221796/1‑149 (MQ=255)
                                                   gCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCtt                                                                                                >  7:236308/1‑149 (MQ=255)
                                                   gCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCtt                                                                                                <  1:221796/149‑1 (MQ=255)
                                                     ttCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTgc                                                                                              <  2:299474/149‑1 (MQ=255)
                                                     ttCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTgc                                                                                              >  8:217412/1‑149 (MQ=255)
                                                     ttCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTgc                                                                                              >  1:63993/1‑149 (MQ=255)
                                                            cgcgcgATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGcc                                                                                       <  1:175753/149‑1 (MQ=255)
                                                               gcgATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTg                                                                                    >  5:152131/1‑149 (MQ=255)
                                                                     cacaCGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGa                                                                               <  6:176195/148‑1 (MQ=255)
                                                                     cacaCGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGAt                                                                              <  2:63993/149‑1 (MQ=255)
                                                                      acacGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTTGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATg                                                                             <  7:216934/149‑1 (MQ=255)
                                                                      acacGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATg                                                                             >  7:244946/1‑149 (MQ=255)
                                                                        acGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGcc                                                                           <  8:236308/149‑1 (MQ=255)
                                                                           tCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTc                                                                        <  7:217412/149‑1 (MQ=255)
                                                                             gcgcATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCaa                                                                      <  3:10284/149‑1 (MQ=255)
                                                                              cgcATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTgc                                                                                              <  8:180443/124‑1 (MQ=255)
                                                                              cgcATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTgc                                                                                              >  7:180443/1‑124 (MQ=255)
                                                                              cgcATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCatg                                                                     >  1:54115/1‑147 (MQ=255)
                                                                              cgcATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAAc                                                                     <  6:280664/149‑1 (MQ=255)
                                                                               gcATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAAc                                                                     >  2:93799/1‑148 (MQ=255)
                                                                               gcATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAAc                                                                     >  6:110194/1‑148 (MQ=255)
                                                                                     ccTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGTCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCATGCTTGATGCCTTCAACGATttct                                                              <  1:185403/149‑1 (MQ=255)
                                                                                      cTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGTCCTTGCGCACGCCGTGCTTGATGCCTTCAACGATttctt                                                             <  4:268689/149‑1 (MQ=255)
                                                                                         gATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGGCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCaa                                                          >  3:187882/1‑149 (MQ=255)
                                                                                         caaggTTTCAACGAAGAACTTGCGTGGGTCGACCTCGCTCGGGTTCTGCGCCATCAGACGGAGCATCGCACCGGTAGAGGCCAGACGCAGGTCGGGATCGATg                                                                                                        >  4:282308/4‑103 (MQ=255)
                                                                                          aTGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCaa                                                          <  1:93799/148‑1 (MQ=255)
                                                                                          aTGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAAc                                                         <  8:244946/149‑1 (MQ=255)
                                                                                                tCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTa                                                    <  2:13430/148‑1 (MQ=255)
                                                                                                     gaagaaCTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCgg                                                      <  4:334788/141‑1 (MQ=255)
                                                                                                     gaagaaCTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTagattg                                               >  3:334788/1‑143 (MQ=255)
                                                                                                     gaagaaCTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCg                                               >  2:148997/1‑148 (MQ=255)
                                                                                                     gaagaaCTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCg                                               >  7:19218/1‑148 (MQ=255)
                                                                                                     gaagaaCTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGt                                              <  2:124261/149‑1 (MQ=255)
                                                                                                      aagaaCTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTa                                             <  7:177342/149‑1 (MQ=255)
                                                                                                      aagaaCTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCATGATGCCTTCAACGATTTCTTCAACCGGTACACCg                                               >  7:182399/1‑147 (MQ=255)
                                                                                                      aagaaCTTGCGTGGGTCGAACTCGCTCGGGTTCTGCACCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTa                                             <  6:299580/149‑1 (MQ=255)
                                                                                                            ttGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTc                                       >  2:125918/1‑149 (MQ=255)
                                                                                                                gTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGa                                                                   >  2:275147/1‑117 (MQ=255)
                                                                                                                gTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGa                                                                   <  1:275147/117‑1 (MQ=255)
                                                                                                                gTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCttt                                    <  5:24895/148‑1 (MQ=255)
                                                                                                                  gggTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGAt                                 >  7:297816/1‑149 (MQ=255)
                                                                                                                  gggTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGAt                                 <  4:187882/149‑1 (MQ=255)
                                                                                                                       gAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTcgc                            <  8:297816/149‑1 (MQ=255)
                                                                                                                           tcgctcgGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCaa                                                                      >  4:68808/1‑103 (MQ=255)
                                                                                                                           tcgctcgGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCaa                                                                      <  3:68808/103‑1 (MQ=255)
                                                                                                                           tcgctcgGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCgccgcc                        >  7:196777/1‑149 (MQ=255)
                                                                                                                             gctcgGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGccgccg                       <  6:315345/148‑1 (MQ=255)
                                                                                                                                    ttCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCCGCCGTACTGGtt               >  7:313599/1‑149 (MQ=255)
                                                                                                                                      atGCGCCGTCTGATGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAAc                                                                     <  7:239760/92‑1 (MQ=255)
                                                                                                                                      atGCGCCGTCTGATGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAAc                                                                     >  8:239760/2‑93 (MQ=255)
                                                                                                                                      atGCGCCGTCTGATGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCCGCCGTACTGGTtga             >  1:17888/2‑149 (MQ=255)
                                                                                                                                       tGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAAc                                                         >  8:36832/1‑104 (MQ=255)
                                                                                                                                       tGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAAc                                                         <  7:36832/104‑1 (MQ=255)
                                                                                                                                       tGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCtt                                     >  6:306031/1‑124 (MQ=255)
                                                                                                                                       tGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCtt                                     <  5:306031/124‑1 (MQ=255)
                                                                                                                                         cgcCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCCGCCGTACTGGTtga             <  1:125918/146‑1 (MQ=255)
                                                                                                                                          gcCGTCTGATGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCCGCCGTACTGGTtga             <  7:3728/145‑1 (MQ=255)
                                                                                                                                          gcCGTCTGATGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCCGCCGTACTGGTtga             >  8:3728/1‑145 (MQ=255)
                                                                                                                                           ccATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCCGCCGTACTGGTTGATGATc        <  6:152131/149‑1 (MQ=255)
                                                                                                                                            cATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCAGACGTACTGGTTGATCATCg       >  2:72102/1‑149 (MQ=255)
                                                                                                                                                 gACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCCGCCGTACTGGTTGATGATCgccagc  <  8:196777/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGATGCACACGTCGCGCATGGCCTTGATGGTTTCACCGAAGAACTTGCGTGGGTCGAACTCGCTCGGGTTCTGCGCCATCAGACGGCGCATGGCACCGGTAGAGGCCAGACGCAGGTCGGTATCGATGTTGACCTTGCGCACGCCGTGCTTGATGCCTTCAACGATTTCTTCAACCGGTACACCGTAGGTTTCTTTGATGTCGCCGCCGTACTGGTTGATGATCGCCAGC  >  NC_002947/5651454‑5651747

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