Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,591,640 C→T 13.6% Y118Y (TAC→TAT mccB → methylcrotonyl‑CoA carboxylase biotin‑containing subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,591,6400CT13.6% 468.9 / 55.1 183Y118Y (TAC→TATmccBmethylcrotonyl‑CoA carboxylase biotin‑containing subunit beta
Reads supporting (aligned to +/- strand):  ref base C (64/94);  new base T (8/17);  total (72/111)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

Only 100 of 200 total aligned reads displayed.

GTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTT  >  NC_002947/4591493‑4591789
                                                                                                                                                   |                                                                                                                                                     
 gttcTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCg                                                                                                                                                    <  2:4866099/147‑1 (MQ=255)
  tCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGc                                                                                                                                                   <  2:3838497/150‑1 (MQ=255)
   ccTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCt                                                                                                                                                  <  2:2114620/150‑1 (MQ=255)
    cTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTg                                                                                                                                                 >  1:2926794/1‑150 (MQ=255)
     tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAATGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                                                                >  2:4399071/1‑150 (MQ=255)
     tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                                                                >  2:1115123/1‑150 (MQ=255)
     tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                                                                >  2:1237762/1‑150 (MQ=255)
     tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                                                                >  2:2292575/1‑150 (MQ=255)
       gagaTCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACg                                                                                                                                              <  2:4156181/150‑1 (MQ=255)
        agaTCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACgg                                                                                                                                             >  1:4281709/1‑150 (MQ=255)
          ggcGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTg                                                                                                                                           >  2:3152132/3‑150 (MQ=255)
          ggcGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTg                                                                                                                                           >  2:595633/3‑150 (MQ=255)
           tCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTga                                                                                                                                          >  2:1951649/1‑150 (MQ=255)
            cGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaa                                                                                                                                         <  1:2975360/150‑1 (MQ=255)
            cGGCCAGCGGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaa                                                                                                                                         <  1:544572/150‑1 (MQ=255)
              gCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaaga                                                                                                                                       >  2:738968/1‑150 (MQ=255)
               ccAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaagaa                                                                                                                                      <  1:2020285/150‑1 (MQ=255)
               ccAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaagaa                                                                                                                                      <  2:873428/150‑1 (MQ=255)
                 aGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGc                                                                                                                                    <  2:5007128/150‑1 (MQ=255)
                     ggCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCt                                                                                                                                <  1:2380621/150‑1 (MQ=255)
                     ggCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCt                                                                                                                                <  2:2010104/150‑1 (MQ=255)
                         gccCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGt                                                                                                                            <  1:4553114/150‑1 (MQ=255)
                           ccATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGc                                                                                                                          >  2:4338254/1‑150 (MQ=255)
                                 ggTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAAc                                                                                                                    >  1:3189251/1‑150 (MQ=255)
                                     tACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATc                                                                                                                >  2:131835/1‑150 (MQ=255)
                                     tACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATc                                                                                                                >  1:4714743/1‑150 (MQ=255)
                                        ggCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGcc                                                                                                             <  2:1201610/150‑1 (MQ=255)
                                                  tGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAAcc                                                                                                   <  1:1985112/150‑1 (MQ=255)
                                                    ccAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGc                                                                                                 <  2:2139955/150‑1 (MQ=255)
                                                     cAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGcc                                                                                                >  1:2422465/1‑150 (MQ=255)
                                                           ccggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAACCCAACGCCCTGCAGAACCGCCTGCCGt                                                                                          <  2:358204/148‑1 (MQ=255)
                                                              gcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCa                                                                                       >  2:1668968/1‑150 (MQ=255)
                                                              gcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCa                                                                                       <  1:1060953/150‑1 (MQ=255)
                                                              gcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCa                                                                                       <  1:1060910/150‑1 (MQ=255)
                                                                gTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGCGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatc                                                                                     >  2:958030/1‑150 (MQ=255)
                                                                     tGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatctatct                                                                                <  1:3875477/150‑1 (MQ=255)
                                                                      gCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtg                                                                               <  1:3746850/150‑1 (MQ=255)
                                                                      accGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCtg                                                                               >  2:928522/2‑150 (MQ=255)
                                                                       ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                                                              <  2:1775460/150‑1 (MQ=255)
                                                                       ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                                                              >  1:3371096/1‑150 (MQ=255)
                                                                       ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCATGCATCTATCtgg                                                                              >  2:4236029/1‑150 (MQ=255)
                                                                       ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCATGCATCTATCtgg                                                                              >  2:3669283/1‑150 (MQ=255)
                                                                        cGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTTCAGAACCGCCTGCCGTGCATCTATCtggt                                                                             >  2:4186479/1‑150 (MQ=255)
                                                                         ggCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCtggtc                                                                            <  2:258614/150‑2 (MQ=255)
                                                                         ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                            <  2:4625471/150‑1 (MQ=255)
                                                                         ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                            <  2:4625405/150‑1 (MQ=255)
                                                                         ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                            <  2:362629/150‑1 (MQ=255)
                                                                         ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                            <  1:1922755/150‑1 (MQ=255)
                                                                         ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                            <  1:3671118/150‑1 (MQ=255)
                                                                         ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                            <  1:3149293/150‑1 (MQ=255)
                                                                             tCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACt                                                                        <  1:1460488/150‑1 (MQ=255)
                                                                               ggCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTcc                                                                      >  1:2522935/1‑150 (MQ=255)
                                                                               ggCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTcc                                                                      <  2:2959663/150‑1 (MQ=255)
                                                                                  cgcgTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCcggc                                                                   <  1:809913/150‑1 (MQ=255)
                                                                                     gTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCggcggc                                                                <  1:5003997/150‑1 (MQ=255)
                                                                                       ggAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGCGGACTCCGGCGgcgc                                                              <  2:370474/150‑1 (MQ=255)
                                                                                           ggCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAAc                                                          <  2:3267628/150‑1 (MQ=255)
                                                                                              gTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCACCTTg                                                       >  1:4378327/1‑150 (MQ=255)
                                                                                               tGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGc                                                      <  2:274989/150‑1 (MQ=255)
                                                                                                 gAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTCTCTGGTGGACTCCGGCGGCGCCAACTTGccc                                                    >  1:4978863/1‑150 (MQ=255)
                                                                                                 gAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGccc                                                    >  1:680551/1‑150 (MQ=255)
                                                                                                 gAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGccc                                                    <  2:3134276/150‑1 (MQ=255)
                                                                                                  aaTGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGcccc                                                   <  1:4861324/150‑1 (MQ=255)
                                                                                                     gCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGcc                                                <  1:3158445/150‑1 (MQ=255)
                                                                                                       aTGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAg                                              >  2:2984099/1‑150 (MQ=255)
                                                                                                       aTGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCCATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGcccaatcag                                              >  2:4352127/1‑144 (MQ=255)
                                                                                                         gATCGTGGCCAACGACGCCTCGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCt                                                                                                                                <  2:3400941/66‑1 (MQ=255)
                                                                                                         gATCGTGGCCAACGACGCCCCGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCt                                                                                                                                >  1:3400941/1‑66 (MQ=255)
                                                                                                           tcgtctgcAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACg                                          >  1:2908130/8‑150 (MQ=255)
                                                                                                             gTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGaa                                        <  1:512134/150‑1 (MQ=255)
                                                                                                                 gcAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGgttt                                    >  1:3532100/2‑148 (MQ=255)
                                                                                                                  cAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTCtt                                   >  2:3291950/1‑150 (MQ=255)
                                                                                                                  cAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGtt                                   <  1:2122634/150‑1 (MQ=255)
                                                                                                                  cAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGtt                                   >  1:1730302/1‑150 (MQ=255)
                                                                                                                   aaCGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGttttt                                  <  2:858398/150‑2 (MQ=255)
                                                                                                                   aaCGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTc                                  >  2:1623711/1‑150 (MQ=255)
                                                                                                                   aaCGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTc                                  <  1:4589583/150‑1 (MQ=255)
                                                                                                                        cgCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGa                             <  1:213036/150‑1 (MQ=255)
                                                                                                                         gCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGAc                            >  2:1036860/1‑150 (MQ=255)
                                                                                                                             cGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTTGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATcgcg                        <  2:1657951/150‑1 (MQ=255)
                                                                                                                              ggTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGa                       <  1:3831610/150‑1 (MQ=255)
                                                                                                                              ggTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGa                       >  1:379051/1‑150 (MQ=255)
                                                                                                                              ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                       <  1:4051909/150‑1 (MQ=255)
                                                                                                                              ggTAAAAGGCGGTTCGTACTACCCGCTGACGGCGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                       <  1:958030/150‑1 (MQ=255)
                                                                                                                                 aaaaGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGcg                    >  1:3587916/1‑149 (MQ=255)
                                                                                                                                  aaaGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCcgg                                                                    >  1:2779206/1‑101 (MQ=255)
                                                                                                                                       cGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCgg              <  2:40776/150‑1 (MQ=255)
                                                                                                                                       cGGTTCGTACGACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCgg              <  2:2990143/150‑1 (MQ=255)
                                                                                                                                        ggTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCggg             >  2:2780863/1‑150 (MQ=255)
                                                                                                                                         gTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGgc            >  1:872426/1‑150 (MQ=255)
                                                                                                                                         gTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGgc            >  1:22833/1‑150 (MQ=255)
                                                                                                                                         gTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGgc            <  2:527189/150‑1 (MQ=255)
                                                                                                                                         gTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGgc            >  1:3002916/1‑150 (MQ=255)
                                                                                                                                              tactacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGcc                                                                                            >  2:1430856/1‑65 (MQ=255)
                                                                                                                                                 tacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCa                    >  1:4650224/1‑134 (MQ=255)
                                                                                                                                                 tacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttc    >  1:1383363/1‑150 (MQ=255)
                                                                                                                                                   cccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttctt  <  1:917125/150‑1 (MQ=255)
                                                                                                                                                   cccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttctt  >  2:3108134/1‑150 (MQ=255)
                                                                                                                                                   cccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttctt  >  2:553543/1‑150 (MQ=255)
                                                                                                                                                   cccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttctt  <  1:3489728/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTT  >  NC_002947/4591493‑4591789

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