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.6% 1799.6 / NA 471L37L (TTG→CTG) cobAuroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (254/215);  minor base A (2/0);  total (256/215)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 471 total aligned reads displayed.

CCCCCGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATTTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGGTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTG  >  NC_008463/3711602‑3711895
                                                                                                                                                   |                                                                                                                                                  
 ccccGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTgg                                                                                                                                                 <  1:3571826/150‑1 (MQ=255)
    cGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCa                                                                                                                                              >  2:788045/1‑150 (MQ=255)
        ccATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAAccc                                                                                                                                          <  2:2480761/150‑1 (MQ=255)
         cATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCg                                                                                                                                         >  1:3724382/1‑150 (MQ=255)
           tGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAg                                                                                                                                       >  2:3935811/1‑150 (MQ=255)
           tGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAg                                                                                                                                       >  1:8206184/1‑150 (MQ=255)
             gcggcgAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCa                                                                                                                                     <  2:5551547/150‑1 (MQ=255)
             gcggcgAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCGGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCa                                                                                                                                     <  1:3126534/150‑1 (MQ=255)
              cggcgAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCAt                                                                                                                                    <  2:1101884/150‑1 (MQ=255)
               ggcgAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATa                                                                                                                                   <  2:8396763/150‑1 (MQ=255)
               ggcgAATGCCGGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATa                                                                                                                                   >  2:7507100/1‑150 (MQ=255)
                   aaTGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTgg                                                                                                                               <  2:7350415/150‑1 (MQ=255)
                      gCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGc                                                                                                                            <  1:7552086/150‑1 (MQ=255)
                        cTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCAc                                                                                                                          >  2:4108081/1‑150 (MQ=255)
                         tGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACt                                                                                                                         <  1:4746089/150‑1 (MQ=255)
                         tGAGGAGAACGACGAGGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACt                                                                                                                         <  2:643636/150‑1 (MQ=255)
                          gaggagAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTg                                                                                                                        >  2:2097147/1‑150 (MQ=255)
                           aggagAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTgt                                                                                                                       >  2:963850/1‑150 (MQ=255)
                            ggagAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTgtc                                                                                                                      <  1:6007587/150‑1 (MQ=255)
                              agaACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTgtccg                                                                                                                    <  2:4061038/150‑1 (MQ=255)
                                  cgacgaTGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGtccgtccg                                                                                                                >  1:3864023/1‑150 (MQ=255)
                                    acgaTGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGcc                                                                                                              >  1:3224084/1‑150 (MQ=255)
                                    acgaTGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGcc                                                                                                              >  1:9849461/1‑150 (MQ=255)
                                      gaTGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCg                                                                                                            >  1:9228222/1‑150 (MQ=255)
                                       aTGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCgg                                                                                                           >  1:8922421/1‑150 (MQ=255)
                                         gAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCggct                                                                                                         <  2:10275396/150‑3 (MQ=255)
                                         gAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCggct                                                                                                         <  1:5335053/150‑3 (MQ=255)
                                           gCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGgctgg                                                                                                       >  1:2659688/1‑150 (MQ=255)
                                             ggTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTc                                                                                                     >  1:4038880/1‑150 (MQ=255)
                                                aaGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCgg                                                                                                  <  1:942579/150‑1 (MQ=255)
                                                   gTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTc                                                                                               <  2:1424825/150‑1 (MQ=255)
                                                   gTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTc                                                                                               <  1:1328186/150‑1 (MQ=255)
                                                     ctggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCgg                                                                                             <  1:7778246/150‑1 (MQ=255)
                                                      tggctggTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGc                                                                                            <  2:6563387/150‑1 (MQ=255)
                                                          tggtggGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAc                                                                                        >  1:2934058/1‑150 (MQ=255)
                                                           ggtgggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCaaaag                                                                                       <  2:8785834/150‑3 (MQ=255)
                                                            gtgggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAAcgc                                                                                      <  1:10068349/150‑1 (MQ=255)
                                                              gggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcgg                                                                                    <  1:802098/150‑1 (MQ=255)
                                                              gggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcgg                                                                                    <  1:8719996/150‑1 (MQ=255)
                                                               ggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGCAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTTGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcggc                                                                                   <  2:940194/150‑1 (MQ=255)
                                                               ggtgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACgcggc                                                                                   <  2:9876497/150‑1 (MQ=255)
                                                                 tgCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGcggcgg                                                                                 >  2:5658114/1‑150 (MQ=255)
                                                                     ggTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGTTGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGc                                                                             <  2:10133064/150‑1 (MQ=255)
                                                                         ccGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCt                                                                         <  1:7066053/150‑1 (MQ=255)
                                                                               aTCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCacac                                                                   >  1:4241483/1‑150 (MQ=255)
                                                                                 ccGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCg                                                                 >  1:4418721/1‑150 (MQ=255)
                                                                                   ggAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCa                                                               <  1:6788894/150‑1 (MQ=255)
                                                                                   ggAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCa                                                               <  2:10792401/150‑1 (MQ=255)
                                                                                    gAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAATCGCGGCGGCTGCCGCTCCACACCGCAg                                                              >  2:1564702/1‑150 (MQ=255)
                                                                                    gAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAg                                                              <  2:5214821/150‑1 (MQ=255)
                                                                                    gAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAg                                                              >  2:3690393/1‑150 (MQ=255)
                                                                                      aCTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGa                                                            <  2:9054547/150‑1 (MQ=255)
                                                                                       cTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGAc                                                           >  1:3946718/1‑150 (MQ=255)
                                                                                       cTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGAc                                                           <  2:8672102/150‑1 (MQ=255)
                                                                                       cTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGAc                                                           >  2:2341311/1‑150 (MQ=255)
                                                                                        tGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACt                                                          <  2:2614488/150‑1 (MQ=255)
                                                                                         gTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACtt                                                         <  1:6554545/150‑1 (MQ=255)
                                                                                         gTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATATGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACtt                                                         >  2:10313817/1‑150 (MQ=255)
                                                                                               ccTGACGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCCACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCa                                                   <  1:146288/150‑1 (MQ=255)
                                                                                                 tGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAgc                                                 <  2:5877481/150‑1 (MQ=255)
                                                                                                 tGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAgc                                                 <  1:5841323/150‑1 (MQ=255)
                                                                                                    aGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGcc                                              <  1:5190466/150‑1 (MQ=255)
                                                                                                    aGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGcc                                              >  2:5369016/1‑150 (MQ=255)
                                                                                                     ggCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCt                                             <  1:10971845/150‑1 (MQ=255)
                                                                                                       cAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGa                                           >  2:7061297/1‑150 (MQ=255)
                                                                                                           tCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAATCGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgct                                       >  2:3416223/1‑150 (MQ=255)
                                                                                                            cgcgCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgctg                                      <  2:9235311/150‑1 (MQ=255)
                                                                                                            cgcgCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgctg                                      >  2:7727369/1‑150 (MQ=255)
                                                                                                             gcgcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgctgc                                     <  1:1742033/150‑1 (MQ=255)
                                                                                                             gcgcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGAtgctgc                                     >  1:6297804/1‑150 (MQ=255)
                                                                                                              cgcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCg                                    >  1:10198621/1‑150 (MQ=255)
                                                                                                               gcCCTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCgg                                   <  2:8943192/150‑1 (MQ=255)
                                                                                                                cccTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGc                                  >  1:9552999/1‑150 (MQ=255)
                                                                                                                 ccTGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCa                                 <  1:2215490/150‑1 (MQ=255)
                                                                                                                  cTGCACTACGAGGATGTTGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCAGACCGACCGGCTGGTCCGTGTCGGCAAACTCGTCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCAt                                >  2:6951898/1‑150 (MQ=255)
                                                                                                                   tGCAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATg                               >  1:1560063/1‑150 (MQ=255)
                                                                                                                     cAGGACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGcc                             <  1:5575375/150‑1 (MQ=255)
                                                                                                                       ggACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCg                           >  2:5127401/1‑150 (MQ=255)
                                                                                                                        gACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGc                          >  2:4343379/1‑150 (MQ=255)
                                                                                                                        gACGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGc                          <  1:4930900/150‑1 (MQ=255)
                                                                                                                         aCGCGGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGcc                         >  1:10450114/1‑150 (MQ=255)
                                                                                                                           gcgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAg                       >  1:1856264/1‑150 (MQ=255)
                                                                                                                            cgGATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAgg                      >  2:5956573/1‑150 (MQ=255)
                                                                                                                              gATGTGGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGa                    >  2:3196378/1‑150 (MQ=255)
                                                                                                                                tgtgGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACg                  >  2:4707335/1‑150 (MQ=255)
                                                                                                                                tgtgGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACg                  >  2:9783483/1‑150 (MQ=255)
                                                                                                                                tgtgGTGATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACg                  >  1:7221953/1‑150 (MQ=255)
                                                                                                                                 gtggtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGc                 >  2:1142134/1‑150 (MQ=255)
                                                                                                                                    gtgATGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGc              <  1:8398383/150‑1 (MQ=255)
                                                                                                                                    gtgATGGTCGACGATCGGGTCCACCCGAGCATACGGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCCGCCGCTCCACACCGCCGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGc              <  2:1854684/150‑1 (MQ=255)
                                                                                                                                      gttGGTCGTCGATCTGTTCAACCCGATCCTACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTTCCGCTCCACTCGGCAGGACTTCATCCATCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgt            <  2:105823/148‑1 (MQ=255)
                                                                                                                                        tGGTCGACGATCTGGTCAACCCGAGCATACTGGAGCACGGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtgg          <  2:10406488/150‑1 (MQ=255)
                                                                                                                                         ggTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggt         <  1:7978377/150‑1 (MQ=255)
                                                                                                                                         ggTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggt         >  2:6674837/1‑150 (MQ=255)
                                                                                                                                          gTCGACGATCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCgtggtg        <  2:3576817/150‑1 (MQ=255)
                                                                                                                                             gacgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTgcgc     <  1:10887746/150‑1 (MQ=255)
                                                                                                                                              acgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCTGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGcc    >  1:7429666/1‑150 (MQ=255)
                                                                                                                                              acgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGcc    >  2:2988128/1‑150 (MQ=255)
                                                                                                                                              acgaTCTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGcc    >  2:9594886/1‑150 (MQ=255)
                                                                                                                                               cgaTCTGGTCAACCCGAGCATACTGGCGCACTGTCCGTCCGCCCGGCTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCt   <  2:8224148/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCCCCGGCCCATGGCGGCGAATGCCTGAGGAGAACGACGATGAGCGGTAAGGTCTGGCTGGTGGGTGCGGGTCCCGGCGATCCGGAACTGTTGACCCTGAAGGCAGTTCGCGCCCTGCAGGACGCGGATGTGGTGATGGTCGACGATTTGGTCAACCCGAGCATACTGGAGCACTGTCCGTCCGCCCGGGTGGTCCGGGTCGGCAAACGCGGCGGCTGCCGCTCCACACCGCAGGACTTCATCCAGCGCCTGATGCTGCGGCATGCCCGCCAGGGACGCAGCGTGGTGCGCCTG  >  NC_008463/3711602‑3711895

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: