Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,241,531 C→T 9.3% intergenic (+23/‑69) 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,5310CT 9.3% 178.0 / 7.0 65intergenic (+23/‑69)PP_1080/grxDmolybdopterin‑binding oxidoreductase/monothiol glutaredoxin
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base T (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 = 8.40e-01

CTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGCC  >  NC_002947/1241386‑1241670
                                                                                                                                                 |                                                                                                                                             
cTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCt                                                                                                                                            >  2:171908/1‑149 (MQ=255)
  gCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAAGCCGAGCGCCGTGCTCg                                                                                                                                          >  3:268165/1‑149 (MQ=255)
   cAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTc                                                                                                                                           <  8:22015/147‑1 (MQ=255)
    agagCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGG‑c                                                                                                                                       <  7:138428/149‑1 (MQ=255)
       tCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAAGCCGAGCGCCGTGCTCGG‑Cctt                                                                                                                                    <  4:268165/149‑1 (MQ=255)
         gCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTc                                                                                                                                   <  4:61207/148‑1 (MQ=255)
         gCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTc                                                                                                                                   <  3:220150/148‑1 (MQ=255)
          cccTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGGAACGGGGCGGGGCATGGCGTAACCGGGCAGGTCGAGGGCGGGTGAGGGACGGCAAGGCCCGGCACGGCGCACCG‑CATTCCg                                                                                                                                 >  8:134681/1‑149 (MQ=255)
          cccTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGG‑Ccttgcc                                                                                                                                 <  3:238592/149‑4 (MQ=255)
          cccTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTcc                                                                                                                                  >  5:97452/1‑148 (MQ=255)
               gcgcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATa                                                                                                                              <  6:39779/147‑1 (MQ=255)
               gcgcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATa                                                                                                                              <  2:145514/147‑1 (MQ=255)
               ccccAACCCGGCGTCCGTCCCAATGCCCTCACCGATCCGCCCTTTCGGCACCTGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATa                                                                                                                              <  7:134681/144‑1 (MQ=255)
                cgcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGTCCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACaa                                                                                                                           <  6:134907/149‑1 (MQ=255)
                   aaCCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGc                                                                                                                        >  2:207534/1‑149 (MQ=255)
                         gCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACAcc                                                                                                                  >  3:37098/1‑149 (MQ=255)
                             gAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGc                                                                                                              >  5:307536/1‑149 (MQ=255)
                               gTGCCATTGACCTGACCGATGAGCTCTTTCGGGAGCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGAAGGTCGAGGCCGCTTGAGTAACGGCAAGGCCGAGCACGGCGGTCGG‑CTTTCCGATACAATGCGACACCGTGcc                                                                                                             <  6:6640/148‑1 (MQ=255)
                                   cAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCcgacga                                                                                                        >  6:17034/1‑149 (MQ=255)
                                   cAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCcgacga                                                                                                        >  1:213930/1‑149 (MQ=255)
                                          ctgaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACACGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTc                                                                                                 <  2:213930/149‑1 (MQ=255)
                                            gaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGcc                                                                                                             >  3:72525/1‑135 (MQ=255)
                                            gaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGcc                                                                                                             <  4:72525/135‑1 (MQ=255)
                                                        cTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGcc                                                                                    >  4:59626/1‑148 (MQ=255)
                                                           tCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGa                                                                                  <  3:308220/147‑1 (MQ=255)
                                                           tCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAg                                                                                 <  5:17034/148‑1 (MQ=255)
                                                           tCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAg                                                                                 <  1:207534/148‑1 (MQ=255)
                                                              ggACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGc                                                                             >  4:69856/1‑149 (MQ=255)
                                                               gACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGc                                                                             >  7:291846/1‑148 (MQ=255)
                                                                aCCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTa                                                                           <  6:307536/149‑1 (MQ=255)
                                                                   agctttCAGGAAAAGCTGCGCTCAATGGAGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGACCACGGCGCTCGG‑CTTTCCGATACAATGCGACACAGTCCCGACGACAGTGTCGGATAGTTCAGACGAGGTGCTACAt                                                                        >  8:180189/4‑149 (MQ=255)
                                                                       ttCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGa                                                                     >  1:316391/1‑148 (MQ=255)
                                                                        tCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                    <  2:316391/148‑1 (MQ=255)
                                                                        tCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                    <  7:301218/148‑1 (MQ=255)
                                                                           ggAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCTGATAGTTCAGCCGAg                                                                                 <  7:288929/132‑1 (MQ=255)
                                                                           ggAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAg                                                                                 >  8:288929/1‑132 (MQ=255)
                                                                            gAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATat‑cat                                                              <  3:69856/149‑1 (MQ=255)
                                                                                     cgcTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGa                                                      <  4:37098/148‑1 (MQ=255)
                                                                                      gcTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                    <  4:87343/134‑1 (MQ=255)
                                                                                      gcTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                    >  3:87343/1‑134 (MQ=255)
                                                                                        tCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCa                                                   >  6:309389/1‑148 (MQ=255)
                                                                                          aaTGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCa                                                   >  4:88386/1‑146 (MQ=255)
                                                                                          aaTGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCa                                                   >  4:274267/1‑146 (MQ=255)
                                                                                          aaTGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCa                                                   >  2:94072/1‑146 (MQ=255)
                                                                                                gTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAagag                                              >  7:173958/1‑145 (MQ=255)
                                                                                                   cccGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGAtt                                        <  5:156334/148‑1 (MQ=255)
                                                                                                   cccGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTg                                       <  3:88386/149‑1 (MQ=255)
                                                                                                    ccGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGGGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGc                                      >  1:229426/1‑149 (MQ=255)
                                                                                                      gTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGcc                                                                                    >  4:8069/1‑102 (MQ=255)
                                                                                                      gTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGcc                                                                                    <  3:8069/102‑1 (MQ=255)
                                                                                                      gTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGcc                                     >  2:67820/1‑148 (MQ=255)
                                                                                                       tGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGcc                                     <  1:160461/147‑1 (MQ=255)
                                                                                                       tGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGcc                                     >  2:160461/1‑147 (MQ=255)
                                                                                                             tCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAacac                             >  2:15065/1‑149 (MQ=255)
                                                                                                                   ccGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCt                       <  3:59626/149‑1 (MQ=255)
                                                                                                                       ttGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTATGCtt                   >  1:124570/1‑149 (MQ=255)
                                                                                                                         gatgAACGGCAAGGCCGATCACGGCGCTCGG‑CTTTACGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTAAATGGATATCCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCtt                   <  1:67820/145‑1 (MQ=255)
                                                                                                                         gAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTa                 >  7:148713/1‑149 (MQ=255)
                                                                                                                            gAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCAGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGAACAAAGAGCAGAATGCCAACAACACCATTCTCCTTTACAt              >  1:219252/1‑149 (MQ=255)
                                                                                                                             aaCGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACAt              >  4:147101/1‑148 (MQ=255)
                                                                                                                               tcggaaagCCGAGCGCCGTGCTCGGCCTTGCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATg             >  7:183731/8‑148 (MQ=255)
                                                                                                                                ggCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGaaa          <  8:291846/149‑1 (MQ=255)
                                                                                                                                ggCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGaaa          <  6:97452/149‑1 (MQ=255)
                                                                                                                                cggaaagCCGAGCGCCGTGCTCGG‑CTTTCCGATACAATGTGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGGGCTACATGGATAG‑CATCGAAACGATCAAAGAGCCGATTGCCAACAACACCATTCTGCTGTACAGGaac          >  4:55501/7‑148 (MQ=255)
                                                                                                                                        cGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGTTAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAACAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGcc  <  7:180189/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGG‑CTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATAT‑CATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGCC  >  NC_002947/1241386‑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: ‑