Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,235,841 G→C 7.1% G318R (GGG→CGG)  uvrB → uvrABC excinuclease protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,235,8410GC 7.1% 174.1 / 4.0 56G318R (GGG→CGG) uvrBuvrABC excinuclease protein B
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base C (2/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

CGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTCCG  >  NC_002947/2235694‑2235986
                                                                                                                                                     |                                                                                                                                                  
cGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCccggccgg                                                                                                                                                   <  2:222710/149‑1 (MQ=255)
 gTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCggg                                                                                                                                                  <  8:94040/149‑1 (MQ=255)
 gTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCggg                                                                                                                                                  <  2:7782/149‑1 (MQ=255)
  tCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCgggg                                                                                                                                                 >  2:99427/1‑149 (MQ=255)
      gAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCCGGCCGGGcgcc                                                                                                                                             <  7:320799/149‑1 (MQ=255)
         tACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGcgccgcc                                                                                                                                          >  4:39523/1‑149 (MQ=255)
              gCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGgc                                                                                                                                                >  1:59506/1‑138 (MQ=255)
              gCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGgc                                                                                                                                                <  2:59506/138‑1 (MQ=255)
              gCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCa                                                                                                                                     >  2:237865/1‑149 (MQ=255)
               cacaAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGccc                                                                                                                                      <  5:126245/147‑1 (MQ=255)
                  aaGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCt                                                                                                                                 >  8:29711/1‑149 (MQ=255)
                    ggCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTg                                                                                                                                >  2:259464/1‑148 (MQ=255)
                        aacaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTAc                                                                                                                             >  1:90531/1‑148 (MQ=255)
                        aacaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGa                                                                                                                           >  2:229440/1‑149 (MQ=255)
                         acaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGcgccgc                                                                                                                                           <  4:262329/132‑1 (MQ=255)
                         acaaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGcgccgc                                                                                                                                           >  3:262329/1‑132 (MQ=255)
                          caaGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGAtt                                                                                                                         <  8:16505/149‑1 (MQ=255)
                            aGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTAc                                                                                                                       <  5:11555/149‑1 (MQ=255)
                               tGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTg                                                                                                                    >  7:6142/1‑149 (MQ=255)
                                    gAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTgccgc                                                                                                                <  5:321704/148‑1 (MQ=255)
                                       gcgcAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGa                                                                                                            <  1:202951/149‑1 (MQ=255)
                                         gcAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATg                                                                                                          <  6:150787/149‑1 (MQ=255)
                                           aGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCATGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATgc                                                                                                         <  1:99427/149‑1 (MQ=255)
                                                   caGCAACGCTCACGCTTCGATCTGGAGATGATCCTCGATCTGGGGTAGTGCAACGGCATCGAAAACTACTAGCGTTACCTGTCCGGGCTCCCGGCCGGGGCGCCGCCGCCCACCCTGTTCGATTACCTGCCGCCCGATGCGctgctg                                                                                                 <  4:252978/146‑1 (MQ=255)
                                                    aGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTACAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGt                                                                                               >  8:177261/1‑148 (MQ=255)
                                                         cGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATcga                                                                                         <  2:90531/149‑1 (MQ=255)
                                                              ccGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATgcgc                                                                                                       <  6:177792/131‑1 (MQ=255)
                                                              ccGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATgcgc                                                                                                       >  5:177792/1‑131 (MQ=255)
                                                                     gATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGt                                                                             <  8:6142/149‑1 (MQ=255)
                                                                           gagaTGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCa                                                                           <  1:237865/145‑1 (MQ=255)
                                                                                    ccTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCGTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAg                                                                >  1:76562/1‑149 (MQ=255)
                                                                                    ccTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAg                                                                <  2:76562/148‑1 (MQ=255)
                                                                                     cTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGt                                                              <  1:229440/149‑1 (MQ=255)
                                                                                      tCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTc                                                             >  5:174438/1‑149 (MQ=255)
                                                                                      tCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCCGGCCGGGcgccc                                                                                                                                            <  3:199053/71‑2 (MQ=39)
                                                                                      tCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCCGGCCGGGcgccc                                                                                                                                            >  4:199053/1‑70 (MQ=39)
                                                                                        gAGCTGGGCTACTACAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTc                                                             <  7:177261/147‑1 (MQ=255)
                                                                                               gCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGctgctg                                                                                                 <  1:304154/104‑1 (MQ=255)
                                                                                               gCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGctgctg                                                                                                 >  2:304154/1‑104 (MQ=255)
                                                                                                cTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCAGCGCCATGt                                                   <  7:63199/149‑1 (MQ=255)
                                                                                                    tGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTAc                                                 >  8:240886/1‑147 (MQ=255)
                                                                                                      cAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAgg                                             <  7:240886/149‑1 (MQ=255)
                                                                                                           gCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGAc                                        >  5:133741/1‑149 (MQ=255)
                                                                                                           gCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGAc                                        >  4:294566/1‑149 (MQ=255)
                                                                                                              tCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGt                                     >  3:149663/1‑149 (MQ=255)
                                                                                                              tCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGt                                     >  4:30147/1‑149 (MQ=255)
                                                                                                                 aaaaCTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCa                                  <  4:149663/149‑1 (MQ=255)
                                                                                                                 aaaaCTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCa                                  >  7:43393/1‑149 (MQ=255)
                                                                                                                        ctCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAg                            <  4:254420/148‑1 (MQ=255)
                                                                                                                        ctCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAg                            <  4:258837/148‑1 (MQ=255)
                                                                                                                           gcgTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGgaga                        <  3:30147/149‑1 (MQ=255)
                                                                                                                                aCCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTg                   >  1:25125/1‑149 (MQ=255)
                                                                                                                                     gcggcggCGCCCCGGCCGGGCGCCGCCGCCCCCCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGCCGAATCACACGTCAGCGGTCCGCCCGTCGGGGCCATGTACAAGGGCGCCCGTTCACGCAAGGAGACCCTGGtgg               >  4:331037/6‑147 (MQ=255)
                                                                                                                                      tCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTc                >  7:300674/1‑147 (MQ=255)
                                                                                                                                      tCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGa              <  3:39523/149‑1 (MQ=255)
                                                                                                                                         gggCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTa           >  6:252998/1‑149 (MQ=255)
                                                                                                                                          ggCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTAc          <  6:174438/149‑1 (MQ=255)
                                                                                                                                          ggCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTAc          >  4:267894/1‑149 (MQ=255)
                                                                                                                                             gcCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACgg        >  1:95044/1‑149 (MQ=255)
                                                                                                                                             gcCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAg                            >  3:109366/1‑128 (MQ=255)
                                                                                                                                                cggccggGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTc    <  7:29711/149‑1 (MQ=255)
                                                                                                                                                 ggccggGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTcc   >  3:336359/1‑149 (MQ=255)
                                                                                                                                                  gccggGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTCCg  >  4:170149/1‑149 (MQ=255)
                                                                                                                                                  gccggGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTCCg  <  4:336359/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTCTGGAGTACCTGCACAAGGCCAACAAGCTGGTCGAAGCGCAGCGTCTGGAGCAACGCACCCGCTTCGATCTGGAGATGATCCTCGAGCTGGGCTACTGCAACGGCATCGAAAACTACTCGCGTTACCTGTCCGGGCGCCCGGCCGGGGCGCCGCCGCCCACCCTGTACGATTACCTGCCGCCCGATGCGCTGCTGGTGATCGACGAATCACACGTCAGCGTTCCGCAGGTCGGCGCCATGTACAAGGGCGACCGTTCACGCAAGGAGACCCTGGTCGAGTACGGCTTCCG  >  NC_002947/2235694‑2235986

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