Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,721,633 C→T 11.1% intergenic (‑38/+26) PP_4178 ← / ← htpG dienelactone hydrolase family protein/chaperone protein HtpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,721,6330CT11.1% 165.3 / 9.3 63intergenic (‑38/+26)PP_4178/htpGdienelactone hydrolase family protein/chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base T (3/4);  total (31/32)
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.96e-01

ATCACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGAAACG  >  NC_002947/4721489‑4721781
                                                                                                                                                |                                                                                                                                                    
aTCACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAg                                                                                                                                                  >  1:116365/1‑149 (MQ=255)
    ccAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCggg                                                                                                                                              <  5:197166/149‑1 (MQ=255)
    ccAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCggg                                                                                                                                              >  6:222532/1‑149 (MQ=255)
    ccAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCggg                                                                                                                                              <  2:348043/149‑1 (MQ=255)
     cAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGc                                                                                                                                             >  1:234271/1‑149 (MQ=255)
      aGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCt                                                                                                                                            <  5:38735/149‑1 (MQ=255)
       gCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCtt                                                                                                                                           <  7:297927/149‑1 (MQ=255)
        cccTGGGCGGCCCAGCGCACCCGTCTCATAGACCAAACGGCTTTCGTAGTTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCttt                                                                                                                                          <  8:75994/149‑1 (MQ=255)
         ccTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCtttt                                                                                                                                         <  4:50753/149‑2 (MQ=255)
         ccTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTg                                                                                                                                         <  2:317990/149‑1 (MQ=255)
                  cccAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCttttttat                                                                                                                                     <  8:205136/144‑6 (MQ=255)
                  cccAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCttttttat                                                                                                                                     >  7:205136/1‑139 (MQ=255)
                  cccAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATc                                                                                                                                <  2:234271/149‑1 (MQ=255)
                    cAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGc                                                                                                                                                 >  7:251373/1‑130 (MQ=255)
                    cAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGc                                                                                                                                                 <  8:251373/130‑1 (MQ=255)
                    cAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCa                                                                                                                               >  5:295776/1‑148 (MQ=255)
                      gcgcACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACt                                                                                                                            <  7:63343/149‑1 (MQ=255)
                       cgcACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACt                                                                                                                            >  5:130803/1‑148 (MQ=255)
                         cACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCaa                                                                                                                         <  6:295776/149‑1 (MQ=255)
                            ccGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCggg                                                                                                                                              <  8:267532/125‑1 (MQ=255)
                            ccGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCggg                                                                                                                                              >  7:267532/1‑125 (MQ=255)
                            ccGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCa                                                                                                                      <  2:116365/149‑1 (MQ=255)
                              ggCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGa                                                                                                                    <  5:222532/149‑1 (MQ=255)
                                     aGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGAACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTc                                                                                                              >  2:144269/1‑148 (MQ=255)
                                       aCCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGAAAGCTCGAc                                                                                                           >  7:319959/1‑149 (MQ=255)
                                          aaaCGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAg                                                                                                        >  4:106156/1‑149 (MQ=255)
                                          aaaCGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAg                                                                                                        <  3:106156/149‑1 (MQ=255)
                                                cTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCtttttta                                                                                                                                      >  6:88543/1‑109 (MQ=255)
                                                cTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCtttttta                                                                                                                                      <  5:88543/113‑5 (MQ=255)
                                                     gTAGGTCTTGCCAGACAGGTGATAAACCATCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCa                                                                                             <  4:158051/149‑1 (MQ=255)
                                                     gTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGtt                                                                                               >  7:319621/1‑147 (MQ=255)
                                                     gTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTc                                                                                              <  1:110768/148‑1 (MQ=255)
                                                               ccAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAAc                                                                                   <  6:130803/149‑1 (MQ=255)
                                                                     aGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAgc                                                                             >  5:87891/1‑149 (MQ=255)
                                                                       gTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAg                                                                                            >  1:330247/1‑132 (MQ=255)
                                                                       gTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAg                                                                                            <  2:330247/132‑1 (MQ=255)
                                                                          aTAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGccgccg                                                                        >  1:57124/1‑149 (MQ=255)
                                                                              aCCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAg                                                                                                                        <  5:164123/97‑1 (MQ=255)
                                                                              aCCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAg                                                                                                                        >  6:164123/1‑97 (MQ=255)
                                                                                 aGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCtt                                                                 <  8:319959/149‑1 (MQ=255)
                                                                                      ttCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGc                                                            <  1:144269/149‑1 (MQ=255)
                                                                                         gacgaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGt                                                         >  8:97476/1‑149 (MQ=255)
                                                                                         gacgaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGt                                                         >  6:300523/1‑149 (MQ=255)
                                                                                                  tttGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCggg                                                                      <  8:282593/127‑1 (MQ=255)
                                                                                                  tttGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCggg                                                                      >  7:282593/1‑127 (MQ=255)
                                                                                                       tcatcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAgg                                           >  3:276446/1‑149 (MQ=255)
                                                                                                         atcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGgc                                         >  4:179258/1‑149 (MQ=255)
                                                                                                         atcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCAAAAAATAAAGCAGAAAGGTGGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGGCCGGTACTTACGGCTTTAACCGGCCCCCAGGgc                                         >  1:103551/1‑149 (MQ=255)
                                                                                                          tcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGgcc                                        >  1:62667/1‑149 (MQ=255)
                                                                                                                 cTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGAt                                 <  7:231832/149‑1 (MQ=255)
                                                                                                                 cTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGAt                                 <  8:319621/149‑1 (MQ=255)
                                                                                                                  tCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGAt                                 >  3:348284/1‑148 (MQ=255)
                                                                                                                  tCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGAt                                 <  4:348284/148‑1 (MQ=255)
                                                                                                                    cTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCg                               >  3:137269/1‑148 (MQ=255)
                                                                                                                           aCTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGAt                       <  7:97476/149‑1 (MQ=255)
                                                                                                                               gTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGAtgtg                    <  2:57124/148‑1 (MQ=255)
                                                                                                                               gTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGAtgtg                    <  1:143202/148‑1 (MQ=255)
                                                                                                                                  aTAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGa                 >  8:317325/1‑148 (MQ=255)
                                                                                                                                   tAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGa                 >  6:272463/1‑147 (MQ=255)
                                                                                                                                   cagcaaaGCCCGCTTCGGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGa                 >  7:104771/5‑147 (MQ=255)
                                                                                                                                    aaaaaaGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAg              >  4:117835/1‑149 (MQ=255)
                                                                                                                                     aaaaaGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGgccgcc                                              <  2:208942/116‑1 (MQ=255)
                                                                                                                                     aaaaaGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGgccgcc                                              >  1:208942/1‑116 (MQ=255)
                                                                                                                                              gccgAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGc                                                                                          >  7:93782/1‑63 (MQ=255)
                                                                                                                                              gccgAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGc                                                                                          <  8:93782/63‑1 (MQ=255)
                                                                                                                                                cgAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGAAACg  <  5:272463/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ATCACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGAAACG  >  NC_002947/4721489‑4721781

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