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% 1913.5 / NA 488V51L (GTG→CTG) cobAuroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (284/203);  minor base A (0/1);  total (284/204)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 494 total aligned reads displayed.

GAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATTTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGGTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCGAAGAAG  >  NC_008463/3711643‑3711938
                                                                                                                                                    |                                                                                                                                                   
 aGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctgg                                                                                                                                                  >  2:7528887/1‑151 (MQ=255)
   cGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTc                                                                                                                                                >  1:10578047/1‑151 (MQ=255)
   cGGTAAGGTCTGGCTGGTGGGTGCGGGTCACTGCGATCCGGAACTGTTGACCCTGAAGGCAGTGCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGTTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTc                                                                                                                                                <  1:5355288/151‑1 (MQ=255)
     gTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCg                                                                                                                                              >  2:8659523/1‑151 (MQ=255)
      tAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCgg                                                                                                                                             >  2:8041174/1‑151 (MQ=255)
            ctggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGc                                                                                                                                       <  2:5815205/151‑1 (MQ=255)
            ctggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGc                                                                                                                                       <  1:10864349/151‑1 (MQ=255)
             tggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCa                                                                                                                                      >  1:8502678/1‑151 (MQ=255)
             tggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCa                                                                                                                                      <  1:2040982/151‑1 (MQ=255)
               gctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCaaa                                                                                                                                    >  1:6780543/1‑151 (MQ=255)
                ctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAc                                                                                                                                   >  1:7356340/1‑151 (MQ=255)
                 tggtggGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcg                                                                                                                                  <  1:7883273/151‑1 (MQ=255)
                 tggtggGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcg                                                                                                                                  >  1:3417539/1‑151 (MQ=255)
                  ggtgggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcgc                                                                                                                                 >  1:4097931/1‑151 (MQ=255)
                     gggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcggc                                                                                                                              >  2:3925048/1‑151 (MQ=255)
                        tgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGTTGTGGTGTTGGGCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCggcggc                                                                                                                           >  1:9977111/1‑151 (MQ=255)
                         gCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCt                                                                                                                          <  1:6088193/151‑1 (MQ=255)
                          cGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTg                                                                                                                         <  1:1843384/151‑1 (MQ=255)
                           gggTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGc                                                                                                                        <  2:4617684/151‑1 (MQ=255)
                                  ggCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCca                                                                                                                 >  1:10699982/1‑151 (MQ=255)
                                    cGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCcaca                                                                                                               <  2:204159/151‑1 (MQ=255)
                                    cGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCcaca                                                                                                               <  2:5963513/151‑1 (MQ=255)
                                     gATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCacac                                                                                                              >  2:4849013/1‑151 (MQ=255)
                                     gATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCacac                                                                                                              <  1:9756722/151‑1 (MQ=255)
                                     gATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCacac                                                                                                              <  1:6250576/151‑1 (MQ=255)
                                      aTCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACAcc                                                                                                             <  2:2597438/151‑1 (MQ=255)
                                       tCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCg                                                                                                            >  2:4277075/1‑151 (MQ=255)
                                         cGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCa                                                                                                          >  2:782250/1‑151 (MQ=255)
                                         cGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCa                                                                                                          <  2:5437447/151‑1 (MQ=255)
                                          ggAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAg                                                                                                         >  2:8950477/1‑151 (MQ=255)
                                          ggAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAg                                                                                                         >  2:4104065/1‑151 (MQ=255)
                                           gAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGAGGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAgg                                                                                                        >  2:6733210/1‑151 (MQ=255)
                                              cTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACt                                                                                                     >  2:849009/1‑151 (MQ=255)
                                                gTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAATCGCGGCGGCTGCCGCTCCACACCGCAGGACTTc                                                                                                   >  1:10326897/1‑151 (MQ=255)
                                                gTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTc                                                                                                   >  1:5533482/1‑151 (MQ=255)
                                                gTTGACCCGGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTc                                                                                                   <  2:9750633/151‑1 (MQ=255)
                                                    aCCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATcc                                                                                               <  1:7499503/151‑1 (MQ=255)
                                                          aaGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCCATCGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGcc                                                                                         >  2:5604003/1‑151 (MQ=255)
                                                             gCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGa                                                                                      >  1:1685402/1‑151 (MQ=255)
                                                              cAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAt                                                                                     <  2:2778342/151‑1 (MQ=255)
                                                               aGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtg                                                                                    <  1:8807511/151‑1 (MQ=255)
                                                               aGTTCGCGCCCTGCAGGACGCGGAGGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtg                                                                                    >  2:5024070/1‑151 (MQ=255)
                                                                         cTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATg                                                                          >  2:3761433/1‑151 (MQ=255)
                                                                         cTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACACGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATg                                                                          <  1:10210852/151‑1 (MQ=255)
                                                                           gCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGcc                                                                        >  2:3645377/1‑151 (MQ=255)
                                                                           gCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGcc                                                                        >  1:3623925/1‑151 (MQ=255)
                                                                                aCGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCa                                                                   >  1:9348800/1‑151 (MQ=255)
                                                                                  gcgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAgg                                                                 <  1:7135991/151‑1 (MQ=255)
                                                                                  gcgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAgg                                                                 >  1:3063515/1‑151 (MQ=255)
                                                                                  gcgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAgg                                                                 <  2:9297093/151‑1 (MQ=255)
                                                                                  gcgGATGTGGTGATGGTCGACGAACTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAgg                                                                 >  2:157736/1‑151 (MQ=255)
                                                                                   cgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAggg                                                                >  2:10578929/1‑151 (MQ=255)
                                                                                   cgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAggg                                                                >  2:1279400/1‑151 (MQ=255)
                                                                                     gATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGAc                                                              >  2:4889169/1‑151 (MQ=255)
                                                                                      aTGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACg                                                             <  2:5103661/151‑1 (MQ=255)
                                                                                       tgtgGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGc                                                            >  1:9205212/1‑151 (MQ=255)
                                                                                       tgtgGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGc                                                            >  2:10407628/1‑151 (MQ=255)
                                                                                        gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCa                                                           <  2:625773/151‑1 (MQ=255)
                                                                                        gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCa                                                           <  2:4401114/151‑1 (MQ=255)
                                                                                        gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCCCCTGATGCTGCGGCATGCCCGCCAGGGACGCa                                                           >  2:4786573/1‑151 (MQ=255)
                                                                                         tggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAg                                                          >  2:1156720/1‑151 (MQ=255)
                                                                                         tggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAg                                                          >  2:9565855/1‑151 (MQ=255)
                                                                                           gtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCg                                                        >  1:4188937/1‑151 (MQ=255)
                                                                                           gtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCg                                                        >  1:970567/1‑151 (MQ=255)
                                                                                           gtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCg                                                        >  1:3214452/1‑151 (MQ=255)
                                                                                             gATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtg                                                      <  2:1886050/151‑1 (MQ=255)
                                                                                             gATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtg                                                      >  1:8682307/1‑151 (MQ=255)
                                                                                              aTGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtgg                                                     >  1:11026790/1‑151 (MQ=255)
                                                                                              aTGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtgg                                                     >  1:1292163/1‑151 (MQ=255)
                                                                                                ggTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggtg                                                   >  1:2687387/1‑151 (MQ=255)
                                                                                                ggTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggtg                                                   >  2:9455294/1‑151 (MQ=255)
                                                                                                    gacgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGcc                                               >  2:10510043/1‑151 (MQ=255)
                                                                                                    gacgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGcc                                               >  2:10129317/1‑151 (MQ=255)
                                                                                                      cgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTg                                             >  1:5285976/1‑151 (MQ=255)
                                                                                                          ctGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAg                                         >  1:3742768/2‑151 (MQ=255)
                                                                                                           tGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAgg                                        >  1:5143166/1‑151 (MQ=255)
                                                                                                            ggTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGGCGCAGCGTGGTGCGCCTGAAAggc                                       <  1:6596613/151‑1 (MQ=255)
                                                                                                            ggTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggc                                       <  1:10004248/151‑1 (MQ=255)
                                                                                                             gTCAACCCGAGCAGACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggcg                                      <  2:649124/151‑1 (MQ=255)
                                                                                                               cAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggcggc                                    >  2:3092578/1‑151 (MQ=255)
                                                                                                               cAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAggcggc                                    <  2:8806454/151‑1 (MQ=255)
                                                                                                                  cccGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGAt                                 >  1:7821425/1‑151 (MQ=255)
                                                                                                                  cccGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGAt                                 <  2:7129836/151‑1 (MQ=255)
                                                                                                                   ccGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATc                                >  2:3316745/1‑151 (MQ=255)
                                                                                                                   ccGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATc                                >  1:5713104/1‑151 (MQ=255)
                                                                                                                    cGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAAGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATcc                               <  1:9019558/151‑1 (MQ=255)
                                                                                                                     gAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCa                              >  2:9847903/1‑151 (MQ=255)
                                                                                                                     gAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCa                              <  2:194653/151‑1 (MQ=255)
                                                                                                                           aCTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATc                        <  1:4893338/151‑1 (MQ=255)
                                                                                                                               gAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCg                    <  2:8239678/151‑1 (MQ=255)
                                                                                                                                  cACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGTTCCATGCATCTTCGGcc                 >  2:6111199/1‑151 (MQ=255)
                                                                                                                                   aCTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCcg                <  1:948712/151‑1 (MQ=255)
                                                                                                                                   aCTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCcg                >  2:589257/1‑151 (MQ=255)
                                                                                                                                   aCTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCcg                >  1:2889069/1‑151 (MQ=255)
                                                                                                                                    cTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGAATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCcgc               <  2:10924696/151‑1 (MQ=255)
                                                                                                                                     tGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCcgcg              >  1:2348016/1‑151 (MQ=255)
                                                                                                                                         cgtccgCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCgg          <  2:10250414/151‑1 (MQ=255)
                                                                                                                                          gtccgCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGc         >  1:451651/1‑151 (MQ=255)
                                                                                                                                             cgCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCgaa      <  2:3145926/151‑1 (MQ=255)
                                                                                                                                                ccGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGAAAGGCGGCGATCCATGCATCTTCGGCCGCGCCGGCgaagaa   <  2:2398248/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: