Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,288,074 G→C 9.1% A13A (GCC→GCG xthA ← exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,288,0740GC 9.1% 229.3 / 12.2 77A13A (GCC→GCGxthAexodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base C (3/4);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.76e-01

GGCCTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTGGGGGTGAACGAA  >  NC_002947/3287929‑3288219
                                                                                                                                                  |                                                                                                                                                 
ggCCTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATgcgggcgg                                                                                                                                                   <  8:200737/146‑1 (MQ=255)
ggCCTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATgcgggcgg                                                                                                                                                   >  7:200737/1‑146 (MQ=255)
ggCCTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGcc                                                                                                                                                >  3:200852/1‑149 (MQ=255)
 gCCTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGccc                                                                                                                                               >  2:98947/1‑149 (MQ=255)
   cTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCg                                                                                                                                              >  5:265686/1‑148 (MQ=255)
      tCTGGCCATGGTAGTGCACGTGGTAACGCTTTGCCTGCACGTCGGCCATCGGGTACTGATCGTCGGTGACCTTGGTTTCCTGCAGGGCGTTCACGTGCGGCTGGTTCTTTTCGATCAGCGCCGCCATCTGATGCGGGCGGGCCCGCAg                                                                                                                                           <  8:181536/148‑1 (MQ=255)
      tCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAgg                                                                                                                                          >  3:157965/1‑149 (MQ=255)
      tCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAgg                                                                                                                                          >  3:152397/1‑149 (MQ=255)
       cTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAgg                                                                                                                                          <  6:128247/148‑1 (MQ=255)
          gCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTTGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGt                                                                                                                                      <  8:149345/149‑1 (MQ=255)
          gCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCCCGCCCGCatca                                                                                                                                         <  3:178958/146‑4 (MQ=255)
           ccATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGGGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGtt                                                                                                                                     <  5:110400/149‑1 (MQ=255)
            cATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTg                                                                                                                                    <  3:105379/149‑1 (MQ=255)
                gTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGt                                                                                                                                <  2:235539/149‑1 (MQ=255)
                      cACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGaaa                                                                                                                           >  7:96687/1‑148 (MQ=255)
                      cACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAAc                                                                                                                          <  4:200852/149‑1 (MQ=255)
                       aCGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGaaa                                                                                                                           >  1:45651/1‑148 (MQ=255)
                         gTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGa                                                                                                                       >  5:191394/1‑149 (MQ=255)
                                 ccAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGa                                                                                                                       >  6:265673/1‑141 (MQ=255)
                                 ccAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGa                                                                                                                       <  5:265673/141‑1 (MQ=255)
                                 ccAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCtt                                                                                                                <  2:28065/148‑1 (MQ=255)
                                 ccAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCtt                                                                                                                >  1:28065/1‑148 (MQ=255)
                                     tGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATg                                                                                                            >  7:123670/1‑148 (MQ=255)
                                       ccTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGa                                                                                                          >  5:263372/1‑148 (MQ=255)
                                          gCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGa                                                                                                      >  3:32984/1‑149 (MQ=255)
                                              gTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAAt                                                                                                  <  4:216429/149‑1 (MQ=255)
                                              gTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAAt                                                                                                  <  4:157965/149‑1 (MQ=255)
                                              gTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAAt                                                                                                  <  4:152397/149‑1 (MQ=255)
                                                   ccAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCgg                                                                                             <  5:289556/149‑1 (MQ=255)
                                                   ccAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCCCGCCCGCatca                                                                                                                                         <  4:154582/105‑4 (MQ=255)
                                                   ccAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCCCGCCCGCatca                                                                                                                                         <  4:149014/105‑4 (MQ=255)
                                                   ccAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCCCGCCCGCatca                                                                                                                                         >  3:154582/1‑102 (MQ=255)
                                                   ccAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCCCGCCCGCatca                                                                                                                                         >  3:149014/1‑102 (MQ=255)
                                                    cAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCgg                                                                                             >  3:218175/1‑148 (MQ=255)
                                                    cAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGt                                                                                            <  8:123670/149‑1 (MQ=255)
                                                       cGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTaa                                                                                          >  5:159608/1‑148 (MQ=255)
                                                               tGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGgcgc                                                                                   >  6:110199/1‑148 (MQ=255)
                                                                 aTCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGAt                                                                                <  6:191394/149‑1 (MQ=255)
                                                                 aTCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCATTTCGATCAGCGACGCCAGCTGATGAGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCCCTAAAAGGCGCGAt                                                                                >  5:94050/1‑149 (MQ=255)
                                                                  tcgtcgCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATg                                                                               >  7:98473/1‑149 (MQ=255)
                                                                    gtcgCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGCTGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGc                                                                              <  1:98947/148‑1 (MQ=255)
                                                                            ccTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCTAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACa                                                                     >  1:229983/1‑149 (MQ=255)
                                                                            ccTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACa                                                                     <  8:96687/149‑1 (MQ=255)
                                                                             cTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTAGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGCCAt                                                                    >  5:102750/1‑149 (MQ=255)
                                                                                  tttCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACATACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGAc                                                                <  2:313614/148‑1 (MQ=255)
                                                                                  tttCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGAc                                                                >  1:313614/1‑148 (MQ=255)
                                                                                    tCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCt                                                              <  2:45651/148‑1 (MQ=255)
                                                                                      cTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCa                                                           >  3:271044/1‑149 (MQ=255)
                                                                                       tGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCa                                                           >  4:55092/1‑148 (MQ=255)
                                                                                       tGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCAt                                                          >  2:329206/1‑149 (MQ=255)
                                                                                        gaaGGACGATCACGTCCGCCTTGAGATTATCGATCAGCGCCGCCACCTCATCCGCCCCGGCCCGCAGGCCGTTGATGATTAAAAAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGCCGCGATCCTAGACGACATCGACCTTCAt                                                          >  2:288387/3‑148 (MQ=255)
                                                                                          aGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCAc                                                       >  3:275394/1‑149 (MQ=255)
                                                                                              cGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGa                                                                                                                    <  2:30470/84‑1 (MQ=255)
                                                                                              cGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCAGCAGGCCGTTGATGTTGAAACAGACGa                                                                                                                    >  1:30470/1‑84 (MQ=255)
                                                                                               gATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCa                                                  >  6:262021/1‑149 (MQ=255)
                                                                                               gATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAAGCCCGGTAAACGGCGCGATGCTAGCCGACAGCGACGTGCATCCCCAGCa                                                  >  8:158195/1‑149 (MQ=255)
                                                                                                 tCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATg                                                <  6:265686/149‑1 (MQ=255)
                                                                                                    cGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATggcg                                             >  3:287877/1‑149 (MQ=255)
                                                                                                    cGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATggcg                                             <  4:271044/149‑1 (MQ=255)
                                                                                                               tGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGc                                  >  5:24085/1‑149 (MQ=255)
                                                                                                                  ttttCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATgg                                               >  4:200295/1‑133 (MQ=255)
                                                                                                                  ttttCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATgg                                               <  3:200295/133‑1 (MQ=255)
                                                                                                                       gATCATCGCCGCCTGGTGATGCGGGCGGGCCCGCAGGCGGTTGATGTTGAATCAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTgctgct                          <  6:94050/149‑1 (MQ=255)
                                                                                                                              gccgccAGCTGATGCGGTCGTGCCCGCAGGCCGTTGATGTTGAAACAGACTATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGc                   <  1:288387/149‑1 (MQ=255)
                                                                                                                                     gCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTggg            <  3:55092/149‑1 (MQ=255)
                                                                                                                                      cTGATGCGGGCGGGGCCGCTGGCCGTAGATGTTTAACCAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTgggg           <  6:102750/149‑1 (MQ=255)
                                                                                                                                      cTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGcgcg                                                                                  >  1:282887/1‑78 (MQ=255)
                                                                                                                                      cTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGcgcg                                                                                  <  2:282887/78‑1 (MQ=255)
                                                                                                                                      cTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTgggg           <  8:98473/149‑1 (MQ=255)
                                                                                                                                      cTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTgggg           <  2:229983/149‑1 (MQ=255)
                                                                                                                                      cTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTgggg           >  6:191196/1‑149 (MQ=255)
                                                                                                                                       tGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAg                                        <  3:257580/119‑1 (MQ=255)
                                                                                                                                       tGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAg                                        >  4:257580/1‑119 (MQ=255)
                                                                                                                                        gATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGt                     <  6:55760/137‑1 (MQ=255)
                                                                                                                                        gATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGt                     >  5:55760/1‑137 (MQ=255)
                                                                                                                                        cctGCGGGCCCGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTggggg          >  2:176388/3‑148 (MQ=255)
                                                                                                                                        cctGCGGGCCCGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTggggg          <  1:176388/146‑1 (MQ=255)
                                                                                                                                         aTGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTggggg          <  1:329206/147‑1 (MQ=255)
                                                                                                                                         aTGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTggggg          <  4:32984/147‑1 (MQ=255)
                                                                                                                                               gcgggcCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTGGGGGTGAACGaa  <  5:110199/149‑1 (MQ=255)
                                                                                                                                                 gggcCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATcc                                                                                                >  6:370038/1‑53 (MQ=255)
                                                                                                                                                 gggcCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATcc                                                                                                <  5:370038/53‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGCCTTTCTGGCCATGGTAGTGCACGTGGTAACCCAGTGCCTGCACGTCGGCCAGCGGGAACTGATCGTCGCTGACCTTGGTTTCCTGCAGGCCGATCACGTCCGGCTGGTGCTTTTCGATCAGCGCCGCCAGCTGATGCGGGCGGGCCCGCAGGCCGTTGATGTTGAAACAGACGATCTTCATGGAAAGACAATCCCGGTAAAAGGCGCGATGCTAGCCGACATCGACCTGCATCACCAGCATGGCGGCCAGGCGAGCCTGCTGCTAACGTGCTGTTGGGGGTGAACGAA  >  NC_002947/3287929‑3288219

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: