Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,297,143 G→C 6.0% D993E (GAC→GAG) ‡ sbcC ← exonuclease subunit SbcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,297,1430GC 6.0% 271.8 / 6.4 83D993E (GAC→GAG) ‡sbcCexonuclease subunit SbcC
Reads supporting (aligned to +/- strand):  ref base G (36/42);  new base C (3/2);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

TTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAGGCGATC  >  NC_002947/2296996‑2297289
                                                                                                                                                   |                                                                                                                                                  
ttCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGgtgt                                                                                                                                                   >  4:139006/1‑149 (MQ=255)
       tCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACgg                                                                                                                                             <  2:303037/148‑1 (MQ=255)
         tGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCTGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCGACGGTCa                                                                                                                                          >  4:24775/1‑149 (MQ=255)
           ccgccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATc                                                                                                                                        <  3:95770/149‑1 (MQ=255)
           ccgccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATc                                                                                                                                        <  2:151449/149‑1 (MQ=255)
           ccgccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGACACCACGGTcc                                                                                                                                          <  5:207551/147‑2 (MQ=255)
           ccgccAGCGCCTGGTGGGCCTGCGGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGACACCACGGTccaa                                                                                                                                        <  8:303453/149‑4 (MQ=255)
             gccAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCtg                                                                                                                                      >  5:309607/1‑149 (MQ=255)
                   gcCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCcgc                                                                                                                                <  7:261230/149‑1 (MQ=255)
                   gcCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCcgc                                                                                                                                >  3:1134/1‑149 (MQ=255)
                   gcCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCcgc                                                                                                                                <  2:154965/149‑1 (MQ=255)
                     cTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTATGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGt                                                                                                                              <  1:104456/149‑1 (MQ=255)
                        gtggGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAt                                                                                                                           <  2:317000/149‑1 (MQ=255)
                          gggCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAtgc                                                                                                                         >  1:298429/1‑149 (MQ=255)
                           ggCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGAGTCCACGGTCATCTGTGCCGCGTGAtgc                                                                                                                         <  6:146680/148‑1 (MQ=255)
                            gCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAtgctg                                                                                                                       >  3:25614/1‑149 (MQ=255)
                             ccTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCATGGTCATCTGTGCCGCGTGAtgctgc                                                                                                                      <  4:321751/149‑1 (MQ=255)
                             ccTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAtgctg                                                                                                                       >  8:293799/1‑148 (MQ=255)
                             ccTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGAtgctg                                                                                                                       >  3:118712/1‑148 (MQ=255)
                                 ctgGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAg                                                                                                                    <  5:53330/147‑1 (MQ=255)
                                  tgGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGgcg                                                                                                                 >  4:25940/1‑149 (MQ=255)
                                       gcgcTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGtt                                                                                                            <  3:139006/149‑1 (MQ=255)
                                       gcgcTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGt                                                                                                             >  7:286187/1‑148 (MQ=255)
                                        cgcTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTc                                                                                                           <  4:1134/149‑1 (MQ=255)
                                             gtcgtcGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCgcgc                                                                                                       <  3:177191/148‑1 (MQ=255)
                                             gtcgtcGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCgcgc                                                                                                       <  5:159239/148‑1 (MQ=255)
                                                   ggCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTgc                                                                                                 <  3:25940/148‑1 (MQ=255)
                                                       tgttgtGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTcg                                                                                                          >  4:85705/1‑135 (MQ=255)
                                                       tgttgtGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTcg                                                                                                          <  3:85705/135‑1 (MQ=255)
                                                        gttgtGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGcagcag                                                                                            >  8:227416/1‑148 (MQ=255)
                                                         ttgtGCACGCAATTCGGCACATTGCTGCTCGTGGTTTGCCGGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGcagcag                                                                                            <  8:252037/147‑1 (MQ=255)
                                                            tgCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGAcgc                                                                                       >  6:43620/1‑149 (MQ=255)
                                                              cACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCgg                                                                                     <  1:175356/149‑1 (MQ=255)
                                                              cACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCgg                                                                                     <  3:91724/149‑1 (MQ=255)
                                                               aCGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGGCCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGc                                                                                    >  3:300155/1‑149 (MQ=255)
                                                               aCGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGc                                                                                    <  6:153215/149‑1 (MQ=255)
                                                                cGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGcc                                                                                   >  8:135078/1‑149 (MQ=255)
                                                                     ttCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCt                                                                              >  3:207150/1‑149 (MQ=255)
                                                                     ttCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCt                                                                              >  4:29742/1‑149 (MQ=255)
                                                                          cgcATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATg                                                                         <  7:293799/149‑1 (MQ=255)
                                                                          cgcATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCACCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATg                                                                         <  1:269932/149‑1 (MQ=255)
                                                                             aTTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCt                                                                      >  1:84737/1‑149 (MQ=255)
                                                                                 ctgctCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGATGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCttttt                                                                  >  8:204292/1‑149 (MQ=255)
                                                                                    ctCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTcc                                                                <  7:227416/148‑1 (MQ=255)
                                                                                     tCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCcgc                                                              >  8:4292/1‑149 (MQ=255)
                                                                                     tCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCcgc                                                              >  3:194241/1‑149 (MQ=255)
                                                                                        tGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGcc                                                           >  3:170111/1‑149 (MQ=255)
                                                                                          gTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATTCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCtt                                                         <  5:43620/149‑1 (MQ=255)
                                                                                               ccAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGAGCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCaaa                                                    >  1:231838/1‑149 (MQ=255)
                                                                                               ccAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTTATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATCGCTTTTTCCGCGCCTTGCaaa                                                    <  4:25614/149‑1 (MQ=255)
                                                                                                cAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAAc                                                   >  2:112731/1‑149 (MQ=255)
                                                                                                cAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAAc                                                   >  1:96338/1‑149 (MQ=255)
                                                                                                    cGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTg                                               <  8:286187/149‑1 (MQ=255)
                                                                                                     gTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTgg                                              <  4:170111/149‑1 (MQ=255)
                                                                                                     gTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTgg                                              >  8:134564/1‑149 (MQ=255)
                                                                                                      ttCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGgc                                             <  4:207150/149‑1 (MQ=255)
                                                                                                           gcAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAAt                                        >  4:301934/1‑149 (MQ=255)
                                                                                                           gcAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAAt                                        >  1:245897/1‑149 (MQ=255)
                                                                                                            cAGGTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAAtt                                       <  6:303198/149‑1 (MQ=255)
                                                                                                             aGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTc                                      <  2:84737/149‑1 (MQ=255)
                                                                                                              gTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCa                                     <  1:112731/149‑1 (MQ=255)
                                                                                                              gTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCa                                     <  5:320304/149‑1 (MQ=255)
                                                                                                              gTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCa                                     <  4:118712/149‑1 (MQ=255)
                                                                                                                 cGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCAtt                                   <  7:134564/148‑1 (MQ=255)
                                                                                                                 cGGCCAAGCCTTGCTCCAGGGCTTGGACCGGGGTGTCCACGGTCAGCTGTGCCGCGTGATGCTGCACGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCAtt                                   <  4:57591/148‑1 (MQ=255)
                                                                                                                  ggCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGa                                 <  5:329411/149‑1 (MQ=255)
                                                                                                                   gCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGAc                                >  8:126100/1‑149 (MQ=255)
                                                                                                                    ccAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACt                               >  2:151776/1‑149 (MQ=255)
                                                                                                                     cAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACtt                              >  7:251927/1‑149 (MQ=255)
                                                                                                                           ttGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGcgcg                         >  6:300551/1‑148 (MQ=255)
                                                                                                                              cTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCa                      <  2:96338/148‑1 (MQ=255)
                                                                                                                              cTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCAt                     <  2:298429/149‑1 (MQ=255)
                                                                                                                              cTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAAGTCAATGAATTgcgccagc                      <  4:194241/148‑5 (MQ=255)
                                                                                                                               tCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATc                    >  7:21020/1‑149 (MQ=255)
                                                                                                                               tCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATc                    >  4:40054/1‑149 (MQ=255)
                                                                                                                                 cAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCa                  >  6:13564/1‑149 (MQ=255)
                                                                                                                                 cAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCa                  >  2:165383/1‑149 (MQ=255)
                                                                                                                                  aGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCAt                 >  8:11813/1‑149 (MQ=255)
                                                                                                                                  aGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCAt                 <  3:40054/149‑1 (MQ=255)
                                                                                                                                   gggCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATg                <  2:231838/149‑1 (MQ=255)
                                                                                                                                       gatgACCGTGGACACCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGc              >  6:105123/4‑147 (MQ=255)
                                                                                                                                          gACCGTGGACACCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCcag           >  7:267185/1‑147 (MQ=255)
                                                                                                                                          gACCGTGGACACCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCcag           >  2:100818/1‑147 (MQ=255)
                                                                                                                                            ccGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAgg       <  1:134937/149‑1 (MQ=255)
                                                                                                                                                 gtgtCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAGGCGATc  >  4:227638/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCGATTTCTGCCGCCAGCGCCTGGTGGGCCTGCTGGCGGCGCTGGTCGTCGGCCTGTTGTGCACGCAATTCGGCGCATTGCTGCTCGTGGTTTGCCAGGCGTTCGCGCAGTTCGGCCAAGGCTTGCTCCAGGGCTTGGACCGTGGTGTCCACGGTCATCTGTGCCGCGTGATGCTGCAGGCGCTGTTCGCGCTCCTGCAGCAGGACGCGGCCCTGCTCGATGGCTTTTTCCGCGCCTTGCAAACGCTGGCGCAATTCATTGACTTGCGCGTCATCCATGGCCAGCAGGCGATC  >  NC_002947/2296996‑2297289

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: