Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,825,091 G→C 6.3% R55G (CGC→GGC)  mutS ← DNA mismatch repair protein MutS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,825,0910GC 6.3% 196.4 / 3.7 63R55G (CGC→GGC) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base C (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

CACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGC  >  NC_002947/1824945‑1825238
                                                                                                                                                  |                                                                                                                                                     
cACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgg                                                                                                                                                      <  1:362389/148‑1 (MQ=255)
 aCCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggc                                                                                                                                                    >  5:406102/1‑149 (MQ=255)
 aCCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggc                                                                                                                                                    >  4:243295/1‑149 (MQ=255)
    gggCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGt                                                                                                                                                 <  2:112952/149‑1 (MQ=255)
        ccTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCCCGCGGTCAgt                                                                                                                                             <  7:354019/149‑2 (MQ=255)
         cTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACGAGCTTGGCCAGGTAGCCTTCGAGCGAATGGAACGGAATCGCGGACATGGGGAGGGACGGGCCGGGCGCCGGACCGCGGGCGGGCAggg                                                                                                                                            >  4:299538/1‑149 (MQ=255)
         cTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAggg                                                                                                                                            <  3:243295/149‑1 (MQ=255)
         cTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAggg                                                                                                                                            <  6:406102/149‑1 (MQ=255)
               ggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTAGCCGGCCGACTGACCgcgggcg                                                                                                                                                   >  6:197575/1‑136 (MQ=255)
               ggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTAGCCGGCCGACTGACCgcgggcg                                                                                                                                                   <  5:203364/136‑1 (MQ=255)
               ggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTAGCCGGCCGACTGACCgcgggcg                                                                                                                                                   <  5:197575/136‑1 (MQ=255)
               ggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTAGCCGGCCGACTGACCgcgggcg                                                                                                                                                   >  6:203364/1‑136 (MQ=255)
                gtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATc                                                                                                                                     >  4:296058/1‑149 (MQ=255)
                   gCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCcag                                                                                                                                  <  4:101562/149‑1 (MQ=255)
                     ggggTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCcag                                                                                                                                  >  4:336707/1‑147 (MQ=255)
                       ggTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGCTATCACGCAGt                                                                                                                              >  4:210372/1‑149 (MQ=255)
                           gCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTc                                                                                                                           <  1:271142/148‑1 (MQ=255)
                            ccGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAATGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgc                                                                                                                         >  2:80548/1‑149 (MQ=255)
                             cGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgcg                                                                                                                        >  6:30014/1‑149 (MQ=255)
                                  tGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCtt                                                                                                                  >  3:189440/1‑149 (MQ=255)
                                   gCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCttt                                                                                                                 >  6:138823/1‑149 (MQ=255)
                                    cTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgcg                                                                                                                        <  4:327812/142‑1 (MQ=255)
                                       acaGATCACCACCGACTCGCCGAGCTTGACGAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                              <  3:396419/148‑1 (MQ=255)
                                       acaGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCttttttc                                                                                                             >  5:160194/1‑147 (MQ=255)
                                                caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCg                                                                                                    >  2:177718/1‑148 (MQ=255)
                                                 accGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                  <  4:124904/149‑1 (MQ=255)
                                                  ccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                              <  5:131574/137‑1 (MQ=255)
                                                  ccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                              >  6:131574/1‑137 (MQ=255)
                                                         gCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGAtct                                                                                          <  4:189440/149‑1 (MQ=255)
                                                         gCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCCCGCGGTCAgtc                                                                                                                                            >  8:80295/1‑99 (MQ=255)
                                                         gCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCCCGCGGTCAgtc                                                                                                                                            <  7:80295/101‑3 (MQ=255)
                                                          ccGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCTACGAAGAtgcg                                                                                         <  5:30014/149‑4 (MQ=255)
                                                          ccGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGAtctc                                                                                         >  4:229086/1‑149 (MQ=255)
                                                           cGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCg                                                                                        >  7:86903/1‑149 (MQ=255)
                                                             aGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                      <  3:336707/149‑1 (MQ=255)
                                                                  gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGaa                                                                                   <  5:138823/147‑1 (MQ=255)
                                                                    ccAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCg                                                                               <  6:45065/149‑1 (MQ=255)
                                                                     cAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGc                                                                              <  3:296058/149‑1 (MQ=255)
                                                                     cAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGc                                                                              >  2:341114/1‑149 (MQ=255)
                                                                       gCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGcc                                                                             >  8:148830/1‑148 (MQ=255)
                                                                            gCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGc                                                                        <  6:160194/149‑1 (MQ=255)
                                                                            gCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCg                                                                       >  5:210331/1‑149 (MQ=255)
                                                                                 gTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAg                                                                   >  6:93486/1‑148 (MQ=255)
                                                                                      cTTCCAGCGAATGGAATGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATc                                                             <  1:80548/149‑1 (MQ=255)
                                                                                       ttCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACAt                                                              >  1:185879/1‑147 (MQ=255)
                                                                                         ccAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCGATGGGAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACAt                                                              >  1:339639/1‑146 (MQ=255)
                                                                                                 tGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGTACATCAGCTGGTCCgg                                                  <  8:240157/149‑1 (MQ=255)
                                                                                                 tGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCgg                                                  <  6:210331/149‑1 (MQ=255)
                                                                                                   gAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCTACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCGTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGt                                                >  1:276547/1‑149 (MQ=255)
                                                                                                         aaTCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGt                                          <  7:148830/149‑1 (MQ=255)
                                                                                                             ccGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCGGGtttt                                       >  7:108066/1‑148 (MQ=255)
                                                                                                               catgtgcgggATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCa                                    <  4:284285/142‑1 (MQ=255)
                                                                                                                 acaTGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGc                                  <  1:177718/149‑1 (MQ=255)
                                                                                                                  caTGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGGTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCt                                 <  8:86903/149‑1 (MQ=255)
                                                                                                                    tGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTc                               <  3:229086/149‑1 (MQ=255)
                                                                                                                     ggggATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTc                               >  6:151640/1‑148 (MQ=255)
                                                                                                                      gggATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTcc                              <  5:93486/148‑1 (MQ=255)
                                                                                                                       gtaTGGACTGGCCGGCCGACTGACCGCTGGCGGTCAGCGTGATATCAAGCAGTTTCGCGGCCTTTTTCGAATATACGTAGAATATCTCGTAGAAGTCGACCATCAGGTAGAAAATCAGCTGGACCGGTTGCTGATTTTTCAGCTTCCAg                            >  6:410966/3‑149 (MQ=255)
                                                                                                                       ggATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAg                            <  3:299538/149‑1 (MQ=255)
                                                                                                                              cTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActgct                     <  3:210372/149‑1 (MQ=255)
                                                                                                                               tGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGGCCGGGTGCTGGTTTTTCAGCTTCCAGTActgctg                    >  6:307576/1‑149 (MQ=255)
                                                                                                                                    ggccgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCa                                    <  4:286455/128‑1 (MQ=255)
                                                                                                                                    ggccgACTGACAGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCa                                    >  3:286455/1‑128 (MQ=255)
                                                                                                                                         cggtaCCGCGGGCGGTCACTGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCggg          <  6:298600/145‑1 (MQ=255)
                                                                                                                                         ccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTTTTAGAAGATATCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTTCTTcataatagtt          >  6:286451/2‑145 (MQ=255)
                                                                                                                                          cTGACCGCGGGCGTTCTGGGTTATATGCAGCAGTTTCGCGGCCTTTTTCGCTTCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCGGGGTGCTg                                            <  4:126152/114‑1 (MQ=255)
                                                                                                                                          cTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTg                                            >  3:126152/1‑114 (MQ=255)
                                                                                                                                            gACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtg       >  2:322244/1‑149 (MQ=255)
                                                                                                                                             aCCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtgt      >  2:163358/1‑149 (MQ=255)
                                                                                                                                             aCCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtgt      >  3:354322/1‑149 (MQ=255)
                                                                                                                                                 cgggaggTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGc  <  5:410966/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGC  >  NC_002947/1824945‑1825238

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: