Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,850,159 C→G 9.2% intergenic (+76/‑111) hemN → / → fnrA oxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,850,1590CG 9.2% 187.5 / 10.1 65intergenic (+76/‑111)hemN/fnrAoxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base G (3/3);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

GGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGC  >  NC_002947/4850017‑4850307
                                                                                                                                               |                                                                                                                                                     
ggTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGc                                                                                                                                                 <  2:280479/149‑1 (MQ=255)
    tGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCgg                                                                                                                                              <  3:346117/148‑1 (MQ=255)
      cATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            <  4:274168/148‑2 (MQ=255)
      cATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            <  2:200576/148‑2 (MQ=255)
      cATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCgggg                                                                                                                                            <  7:212310/148‑1 (MQ=255)
         ggTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTc                                                                                                                                        >  4:141748/1‑149 (MQ=255)
                   gcgGATCTGGCGACGCCCAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTTTGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAgcgcgc                                                                                                                             <  2:169212/149‑1 (MQ=255)
                         cTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACcat                                                                                                                        >  1:330864/1‑148 (MQ=255)
                            gCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACcat                                                                                                                        <  2:217179/145‑1 (MQ=255)
                              gATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATcgc                                                                                                                  >  2:270433/1‑149 (MQ=255)
                              gATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATcgc                                                                                                                  >  5:29754/1‑149 (MQ=255)
                               aTGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATcgcg                                                                                                                 <  2:31087/149‑1 (MQ=255)
                                   acaACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGcaca                                                                                                             >  6:347862/1‑149 (MQ=255)
                                   acaACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGcaca                                                                                                             >  4:184564/1‑149 (MQ=255)
                                        cGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGc                                                                                                        <  3:184564/149‑1 (MQ=255)
                                         gCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGcc                                                                                                       >  7:54879/1‑149 (MQ=255)
                                         gCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGcc                                                                                                       >  4:193877/1‑149 (MQ=255)
                                          ccAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCt                                                                                                      <  3:141748/149‑1 (MQ=255)
                                            aGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTa                                                                                                     >  8:197114/1‑148 (MQ=255)
                                                    tcGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGACACAATATGAGCCAACAAGGAAACCCGGCAGGCCACATCCCGGGGTTACGCCACCGCGCACCAACATCGGGCACAATGGCCTACACCAtga                                                                                             >  1:101475/1‑148 (MQ=255)
                                                    tcGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCtgtg                                                                                            >  3:250988/1‑149 (MQ=255)
                                                        gAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGt                                                                                        <  5:154169/149‑1 (MQ=255)
                                                        gAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCACACCATCATCGCGCACACTGGCCTACACCCTGTGTCGt                                                                                        >  2:304520/1‑149 (MQ=255)
                                                         aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGcaca                                                                                                             >  5:222149/1‑127 (MQ=255)
                                                         aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGcaca                                                                                                             <  6:222149/127‑1 (MQ=255)
                                                         aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTg                                                                                       <  8:54879/149‑1 (MQ=255)
                                                          gTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                      <  7:202100/149‑1 (MQ=255)
                                                          gTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                      <  2:330864/149‑1 (MQ=255)
                                                          gTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                      <  5:347862/149‑1 (MQ=255)
                                                           tGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCGAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCTTGa                                                                                      <  2:6175/148‑1 (MQ=255)
                                                             aTCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGt                                                                                    <  7:197114/148‑1 (MQ=255)
                                                                     aaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCACGCCACCTGCCGGGGTTTCGCCAGCGCTCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTAc                                                                           >  1:19649/1‑149 (MQ=255)
                                                                       cAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGa                                                                         <  6:29754/149‑1 (MQ=255)
                                                                       cAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGa                                                                         >  2:303159/1‑149 (MQ=255)
                                                                                  aCCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGc                                                              >  3:275426/1‑149 (MQ=255)
                                                                                      cATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTc                                                          <  1:77910/149‑1 (MQ=255)
                                                                                          gACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGAGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCcac                                                      <  2:101475/149‑1 (MQ=255)
                                                                                           aCAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAt                                                                    >  6:369267/1‑134 (MQ=255)
                                                                                           aCAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAt                                                                    <  5:369267/134‑1 (MQ=255)
                                                                                            cAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCcac                                                      >  3:19390/1‑147 (MQ=255)
                                                                                            cAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACaa                                                    <  5:129417/149‑1 (MQ=255)
                                                                                             aaGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAg                                                   >  8:32775/1‑149 (MQ=255)
                                                                                                 ttCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGCCTTATGTGTGCTTTCCCACAAGGATc                                               >  8:186357/1‑149 (MQ=255)
                                                                                                   ccTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCg                                              <  1:303159/148‑1 (MQ=255)
                                                                                                     ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGACCAGCGCGCACCATCATAGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTTTGCTTTCCCACAAGGATCGAt                                            >  1:285755/1‑149 (MQ=255)
                                                                                                       gTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGAtt                                           <  7:32775/147‑1 (MQ=255)
                                                                                                       gTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                          <  3:193877/148‑1 (MQ=255)
                                                                                                        tCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCa                                         >  8:262013/1‑148 (MQ=255)
                                                                                                        tCAGCCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTa                                                                            <  2:199233/113‑1 (MQ=255)
                                                                                                           ggCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATg                                     >  5:239599/1‑149 (MQ=255)
                                                                                                                 acTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTcc                                  >  1:50214/1‑147 (MQ=255)
                                                                                                                        gCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTc                         <  6:239599/149‑1 (MQ=255)
                                                                                                                        gCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGCCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCtatt                                                          >  6:283813/1‑113 (MQ=255)
                                                                                                                           cacaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCaaa                      <  4:250988/149‑1 (MQ=255)
                                                                                                                            acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAAc                     >  5:210738/1‑149 (MQ=255)
                                                                                                                            acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAAc                     >  1:245635/1‑149 (MQ=255)
                                                                                                                            acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAAc                     >  6:44071/1‑149 (MQ=255)
                                                                                                                             caAGGACTCCCGGCAGGCCAGCTGCCGGGGTTTCGCCAGCGCGCACCATCTTTCTGCACTCTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGAGTTATGTGTGCTTTCCCTCAAGGATCGATTCAAATGTCCGAGCCAGTCAAACt                    <  1:302171/149‑1 (MQ=255)
                                                                                                                                 aaaccccGGCAGGTGGCCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCAACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTgcgc                >  3:313459/5‑149 (MQ=255)
                                                                                                                                    tCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCca             >  6:288182/1‑149 (MQ=255)
                                                                                                                                     cccGGCAGGTGGCCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCca             >  5:40634/1‑148 (MQ=255)
                                                                                                                                     cccGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCcac            >  8:151558/1‑149 (MQ=255)
                                                                                                                                      ccGGCAGGTGGCCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCcaca           >  2:309665/1‑149 (MQ=255)
                                                                                                                                        ggCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGcc               >  7:64542/1‑143 (MQ=255)
                                                                                                                                        ggCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGcc               <  8:64542/143‑1 (MQ=255)
                                                                                                                                               cACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGc  <  2:245635/149‑1 (MQ=255)
                                                                                                                                               cACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGc  <  4:275426/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGC  >  NC_002947/4850017‑4850307

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: