Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,449,744 G→C 6.3% intergenic (+259/‑183) nicR → / → nicA HTH‑type transcriptional repressor NicR/nicotinate dehydrogenase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,449,7440GC 6.3% 196.2 / 3.2 64intergenic (+259/‑183)nicR/nicAHTH‑type transcriptional repressor NicR/nicotinate dehydrogenase subunit A
Reads supporting (aligned to +/- strand):  ref base G (24/35);  new base C (2/2);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

CGTAGGAGCGGGTTCACCCGCGAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTT  >  NC_002947/4449599‑4449891
                                                                                                                                                  |                                                                                                                                                   
cgTAGGAGCGGGTTCACCCGCGAATGGTCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCgcgg                                                                                                                                                  >  8:228456/1‑149 (MQ=255)
    ggAGCGGGTTCACCCGCGAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCgcgggcgg                                                                                                                                              >  8:159870/1‑149 (MQ=255)
     gAGCGGGTTCACCCGCGAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGggcggc                                                                                                                                             >  1:101360/1‑149 (MQ=255)
      aGCGGGTTCACCCGCGAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGgcggcg                                                                                                                                            >  2:338531/1‑149 (MQ=255)
         gggTCCACCCGCGAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTc                                                                                                                                         >  1:201463/1‑149 (MQ=255)
          ggTTCACCCGCGAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCg                                                                                                                                        <  8:229027/149‑1 (MQ=255)
          ggTTCACCCGCGAATGGGCCGGCCCTGCTGGCCAATATCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCCCGCGCGgcgcatc                                                                                                                                        <  1:34580/149‑4 (MQ=255)
           gTTCACCCGCGAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCgg                                                                                                                                       <  4:287874/149‑1 (MQ=255)
                 ccGCGAATGGGCCGGCCCTGCTGGCCAGATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTcc                                                                                                                                  >  1:16747/1‑149 (MQ=255)
                 ccGCGAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGAGGCGCCGCGCGGGCGGCGCTCGGTTTccc                                                                                                                                 <  7:159870/149‑1 (MQ=255)
                  cgcgAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGTGGCTCGGCACAGCGCGATGCGCCGCTCGGGCGGCGCTCGGTTTCCCg                                                                                                                                <  2:267754/149‑1 (MQ=255)
                   gcgAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCgg                                                                                                                               <  4:232345/149‑1 (MQ=255)
                     gAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGgc                                                                                                                             <  7:52674/149‑1 (MQ=255)
                            ccggccCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGAc                                                                                                                       >  8:53538/1‑148 (MQ=255)
                               gccCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCAGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATc                                                                                                                   <  3:262222/149‑1 (MQ=255)
                                 ccTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCtt                                                                                                                 <  2:2681/149‑1 (MQ=255)
                                     ctgGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCa                                                                                                               <  3:347591/147‑1 (MQ=255)
                                        gCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCt                                                                                                          >  5:191825/1‑149 (MQ=255)
                                          cAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGc                                                                                                        >  6:67119/1‑149 (MQ=255)
                                          cAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGc                                                                                                        <  4:7646/149‑1 (MQ=255)
                                             aTCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGccc                                                                                                      <  4:219850/149‑1 (MQ=255)
                                                 ccACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCccac                                                                                                   >  4:344830/1‑148 (MQ=255)
                                                  cACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCccacca                                                                                                 <  3:344830/149‑1 (MQ=255)
                                                  cACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCccacca                                                                                                 <  5:356340/149‑1 (MQ=255)
                                                  cACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCGACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGCCAATCTTCATCCCTGCCccacca                                                                                                 >  3:252154/1‑149 (MQ=255)
                                                   aCCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCccacca                                                                                                 >  8:11411/1‑148 (MQ=255)
                                                     cTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACcatca                                                                                              <  7:11411/149‑1 (MQ=255)
                                                     cTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACcatca                                                                                              <  2:16747/149‑1 (MQ=255)
                                                     cTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACcatca                                                                                              <  6:191825/149‑1 (MQ=255)
                                                      tGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGAGCGCCGCCCGCGCGgcgca                                                                                                                                          >  4:41588/1‑102 (MQ=25)
                                                      tGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGAGCGCCGCCCGCGCGgcgca                                                                                                                                          <  3:41588/103‑2 (MQ=25)
                                                      tGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGCGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCCGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACcatca                                                                                              >  8:34663/1‑148 (MQ=255)
                                                        ccGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATa                                                                                            <  4:312251/148‑1 (MQ=255)
                                                          gATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGc                                                                                                        >  6:103891/1‑134 (MQ=255)
                                                          gATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGc                                                                                                        <  5:103891/134‑1 (MQ=255)
                                                                        ccaccaCCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCAt                                                                           <  2:129404/149‑1 (MQ=255)
                                                                            caccGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATctct                                                                       >  7:264390/1‑149 (MQ=255)
                                                                            caccGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATctct                                                                       >  1:158011/1‑149 (MQ=255)
                                                                             accGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTg                                                                      <  1:338531/149‑1 (MQ=255)
                                                                             accGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTg                                                                      <  2:101360/149‑1 (MQ=255)
                                                                              ccGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCGCGGTTTACCGGGCGCCGACAATCTTCATCCCTGCCCCACCCACATAGTTGTGTGGAAATTAATATATgg                                                                     >  1:188052/1‑149 (MQ=255)
                                                                               cGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGc                                                                    <  5:337077/149‑1 (MQ=255)
                                                                                           tAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGct                                                        >  3:348101/1‑149 (MQ=255)
                                                                                           tAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGAGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTATCTGCAAATTCATATCTGGCCTGCTTATTGc                                                         >  6:213891/1‑148 (MQ=255)
                                                                                               cccGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGctcact                                                    >  2:87518/1‑149 (MQ=255)
                                                                                               cccGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGctcact                                                    <  5:67119/149‑1 (MQ=255)
                                                                                                      cTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTc                                              <  7:228456/148‑1 (MQ=255)
                                                                                                       tACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCAt                                            >  4:319657/1‑149 (MQ=255)
                                                                                                         cacacacaAACTCTGTGGCTCGGCACAGCGCGATTCGCCGCGCGGTCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCAt                                            <  5:213891/147‑1 (MQ=255)
                                                                                                         cacacacaAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCAt                                            <  3:79528/147‑1 (MQ=255)
                                                                                                          acacacaAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTa                                         >  4:186417/1‑149 (MQ=255)
                                                                                                            acacaAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTg                                                                                    <  6:231551/104‑1 (MQ=255)
                                                                                                            acacaAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTg                                                                                    >  5:231551/1‑104 (MQ=255)
                                                                                                              acaACCTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCGTCATCCATGCCCCACCATAATAGTTCTCGGCAAATTAATCTCTGGCCTGCTGCTGGCTCACTCATTTCAGGTAGt                                       >  8:341212/1‑147 (MQ=255)
                                                                                                              acaAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGt                                       >  3:352281/1‑147 (MQ=255)
                                                                                                                aaaCTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGt                                   <  4:348101/149‑1 (MQ=255)
                                                                                                                           cTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTGTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGaaa                        <  4:18557/149‑1 (MQ=255)
                                                                                                                              ggCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAAc                       <  8:264390/147‑1 (MQ=255)
                                                                                                                                  caGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTcgc                 <  4:173048/149‑1 (MQ=255)
                                                                                                                                       cgATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGTGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCAt            <  1:87518/149‑1 (MQ=255)
                                                                                                                                       cgATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCAt            <  7:34663/149‑1 (MQ=255)
                                                                                                                                        gATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCAtt           >  3:269449/1‑149 (MQ=255)
                                                                                                                                         gagcgcCGCCCGCGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCAttt          >  8:319739/3‑149 (MQ=255)
                                                                                                                                         aTGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCAttt          <  2:158011/149‑1 (MQ=255)
                                                                                                                                            cgccgcGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGa        <  7:53538/148‑1 (MQ=255)
                                                                                                                                             gccgaGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTCCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGa        <  1:173930/147‑1 (MQ=255)
                                                                                                                                               cgcgcgGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTcc    <  4:269449/149‑1 (MQ=255)
                                                                                                                                                 cgcgGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCtt  <  5:148904/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGTAGGAGCGGGTTCACCCGCGAATGGGCCGGCCCTGCTGGCCAATCTCCACCTGCCGATCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACACAAACTCGGGGGCTCGGCACAGCGCGATGCGCCGCGCGGGCGGCGCTCGGTTTCCCGGGCGCCGACAATCTTCATCCCTGCCCCACCATCATAGTTCTCTGCAAATTCATCTCTGGCCTGCTTCTTGCTCACTCATTTCATGTAGTGAGTACTTACCGAAACGCTCGCGTCATTTGAGTCCTT  >  NC_002947/4449599‑4449891

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: