Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,949,837 A→T 7.1% intergenic (+59/+54) 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,8370AT 7.1% 311.3 / 10.1 100intergenic (+59/+54)PP_2583/oguAacetyltransferase/8‑oxoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (45/46);  new base T (3/4);  total (49/51)
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.86e-01

Only 100 of 108 total aligned reads displayed.

CCTGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTCACCACGAC  >  NC_002947/2949688‑2949986
                                                                                                                                                      |                                                                                                                                                     
ccTGTACGAGGAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCa                                                                                                                                                       >  2:2401428/1‑150 (MQ=255)
ccTGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCa                                                                                                                                                       >  1:181372/1‑150 (MQ=255)
ccTGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCa                                                                                                                                                       >  2:449369/1‑150 (MQ=255)
 cTGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCaa                                                                                                                                                      >  1:2230413/1‑150 (MQ=255)
 cTGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTCCCTAGGTGACAGCGCAGCCCaa                                                                                                                                                      >  2:20814/1‑150 (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)
         gAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATTGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTgc                                                                                                                                              >  1:1666365/1‑150 (MQ=255)
              tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCTTTGGGCTgcgctgc                                                                                                                                         <  1:626890/150‑4 (MQ=255)
              tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCAt                                                                                                                                          >  1:763757/1‑149 (MQ=255)
              tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCAt                                                                                                                                          <  2:763757/149‑1 (MQ=255)
              tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATg                                                                                                                                         <  2:1677283/150‑1 (MQ=255)
              tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATg                                                                                                                                         <  1:582800/150‑1 (MQ=255)
                       gATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCtt                                                                                                                                    >  2:1218174/1‑146 (MQ=255)
                       gATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCtt                                                                                                                                    <  1:1218174/146‑1 (MQ=255)
                           cGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCAcc                                                                                                                            <  1:1661538/150‑1 (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:1507098/150‑1 (MQ=255)
                                                 cACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCa                                                                                                      <  1:2351998/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:318561/150‑1 (MQ=255)
                                                                           ggTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGa                                                                            <  1:2002678/150‑1 (MQ=255)
                                                                             tCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATg                                                                          >  2:2275800/1‑150 (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                                                            >  2:2233878/1‑150 (MQ=255)
                                                                                           gTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtgt                                                            >  1:1247775/1‑150 (MQ=255)
                                                                                            tGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtgtg                                                           >  2:1518187/1‑150 (MQ=255)
                                                                                               acatgacatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGtt                                                              >  2:507114/1‑144 (MQ=255)
                                                                                               acatgacatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGtt                                                              <  1:507114/144‑1 (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:225415/150‑1 (MQ=255)
                                                                                                    acatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCa                                                   <  2:1447386/150‑1 (MQ=255)
                                                                                                      atCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTgcgct                                                                                                                                           >  1:1522421/1‑59 (MQ=255)
                                                                                                      atCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTgcgct                                                                                                                                           <  2:1522421/60‑2 (MQ=255)
                                                                                                      atCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGc                                                                       >  1:49887/1‑128 (MQ=255)
                                                                                                       tCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTa                                                <  1:967867/150‑1 (MQ=255)
                                                                                                        cGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCCCCCGACGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTAc                                               <  2:1666365/150‑1 (MQ=255)
                                                                                                              gtgCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCgg                                         <  1:2389936/150‑1 (MQ=255)
                                                                                                                 cgcTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGcc                                                                                         >  2:133824/1‑99 (MQ=255)
                                                                                                                 cgcTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGcc                                                                                         <  1:133824/99‑1 (MQ=255)
                                                                                                                  gcTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATgg                                                                         >  2:566633/1‑114 (MQ=255)
                                                                                                                  gcTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATgg                                                                         <  1:566633/114‑1 (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                                                                                          >  2:1463842/1‑92 (MQ=255)
                                                                                                                       ccGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGc                                                                                          <  1:1463842/92‑1 (MQ=255)
                                                                                                                       ccGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGAc                                >  1:1178612/1‑150 (MQ=255)
                                                                                                                       ccGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGAc                                >  2:1387151/1‑150 (MQ=255)
                                                                                                                       ccGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCCCGGCAGGTCGAc                                >  2:623087/1‑150 (MQ=255)
                                                                                                                        cGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACg                               >  1:1313998/1‑150 (MQ=255)
                                                                                                                        cGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACg                               >  2:1457579/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                                                                                                   >  1:1584221/1‑71 (MQ=255)
                                                                                                                                   tAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGc                                                                                                   <  2:1584221/71‑1 (MQ=255)
                                                                                                                                   tAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTg                    <  1:685187/150‑1 (MQ=255)
                                                                                                                                        ggCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGt                                    <  1:2175242/129‑1 (MQ=255)
                                                                                                                                        ggCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGt                                    >  2:2175242/1‑129 (MQ=255)
                                                                                                                                         gcagcgcagcCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTc              >  1:2153190/1‑150 (MQ=255)
                                                                                                                                          catgcgcagcCCTTTGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCa                                                                                 >  1:1459830/4‑83 (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)
                                                                                                                                                      |                                                                                                                                                     
CCTGTACCAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTCACCACGAC  >  NC_002947/2949688‑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: