Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,681,711 C→G 6.1% A63G (GCC→GGC)  recJ → ssDNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,681,7110CG 6.1% 300.9 / 9.7 99A63G (GCC→GGC) recJssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base C (50/42);  new base G (3/3);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

TTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTACG  >  NC_002947/1681563‑1681857
                                                                                                                                                      |                                                                                                                                                  
tttCGGGGTACCCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGcc                                                                                                                                                   <  2:1171545/150‑1 (MQ=255)
 ttCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCa                                                                                                                                                  <  1:137267/150‑1 (MQ=255)
  tCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCAt                                                                                                                                                 <  1:2089054/150‑1 (MQ=255)
   cTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATc                                                                                                                                                <  1:1120908/150‑1 (MQ=255)
      gggcaCCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAc                                                                                                                                             >  2:1162425/5‑150 (MQ=255)
      ggTAACCCGCCCCCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGCGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAc                                                                                                                                             <  1:1415868/150‑1 (MQ=255)
       gTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGAGAGGAGGCGCTGTACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAcc                                                                                                                                            >  2:2069615/1‑150 (MQ=255)
          aCCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagc                                                                                                                                         <  1:1284974/150‑1 (MQ=255)
           ccTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcg                                                                                                                                        >  2:2202101/1‑150 (MQ=255)
             tGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgcc                                                                                                                                      >  2:1524560/1‑150 (MQ=255)
                   ccTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgcc                                                                                                                                      >  1:826556/1‑144 (MQ=255)
                   ccTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgcc                                                                                                                                      <  2:826556/144‑1 (MQ=255)
                     ttGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCAGGCCTGCTGCCCTACCAGCAGCTCAAGGGCAGCGAAGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATc                                                                                                                              >  2:851498/1‑150 (MQ=255)
                     ttGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATc                                                                                                                              >  1:1260458/1‑150 (MQ=255)
                       gCTGACCCGCCTTTACGCCGCGCGCGGNGTGCAGAGCGAGNCCGAGCNGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCt                                                                                                                            >  1:1872248/1‑150 (MQ=255)
                          gACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGa                                                                                                                                              >  2:2208517/1‑129 (MQ=255)
                          gACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGa                                                                                                                                              <  1:2208517/129‑1 (MQ=255)
                          gACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAt                                                                                                                         <  2:1472072/150‑1 (MQ=255)
                             ccGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCCGCTCAAGGGCATCACAGCCGGTGCGGGACCAGCGGCAGAATCCAGCGACCGTCGCCATCGCCTcccacccga                                                                                                                      >  2:23075/1‑143 (MQ=255)
                             ccGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcgt                                                                                                                      >  1:118291/1‑150 (MQ=255)
                             ccGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcgt                                                                                                                      >  1:2169564/1‑150 (MQ=255)
                               gCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcgtcc                                                                                                                    >  2:531355/1‑149 (MQ=255)
                               gCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGACGCCCATGTGGCCCTCCTGGTCGAAGCCATCGACCAGCGCCAGTGCTTCCTCAgcgtcg                                                                                                                    <  2:615856/150‑1 (MQ=255)
                                tctgcaCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCgg                                                                                                                   >  1:536751/6‑150 (MQ=255)
                                tctgcaCGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCgg                                                                                                                   >  1:1296206/6‑150 (MQ=255)
                                   ttACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGa                                                                                                                <  1:2202101/150‑1 (MQ=255)
                                       gccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccagcaggtc                                                                                                                                    <  2:1700261/126‑6 (MQ=255)
                                       gccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccagcaggtc                                                                                                                                    >  1:1700261/1‑121 (MQ=255)
                                       gccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCa                                                                                                                                                  <  1:867411/112‑1 (MQ=255)
                                       gccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCa                                                                                                                                                  >  2:867411/1‑112 (MQ=255)
                                       gccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTc                                                                                                            >  2:1327983/1‑150 (MQ=255)
                                       gccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCGCTAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTc                                                                                                            >  1:1330512/1‑150 (MQ=255)
                                         cgcgcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGa                                                                                                          <  1:495322/150‑1 (MQ=255)
                                          gcgcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         <  2:1260458/150‑1 (MQ=255)
                                            gcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTCGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGc                                                                                                       >  1:2223607/1‑150 (MQ=255)
                                              gcggcgTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTc                                                                                                                           >  1:357633/1‑128 (MQ=255)
                                              gcggcgTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTc                                                                                                                           <  2:357633/128‑1 (MQ=255)
                                                    tGCAGAGCGAGGCCGAGCTGGACAAGAGCCTCCCGTGCCTTCTGCTCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTg                                                                                               <  2:775928/150‑1 (MQ=255)
                                                     gCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTTTTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGc                                                                                              <  1:23075/150‑1 (MQ=255)
                                                        gagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCCGCTGGTCGATGGCTTCGAccagcag                                                                                                                                       <  2:955402/106‑3 (MQ=255)
                                                        gagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCCGCTGGTCGATGGCTTCGAccagcag                                                                                                                                       >  1:955402/1‑104 (MQ=255)
                                                               gCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGc                                                                                    >  2:1179660/1‑150 (MQ=255)
                                                                     cTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTc                                                                              <  1:2283456/150‑1 (MQ=255)
                                                                     cTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTc                                                                              <  2:473361/150‑1 (MQ=255)
                                                                     cTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATACCCTCGGTGCCACCGCCAGCACCGTc                                                                              <  1:1484148/150‑1 (MQ=255)
                                                                       ggACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccag                                                                                                                                          <  1:393028/88‑1 (MQ=255)
                                                                       ggACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGATGGCTTCGAccag                                                                                                                                          >  2:393028/1‑88 (MQ=255)
                                                                        gACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGCCGCCTCCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCggg                                                                           >  1:744332/1‑150 (MQ=255)
                                                                           aaGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGt                                                                         <  1:569100/149‑1 (MQ=255)
                                                                           aaGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGt                                                                         >  2:569100/1‑149 (MQ=255)
                                                                           aaGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTg                                                                        >  1:134522/1‑150 (MQ=255)
                                                                            agagCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGc                                                                       <  1:1524560/150‑1 (MQ=255)
                                                                               gCCTGGCACGCCTGCTGCCCTACCAGCCGATCAAAGGGACTGGAGTTCCGGTGTATCTGCTGGTCGAAGCCCAAGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTgg                                                                    <  1:1685485/150‑1 (MQ=255)
                                                                                ccTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTCCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTggg                                                                   >  2:1092271/1‑150 (MQ=255)
                                                                                ccTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGCTGCATCGCTCTTGCAGCGCCAGCGCATCCTCATCGTTGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTggg                                                                   <  1:1571890/150‑1 (MQ=255)
                                                                                  tGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCa                                                                                      <  2:684381/129‑1 (MQ=255)
                                                                                  tGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCa                                                                                      >  1:684381/1‑129 (MQ=255)
                                                                                         ccTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgcgcct                                                          <  2:1872248/150‑1 (MQ=255)
                                                                                         ccTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcct                                                                <  2:1016749/144‑1 (MQ=255)
                                                                                         ccTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGCGCTGGgcct                                                                >  1:1016749/1‑144 (MQ=255)
                                                                                            gctgcCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCctgct                                                       >  2:869966/1‑150 (MQ=255)
                                                                                            gctgcCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAACACCGTCGGGGTGCTGGGCCTGCGCctgct                                                       <  2:2169564/150‑1 (MQ=255)
                                                                                             ctgcCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCctgctg                                                      <  1:647492/150‑1 (MQ=255)
                                                                                                cccTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                   <  1:1179660/150‑1 (MQ=255)
                                                                                                 ccTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTg                                                  >  1:1692161/1‑150 (MQ=255)
                                                                                                   tACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCg                                                >  1:1183347/1‑150 (MQ=255)
                                                                                                   tACCAGCAGCCCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCg                                                >  2:718649/1‑150 (MQ=255)
                                                                                                    aCCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCgg                                               <  2:536751/150‑1 (MQ=255)
                                                                                                      cagcagCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCt                                                                      >  2:973394/1‑125 (MQ=255)
                                                                                                      cagcagCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCt                                                                      <  1:973394/125‑1 (MQ=255)
                                                                                                             tCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCg                                      >  2:102817/1‑150 (MQ=255)
                                                                                                             tCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCg                                      >  2:374032/1‑150 (MQ=255)
                                                                                                                  gggCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTa                                  >  1:934393/1‑150 (MQ=255)
                                                                                                                  gggCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTa                                  >  2:1000142/1‑150 (MQ=255)
                                                                                                                    gCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTAcc                                >  2:1872263/1‑150 (MQ=255)
                                                                                                                     cagcGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCt                               <  1:851498/147‑1 (MQ=255)
                                                                                                                      aTCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTg                              >  1:1706190/1‑150 (MQ=255)
                                                                                                                          aGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         <  1:2265811/71‑1 (MQ=255)
                                                                                                                          aGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                         >  2:2265811/1‑71 (MQ=255)
                                                                                                                              cGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCaa                      >  1:1242663/1‑150 (MQ=255)
                                                                                                                               ggTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAAc                     <  2:1183347/150‑1 (MQ=255)
                                                                                                                                      gACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTAAAGGCCGGCGGTGGCACGGCGGGCACCGTCGGGGTGCTGGGACTGCGCCTTCTGGGTGCTGGTCATGTCGATTACCTGGTGCCCAACCGCTTc               <  1:531355/150‑1 (MQ=255)
                                                                                                                                      gACCTGCTGGTCGAAGACATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTc               <  1:1872263/150‑1 (MQ=255)
                                                                                                                                        ccTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCCACCTCGATGCCGACGGTGCCACCGCCAGCACCGTCGCGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGa             <  2:118291/150‑1 (MQ=255)
                                                                                                                                         ctgctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCAcc                                                                                 >  2:1268165/1‑81 (MQ=255)
                                                                                                                                         ctgctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCAcc                                                                                 <  1:1268165/81‑1 (MQ=255)
                                                                                                                                         ctgctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAg            >  1:2132060/1‑150 (MQ=255)
                                                                                                                                         ctgctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAg            >  1:1412818/1‑150 (MQ=255)
                                                                                                                                           gctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTa          >  2:1534442/1‑150 (MQ=255)
                                                                                                                                           gctgGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTa          >  2:1457818/1‑150 (MQ=255)
                                                                                                                                           gctgGGCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTa          >  1:98150/1‑150 (MQ=255)
                                                                                                                                                tCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATg        <  1:1697472/147‑1 (MQ=255)
                                                                                                                                                tCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATg        >  2:1697472/1‑147 (MQ=255)
                                                                                                                                                  gAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCctgct                                                       <  2:217772/98‑1 (MQ=255)
                                                                                                                                                  gAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCctgct                                                       >  1:217772/1‑98 (MQ=255)
                                                                                                                                                  gAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTAcc                                <  2:2062299/121‑1 (MQ=255)
                                                                                                                                                  gAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTAcc                                >  1:2062299/1‑121 (MQ=255)
                                                                                                                                                   aaGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTACg  <  2:1412818/150‑1 (MQ=255)
                                                                                                                                                    aGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGt                                                                                                <  2:1440821/55‑1 (MQ=255)
                                                                                                                                                    aGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGt                                                                                                >  1:1440821/1‑55 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCTACG  >  NC_002947/1681563‑1681857

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