Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,949,839 A→T 6.9% intergenic (+61/+52) PP_2583 → / ← oguA acetyltransferase/8‑oxoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,949,8390AT 6.9% 313.5 / 9.4 101intergenic (+61/+52)PP_2583/oguAacetyltransferase/8‑oxoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (46/48);  new base T (3/4);  total (49/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 103 total aligned reads displayed.

TGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTCACCACGAC  >  NC_002947/2949690‑2949986
                                                                                                                                                      |                                                                                                                                                   
tGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACTGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCaaa                                                                                                                                                     <  1:1967805/150‑1 (MQ=255)
   aCCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCGTGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAggg                                                                                                                                                  <  2:729151/150‑1 (MQ=255)
   aCCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAggg                                                                                                                                                  <  1:1819589/150‑1 (MQ=255)
    ccAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGc                                                                                                                                                 >  2:819318/1‑150 (MQ=255)
     cAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCt                                                                                                                                                >  1:1390468/1‑150 (MQ=255)
     cAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCt                                                                                                                                                >  2:1603432/1‑150 (MQ=255)
      aGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTg                                                                                                                                               >  1:380854/1‑150 (MQ=255)
       gAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATTGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTgc                                                                                                                                              >  1:1666365/1‑150 (MQ=255)
            tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCTTTGGGCTgcgctgc                                                                                                                                         <  1:626890/150‑4 (MQ=255)
            tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCAt                                                                                                                                          <  2:763757/149‑1 (MQ=255)
            tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCAt                                                                                                                                          >  1:763757/1‑149 (MQ=255)
            tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATg                                                                                                                                         <  1:582800/150‑1 (MQ=255)
            tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATg                                                                                                                                         <  2:1677283/150‑1 (MQ=255)
              ggCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGcc                                                                                                                                       <  1:1534171/150‑1 (MQ=255)
                     gATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCtt                                                                                                                                    <  1:1218174/146‑1 (MQ=255)
                     gATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCtt                                                                                                                                    >  2:1218174/1‑146 (MQ=255)
                         cGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCAcc                                                                                                                            <  1:1661538/150‑1 (MQ=255)
                          gTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCAccc                                                                                                                           >  2:2002678/1‑150 (MQ=255)
                               aCCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCAGCAGCTTGCCTAGGTGGCAGGGCCGCCCAAAGGGCTGCGCATGCCCTTCAGTCACCCGGCGa                                                                                                                      <  2:398576/150‑1 (MQ=255)
                                ccgtcgccgtcgGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGAt                                                                                                                     <  2:380854/150‑1 (MQ=255)
                                 cgtcgccgtcgGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATg                                                                                                                    <  1:1603432/150‑1 (MQ=255)
                                          cgGACCACGGCGGCATGGCCTGGGATGAGATGGTCATGGACATTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCNGCCCNAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGa                                                                                                           <  1:2518/150‑1 (MQ=255)
                                            gACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGAGGACATCACACCGTAGTGTGCGCTACACGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATCCCGAGCTGATg                                                                                                         >  1:1807149/1‑150 (MQ=255)
                                               cACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCa                                                                                                      <  1:2351998/150‑1 (MQ=255)
                                               cACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCa                                                                                                      <  1:1507098/150‑1 (MQ=255)
                                                  ggcggcATGGCCTGGCATGAGATGGGCATGGAGCTTTGAGTGGTCATGACATCGTAGTGTGCGCTACCCGCAGCTTGTCTAGGTGGGAGCGCAGCCCAAAGGNCTGCGCNTGCCCTTCATTCNCCCGGCGATGCCGAGCTGATGCCAGGc                                                                                                   <  1:1628964/150‑1 (MQ=255)
                                                       cATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCaa                                                                                              <  1:1337269/150‑1 (MQ=255)
                                                        aTGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAg                                                                                             >  2:318561/1‑150 (MQ=255)
                                                             cTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCa                                                                                        >  2:2205476/1‑150 (MQ=255)
                                                              tGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGa                                                                                                           <  2:638875/130‑1 (MQ=255)
                                                              tGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGa                                                                                                           >  1:638875/1‑130 (MQ=255)
                                                                     agaTGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAg                                                                                >  1:723855/1‑150 (MQ=255)
                                                                       aTGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCATGCGCAGCCCTTTGGGCTgcgctgc                                                                                                                                         <  1:826430/92‑4 (MQ=37)
                                                                       aTGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCATGCGCAGCCCTTTGGGCTgcgctgc                                                                                                                                         >  2:826430/1‑89 (MQ=37)
                                                                       aTGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCt                                                                              >  1:1453683/1‑150 (MQ=255)
                                                                       aTGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCt                                                                              <  2:2230413/150‑1 (MQ=255)
                                                                         ggTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGa                                                                            <  1:2002678/150‑1 (MQ=255)
                                                                         ggTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGa                                                                            <  1:318561/150‑1 (MQ=255)
                                                                           tCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATg                                                                          >  2:2275800/1‑150 (MQ=255)
                                                                            cATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATgg                                                                         <  1:819318/150‑1 (MQ=255)
                                                                             aTGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATggg                                                                        >  2:685187/1‑150 (MQ=255)
                                                                              tGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGc                                                                       >  2:412136/1‑150 (MQ=255)
                                                                              gggACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGc                                                                       <  1:1785786/149‑1 (MQ=255)
                                                                                   cTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCATGCGCAGCCCTTTGGGCTgcgctgc                                                                                                                                         <  2:892619/80‑4 (MQ=25)
                                                                                   cTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCATGCGCAGCCCTTTGGGCTgcgctgc                                                                                                                                         >  1:892619/1‑77 (MQ=25)
                                                                                    tttGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTg                                                                 >  2:1970790/1‑150 (MQ=255)
                                                                                    tttGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTg                                                                 >  1:225415/1‑150 (MQ=255)
                                                                                        aGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCGGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtg                                                             <  1:16018/150‑1 (MQ=255)
                                                                                         gTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCa                                                                                                      <  2:190428/108‑1 (MQ=255)
                                                                                         gTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCa                                                                                                      >  1:190428/1‑108 (MQ=255)
                                                                                         gTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtgt                                                            >  1:1247775/1‑150 (MQ=255)
                                                                                         gTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtgt                                                            >  2:2233878/1‑150 (MQ=255)
                                                                                          tGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtgtg                                                           >  2:1518187/1‑150 (MQ=255)
                                                                                             acatgacatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGtt                                                              <  1:507114/144‑1 (MQ=255)
                                                                                             acatgacatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGtt                                                              >  2:507114/1‑144 (MQ=255)
                                                                                               atgacatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAGGGCTGCGCATg                                                                                                                                         >  1:1305955/1‑66 (MQ=255)
                                                                                               atgacatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAGGGCTGCGCATg                                                                                                                                         <  2:1305955/66‑1 (MQ=255)
                                                                                                tgacatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTc                                                     <  1:2233878/150‑1 (MQ=255)
                                                                                                  acatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCa                                                   <  1:290196/150‑1 (MQ=255)
                                                                                                  acatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCa                                                   <  2:1447386/150‑1 (MQ=255)
                                                                                                  acatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCa                                                   <  2:225415/150‑1 (MQ=255)
                                                                                                    atCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTgcgct                                                                                                                                           >  1:1522421/1‑59 (MQ=255)
                                                                                                    atCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGc                                                                       <  2:49887/128‑1 (MQ=255)
                                                                                                    atCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGc                                                                       >  1:49887/1‑128 (MQ=255)
                                                                                                     tCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTa                                                <  1:967867/150‑1 (MQ=255)
                                                                                                      cGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCCCCCGACGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTAc                                               <  2:1666365/150‑1 (MQ=255)
                                                                                                               cgcTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGcc                                                                                         >  2:133824/1‑99 (MQ=255)
                                                                                                               cgcTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGcc                                                                                         <  1:133824/99‑1 (MQ=255)
                                                                                                                gcTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATgg                                                                         <  1:566633/114‑1 (MQ=255)
                                                                                                                gcTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATgg                                                                         >  2:566633/1‑114 (MQ=255)
                                                                                                                gcTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAgg                                     >  2:742440/1‑150 (MQ=255)
                                                                                                                  tACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTc                                   <  2:181372/150‑1 (MQ=255)
                                                                                                                   aCCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCaa                                                                                              <  2:1691616/90‑1 (MQ=255)
                                                                                                                   aCCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCaa                                                                                              >  1:1691616/1‑90 (MQ=255)
                                                                                                                     ccGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGc                                                                                          <  1:1463842/92‑1 (MQ=255)
                                                                                                                     ccGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGc                                                                                          >  2:1463842/1‑92 (MQ=255)
                                                                                                                     ccGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGAc                                >  2:1387151/1‑150 (MQ=255)
                                                                                                                     ccGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGAc                                >  1:1178612/1‑150 (MQ=255)
                                                                                                                     ccGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCCCGGCAGGTCGAc                                >  2:623087/1‑150 (MQ=255)
                                                                                                                      cGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACg                               >  2:1457579/1‑150 (MQ=255)
                                                                                                                      cGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACg                               >  1:1313998/1‑150 (MQ=255)
                                                                                                                       gCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACgg                              >  2:1219609/1‑150 (MQ=255)
                                                                                                                         aGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTg                            <  1:2275800/150‑1 (MQ=255)
                                                                                                                           cTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGc                          <  1:449369/150‑1 (MQ=255)
                                                                                                                            ttGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGc                                        >  2:764746/1‑135 (MQ=255)
                                                                                                                            ttGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGc                                        <  1:764746/135‑1 (MQ=255)
                                                                                                                                caaGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCt                     >  2:1801891/3‑150 (MQ=255)
                                                                                                                                cTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCt                     <  1:1970790/150‑1 (MQ=255)
                                                                                                                                 tAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGc                                                                                                   <  2:1584221/71‑1 (MQ=255)
                                                                                                                                 tAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGc                                                                                                   >  1:1584221/1‑71 (MQ=255)
                                                                                                                                 tAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTg                    <  1:685187/150‑1 (MQ=255)
                                                                                                                                      ggCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGt                                    >  2:2175242/1‑129 (MQ=255)
                                                                                                                                      ggCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGt                                    <  1:2175242/129‑1 (MQ=255)
                                                                                                                                       gcagcgcagcCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTc              >  1:2153190/1‑150 (MQ=255)
                                                                                                                                        catgcgcagcCCTTTGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCa                                                                                 >  1:1459830/4‑83 (MQ=38)
                                                                                                                                        catgcgcagcCCTTTGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCa                                                                                 <  2:1459830/80‑1 (MQ=38)
                                                                                                                                         agcgcagcCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCtt            >  1:905228/1‑150 (MQ=255)
                                                                                                                                                gcCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTcacca      <  2:1390468/150‑1 (MQ=255)
                                                                                                                                                  ccAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTCACCacg    >  1:690744/1‑150 (MQ=255)
                                                                                                                                                    aaaGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTCACCacgac  >  2:1453652/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTCACCACGAC  >  NC_002947/2949690‑2949986

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