Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,591,604 G→C 6.6% V106V (GTG→GTC 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,6040GC 6.6% 602.7 / 24.2 182V106V (GTG→GTCmccBmethylcrotonyl‑CoA carboxylase biotin‑containing subunit beta
Reads supporting (aligned to +/- strand):  ref base G (68/102);  new base C (4/8);  total (72/110)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 185 total aligned reads displayed.

CGCGAGCGTATCGACCGCCTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGT  >  NC_002947/4591455‑4591753
                                                                                                                                                     |                                                                                                                                                     
 gcgAGCGTATCGACCGCCTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTgg                                                                                                                                                      <  1:3864469/150‑1 (MQ=255)
     gCGTATCGACCGCCTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCaa                                                                                                                                                  <  2:1606909/150‑1 (MQ=255)
          tCGACCGCCTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAacgacg                                                                                                                                             >  2:4051909/1‑150 (MQ=255)
            gACCGCCTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGcc                                                                                                                                           >  1:3522570/1‑150 (MQ=255)
             aCCGCCTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCa                                                                                                                                          >  2:1922755/1‑150 (MQ=255)
               cGCCTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACg                                                                                                                                        <  1:2420232/150‑1 (MQ=255)
               cGCCTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCAGCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACg                                                                                                                                        >  1:2990054/1‑150 (MQ=255)
                 ccTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGt                                                                                                                                      <  1:3201680/150‑1 (MQ=255)
                 ccTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGt                                                                                                                                      >  1:2959663/1‑150 (MQ=255)
                     ctgGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTaaaa                                                                                                                                  <  2:2594237/150‑1 (MQ=255)
                     ctgGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTaaaa                                                                                                                                  <  1:4787847/150‑1 (MQ=255)
                     ctgGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTaaaa                                                                                                                                  <  2:425423/150‑1 (MQ=255)
                     ctgGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTaaaa                                                                                                                                  <  2:3086526/150‑1 (MQ=255)
                       ggACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAgg                                                                                                                                <  2:4369160/150‑1 (MQ=255)
                        gACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGc                                                                                                                               >  1:1337453/1‑150 (MQ=255)
                         aCCCGGGCTCGCCGTTCCTCGAGCTCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCg                                                                                                                              <  2:4654902/150‑1 (MQ=255)
                         aCCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCg                                                                                                                              <  2:5054802/150‑1 (MQ=255)
                         aCCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCg                                                                                                                              <  1:3690077/150‑1 (MQ=255)
                           ccGGGCTCGCCGTTCCTCGAGCTCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGt                                                                                                                            <  2:3134262/150‑1 (MQ=255)
                           ccGGGCTCGCCGATCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGt                                                                                                                            >  2:3620683/1‑150 (MQ=255)
                            cGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGtt                                                                                                                           <  2:4025084/150‑1 (MQ=255)
                            cGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGtt                                                                                                                           <  2:885167/150‑1 (MQ=255)
                                 tCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGtac                                                                                                                      >  2:3489728/1‑150 (MQ=255)
                                   gCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGtacta                                                                                                                    >  2:213036/1‑150 (MQ=255)
                                       ttCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCg                                                                                                                <  2:4866287/150‑1 (MQ=255)
                                          cTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTg                                                                                                             >  1:2971476/1‑150 (MQ=255)
                                           tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAATGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                            >  2:4399071/1‑150 (MQ=255)
                                           tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                            >  1:3173867/1‑150 (MQ=255)
                                           tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                            >  2:2122634/1‑150 (MQ=255)
                                           tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                            >  2:1115123/1‑150 (MQ=255)
                                           tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                            >  2:1237762/1‑150 (MQ=255)
                                             gagaTCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACg                                                                                                          <  2:4156181/150‑1 (MQ=255)
                                             gagaTCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACg                                                                                                          >  1:3951413/1‑150 (MQ=255)
                                              agaTCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACgg                                                                                                         >  1:4281709/1‑150 (MQ=255)
                                                ggcGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTg                                                                                                       >  2:3152132/3‑150 (MQ=255)
                                                ggcGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTg                                                                                                       >  2:595633/3‑150 (MQ=255)
                                                  cGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaa                                                                                                     <  2:2012609/150‑1 (MQ=255)
                                                  cGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaa                                                                                                     <  2:557425/150‑1 (MQ=255)
                                                  cGGCCAGCGGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaa                                                                                                     <  1:544572/150‑1 (MQ=255)
                                                    gCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaaga                                                                                                   >  2:738968/1‑150 (MQ=255)
                                                     ccAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaagaa                                                                                                  <  2:873428/150‑1 (MQ=255)
                                                     ccAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaagaa                                                                                                  >  1:4052857/1‑150 (MQ=255)
                                                       aGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGc                                                                                                <  2:5007128/150‑1 (MQ=255)
                                                         cTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCAc                                                                                              <  1:2472935/150‑1 (MQ=255)
                                                           ggCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCt                                                                                            <  1:2380621/150‑1 (MQ=255)
                                                           ggCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCt                                                                                            <  1:826512/150‑1 (MQ=255)
                                                                 ccATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGc                                                                                      >  2:4338254/1‑150 (MQ=255)
                                                                 ccATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGc                                                                                      >  1:2590505/1‑150 (MQ=255)
                                                                 ccATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGc                                                                                      >  1:988662/1‑150 (MQ=255)
                                                                  cATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGca                                                                                     >  2:4589583/1‑150 (MQ=255)
                                                                          gtACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCAt                                                                             <  2:2583985/150‑1 (MQ=255)
                                                                           tACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATc                                                                            >  1:4714743/1‑150 (MQ=255)
                                                                           tACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATc                                                                            >  2:131835/1‑150 (MQ=255)
                                                                            aCGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCg                                                                           <  1:3829265/150‑1 (MQ=255)
                                                                              ggCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGcc                                                                         <  1:343482/150‑1 (MQ=255)
                                                                                 gAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTg                                                                      >  2:1822354/1‑150 (MQ=255)
                                                                                  aaGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGc                                                                     >  1:4629442/1‑150 (MQ=255)
                                                                                        tGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAAcc                                                               <  1:1834160/150‑1 (MQ=255)
                                                                                        tGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAAcc                                                               <  1:1985112/150‑1 (MQ=255)
                                                                                                gcgggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCg                                                       >  1:52403/1‑150 (MQ=255)
                                                                                                 ccggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAACCCAACGCCCTGCAGAACCGCCTGCCGt                                                      <  2:358204/148‑1 (MQ=255)
                                                                                                    gcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCa                                                   <  1:1060923/150‑1 (MQ=255)
                                                                                                    gcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCa                                                   <  1:1060910/150‑1 (MQ=255)
                                                                                                    gcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCa                                                   >  2:1668968/1‑150 (MQ=255)
                                                                                                      gTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGCGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatc                                                 >  2:958030/1‑150 (MQ=255)
                                                                                                         atcaccGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCatttat                                              >  2:1166905/5‑150 (MQ=255)
                                                                                                         aTTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatctat                                              <  1:4491634/150‑1 (MQ=255)
                                                                                                          tcaccGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCatttatc                                             <  1:2913814/147‑1 (MQ=255)
                                                                                                           caccGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCatttatct                                            <  1:3569528/148‑1 (MQ=255)
                                                                                                            gCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtg                                           <  1:3746850/150‑1 (MQ=255)
                                                                                                            accGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCtg                                           >  2:928522/2‑150 (MQ=255)
                                                                                                             ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                          >  1:3371096/1‑150 (MQ=255)
                                                                                                             ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                          <  2:1775460/150‑1 (MQ=255)
                                                                                                             ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                          >  1:3371247/1‑150 (MQ=255)
                                                                                                             ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCATGCATCTATCtgg                                          >  2:4235937/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:362629/150‑1 (MQ=255)
                                                                                                               ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                        <  1:3149293/150‑1 (MQ=255)
                                                                                                               ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                        <  2:4625456/150‑1 (MQ=255)
                                                                                                                   tCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACt                                    <  1:1460488/150‑1 (MQ=255)
                                                                                                                           gTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCggcggc                            <  1:5003997/150‑1 (MQ=255)
                                                                                                                           gTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCggcggc                            >  1:4583568/1‑150 (MQ=255)
                                                                                                                             ggAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGCGGACTCCGGCGgcgc                          <  2:370474/150‑1 (MQ=255)
                                                                                                                                 ggCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAAc                      <  1:2160913/150‑1 (MQ=255)
                                                                                                                                 ggCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAAc                      <  1:2991430/150‑1 (MQ=255)
                                                                                                                                 ggCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAAc                      <  2:3267628/150‑1 (MQ=255)
                                                                                                                                     tGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGc                  <  2:274989/150‑1 (MQ=255)
                                                                                                                                      tgactGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGcc                 <  2:2167146/146‑1 (MQ=255)
                                                                                                                                      ggAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGcc                 <  2:1990270/150‑1 (MQ=255)
                                                                                                                                       gAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTCTCTGGTGGACTCCGGCGGCGCCAACTTGccc                >  1:4978863/1‑150 (MQ=255)
                                                                                                                                       gAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGccc                <  2:3134276/150‑1 (MQ=255)
                                                                                                                                       gAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGccc                >  1:680699/1‑150 (MQ=255)
                                                                                                                                        aaTGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGcccc               <  1:4861324/150‑1 (MQ=255)
                                                                                                                                            cATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGcccaatca           <  1:2039566/150‑6 (MQ=255)
                                                                                                                                             aTGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAg          >  2:2984099/1‑150 (MQ=255)
                                                                                                                                             aTGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCCATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGcccaatcag          >  2:4352127/1‑144 (MQ=255)
                                                                                                                                               gATCGTGGCCAACGACGCCTCGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCt                                                                                            <  2:3400941/66‑1 (MQ=255)
                                                                                                                                                   gTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGaa    <  1:512134/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGAGCGTATCGACCGCCTGCTGGACCCGGGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGT  >  NC_002947/4591455‑4591753

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: