Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,711,791 G→C V51L (GTG→CTG)  cobA → uroporphyrinogen‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,711,7910GC99.8% 2295.7 / NA 582V51L (GTG→CTG) cobAuroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (346/235);  minor base A (0/1);  total (346/236)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 594 total aligned reads displayed.

GAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATTTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGGTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCGAAGAAG  >  NC_008463/3711643‑3711938
                                                                                                                                                    |                                                                                                                                                   
gAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctg                                                                                                                                                   <  2:5448183/151‑1 (MQ=255)
gAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctg                                                                                                                                                   <  1:10724458/151‑1 (MQ=255)
    ggTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTcc                                                                                                                                               <  1:9676113/151‑1 (MQ=255)
      tAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCgg                                                                                                                                             >  1:208396/1‑151 (MQ=255)
       acggTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCggg                                                                                                                                            >  1:10919733/3‑151 (MQ=255)
       aaGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCggg                                                                                                                                            >  2:6636409/1‑151 (MQ=255)
        aGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGt                                                                                                                                           >  1:9365170/1‑151 (MQ=255)
        aGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGt                                                                                                                                           >  1:5571106/1‑151 (MQ=255)
        aGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGt                                                                                                                                           >  1:11400778/1‑151 (MQ=255)
         tgTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTc                                                                                                                                          <  2:7655923/150‑1 (MQ=255)
          gTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCg                                                                                                                                         >  2:3093971/1‑151 (MQ=255)
             tggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCa                                                                                                                                      >  2:11427587/1‑151 (MQ=255)
              ggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCaa                                                                                                                                     <  2:7860740/151‑1 (MQ=255)
              ggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCaa                                                                                                                                     >  1:3003957/1‑151 (MQ=255)
              ggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCaa                                                                                                                                     >  1:2870852/1‑151 (MQ=255)
               gctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCaaa                                                                                                                                    <  1:467661/151‑1 (MQ=255)
                ctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAc                                                                                                                                   <  1:7580131/151‑1 (MQ=255)
                 tggtggGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcg                                                                                                                                  >  1:1619170/1‑151 (MQ=255)
                   gtgggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcgcg                                                                                                                                <  2:8766526/151‑1 (MQ=255)
                    tgggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcgg                                                                                                                               <  2:1726433/151‑1 (MQ=255)
                       gtgCGGGTCCCGGCTATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGAGATGGTCGACGATCTGGTCAACCCGAGCATACTGTAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGcggcgg                                                                                                                            >  2:3512887/1‑151 (MQ=255)
                                ccGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTc                                                                                                                   <  1:1223098/151‑1 (MQ=255)
                                 cGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTcc                                                                                                                  >  1:6649871/1‑151 (MQ=255)
                                  ggCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCca                                                                                                                 <  1:7883104/151‑1 (MQ=255)
                                  ggCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCca                                                                                                                 >  1:3693330/1‑151 (MQ=255)
                                                gTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTc                                                                                                   <  2:10821946/151‑1 (MQ=255)
                                                gTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTc                                                                                                   <  2:1597533/151‑1 (MQ=255)
                                                gTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTc                                                                                                   >  2:1358890/1‑151 (MQ=255)
                                                  tGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCAt                                                                                                 >  1:8531750/1‑151 (MQ=255)
                                                   gACGCTGAAGGCAGTTCGCGCCCTGCAGGACGCTGATGTTGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTGCGCCCGGATGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATc                                                                                                <  2:6378303/151‑1 (MQ=255)
                                                    aCCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATcc                                                                                               <  2:6749839/151‑1 (MQ=255)
                                                     cccTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCa                                                                                              <  1:11173523/151‑1 (MQ=255)
                                                     cccTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCa                                                                                              >  2:268195/1‑151 (MQ=255)
                                                        tGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAgcg                                                                                           <  2:3364910/151‑1 (MQ=255)
                                                         gAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAgcgc                                                                                          <  2:5854735/151‑1 (MQ=255)
                                                             gCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCCCACCGCAGGACTTCATCCAGCGCCTGa                                                                                      <  1:826231/151‑1 (MQ=255)
                                                              cAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAt                                                                                     <  2:4745647/151‑1 (MQ=255)
                                                               aGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtg                                                                                    >  2:4931718/1‑151 (MQ=255)
                                                               aGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtg                                                                                    <  2:10632689/151‑1 (MQ=255)
                                                               aGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtg                                                                                    >  1:3266491/1‑151 (MQ=255)
                                                                  tCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgctg                                                                                 <  2:426089/151‑1 (MQ=255)
                                                                   cgcgCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgctgc                                                                                >  2:5002814/1‑151 (MQ=255)
                                                                    gcgcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCg                                                                               >  1:6124328/1‑151 (MQ=255)
                                                                     cgcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCgg                                                                              >  2:8256255/1‑151 (MQ=255)
                                                                      tcccTTCAGGACGCGGTTTTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGTCTTGTCCGGGTCGTCAAACGCGGCGTCTGCCGCTCCACTCCGCTGTACTTCATCCAGCGCCTGATGCTGCggt                                                                             <  2:6089234/150‑2 (MQ=255)
                                                                      gcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGc                                                                             >  2:7244529/1‑151 (MQ=255)
                                                                      gcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGc                                                                             >  2:8609800/1‑151 (MQ=255)
                                                                            cAGGACGCGGATGTGGTGATGTTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGccc                                                                       <  2:1514219/151‑1 (MQ=255)
                                                                               gACGAGGATGTGGTTATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCAGGCTGGTCCGGGTCGGCAAACGCGCCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGcc                                                                    >  2:7474343/1‑151 (MQ=255)
                                                                                aCGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCa                                                                   >  1:2781407/1‑151 (MQ=255)
                                                                                 cgcgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAg                                                                  <  1:10201448/151‑1 (MQ=255)
                                                                                  gcgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAgg                                                                 >  2:2777766/1‑151 (MQ=255)
                                                                                  gcgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAgg                                                                 >  1:5366142/1‑151 (MQ=255)
                                                                                   cgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAggg                                                                <  1:10476270/151‑1 (MQ=255)
                                                                                    ggATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGa                                                               <  2:10071593/151‑1 (MQ=255)
                                                                                     gATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGAc                                                              >  2:7982699/1‑151 (MQ=255)
                                                                                     gATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGAc                                                              <  1:6420920/151‑1 (MQ=255)
                                                                                        gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCa                                                           <  2:2093105/151‑1 (MQ=255)
                                                                                         tggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAg                                                          >  1:2345175/1‑151 (MQ=255)
                                                                                          ggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGc                                                         <  2:9172262/151‑1 (MQ=255)
                                                                                           gtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCg                                                        >  2:1162869/1‑151 (MQ=255)
                                                                                            tgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgt                                                       >  2:240248/1‑151 (MQ=255)
                                                                                            tgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgt                                                       >  1:3475163/1‑151 (MQ=255)
                                                                                              aTGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtgg                                                     >  1:1812177/1‑151 (MQ=255)
                                                                                               tGGTCGACGATCTGGTCAACCCGAGCATCCTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggt                                                    >  2:10503776/1‑151 (MQ=255)
                                                                                                ggTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggtg                                                   <  2:8700827/151‑1 (MQ=255)
                                                                                                    gacgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGcc                                               >  1:11506296/1‑151 (MQ=255)
                                                                                                    gacgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGcc                                               >  2:755938/1‑151 (MQ=255)
                                                                                                     acgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCt                                              >  2:3819706/1‑151 (MQ=255)
                                                                                                      cgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTg                                             >  2:2205100/1‑151 (MQ=255)
                                                                                                      cgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTg                                             >  2:9595097/1‑151 (MQ=255)
                                                                                                       gaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGa                                            >  1:6926864/1‑151 (MQ=255)
                                                                                                       gaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGa                                            >  1:8826432/1‑151 (MQ=255)
                                                                                                        atctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGaa                                           <  1:8697288/148‑1 (MQ=255)
                                                                                                        atctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGaa                                           >  1:8093184/4‑151 (MQ=255)
                                                                                                        atctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGaa                                           >  2:4449817/4‑151 (MQ=255)
                                                                                                         tctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGaaa                                          <  1:1409826/149‑1 (MQ=255)
                                                                                                         tctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGaaa                                          <  2:9305041/149‑1 (MQ=255)
                                                                                                          ctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAg                                         <  1:4330048/150‑1 (MQ=255)
                                                                                                          ctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAg                                         >  2:7747548/2‑151 (MQ=255)
                                                                                                            ggTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggc                                       >  2:4072357/1‑151 (MQ=255)
                                                                                                             gTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCTCACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGTAAggcg                                      <  2:3644717/151‑1 (MQ=255)
                                                                                                             gTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggcg                                      <  2:5298784/151‑1 (MQ=255)
                                                                                                                aaCCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGgcggcg                                   >  1:5867097/1‑151 (MQ=255)
                                                                                                                    cGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATcc                               >  2:10288141/1‑151 (MQ=255)
                                                                                                                     gAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCa                              >  2:11044461/1‑151 (MQ=255)
                                                                                                                        cATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGc                           >  2:1886375/1‑151 (MQ=255)
                                                                                                                               gAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCg                    <  1:2230223/151‑1 (MQ=255)
                                                                                                                                aGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCgg                   >  1:10037068/1‑151 (MQ=255)
                                                                                                                                  cACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGcc                 <  2:9702731/151‑1 (MQ=255)
                                                                                                                                  cACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGcc                 <  1:5679474/151‑1 (MQ=255)
                                                                                                                                  cACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGcc                 >  2:6995164/1‑151 (MQ=255)
                                                                                                                                   aCTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCcg                >  1:9170379/1‑151 (MQ=255)
                                                                                                                                     tGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCcgcg              >  1:2473661/1‑151 (MQ=255)
                                                                                                                                        ccgtccgCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCg           >  1:4059040/1‑151 (MQ=255)
                                                                                                                                          gtccgCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGc         >  1:9702314/1‑151 (MQ=255)
                                                                                                                                           tccgCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCg        <  1:5431233/151‑1 (MQ=255)
                                                                                                                                              gCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCgaag     <  2:11184751/151‑1 (MQ=255)
                                                                                                                                                ccGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCTTGCACGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCTATCCATGCATCTTCGGCCGCGCCGGCgaagaa   <  2:4651119/151‑1 (MQ=255)
                                                                                                                                                ccGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCgaagaa   <  1:5190277/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATTTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGGTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCGAAGAAG  >  NC_008463/3711643‑3711938

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: