Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,502,089 C→G 6.5% H176Q (CAC→CAG)  tehB → tellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,502,0890CG 6.5% 163.4 / 2.6 58H176Q (CAC→CAG) tehBtellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base G (2/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

GTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGA  >  NC_000913/1501941‑1502238
                                                                                                                                                    |                                                                                                                                                     
gtggtTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTgcacc                                                                                                                                                      >  2:434105/1‑150 (MQ=255)
  ggtTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTgcaccgc                                                                                                                                                    >  1:1152531/1‑150 (MQ=255)
     tACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACTGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATTTTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTgcaccgcacc                                                                                                                                                 >  1:101444/1‑150 (MQ=255)
       cAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGcaccgcaccgc                                                                                                                                               >  1:1044313/1‑149 (MQ=255)
       cAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGa                                                                                                                                               >  1:628979/1‑150 (MQ=255)
            tGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGTGCGGTGCAGCTCGCCGACgtct                                                                                                                                          <  2:1008990/150‑4 (MQ=255)
              tttGTGGCTGCGATGGATACCGCGGATTATCCATGTACCGTCGGCATCCCGTTTGCCTTCAAAGATGGAGAATTACGTCGATATTACGAAGGCTGCGAGTGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAAc                                                                                                                                        <  2:1163620/149‑1 (MQ=255)
                tgtgGCGGCGATTGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGATGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCTACGCCAACgg                                                                                                                                      >  1:802998/1‑150 (MQ=255)
                  tgGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTa                                                                                                                                    >  1:753985/1‑150 (MQ=255)
                   ggcggcGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGGGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTaa                                                                                                                                   <  1:4598/150‑1 (MQ=255)
                    gcggcgATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAAt                                                                                                                                  >  2:1125822/1‑150 (MQ=255)
                    gcggcgATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTAAAAGAGGGAGAATTACGTCGATATTACGAAGGAAGGGAGAGGGGGAAATACAATGACGACGTCCGCGCGCTGCACCGCACCGACGCCAACGGTAAt                                                                                                                                  >  2:621281/1‑150 (MQ=255)
                     cggcgATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATc                                                                                                                                 >  1:73525/1‑150 (MQ=255)
                        cgATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTAGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAATAAGTCGGCGAGCAGCAACGCACCGACGCCAACGGTAATCGTa                                                                                                                              >  2:367547/1‑150 (MQ=255)
                              aTACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAAc                                                                                                                        <  2:1316119/150‑1 (MQ=255)
                                   gCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGt                                                                                                                   >  2:393147/1‑150 (MQ=255)
                                       aTTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCg                                                                                                               >  2:855423/1‑150 (MQ=255)
                                       aTTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCg                                                                                                               <  1:279769/150‑1 (MQ=255)
                                         tATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGcc                                                                                                             <  2:160359/150‑1 (MQ=255)
                                            ccATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGCCGCCCACGGTAATCGTATTAAACTGCGTTTCGCCACg                                                                                                          >  2:96907/1‑150 (MQ=255)
                                                gTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGc                                                                                                      <  2:272761/150‑1 (MQ=255)
                                                 tACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGTGATATACGTCGATATTACGAAGGCTGGGAGAGGGGGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCt                                                                                                     <  1:678297/150‑1 (MQ=255)
                                                 tACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCt                                                                                                     <  2:73525/150‑1 (MQ=255)
                                                   ccGTCGGCTTCCCGTTTCCCTTCAAAGAGGGAGAATTCCGTCGATATTACGCAGGCTCGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTgg                                                                                                   <  2:35190/150‑1 (MQ=255)
                                                       cGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGACGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACg                                                                                               <  1:551292/150‑1 (MQ=255)
                                                             cccGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTaaaaa                                                                                         >  1:879854/1‑150 (MQ=255)
                                                             cccGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCGGGCACGTaaaaa                                                                                         >  2:1259286/1‑150 (MQ=255)
                                                               cGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAAt                                                                                       >  1:1132297/1‑150 (MQ=255)
                                                                    gCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGTGCGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGcccc                                                                                  >  2:146712/1‑150 (MQ=255)
                                                                    gCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATTCTGGCACGTAAAAAATGAccc                                                                                  <  2:982423/150‑1 (MQ=255)
                                                                    gCCTTCAAAGAGGGAGAATTACCTCGATATTACGAAGGCTGGGAGAGTGTGAAATACAATGAAGACGGCGGCGAGCTGCAGCGCACCGACGCCAACGGAAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGAccc                                                                                  <  1:154712/150‑1 (MQ=255)
                                                                        tCAAAGAGGGAGAATTACGTCGATATTACGAAGGATGGGAGAGGGTGAAATACAATGAAGACGTCGGAGAGCTGAACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGACACGATGCTGGCACGTAAAAAATGACCCTGTa                                                                              >  1:219816/1‑150 (MQ=255)
                                                                         cAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTaa                                                                             >  2:354643/1‑150 (MQ=255)
                                                                           aaGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGCCCCGGTAAGc                                                                           >  2:1074480/1‑150 (MQ=255)
                                                                           aaGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGc                                                                           <  2:598858/150‑1 (MQ=255)
                                                                                               tatTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTcc                                                       <  2:73735/150‑1 (MQ=255)
                                                                                                 ttACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTcccc                                                     <  2:525503/150‑1 (MQ=255)
                                                                                                   aCGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCAc                                                   >  1:266302/1‑150 (MQ=255)
                                                                                                    cGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACg                                                  >  1:220029/1‑150 (MQ=255)
                                                                                                      aaGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCt                                                >  1:686462/1‑150 (MQ=255)
                                                                                                       aGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTg                                               <  1:1023949/150‑1 (MQ=255)
                                                                                                        ggatGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGAATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCTTGAAAATTCCCCACGCTga                                              >  2:165826/4‑150 (MQ=255)
                                                                                                          cTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTgaga                                            >  2:744039/1‑150 (MQ=255)
                                                                                                              gagaGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACTGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAgatgat                                        <  1:166143/150‑1 (MQ=255)
                                                                                                               agagGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGAtt                                       >  2:791171/1‑150 (MQ=255)
                                                                                                               agagGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGAtt                                       <  2:613683/150‑1 (MQ=255)
                                                                                                                   ggTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACt                                   <  1:393147/150‑1 (MQ=255)
                                                                                                                    gTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTg                                  >  1:710550/1‑150 (MQ=255)
                                                                                                                     tGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGt                                 <  1:855423/150‑1 (MQ=255)
                                                                                                                        aaTACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGttct                              <  1:1199167/150‑1 (MQ=255)
                                                                                                                        aaTACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGttct                              <  2:628979/150‑1 (MQ=255)
                                                                                                                          tACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGttcttc                            >  1:610941/1‑150 (MQ=255)
                                                                                                                          tACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGttctta                            >  1:1098452/1‑149 (MQ=255)
                                                                                                                               tGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTACACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTc                       <  2:364618/150‑1 (MQ=255)
                                                                                                                               cgttggcGTCGGTGCGGTGCAGCTCGCCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTc                       >  2:1291647/7‑150 (MQ=255)
                                                                                                                                taaGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCg                      <  2:1152531/149‑1 (MQ=255)
                                                                                                                                  aGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGt                    <  2:753985/150‑1 (MQ=255)
                                                                                                                                   ggcGTCGGTGCGGTGCAGCTCGCCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTa                   >  2:1056926/3‑150 (MQ=255)
                                                                                                                                     cGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAg                 <  1:562870/150‑1 (MQ=255)
                                                                                                                                      gTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGc                >  2:905646/1‑150 (MQ=255)
                                                                                                                                       tCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCa               <  1:593113/150‑1 (MQ=255)
                                                                                                                                        cGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGGat              >  1:218007/1‑150 (MQ=255)
                                                                                                                                                ggcaccgcaccGACGCCAACGGTAATCGTATTAAACTCCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTa      <  2:120946/149‑1 (MQ=255)
                                                                                                                                                   accgcacagACGCCAACGTTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTACCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCg   <  1:96907/150‑1 (MQ=255)
                                                                                                                                                    ccgcaccgACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGa  <  1:750146/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTTGCCTTCAAAGAGGGAGAATTACGTCGATATTACGAAGGCTGGGAGAGGGTGAAATACAATGAAGACGTCGGCGAGCTGCACCGCACCGACGCCAACGGTAATCGTATTAAACTGCGTTTCGCCACGATGCTGGCACGTAAAAAATGACCCGGTAAGCACAAAACGCGTGAAAATTCCCCACGCTGAGATGATTTACTGTTCTTCTTTTCGGTAAGCATATTTTTTATCGA  >  NC_000913/1501941‑1502238

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑