Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,241,529 G→C 9.2% intergenic (+21/‑71) PP_1080 → / → grxD molybdopterin‑binding oxidoreductase/monothiol glutaredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,241,5290GC 9.2% 182.1 / 9.6 65intergenic (+21/‑71)PP_1080/grxDmolybdopterin‑binding oxidoreductase/monothiol glutaredoxin
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base C (3/3);  total (33/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 = 5.81e-01

GCACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGCC  >  NC_002947/1241382‑1241670
                                                                                                                                                   |                                                                                                                                               
cagccTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACgg                                                                                                                                                <  7:35982/146‑1 (MQ=255)
    cTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCt                                                                                                                                            >  2:171908/1‑149 (MQ=255)
      gCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAAGCCGAGCGCCGTGCTCg                                                                                                                                          >  3:268165/1‑149 (MQ=255)
       cAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTc                                                                                                                                           <  8:22015/147‑1 (MQ=255)
        agagCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGGc                                                                                                                                       <  7:138428/149‑1 (MQ=255)
           tCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAAGCCGAGCGCCGTGCTCGGCctt                                                                                                                                    <  4:268165/149‑1 (MQ=255)
             gCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTc                                                                                                                                   <  4:61207/148‑1 (MQ=255)
             gCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTc                                                                                                                                   <  3:220150/148‑1 (MQ=255)
              cccTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGGAACGGGGCGGGGCATGGCGTAACCGGGCAGGTCGAGGGCGGGTGAGGGACGGCAAGGCCCGGCACGGCGCACCGCATTCCg                                                                                                                                 >  8:134681/1‑149 (MQ=255)
              cccTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGGCcttgcc                                                                                                                                 <  3:238592/149‑4 (MQ=255)
              cccTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTcc                                                                                                                                  >  5:97452/1‑148 (MQ=255)
                   gcgcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATa                                                                                                                              <  6:39779/147‑1 (MQ=255)
                   gcgcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATa                                                                                                                              <  2:145514/147‑1 (MQ=255)
                   ccccAACCCGGCGTCCGTCCCAATGCCCTCACCGATCCGCCCTTTCGGCACCTGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATa                                                                                                                              <  7:134681/144‑1 (MQ=255)
                    cgcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGTCCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACaa                                                                                                                           <  6:134907/149‑1 (MQ=255)
                       aaCCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGc                                                                                                                        >  2:207534/1‑149 (MQ=255)
                             gCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACAcc                                                                                                                  >  3:37098/1‑149 (MQ=255)
                                 gAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGc                                                                                                              >  5:307536/1‑149 (MQ=255)
                                   gTGCCATTGACCTGACCGATGAGCTCTTTCGGGAGCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGAAGGTCGAGGCCGCTTGAGTAACGGCAAGGCCGAGCACGGCGGTCGGCTTTCCGATACAATGCGACACCGTGcc                                                                                                             <  6:6640/148‑1 (MQ=255)
                                       cAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCcgacga                                                                                                        >  1:213930/1‑149 (MQ=255)
                                       cAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCcgacga                                                                                                        >  6:17034/1‑149 (MQ=255)
                                              ctgaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACACGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTc                                                                                                 <  2:213930/149‑1 (MQ=255)
                                                gaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGcc                                                                                                             <  4:72525/135‑1 (MQ=255)
                                                gaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGcc                                                                                                             >  3:72525/1‑135 (MQ=255)
                                                            cTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGcc                                                                                    >  4:59626/1‑148 (MQ=255)
                                                               tCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGa                                                                                  <  3:308220/147‑1 (MQ=255)
                                                               tCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAg                                                                                 <  5:17034/148‑1 (MQ=255)
                                                               tCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAg                                                                                 <  1:207534/148‑1 (MQ=255)
                                                                  ggACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGc                                                                             >  4:69856/1‑149 (MQ=255)
                                                                   gACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGc                                                                             >  7:291846/1‑148 (MQ=255)
                                                                    aCCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTa                                                                           <  6:307536/149‑1 (MQ=255)
                                                                       agctttCAGGAAAAGCTGCGCTCAATGGAGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGACCACGGCGCTCGGCTTTCCGATACAATGCGACACAGTCCCGACGACAGTGTCGGATAGTTCAGACGAGGTGCTACAt                                                                        >  8:180189/4‑149 (MQ=255)
                                                                           ttCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGa                                                                     >  1:316391/1‑148 (MQ=255)
                                                                            tCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                    <  7:301218/148‑1 (MQ=255)
                                                                            tCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                    <  2:316391/148‑1 (MQ=255)
                                                                               ggAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCTGATAGTTCAGCCGAg                                                                                 <  7:288929/132‑1 (MQ=255)
                                                                               ggAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAg                                                                                 >  8:288929/1‑132 (MQ=255)
                                                                                gAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATatcat                                                              <  3:69856/149‑1 (MQ=255)
                                                                                         cgcTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGa                                                      <  4:37098/148‑1 (MQ=255)
                                                                                          gcTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                    <  4:87343/134‑1 (MQ=255)
                                                                                          gcTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                    >  3:87343/1‑134 (MQ=255)
                                                                                            tCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCa                                                   >  6:309389/1‑148 (MQ=255)
                                                                                              aaTGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCa                                                   >  4:88386/1‑146 (MQ=255)
                                                                                              aaTGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCa                                                   >  4:274267/1‑146 (MQ=255)
                                                                                              aaTGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCa                                                   >  2:94072/1‑146 (MQ=255)
                                                                                                    gTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAagag                                              >  7:173958/1‑145 (MQ=255)
                                                                                                       cccGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGAtt                                        <  5:156334/148‑1 (MQ=255)
                                                                                                       cccGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTg                                       <  3:88386/149‑1 (MQ=255)
                                                                                                        ccGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGGGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGc                                      >  1:229426/1‑149 (MQ=255)
                                                                                                          gTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGcc                                                                                    >  4:8069/1‑102 (MQ=255)
                                                                                                          gTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGcc                                                                                    <  3:8069/102‑1 (MQ=255)
                                                                                                          gTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGcc                                     >  2:67820/1‑148 (MQ=255)
                                                                                                           tGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGcc                                     <  1:160461/147‑1 (MQ=255)
                                                                                                           tGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGcc                                     >  2:160461/1‑147 (MQ=255)
                                                                                                                 tCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAacac                             >  2:15065/1‑149 (MQ=255)
                                                                                                                       ccGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCt                       <  3:59626/149‑1 (MQ=255)
                                                                                                                           ttGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTATGCtt                   >  1:124570/1‑149 (MQ=255)
                                                                                                                             gatgAACGGCAAGGCCGATCACGGCGCTCGGCTTTACGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTAAATGGATATCCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCtt                   <  1:67820/145‑1 (MQ=255)
                                                                                                                             gAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTa                 >  7:148713/1‑149 (MQ=255)
                                                                                                                                gAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCAGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGAACAAAGAGCAGAATGCCAACAACACCATTCTCCTTTACAt              >  1:219252/1‑149 (MQ=255)
                                                                                                                                 aaCGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACAt              >  4:147101/1‑148 (MQ=255)
                                                                                                                                   tcggaaagCCGAGCGCCGTGCTCGGCCTTGCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATg             >  7:183731/8‑148 (MQ=255)
                                                                                                                                    ggCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGaaa          <  6:97452/149‑1 (MQ=255)
                                                                                                                                    ggCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGaaa          <  8:291846/149‑1 (MQ=255)
                                                                                                                                    cggaaagCCGAGCGCCGTGCTCGGCTTTCCGATACAATGTGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGGGCTACATGGATAGCATCGAAACGATCAAAGAGCCGATTGCCAACAACACCATTCTGCTGTACAGGaac          >  4:55501/7‑148 (MQ=255)
                                                                                                                                            cGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGTTAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAACAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGcc  <  7:180189/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGCC  >  NC_002947/1241382‑1241670

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