Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,981,210 G→A 5.1% A179V (GCC→GTC)  cobW ← adenosylcobalamin biosynthesis protein CobW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,981,2100GA 5.1% 261.9 / 2.6 89A179V (GCC→GTC) cobWadenosylcobalamin biosynthesis protein CobW
Reads supporting (aligned to +/- strand):  ref base G (42/42);  new base A (3/2);  total (45/44)
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.70e-01

GCAGCACCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGACG  >  NC_002947/3981064‑3981359
                                                                                                                                                  |                                                                                                                                                     
gcagcaCCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCa                                                                                                                                                    <  2:2050400/150‑1 (MQ=255)
gcagcaCCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCa                                                                                                                                                    >  2:2208272/1‑150 (MQ=255)
   gcaCCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCCGCt                                                                                                                                                 >  2:845260/1‑150 (MQ=255)
    caCCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTg                                                                                                                                                <  1:2080176/150‑1 (MQ=255)
    caCCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTg                                                                                                                                                <  2:856539/150‑1 (MQ=255)
          cAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGAGGGACAGCGGGTGAGc                                                                                                                                          >  2:1869946/1‑150 (MQ=255)
             gggCAGCTTGCCGCTGCTGGCCTCGCTCACTTTCCCAGACGGGGGCAGGTCTTCGGGGGGGGCTGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGCCGGCCTTGTTCAACTCCACCCCGTCGGCGCTGGCCAGCTGGTCATCGaa                                                                                                                                       <  2:1728317/149‑1 (MQ=255)
                cAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAg                                                                                                                                    >  1:730533/1‑150 (MQ=255)
                    ttGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGc                                                                                                                                          <  2:1336020/140‑1 (MQ=255)
                    ttGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGc                                                                                                                                          >  1:1336020/1‑140 (MQ=255)
                      gCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTg                                                                                                                              >  2:1508491/1‑150 (MQ=255)
                           tgctgGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAggg                                                                                                                         >  1:2273233/1‑150 (MQ=255)
                              tgGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCg                                                                                                                      >  2:199830/1‑150 (MQ=255)
                              tgGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCg                                                                                                                      >  2:1700397/1‑150 (MQ=255)
                               ggCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACAACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGa                                                                                                                     >  2:58033/1‑150 (MQ=255)
                               ggCCTCGAGCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTTCACCTGTACATCGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGa                                                                                                                     <  2:695242/150‑1 (MQ=255)
                                  cTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGACCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTc                                                                                                                  >  2:419301/1‑150 (MQ=255)
                                   tCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCa                                                                                                                 <  2:883890/150‑1 (MQ=255)
                                   tCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCACGTCACCCACCCTGGCCACGACCAGGTCGGCGCTGGCAAGATGGCCAGCGAACAGCTCGTGCAGGGGCGATTCa                                                                                                                 >  1:1251448/1‑150 (MQ=255)
                                       tCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGt                                                                                                             >  2:1227481/1‑150 (MQ=255)
                                         aCTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGACCAGCTGGCCAGCgccga                                                                                                                                     >  2:226383/1‑120 (MQ=255)
                                         aCTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGACCAGCTGGCCAGCgccga                                                                                                                                     <  1:226383/124‑5 (MQ=255)
                                                ccgccgGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAg                                                                                                    <  1:1249386/150‑1 (MQ=255)
                                                 cgccgGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAgg                                                                                                   <  1:1554597/150‑1 (MQ=255)
                                                  gccgGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAggg                                                                                                  >  1:1155517/1‑150 (MQ=255)
                                                      gcggcAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGa                                                                                              <  1:517957/150‑1 (MQ=255)
                                                      gcggcAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTAAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGa                                                                                              >  1:915067/1‑150 (MQ=255)
                                                       cggcAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAg                                                                                             <  1:199830/150‑1 (MQ=255)
                                                                ttcGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGCTCAGGTCGTCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTg                                                                                    >  1:1254887/1‑150 (MQ=255)
                                                                ttcGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTg                                                                                    >  2:1253545/1‑150 (MQ=255)
                                                                    gCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTa                                                                                                     <  1:818284/129‑1 (MQ=255)
                                                                    gCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTa                                                                                                     >  2:818284/1‑129 (MQ=255)
                                                                    gCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTgg                                                                                   <  1:2081960/147‑1 (MQ=255)
                                                                    gCCACCTCGGCACGCACCTTGGCCAGGCATTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTgg                                                                                   >  2:2081960/1‑147 (MQ=255)
                                                                     ccACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCAt                                                                               >  2:2181713/1‑150 (MQ=255)
                                                                        ccTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGa                                                                            >  1:549302/1‑150 (MQ=255)
                                                                         cTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGAt                                                                                                                    >  1:1736437/1‑109 (MQ=255)
                                                                         cTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGAt                                                                                                                    >  1:1737758/1‑109 (MQ=255)
                                                                         cTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGAt                                                                                                                    <  2:1737758/109‑1 (MQ=255)
                                                                         cTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCtt                                                                                          <  2:893407/135‑1 (MQ=255)
                                                                         cTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCtt                                                                                          >  1:893407/1‑135 (MQ=255)
                                                                         cTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGAc                                                                           <  1:2219286/150‑1 (MQ=255)
                                                                         cTCCGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGCCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGAt                                                                                                                    <  2:1736437/109‑1 (MQ=255)
                                                                          tCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAgg                                                                                                   >  1:2107805/1‑125 (MQ=255)
                                                                          tCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAgg                                                                                                   <  2:2107805/125‑1 (MQ=255)
                                                                          tCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGAcc                                                                          <  1:871388/150‑1 (MQ=255)
                                                                          tCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGAcc                                                                          <  1:2372495/150‑1 (MQ=255)
                                                                          tCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGAcc                                                                          <  1:1508491/150‑1 (MQ=255)
                                                                                  cacCTTGGGCCGCCCTTAAGCGTCGACTCTGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCATCACGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTGCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCgg                                                                  <  1:419301/150‑1 (MQ=255)
                                                                                  cacCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCgg                                                                  <  1:464557/150‑1 (MQ=255)
                                                                                   acCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCggg                                                                 >  2:668154/1‑150 (MQ=255)
                                                                                   acCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCggg                                                                 >  1:270210/1‑150 (MQ=255)
                                                                                     cTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTa                                                               >  2:2030476/1‑150 (MQ=255)
                                                                                     cTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTa                                                               <  1:2030476/150‑1 (MQ=255)
                                                                                      ttGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAg                                                              >  1:729249/1‑150 (MQ=255)
                                                                                         gCCAGGCCTTCGGCGTCGCTCAGGTCGACCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAggcg                                                           <  2:315478/150‑1 (MQ=255)
                                                                                             ggCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCaa                                                       >  2:500391/1‑150 (MQ=255)
                                                                                             ggCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTCGGGATCGACCTGGTCCGGGTAGGCGGCaa                                                       >  2:1012846/1‑150 (MQ=255)
                                                                                              gCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCaaa                                                      <  1:668154/150‑1 (MQ=255)
                                                                                               ccTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAg                                                     <  1:2208272/150‑1 (MQ=255)
                                                                                                cTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCGGGGTAGGCGGCAAAgg                                                    >  2:1434146/1‑150 (MQ=255)
                                                                                                 ttCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGNTCGAGCNTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGt                                                   <  1:1440295/150‑1 (MQ=255)
                                                                                                 ttCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGt                                                   <  2:549302/150‑1 (MQ=255)
                                                                                                   cGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGc                                                 <  2:270210/150‑1 (MQ=255)
                                                                                                         cGATCAGGTCGGCCTTGTACAGCACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGCTCCTGGTTAGGGTCGAGCTTGCGCTGGGCATCCACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCgg                                           <  1:520822/150‑1 (MQ=255)
                                                                                                            tCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAggc                                                            >  2:475492/1‑130 (MQ=255)
                                                                                                            tCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAggc                                                            <  1:475492/130‑1 (MQ=255)
                                                                                                            tCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCa                                        <  2:1820090/150‑1 (MQ=255)
                                                                                                                gTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGc                                    >  2:209386/1‑150 (MQ=255)
                                                                                                                   ggCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCgg                                 <  1:2181713/150‑1 (MQ=255)
                                                                                                                   ggCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCgg                                 <  1:209386/150‑1 (MQ=255)
                                                                                                                   ggCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGATAGGGTCGAGCTTGCGCGCCAGATCTACGTTGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCgg                                 <  1:1253545/150‑1 (MQ=255)
                                                                                                                      cTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCt                              >  1:1442817/1‑150 (MQ=255)
                                                                                                                          ttCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTgg                                                                                   >  1:1196597/1‑93 (MQ=255)
                                                                                                                          ttCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTgg                                                                                   <  2:1196597/93‑1 (MQ=255)
                                                                                                                          ttCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCg                          >  1:1558170/1‑150 (MQ=255)
                                                                                                                           tCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGAc                                                                           >  2:899892/1‑100 (MQ=255)
                                                                                                                           tCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGAc                                                                           <  1:899892/100‑1 (MQ=255)
                                                                                                                            cAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGAc                        <  1:1227481/150‑1 (MQ=255)
                                                                                                                             aaCACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGAcc                       >  2:330963/1‑150 (MQ=255)
                                                                                                                              acacCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCa                      <  1:330963/150‑1 (MQ=255)
                                                                                                                                accaccGTGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGGCCACg                    >  2:1300220/1‑150 (MQ=255)
                                                                                                                                accaccAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGAATCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACg                    >  2:823508/1‑150 (MQ=255)
                                                                                                                                 ccaccaGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACgg                   <  1:845260/150‑1 (MQ=255)
                                                                                                                                  caccaGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATTGACACAGTCCCGGAAGTCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGt                  >  2:1873384/1‑150 (MQ=255)
                                                                                                                                         tCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGAcc           >  1:1837774/1‑150 (MQ=255)
                                                                                                                                         tCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGAcc           <  2:915067/150‑1 (MQ=255)
                                                                                                                                           ttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                           >  2:1926058/3‑132 (MQ=255)
                                                                                                                                           ttcgcTGACCAGCTGGCCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTc                           <  1:1926058/130‑1 (MQ=255)
                                                                                                                                                  gCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGAcc           <  1:822262/141‑1 (MQ=255)
                                                                                                                                                  gCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGAcc           >  2:822262/1‑141 (MQ=255)
                                                                                                                                                  gCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGACg  <  2:1155517/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCAGCACCTCCAGTGGCAGCTTGCCGCTGCTGGCCTCGATCACTTTCACCGCCGGCGGCAGTTCTTCGGCCACCTCGGCACGCACCTTGGCCAGGCCTTCGGCGTCGATCAGGTCGGCCTTGTTCAACACCACCAGGTCGGCGCTGGCCAGCTGGTCAGCGAACAGCTCGTGCAGGGGCGATTCATGGTCCAGGTTAGGGTCGAGCTTGCGCTGGGCATCGACCTGGTCCGGGTAGGCGGCAAAGGTGCCAGCGGCCACGGCCGGGCTGTCGACCACGGTAATGACCGCATCGACG  >  NC_002947/3981064‑3981359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: