Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,591,700 G→C 15.2% P138P (CCG→CCC 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,7000GC15.2% 575.0 / 85.8 210P138P (CCG→CCCmccBmethylcrotonyl‑CoA carboxylase biotin‑containing subunit beta
Reads supporting (aligned to +/- strand):  ref base G (94/84);  new base C (18/14);  total (112/98)
Fisher's exact test for biased strand distribution p-value = 8.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 212 total aligned reads displayed.

CGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTG  >  NC_002947/4591552‑4591849
                                                                                                                                                    |                                                                                                                                                     
cgggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGt                                                                                                                                                      <  1:5368733/150‑1 (MQ=255)
cgggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGt                                                                                                                                                      <  1:5111815/150‑1 (MQ=255)
  ggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACACACCATCGCCCTGCAGAACCGCCTGCCGTGc                                                                                                                                                    >  2:2941799/1‑150 (MQ=255)
  ggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGc                                                                                                                                                    <  2:3751336/150‑1 (MQ=255)
       gATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatcta                                                                                                                                               >  1:5178404/1‑150 (MQ=255)
        atcaccGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCatttat                                                                                                                                              >  2:742042/5‑150 (MQ=255)
        atcaccGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCatttat                                                                                                                                              >  1:2111165/5‑150 (MQ=255)
            ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                                                                                                                          >  2:3622930/1‑150 (MQ=255)
            ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                                                                                                                          >  2:2034717/1‑150 (MQ=255)
             cGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCtggt                                                                                                                                         <  2:5268315/150‑1 (MQ=255)
              ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                                                                                        <  1:2154498/150‑1 (MQ=255)
               gCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCtggtcg                                                                                                                                       >  1:4384856/1‑148 (MQ=255)
               gCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtgg                                                                                                                                       >  1:32587/1‑150 (MQ=255)
                       cgcgTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                                                                                        >  2:1872437/1‑141 (MQ=255)
                             gAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGcc                                                                                                                         <  2:437543/150‑1 (MQ=255)
                              aaGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCa                                                                                                                        <  2:1334602/150‑1 (MQ=255)
                               aGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCaa                                                                                                                       <  2:4483930/150‑1 (MQ=255)
                               aGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCaa                                                                                                                       <  2:3788477/150‑1 (MQ=255)
                                ggCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAAc                                                                                                                      <  1:4448360/150‑1 (MQ=255)
                                  cGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACtt                                                                                                                    >  2:5123353/1‑150 (MQ=255)
                                  cGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACtt                                                                                                                    >  1:2747064/1‑150 (MQ=255)
                                    tGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatc                                                                                                                                                 >  1:110789/1‑119 (MQ=255)
                                    tGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGc                                                                                                                  <  1:2475313/150‑1 (MQ=255)
                                    tGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGc                                                                                                                  <  1:3182510/150‑1 (MQ=255)
                                    tGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGc                                                                                                                  >  2:3438545/1‑150 (MQ=255)
                                    tGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGc                                                                                                                  <  1:149566/150‑1 (MQ=255)
                                    gtgactGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGATAACCGACTGCCCTGCATTTATCTGGTCGATTCCGGAGGCGCGAACCTGc                                                                                                                  <  1:3848079/145‑1 (MQ=255)
                                         tGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGcccaat                                                                                                             >  1:1978902/1‑147 (MQ=255)
                                           cATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtg                                                                                                                                           <  1:2200138/118‑1 (MQ=255)
                                           cATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtg                                                                                                                                           >  2:2200128/1‑118 (MQ=255)
                                            aTGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGgcgc                                                                                                                          <  2:1995988/134‑1 (MQ=255)
                                                tCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACg                                                                                                      >  2:1197597/1‑150 (MQ=255)
                                                tCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACg                                                                                                      >  2:3955634/1‑150 (MQ=255)
                                                  gTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAACCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGaa                                                                                                    <  2:2428605/150‑1 (MQ=255)
                                                  gTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGCACCGCCTGCCGTGCAGCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGac                                                                                                    >  2:1683359/1‑149 (MQ=255)
                                                   tGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAg                                                                                                   >  2:1926760/1‑150 (MQ=255)
                                                    ggCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                                                                                        >  2:2561405/1‑112 (MQ=255)
                                                      gcAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGCCGGCGCGAACCTGCCCAATCAGGACGAGgttt                                                                                                >  2:3301760/2‑148 (MQ=255)
                                                        aaCGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGttttt                                                                                              <  2:4067895/150‑2 (MQ=255)
                                                        aaCGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTc                                                                                              <  1:5225741/150‑1 (MQ=255)
                                                            atgcaaCGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTtccgg                                                                                          <  2:2329331/145‑3 (MQ=255)
                                                              gCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGAc                                                                                        >  2:2527224/1‑150 (MQ=255)
                                                                aaCGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGTGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATcg                                                                                      >  1:1416305/2‑150 (MQ=255)
                                                                  cGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATcgcg                                                                                    >  1:2547324/1‑150 (MQ=255)
                                                                  cGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACcgcg                                                                                    >  1:2672075/1‑150 (MQ=255)
                                                                   ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGTATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                                                                                   <  2:1428642/150‑1 (MQ=255)
                                                                     tAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGc                                                                                 >  2:3819518/1‑150 (MQ=255)
                                                                        aaGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCAtt                                                                              >  2:4096295/1‑150 (MQ=255)
                                                                          ggCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTc                                                                            <  1:3523912/150‑1 (MQ=255)
                                                                           gCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCg                                                                           <  2:5281155/150‑1 (MQ=255)
                                                                            cGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCgg                                                                          <  2:5080356/150‑1 (MQ=255)
                                                                            cGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCgg                                                                          <  1:3819518/150‑1 (MQ=255)
                                                                                   tactacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAtt                                                                   >  1:4857775/1‑150 (MQ=255)
                                                                                   tactacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAtt                                                                   >  2:4240449/1‑150 (MQ=255)
                                                                                    actacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAttt                                                                  >  1:1510495/1‑150 (MQ=255)
                                                                                      tacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttc                                                                >  1:1232250/1‑150 (MQ=255)
                                                                                       acCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttct                                                               >  2:1526872/1‑150 (MQ=255)
                                                                                        cccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttctt                                                              <  1:3747965/150‑1 (MQ=255)
                                                                                          cGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCa                                                            <  1:912128/150‑1 (MQ=255)
                                                                                           gCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCaa                                                           <  1:3994865/150‑1 (MQ=255)
                                                                                             tcaCGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGATCATTTCGGGCGCATCTTCTACaatc                                                         <  2:280898/148‑3 (MQ=255)
                                                                                             tGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAAcc                                                         <  1:4096295/150‑1 (MQ=255)
                                                                                              gACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCa                                                        <  2:1107333/150‑1 (MQ=255)
                                                                                               aCGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAg                                                       <  2:4851366/150‑1 (MQ=255)
                                                                                                 ggTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGATCATTTCGGGCGCATCTTCTACAATCAGGc                                                     >  2:315423/1‑150 (MQ=255)
                                                                                                 ggTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGATCATTTCGGGCGCATCTTCTACAATCAGGc                                                     >  2:4094760/1‑150 (MQ=255)
                                                                                                  gTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGcc                                                    <  1:3126035/150‑1 (MQ=255)
                                                                                                    gaagaaGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCaa                                                  <  1:1847574/150‑1 (MQ=255)
                                                                                                    gaagaaGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCaa                                                  <  1:3433415/150‑1 (MQ=255)
                                                                                                     aagaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAAc                                                 <  1:5394133/150‑1 (MQ=255)
                                                                                                       gaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACAt                                               >  1:702478/1‑150 (MQ=255)
                                                                                                       gaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACAt                                               >  2:1045480/1‑150 (MQ=255)
                                                                                                       gaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACAt                                               >  1:2371854/1‑150 (MQ=255)
                                                                                                         agCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGa                                             <  1:3916243/150‑1 (MQ=255)
                                                                                                            aCCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAgcg                                          >  2:2626158/1‑150 (MQ=255)
                                                                                                             ccTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAgcgc                                         <  1:1679305/150‑1 (MQ=255)
                                                                                                               tGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAgcgcgc                                       >  2:5185878/1‑150 (MQ=255)
                                                                                                                 cGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCaa                                     <  1:2610695/150‑1 (MQ=255)
                                                                                                                 cGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCaa                                     <  1:3447915/150‑1 (MQ=255)
                                                                                                                  gTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAg                                    <  1:1278093/150‑1 (MQ=255)
                                                                                                                   tGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAgg                                   >  1:4986343/1‑150 (MQ=255)
                                                                                                                    gCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGc                                  <  2:1891967/150‑1 (MQ=255)
                                                                                                                     cacaAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCa                                 >  1:295428/1‑150 (MQ=255)
                                                                                                                      acaAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGTGCGCGCAAGGCAt                                >  2:3704301/1‑150 (MQ=255)
                                                                                                                      acaAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCAt                                >  1:3331020/1‑150 (MQ=255)
                                                                                                                      acaAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCAt                                >  2:3584927/1‑150 (MQ=255)
                                                                                                                        aaaCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCCTTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATcc                              >  1:4799752/1‑150 (MQ=255)
                                                                                                                         aaCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATccc                             <  2:5381509/150‑1 (MQ=255)
                                                                                                                          aCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCg                            >  2:4581914/1‑150 (MQ=255)
                                                                                                                           ccATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGc                           <  1:1197597/150‑1 (MQ=255)
                                                                                                                             aTCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAg                         >  1:1590911/1‑150 (MQ=255)
                                                                                                                              tCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGa                        <  2:1463987/150‑1 (MQ=255)
                                                                                                                                    tGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCg                  <  1:2179829/150‑1 (MQ=255)
                                                                                                                                    tGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCg                  <  1:3052275/150‑1 (MQ=255)
                                                                                                                                      cAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCgtg                >  1:2627127/1‑150 (MQ=255)
                                                                                                                                        gAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCgtggt              >  1:1876800/1‑150 (MQ=255)
                                                                                                                                        gAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCgtggt              <  1:3650784/150‑1 (MQ=255)
                                                                                                                                                tGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATGTTCTGCAACCAGGCCAACATGAGCGCGCCAGGCATCCCGCAGATCGCCGGGGTGATGGGcg      >  2:4520349/1‑149 (MQ=255)
                                                                                                                                                 gCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTc     >  1:2002173/1‑150 (MQ=255)
                                                                                                                                                 gCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTc     >  2:10450/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCGTGGTGATGGGCTCGTG  >  NC_002947/4591552‑4591849

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: