Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,297,648 T→A 6.1% D825V (GAC→GTC)  sbcC ← exonuclease subunit SbcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,297,6480TA 6.1% 192.0 / 2.7 66D825V (GAC→GTC) sbcCexonuclease subunit SbcC
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base A (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGG  >  NC_002947/2297503‑2297790
                                                                                                                                                 |                                                                                                                                              
ctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCa                                                                                                                                             <  3:317139/149‑1 (MQ=255)
ctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCa                                                                                                                                             <  2:331444/149‑1 (MQ=255)
 tgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTTGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGGTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAg                                                                                                                                            >  7:131019/1‑149 (MQ=255)
   ccAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGACCAGCt                                                                                                                                          >  5:312886/1‑149 (MQ=255)
   ccAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGACCAGCt                                                                                                                                          >  5:317881/1‑149 (MQ=255)
             gcgcTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCg                                                                                                                                <  7:218099/149‑1 (MQ=255)
                 tGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCt                                                                                                                            <  7:137274/149‑1 (MQ=255)
                   gCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtg                                                                                                                          >  4:172099/1‑149 (MQ=255)
                         tCGCCGAACTGCTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGctc                                                                                                                    >  8:115605/1‑149 (MQ=255)
                            ccGAGCAGTTCGGCGAGTGAGGCAAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTc                                                                                                                 <  7:192537/149‑1 (MQ=255)
                              gtgCAGTTCGGCGAGTGAGGCCAGGTCTTGCTGGCGTTGCTCGTAGAGGGAGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGGACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtgccg                                                                                                                       <  2:291647/139‑3 (MQ=255)
                              gaactgctCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCa                                                                                                               <  7:115605/142‑1 (MQ=255)
                              gAGCAGTTCGGAGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGAGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtgcagcagatgct                                                                                                               >  1:291647/1‑142 (MQ=255)
                                  aGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGCCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTc                                                                                                           <  1:121389/149‑1 (MQ=255)
                                      cGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTc                                                                                                        <  2:45832/148‑1 (MQ=255)
                                       ggCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTcc                                                                                                       <  3:172099/148‑1 (MQ=255)
                                        gCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCtt                                                                                                     >  1:147696/1‑149 (MQ=255)
                                         cGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTg                                                                                                    <  6:162168/149‑1 (MQ=255)
                                                ggCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGtt                                                                                             >  6:79016/1‑149 (MQ=255)
                                                   cAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCa                                                                                          <  4:181519/149‑1 (MQ=255)
                                                   cAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCa                                                                                          <  2:147696/149‑1 (MQ=255)
                                                            cTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGACCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTtgctgc                                                                                  <  6:317881/148‑1 (MQ=255)
                                                            cTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGACCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTtgctgc                                                                                  <  6:312886/148‑1 (MQ=255)
                                                             tGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAg                                                                                <  3:39331/149‑1 (MQ=255)
                                                                 gTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGgcgc                                                                            >  8:253365/1‑149 (MQ=255)
                                                                    gCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTg                                                                          <  7:106190/148‑1 (MQ=255)
                                                                    gCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGATCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTg                                                                          <  5:79016/148‑1 (MQ=255)
                                                                      tcgtcgAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTc                                                                                            <  5:293060/128‑1 (MQ=255)
                                                                      tcgtcgAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTc                                                                                            >  6:293060/1‑128 (MQ=255)
                                                                      tcgtcgAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTggg                                                                        >  6:95006/1‑148 (MQ=255)
                                                                             gggCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGAGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCAAGTTGCTGCAGGCGATGGGCAATctgc                                                                >  7:73224/1‑149 (MQ=255)
                                                                              ggcggcCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATctgct                                                               <  5:12348/149‑1 (MQ=255)
                                                                                  gcCAGTTTCTGCCGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTc                                                           >  6:72012/1‑149 (MQ=255)
                                                                                   ccAGTTTCTGCTATAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTc                                                           >  8:63955/1‑148 (MQ=255)
                                                                                    cAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTc                                                           >  7:237876/1‑147 (MQ=255)
                                                                                    cAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGa                                                         <  8:237876/149‑1 (MQ=255)
                                                                                      gtgtCTGCTGTACCTGCTGCTGGCCTTGCACGCGCGCCTGCTGGTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAgg                                                       <  8:131019/146‑1 (MQ=255)
                                                                                           tgctgTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAg                                                  >  4:323093/1‑149 (MQ=255)
                                                                                             ctttAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGa                                                   <  7:63955/143‑1 (MQ=255)
                                                                                             ctgTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGa                                                   <  1:295208/146‑1 (MQ=255)
                                                                                                tAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAg                                              <  6:141691/148‑1 (MQ=255)
                                                                                                    tgctgctgGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCAtt                                         >  6:132143/1‑149 (MQ=255)
                                                                                                      ctgctgGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTg                                        <  2:290095/148‑1 (MQ=255)
                                                                                                         ctgGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGc                                      <  7:253365/147‑1 (MQ=255)
                                                                                                         ctgGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCg                                    >  3:293367/1‑149 (MQ=255)
                                                                                                             ccTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGgtc                                >  5:104378/1‑149 (MQ=255)
                                                                                                               ttGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCt                                                                                                      <  5:281352/77‑1 (MQ=255)
                                                                                                               ttGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCt                                                                                                      >  6:281352/1‑77 (MQ=255)
                                                                                                                 gCACGCGCGCCTGCTGTTGGTGACGCAGCTTGTCCAGCTTGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTTGTCGAGGGCGAGGAAGGCATTGGCCGGgtcgt                              <  8:73224/147‑1 (MQ=255)
                                                                                                                         gcCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGc                    >  6:164351/1‑149 (MQ=255)
                                                                                                                          ccTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGcc                   >  1:196337/1‑149 (MQ=255)
                                                                                                                               tgTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAg              >  4:217501/1‑149 (MQ=255)
                                                                                                                               tgTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAg              >  3:269389/1‑149 (MQ=255)
                                                                                                                               tgTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAg              >  6:251325/1‑149 (MQ=255)
                                                                                                                                gTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAgc             >  4:142437/1‑149 (MQ=255)
                                                                                                                                 ttGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGgtc                                <  2:125437/129‑1 (MQ=255)
                                                                                                                                 ttGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGgtc                                >  1:125437/1‑129 (MQ=255)
                                                                                                                                    gTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCtt         >  3:150065/1‑149 (MQ=255)
                                                                                                                                     tCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCtt                 >  4:6993/1‑140 (MQ=255)
                                                                                                                                     tCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCtt                 <  3:6993/140‑1 (MQ=255)
                                                                                                                                       aCGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTcc       <  6:104378/148‑1 (MQ=255)
                                                                                                                                         gCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCctct    >  4:251656/1‑149 (MQ=255)
                                                                                                                                         gCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCctct    >  1:84223/1‑149 (MQ=255)
                                                                                                                                           aGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTgg  <  4:269389/149‑1 (MQ=255)
                                                                                                                                           aGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTgg  <  5:132143/149‑1 (MQ=255)
                                                                                                                                           aGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTgg  <  5:72012/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTCCTTGCGCTGTTCCAGTTGCTGCAGGCGCTGGGCAATCTGCTGGTCGAGGGCGAGGAAGGCATTGGCCGGGTCGTCGTTGAGGGCCTTGAGCGCTTCCTCTGG  >  NC_002947/2297503‑2297790

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: ‑