Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,711,749 T→C L37L (TTG→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,7490TC99.8% 1809.7 / NA 464L37L (TTG→CTG) cobAuroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (269/194);  minor base G (1/0);  total (270/194)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 464 total aligned reads displayed.

CCCCCCGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATTTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGGTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGA  >  NC_008463/3711601‑3711896
                                                                                                                                                    |                                                                                                                                                   
ccccccGGCCCATGGCGGCGAATGCCTGAGGAGAACGCCGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTg                                                                                                                                                   <  1:9116974/151‑1 (MQ=255)
 cccccGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTgg                                                                                                                                                  >  1:733101/1‑151 (MQ=255)
  ccccGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGt                                                                                                                                                 <  2:5676961/151‑1 (MQ=255)
  ccccGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGt                                                                                                                                                 >  1:422576/1‑151 (MQ=255)
  ccccGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGt                                                                                                                                                 <  1:4154838/151‑1 (MQ=255)
    ccGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCa                                                                                                                                               <  1:236923/151‑1 (MQ=255)
         ccATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCg                                                                                                                                          >  1:3991683/1‑151 (MQ=255)
           aTGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAg                                                                                                                                        >  2:3211416/1‑151 (MQ=255)
           aTGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAg                                                                                                                                        <  1:6754913/151‑1 (MQ=255)
           aTGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAg                                                                                                                                        >  1:9338244/1‑151 (MQ=255)
           aTGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAg                                                                                                                                        >  2:1028216/1‑151 (MQ=255)
              gcggcgAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCAt                                                                                                                                     >  1:8044792/1‑151 (MQ=255)
              gcggcgAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCAt                                                                                                                                     <  1:9025470/151‑1 (MQ=255)
                   gAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTgg                                                                                                                                >  2:8775771/1‑151 (MQ=255)
                    aaTGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGa                                                                                                                               <  2:8764447/151‑1 (MQ=255)
                      tGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGc                                                                                                                             >  2:132291/1‑151 (MQ=255)
                        ccTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCAc                                                                                                                           <  1:5967405/151‑1 (MQ=255)
                         cggaggagCACGACGAGGAGCGGTACGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCACCCCGAGCATACTGGAGCACt                                                                                                                          <  2:2415244/149‑1 (MQ=255)
                               agaACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTgtccgt                                                                                                                    <  2:3818778/151‑1 (MQ=255)
                                  acgacgATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGtccgtccg                                                                                                                 <  2:8570494/151‑1 (MQ=255)
                                   cgacgaTGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGc                                                                                                                >  2:817815/1‑151 (MQ=255)
                                   cgacgaTGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGc                                                                                                                >  2:4471854/1‑151 (MQ=255)
                                   cgacgaTGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGc                                                                                                                >  2:7649310/1‑151 (MQ=255)
                                     acgaTGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGccc                                                                                                              <  2:8976825/151‑1 (MQ=255)
                                         tGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCggct                                                                                                          <  2:3337109/151‑3 (MQ=255)
                                         tGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGGTGACCCTGAAGCCAGTTCGCGCCCTGCAGGACCCGGATGTCGAGATGGTCGACGATCTGGTCAACCCGCGCATAATGGAGCACTGTCCGGCCGCCCggct                                                                                                          >  2:4190541/1‑149 (MQ=255)
                                           aGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctgg                                                                                                        >  1:6353523/1‑151 (MQ=255)
                                           aGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctgg                                                                                                        <  2:9743200/151‑1 (MQ=255)
                                              ggTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTcc                                                                                                     <  2:7521903/151‑1 (MQ=255)
                                                tAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCgg                                                                                                   >  1:9700824/1‑151 (MQ=255)
                                                 aaGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCggg                                                                                                  >  2:931666/1‑151 (MQ=255)
                                                 aaGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCggg                                                                                                  >  1:1572416/1‑151 (MQ=255)
                                                 aaGGTCTGGCTGGTGGGTGCGGGGCCCCGCGATCCGGACCTGTTGCCCCAGAAGGCAGGTCGCGCCCCGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCggg                                                                                                  <  2:3951911/151‑1 (MQ=255)
                                                  aGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGt                                                                                                 >  1:4776010/1‑151 (MQ=255)
                                                   ggTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGACCTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTc                                                                                                >  2:5313459/1‑151 (MQ=255)
                                                   ggTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTc                                                                                                >  2:1204146/1‑151 (MQ=255)
                                                   ggTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTc                                                                                                >  1:5510176/1‑151 (MQ=255)
                                                      ctggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGc                                                                                             >  1:9983129/1‑151 (MQ=255)
                                                           tggtggGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcg                                                                                        >  1:2792550/1‑151 (MQ=255)
                                                           tggtggGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcg                                                                                        >  1:10530674/1‑151 (MQ=255)
                                                             gtgggtgCGGGTCCCGGCGATCCGGAACTGATGACCCTGAAGGCAGTTCGCGCCCAGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCAAACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcgcg                                                                                      <  2:1462069/151‑1 (MQ=255)
                                                              tgggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcgg                                                                                     <  1:7558920/151‑1 (MQ=255)
                                                               gggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGACGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcggc                                                                                    <  1:3522595/151‑1 (MQ=255)
                                                                ggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcggcg                                                                                   <  2:8010988/151‑1 (MQ=255)
                                                                 gtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGcggcgg                                                                                  >  1:8225061/1‑151 (MQ=255)
                                                                    cGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTg                                                                               >  1:8748878/1‑151 (MQ=255)
                                                                      ggTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGcc                                                                             <  2:2010694/151‑1 (MQ=255)
                                                                      ggTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGcc                                                                             >  1:8135920/1‑151 (MQ=255)
                                                                         cccGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCt                                                                          >  1:9545234/1‑151 (MQ=255)
                                                                          ccGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTc                                                                         >  2:5841564/1‑151 (MQ=255)
                                                                            ggCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCca                                                                       >  2:10398274/1‑151 (MQ=255)
                                                                              cGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCcaca                                                                     >  1:10012584/1‑151 (MQ=255)
                                                                              cGAGCCGGAACTGTTGACCCTGCAGGCAGTTCGCGCCCGGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCcaca                                                                     <  2:5458288/151‑1 (MQ=255)
                                                                                    ggAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAg                                                               <  1:4295187/151‑1 (MQ=255)
                                                                                    ggAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAg                                                               >  2:4974329/1‑151 (MQ=255)
                                                                                     gAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAgg                                                              >  1:925348/1‑151 (MQ=255)
                                                                                      aaCTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGa                                                             <  2:6861683/151‑1 (MQ=255)
                                                                                        cTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAATCGCGGCGGCTGCCGCTCCACACCGCAGGACt                                                           >  1:4543155/1‑151 (MQ=255)
                                                                                            tGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCAt                                                       >  1:7132904/1‑151 (MQ=255)
                                                                                                ccTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAg                                                   <  1:5808151/151‑1 (MQ=255)
                                                                                                 cTGTAGGCAGTTCGCGCCCTGCCGGAGGCGGCTGTGGTGATGGTCGACGCTCTGGTCAACCCGAGCATAGTGGAGCTCTGTCCGTCGGCGCGGCTGGTCCGGGTCGGCGACCGCGGCGGCTGCCGCTCCACACCGCAGGAGTTCATCCAgc                                                  >  2:6592579/1‑151 (MQ=255)
                                                                                                   gAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAgcgc                                                >  2:6038737/1‑151 (MQ=255)
                                                                                                    aaGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGcc                                               >  2:3495219/1‑151 (MQ=255)
                                                                                                    aaGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGcc                                               >  2:3932909/1‑151 (MQ=255)
                                                                                                        cAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAt                                           >  1:2698811/1‑151 (MQ=255)
                                                                                                        cAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAt                                           >  2:10053937/1‑151 (MQ=255)
                                                                                                           ttCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgct                                        <  1:2989265/151‑1 (MQ=255)
                                                                                                           ttCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgct                                        <  2:7040497/151‑1 (MQ=255)
                                                                                                             cgcgCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgctgc                                      <  1:1782739/151‑1 (MQ=255)
                                                                                                               cgcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCTACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCgg                                    <  2:10491519/151‑1 (MQ=255)
                                                                                                                gcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGc                                   >  1:3056351/1‑151 (MQ=255)
                                                                                                                 cccTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCa                                  >  1:4974671/1‑151 (MQ=255)
                                                                                                                 cccTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCa                                  <  1:1968733/151‑1 (MQ=255)
                                                                                                                  ccTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCAt                                 >  2:3685433/1‑151 (MQ=255)
                                                                                                                  ccTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCAt                                 <  1:6868490/151‑1 (MQ=255)
                                                                                                                  ccTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCAt                                 >  1:3805881/1‑151 (MQ=255)
                                                                                                                  ccTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGATGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCAt                                 <  2:2561963/151‑1 (MQ=255)
                                                                                                                    tGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGc                               <  2:5078135/151‑1 (MQ=255)
                                                                                                                      cAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGCGCATACTGGAGCACTGTCCGGCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGccc                             <  2:1809590/151‑1 (MQ=255)
                                                                                                                            gcgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAgg                       >  2:10620597/1‑151 (MQ=255)
                                                                                                                             cgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAggg                      >  2:1535044/1‑151 (MQ=255)
                                                                                                                              ggATGTGGTGATGGTCGACGATCTGGTCCACCCGAGCATACTGGAGCACTGGCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGTCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGa                     <  2:2195904/151‑1 (MQ=255)
                                                                                                                              ggATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGa                     >  1:7905046/1‑151 (MQ=255)
                                                                                                                               gATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGAc                    <  1:2573425/151‑1 (MQ=255)
                                                                                                                               gATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGAc                    <  1:3632691/151‑1 (MQ=255)
                                                                                                                                aTGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACg                   >  2:2066845/1‑151 (MQ=255)
                                                                                                                                aTGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACg                   >  2:473489/1‑151 (MQ=255)
                                                                                                                                aTGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGACAACGCGGCGGCTGCCGCTCCTCACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACg                   >  2:6514820/1‑151 (MQ=255)
                                                                                                                                 tgtgGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGc                  >  1:5234367/1‑151 (MQ=255)
                                                                                                                                  gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCa                 >  1:508891/1‑151 (MQ=255)
                                                                                                                                  gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCa                 <  2:9545234/151‑1 (MQ=255)
                                                                                                                                    ggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGc               >  2:2926931/1‑151 (MQ=255)
                                                                                                                                      tgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgt             >  1:3394539/1‑151 (MQ=255)
                                                                                                                                         tGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggt          >  1:6698063/1‑151 (MQ=255)
                                                                                                                                           gTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTgc        >  2:3124695/1‑151 (MQ=255)
                                                                                                                                            tCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTgcg       >  2:6726598/1‑151 (MQ=255)
                                                                                                                                                cgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTg   >  2:8357705/1‑151 (MQ=255)
                                                                                                                                                cgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTg   >  1:1342695/1‑151 (MQ=255)
                                                                                                                                                cgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTg   <  1:10274137/151‑1 (MQ=255)
                                                                                                                                                 gaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGa  >  2:6345679/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCCCCGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATTTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGGTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTGA  >  NC_008463/3711601‑3711896

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: