Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,301,752 T→C 5.1% Q674R (CAG→CGG) ‡ recB ← Chi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,301,7520TC 5.1% 235.0 / 2.0 78Q674R (CAG→CGG) ‡recBChi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base C (2/2);  total (40/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 = 9.66e-01

CAGTTCGCCGGCGGCCTGCTGCAGCAGCTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCG  >  NC_002947/5301607‑5301898
                                                                                                                                                 |                                                                                                                                                  
cAGTTCGCCGGCGGCCTGCTGCAGCAGCTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGcc                                                                                                                                                 <  5:153597/149‑1 (MQ=255)
  gTTCGCCGGCGGCCTGCTGCAGCAGCTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAg                                                                                                                                               <  3:78140/149‑1 (MQ=255)
          gcggcCTGCTGCAGCAGCTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTcgcg                                                                                                                                       <  5:181424/149‑1 (MQ=255)
             gcCTGCTGCAGCAGCTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTa                                                                                                                                    <  2:111827/149‑1 (MQ=255)
                  ctgcAGCAGCTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCAGCGCCAgggggt                                                                                                                                          <  6:395401/138‑5 (MQ=255)
                  ctgcAGCAGCTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCAGCGCCAgggggt                                                                                                                                          >  5:395401/1‑134 (MQ=255)
                  ctgcAGCAGCTCGGCCAGGGGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAACTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAgg                                                                                                                                <  4:28880/148‑1 (MQ=255)
                         agcTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACg                                                                                                                        >  8:192601/1‑149 (MQ=255)
                           cTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACa                                                                                                                                  >  4:8016/1‑137 (MQ=255)
                           cTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACa                                                                                                                                  <  3:8016/137‑1 (MQ=255)
                           cTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCa                                                                                                                      >  2:144495/1‑149 (MQ=255)
                           cTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCa                                                                                                                      <  4:268787/149‑1 (MQ=255)
                            tCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCAt                                                                                                                     >  8:94546/1‑149 (MQ=255)
                            gcGGCCAGGGGCAGCAGGTTGGGCAGCACCCGCTCGCCAGCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCTTGCTGCAGGTGGCGCAACAGCGGCAACACCCCCTGGCGCTGCTAGGTGTCGCGTTACAGGCCAAAACGCac                                                                                                                     <  4:211684/148‑2 (MQ=255)
                              ggCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCa                                                                                                                   >  8:373825/1‑149 (MQ=255)
                              ggCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCa                                                                                                                   <  2:184061/149‑1 (MQ=255)
                                   ggTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAg                                                                                                              <  7:192601/149‑1 (MQ=255)
                                       cagcagGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGctct                                                                                                          <  7:284784/149‑1 (MQ=255)
                                         gcagGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCAt                                                                                                                     >  2:195896/1‑136 (MQ=255)
                                         gcagGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCAt                                                                                                                     <  1:195896/136‑1 (MQ=255)
                                            ggTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTcc                                                                                                        >  4:247714/1‑146 (MQ=255)
                                             gTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAACCGCATCGCCCAGCTCTCCCAgc                                                                                                    >  8:137102/1‑148 (MQ=255)
                                                  tCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCAGCGCCAgggg                                                                                                                                            <  5:371334/104‑3 (MQ=255)
                                                  tCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCAGCGCCAgggg                                                                                                                                            >  6:371334/1‑102 (MQ=255)
                                                    aGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTcc                                                                                                        >  3:189280/1‑138 (MQ=255)
                                                    aGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTcc                                                                                                        <  4:189280/138‑1 (MQ=255)
                                                       aCCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCTTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTc                                                                                          <  8:57905/149‑1 (MQ=255)
                                                       aCCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTc                                                                                          <  3:247714/149‑1 (MQ=255)
                                                        cccGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTcc                                                                                         <  3:187292/149‑1 (MQ=255)
                                                         ccGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGc                                                                                               >  3:399129/1‑142 (MQ=255)
                                                         ccGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGc                                                                                               <  4:399129/142‑1 (MQ=255)
                                                           gCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTgg                                                                                      <  1:57610/149‑1 (MQ=255)
                                                               gCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTaa                                                                                  <  7:373825/149‑1 (MQ=255)
                                                                agaTCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAg                                                                                 <  1:144495/147‑1 (MQ=255)
                                                                      cTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCgg                                                                             <  7:243735/147‑1 (MQ=255)
                                                                      cTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCgg                                                                             >  8:243735/1‑147 (MQ=255)
                                                                       tGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTc                                                                           >  3:106816/1‑148 (MQ=255)
                                                                         cggcggATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAg                                                                        >  7:57730/1‑149 (MQ=255)
                                                                          gggggATCAGTGTGCGCGGTTGCTGGTAGTGATGCAGGAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCTCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAgc                                                                       <  8:57730/146‑1 (MQ=255)
                                                                          ggcgtATCTGTGTTCGCGGTAGCTGGAAGTGATTAAGCAGGTGGCGCAACATCGGCAACACCCCCTTGCGCTGCCATGTGTCGCGGTACAGTCGAAAACGCATCACCCAGCTCTCCCAGAACCTCTCGTCCTGGTTAAGCCGGTCCAg                                                                        <  6:267383/148‑1 (MQ=255)
                                                                          ggcggATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAgc                                                                       <  8:211556/149‑1 (MQ=255)
                                                                          ggcggATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAgc                                                                       <  1:7821/149‑1 (MQ=255)
                                                                           gcggATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCAAgcg                                                                      >  8:192013/1‑149 (MQ=255)
                                                                              gATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGcag                                                                   >  8:253102/1‑149 (MQ=255)
                                                                                   gtgtGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCACGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGcagccagc                                                              >  2:166782/1‑149 (MQ=255)
                                                                                     gtgCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGa                                                            >  3:110791/1‑149 (MQ=255)
                                                                                       gcgcGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGcagcca                                                                >  3:115322/1‑143 (MQ=255)
                                                                                       gcgcGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGcagcca                                                                <  4:115322/143‑1 (MQ=255)
                                                                                         gcgGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCACGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAgg                                                        <  1:166782/149‑1 (MQ=255)
                                                                                                cTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAgtg                                                 <  7:162699/149‑1 (MQ=255)
                                                                                                      gTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGc                                           >  6:144200/1‑149 (MQ=255)
                                                                                                       tCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGcagccagc                                                              <  1:195013/129‑1 (MQ=255)
                                                                                                       tCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGcagccagc                                                              >  2:195013/1‑129 (MQ=255)
                                                                                                         aTGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAg                                                                        <  1:92193/117‑1 (MQ=255)
                                                                                                         aTGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAg                                                                        >  2:92193/1‑117 (MQ=255)
                                                                                                         aTGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTg                                        >  4:356559/1‑149 (MQ=255)
                                                                                                          tGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGc                                       <  5:144200/149‑1 (MQ=255)
                                                                                                          tGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGc                                       >  4:186339/1‑149 (MQ=255)
                                                                                                          tGCAGCAGGGGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGc                                       <  7:94546/149‑1 (MQ=255)
                                                                                                              gcagGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCa                                     >  4:180794/1‑147 (MQ=255)
                                                                                                                       gcAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCgg                                                                             <  6:32617/98‑1 (MQ=255)
                                                                                                                       gcAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCgg                                                                             >  5:32617/1‑98 (MQ=255)
                                                                                                                       gcAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCtt                           >  6:166139/1‑148 (MQ=255)
                                                                                                                       gcAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCtt                           >  2:80737/1‑148 (MQ=255)
                                                                                                                            aTCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAg                     <  7:253102/149‑1 (MQ=255)
                                                                                                                             tCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCtt                           >  3:299791/1‑142 (MQ=255)
                                                                                                                             tCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCtt                           <  4:299791/142‑1 (MQ=255)
                                                                                                                                gCAACACCCCCTGGCGCTGCCAGGAGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCACACCCGCTCGTACTGGTTAAGCAGGTCCAGCGCAGCCCGCGACAGGCCCAGTGTAAGGCTTTCAAGAGCAGCCTTGAGCAcgcg                  >  2:190723/1‑148 (MQ=255)
                                                                                                                                  aaCACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTc               >  6:384149/1‑149 (MQ=255)
                                                                                                                                      cccccTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCCGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAg           >  8:147157/1‑149 (MQ=255)
                                                                                                                                       ccccTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCCGCCCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGCCCGGCCCAGTGTAAGGCTTGCCCGAGCAGCCTTGAGCAGGCGCTCCGAGt          >  5:412337/1‑149 (MQ=255)
                                                                                                                                       ccccTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTc                                                                           >  6:298324/1‑84 (MQ=255)
                                                                                                                                       ccccTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTc                                                                           <  5:298324/84‑1 (MQ=255)
                                                                                                                                       ccccTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGt          >  7:71187/1‑149 (MQ=255)
                                                                                                                                       ccccTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGt          >  2:157437/1‑149 (MQ=255)
                                                                                                                                       ccccTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGt          >  4:252568/1‑149 (MQ=255)
                                                                                                                                             gcgcTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTCCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCt    >  8:365347/1‑149 (MQ=255)
                                                                                                                                               gcTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCg  >  8:130110/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAGTTCGCCGGCGGCCTGCTGCAGCAGCTCGGCCAGGTGCAGCAGGTTGGTCAGCACCCGCTCGCCATCACTGCGGCGGATCAGTGTGCGCGGTAGCTGGAAGTCATGCAGCAGGTGGCGCAACATCGGCAACACCCCCTGGCGCTGCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCG  >  NC_002947/5301607‑5301898

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: