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,7910GC100.0% 1933.8 / NA 491V51L (GTG→CTG) cobAuroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (279/212);  total (279/212)

Only 100 of 506 total aligned reads displayed.

AGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATTTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGGTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCGAAGAA  >  NC_008463/3711644‑3711937
                                                                                                                                                   |                                                                                                                                                  
aGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctg                                                                                                                                                  >  1:4780967/1‑150 (MQ=255)
 gCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctgg                                                                                                                                                 >  2:5857992/1‑150 (MQ=255)
 gCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctgg                                                                                                                                                 >  1:3114846/1‑150 (MQ=255)
  cGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctggt                                                                                                                                                <  2:4990389/150‑1 (MQ=255)
   ggTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTc                                                                                                                                               >  2:2304395/1‑150 (MQ=255)
   ggTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTc                                                                                                                                               >  1:4532644/1‑150 (MQ=255)
    gTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTcc                                                                                                                                              <  2:3646408/150‑1 (MQ=255)
    gTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTcc                                                                                                                                              <  1:2625358/150‑1 (MQ=255)
     tAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCg                                                                                                                                             <  1:8793982/150‑1 (MQ=255)
     tAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCg                                                                                                                                             <  1:7936120/150‑1 (MQ=255)
       aGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCggg                                                                                                                                           >  1:4934281/1‑150 (MQ=255)
       aGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCggg                                                                                                                                           >  2:7727920/1‑150 (MQ=255)
             ggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCa                                                                                                                                     >  2:961354/1‑150 (MQ=255)
                tggtggGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAc                                                                                                                                  <  2:10610979/150‑1 (MQ=255)
                 ggtgggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcg                                                                                                                                 >  1:5234323/1‑150 (MQ=255)
                    gggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcgg                                                                                                                              >  2:6186779/1‑150 (MQ=255)
                      gtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcggcg                                                                                                                            >  1:6960693/1‑150 (MQ=255)
                        gCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCggcggc                                                                                                                          <  1:10000485/150‑1 (MQ=255)
                         cGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCt                                                                                                                         <  1:3751048/150‑1 (MQ=255)
                           ggTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGc                                                                                                                       >  1:3260776/1‑150 (MQ=255)
                              cccGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGACCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGc                                                                                                                    <  2:6374681/150‑1 (MQ=255)
                                  gCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCca                                                                                                                >  1:1733663/1‑150 (MQ=255)
                                    gATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCcaca                                                                                                              >  1:4457766/1‑150 (MQ=255)
                                     aTCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCacac                                                                                                             >  1:6780805/1‑150 (MQ=255)
                                     aTCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCacac                                                                                                             <  1:2911850/150‑1 (MQ=255)
                                     aTCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGGGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGGCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCacac                                                                                                             >  2:2634973/1‑150 (MQ=255)
                                        cGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGc                                                                                                          >  1:8540502/1‑150 (MQ=255)
                                         ggAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCa                                                                                                         <  2:2834479/150‑1 (MQ=255)
                                         ggAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCa                                                                                                         >  1:10527111/1‑150 (MQ=255)
                                         ggAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCa                                                                                                         <  1:907679/150‑1 (MQ=255)
                                          gAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAg                                                                                                        >  1:7269108/1‑150 (MQ=255)
                                               gTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACtt                                                                                                   >  1:5471259/1‑150 (MQ=255)
                                                ttGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCTCACCGCAGGACTTc                                                                                                  >  2:5376402/1‑150 (MQ=255)
                                                 tGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCa                                                                                                 <  1:10464408/150‑1 (MQ=255)
                                                   aCCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATc                                                                                               <  2:978924/150‑1 (MQ=255)
                                                       tGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAgc                                                                                           <  2:6980428/150‑1 (MQ=255)
                                                         aaGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAgcgc                                                                                         >  1:4135510/1‑150 (MQ=255)
                                                          aGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGcc                                                                                        >  2:3567661/1‑150 (MQ=255)
                                                           ggCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCt                                                                                       >  1:704240/1‑150 (MQ=255)
                                                               gTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtg                                                                                   <  2:6096496/150‑1 (MQ=255)
                                                               gTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtg                                                                                   <  2:560171/150‑1 (MQ=255)
                                                               gTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtg                                                                                   >  1:9403921/1‑150 (MQ=255)
                                                                   gcgcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgctgc                                                                               <  2:8993156/150‑1 (MQ=255)
                                                                     gcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCgg                                                                             >  2:836275/1‑150 (MQ=255)
                                                                       ccTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCa                                                                           >  1:97478/1‑150 (MQ=255)
                                                                        cTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCAt                                                                          >  1:6615933/1‑150 (MQ=255)
                                                                         tGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATg                                                                         <  1:1034633/150‑1 (MQ=255)
                                                                           cAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGcc                                                                       <  1:5729598/150‑1 (MQ=255)
                                                                            aGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGccc                                                                      >  1:2370253/1‑150 (MQ=255)
                                                                            aGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGccc                                                                      >  2:4073088/1‑150 (MQ=255)
                                                                             ggACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCg                                                                     >  2:5679192/1‑150 (MQ=255)
                                                                                  cgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAgg                                                                <  2:10409443/150‑1 (MQ=255)
                                                                                   ggATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAggg                                                               >  1:132469/1‑150 (MQ=255)
                                                                                     aTGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGAc                                                             >  2:8782402/1‑150 (MQ=255)
                                                                                       gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGc                                                           >  1:1607459/1‑150 (MQ=255)
                                                                                       gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGc                                                           >  2:881216/1‑150 (MQ=255)
                                                                                       gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGc                                                           <  2:6562841/150‑1 (MQ=255)
                                                                                        tggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCa                                                          <  1:3984511/150‑1 (MQ=255)
                                                                                          gtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGc                                                        >  2:10692429/1‑150 (MQ=255)
                                                                                              tGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtgg                                                    <  2:3114606/150‑1 (MQ=255)
                                                                                               ggTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggt                                                   <  2:3094735/150‑1 (MQ=255)
                                                                                               ggTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggt                                                   >  2:789145/1‑150 (MQ=255)
                                                                                                gTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggtg                                                  <  2:3456148/150‑1 (MQ=255)
                                                                                                gTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggtg                                                  >  2:4684620/1‑150 (MQ=255)
                                                                                                  cgacgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTgcg                                                >  2:1344081/1‑150 (MQ=255)
                                                                                                   gacgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTgcgc                                               >  2:8699957/1‑150 (MQ=255)
                                                                                                    acgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGcc                                              >  1:1484367/1‑150 (MQ=255)
                                                                                                    acgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGcc                                              <  1:3419615/150‑1 (MQ=255)
                                                                                                     cgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCt                                             <  1:6035274/150‑1 (MQ=255)
                                                                                                      gaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTg                                            >  2:5146736/1‑150 (MQ=255)
                                                                                                       atctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGAGGTGCGCCTGa                                           >  1:9652030/4‑150 (MQ=255)
                                                                                                       atctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGGCGGCAAGCGGGGCGGCTGGCGCTCCACACCGCTGGACTTCATCCAGCGCGTGATGCTGCGGCTTGCCCGCCTGGGACGCAGCGTGGTGCGCCTgt                                           >  2:7919476/4‑149 (MQ=255)
                                                                                                        tctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGaa                                          >  2:4320175/3‑150 (MQ=255)
                                                                                                        tctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGaa                                          >  2:8210859/3‑150 (MQ=255)
                                                                                                         ctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGaaa                                         >  2:7212783/2‑150 (MQ=255)
                                                                                                            gTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggc                                      <  1:2159693/150‑1 (MQ=255)
                                                                                                              cAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggcgg                                    >  2:8595150/1‑150 (MQ=255)
                                                                                                               aaCCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggcggc                                   <  2:3326018/150‑1 (MQ=255)
                                                                                                               aaCCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggcggc                                   <  2:6848343/150‑1 (MQ=255)
                                                                                                                       cATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATg                           <  1:6294610/150‑1 (MQ=255)
                                                                                                                        aTACTGGTGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCATGACTTCATACTGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTTCGCCTGAAAGGCGTCGATCCTTGc                          <  2:9249575/150‑1 (MQ=255)
                                                                                                                         tACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCTCACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCa                         >  2:1925764/1‑150 (MQ=255)
                                                                                                                         tACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCa                         <  2:9478517/150‑1 (MQ=255)
                                                                                                                           cTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACTCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATc                       >  2:1563127/1‑150 (MQ=255)
                                                                                                                           cTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATc                       <  2:4197977/150‑1 (MQ=255)
                                                                                                                           cTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATc                       >  1:8271323/1‑150 (MQ=255)
                                                                                                                            tGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCt                      <  2:4827063/150‑1 (MQ=255)
                                                                                                                               aGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCg                   <  2:9083657/150‑1 (MQ=255)
                                                                                                                                gCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCgg                  >  1:7665044/1‑150 (MQ=255)
                                                                                                                                 cACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGc                 <  2:7376169/150‑1 (MQ=255)
                                                                                                                                 cACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGc                 <  2:1194079/150‑1 (MQ=255)
                                                                                                                                  aCTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGcc                >  1:7393344/1‑150 (MQ=255)
                                                                                                                                   cTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCcg               >  1:5838133/1‑150 (MQ=255)
                                                                                                                                    tGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCcgc              <  2:10080684/150‑1 (MQ=255)
                                                                                                                                       ccgtccgCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGcc           <  2:8019872/150‑1 (MQ=255)
                                                                                                                                             gCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCgaa     >  2:2135676/1‑150 (MQ=255)
                                                                                                                                              cccGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCgaag    >  2:2810293/1‑150 (MQ=255)
                                                                                                                                              cccGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCgaag    >  2:2238105/1‑150 (MQ=255)
                                                                                                                                               ccGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCgaaga   >  1:1954044/1‑150 (MQ=255)
                                                                                                                                                cGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCgaagaa  <  1:10649862/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATTTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGGTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCGAAGAA  >  NC_008463/3711644‑3711937

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: