Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,825,091 G→C 7.8% 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 7.8% 192.3 / 6.4 64R55G (CGC→GGC) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base C (3/2);  total (32/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 = 3.72e-01

CACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTG  >  NC_002947/1824945‑1825237
                                                                                                                                                  |                                                                                                                                                    
cACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcggg                                                                                                                                                    <  2:228359/149‑1 (MQ=255)
 aCCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggc                                                                                                                                                   >  2:232160/1‑149 (MQ=255)
     tgCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGTTCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGTTGGACTGGCCGGCCGACTGACCGCTGGCGGTc                                                                                                                                               <  1:227246/148‑1 (MQ=255)
        ccTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAgg                                                                                                                                            <  1:232160/149‑1 (MQ=255)
        ccTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAgg                                                                                                                                            >  6:14239/1‑149 (MQ=255)
        ccTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAgg                                                                                                                                            <  5:14239/149‑1 (MQ=255)
             cTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGat                                                                                                                                       <  1:363635/149‑1 (MQ=255)
              tggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGCGTGata                                                                                                                                      >  8:104097/1‑149 (MQ=255)
                    cGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCcagc                                                                                                                                <  4:155545/149‑1 (MQ=255)
                           gCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATcc                                                                                                                                   >  3:173304/1‑139 (MQ=255)
                           gCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATcc                                                                                                                                   <  4:173304/139‑1 (MQ=255)
                           gCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcg                                                                                                                         >  4:262510/1‑149 (MQ=255)
                                   gCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGACCGGAATCCCGCACATGGGGATGGACTGGACGGCCGCCTGACCGCGGGCGGTGAGGGGGATTTGCAGCAGTTTGGCGGCCttt                                                                                                                 >  7:222174/1‑149 (MQ=255)
                                    cTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCtttt                                                                                                                <  2:161023/149‑1 (MQ=255)
                                       acaGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                              <  3:152989/148‑1 (MQ=255)
                                         aGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCa                                                                                                           >  2:87959/1‑149 (MQ=255)
                                           aTCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATGCAGCAGTTTCGCGGCCTTTTTCGCATc                                                                                                         >  6:17563/1‑149 (MQ=255)
                                                caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCg                                                                                                    >  5:168427/1‑148 (MQ=255)
                                                 accGACTCGCCGAGTATGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCg                                                                                                    >  2:170257/1‑147 (MQ=255)
                                                       tCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGAt                                                                                            >  5:360958/1‑149 (MQ=255)
                                                         gCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGAtct                                                                                          >  2:73337/1‑149 (MQ=255)
                                                           cGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCg                                                                                        <  3:273939/149‑1 (MQ=255)
                                                             aGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTa                                                                                      <  6:168427/149‑1 (MQ=255)
                                                                  gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgcg                                                                                                                       >  2:134773/1‑112 (MQ=255)
                                                                  gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgcg                                                                                                                       <  1:134773/112‑1 (MQ=255)
                                                                  gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTAAGGGTGATATCCAGCAGTTTAGCGGCCTTTTTCGCATATTCGTAGAAGATATCGTAGAAg                                                                                  >  6:232244/1‑148 (MQ=255)
                                                                    ccAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTTATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTc                                                                                <  5:232244/148‑1 (MQ=255)
                                                                     cAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGc                                                                              >  8:273749/1‑149 (MQ=255)
                                                                            gCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGc                                                                        <  7:104097/148‑1 (MQ=255)
                                                                             ccAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCgg                                                                      <  1:170257/149‑1 (MQ=255)
                                                                               aGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTa                                                                    <  6:360958/149‑1 (MQ=255)
                                                                                    gCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAAc                                                                <  4:315123/148‑1 (MQ=255)
                                                                                    gCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACa                                                               >  6:87962/1‑149 (MQ=255)
                                                                                       ttCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACAt                                                              >  2:305656/1‑147 (MQ=255)
                                                                                         ccAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAg                                                          <  7:273749/148‑1 (MQ=255)
                                                                                          cAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTTCGTAGACGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGc                                                         <  1:73337/149‑1 (MQ=255)
                                                                                           aGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGGAGTCGCCCATGCGGTAGAACATCAGCTg                                                       >  6:164547/1‑149 (MQ=255)
                                                                                                 tGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCATCTACGAAGAtgcg                                                                                         <  7:85062/110‑4 (MQ=255)
                                                                                                 tGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCATCTACGAAGAtgcg                                                                                         >  8:85062/1‑107 (MQ=255)
                                                                                                   gAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTTATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGt                                               <  5:339161/149‑1 (MQ=255)
                                                                                                   gAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGt                                               >  6:302984/1‑149 (MQ=255)
                                                                                                    aaCGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCt                                                        <  7:164992/139‑1 (MQ=255)
                                                                                                    aaCGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTg                                              >  5:147183/1‑149 (MQ=255)
                                                                                                     aCGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGc                                             >  1:104773/1‑149 (MQ=255)
                                                                                                             ccGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGCTATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGttttt                                     >  2:199228/1‑149 (MQ=255)
                                                                                                              cGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTc                                    <  1:87959/149‑1 (MQ=255)
                                                                                                                    tGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGGCAGGGTGATATCCCGCAGTTTCGCGGCCTTTTTAGCATCTTCGTTGAAGATATCGTAGAAGTAGCCCATGCGGGAGAACATCAGCTGGTCCGGGGGGTGGGTTTTGAGCtga                              >  6:224755/1‑147 (MQ=255)
                                                                                                                      gggATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGGAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCa                            <  5:164547/149‑1 (MQ=255)
                                                                                                                      gggATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCa                            <  3:262510/149‑1 (MQ=255)
                                                                                                                         aTGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTa                         <  5:87962/149‑1 (MQ=255)
                                                                                                                         aTGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTa                         <  5:224755/149‑1 (MQ=255)
                                                                                                                                  ccggccgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCgg                                                                      <  8:293195/96‑1 (MQ=255)
                                                                                                                                  ccggccgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCgg                                                                      >  7:293195/1‑96 (MQ=255)
                                                                                                                                   cggccgACTGACGGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCGTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGcatc               <  5:17563/149‑1 (MQ=255)
                                                                                                                                      cacccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCatcatc            >  8:206869/5‑149 (MQ=255)
                                                                                                                                       cgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCg           >  7:344262/1‑149 (MQ=255)
                                                                                                                                       acccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGc                                 <  5:26255/124‑1 (MQ=255)
                                                                                                                                       acccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGc                                 >  6:26255/4‑127 (MQ=255)
                                                                                                                                        gACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTTGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCgg          >  2:259990/1‑149 (MQ=255)
                                                                                                                                        gACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCgg          >  1:359573/1‑149 (MQ=255)
                                                                                                                                           tGACCGCGGGCGGTCATGGTGATATACATGATTTTCGCGGCCTTTTTCGGATCTTAGTATAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCCAGCTGGTCCGGGTGCTGGTTTTTCAGCtt                               <  1:27255/126‑1 (MQ=255)
                                                                                                                                           tGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCtt                               >  2:27255/1‑125 (MQ=255)
                                                                                                                                           tGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAg                                                                   >  1:64862/1‑90 (MQ=255)
                                                                                                                                           tGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAg                                                                   <  2:64862/90‑1 (MQ=255)
                                                                                                                                             aCCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCa                            <  7:161015/126‑1 (MQ=255)
                                                                                                                                             aCCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCa                            >  8:161015/1‑126 (MQ=255)
                                                                                                                                               cgcgGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCt                                                        <  2:120297/96‑1 (MQ=255)
                                                                                                                                                gcgggcggTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGtgtgtg  >  3:141949/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTG  >  NC_002947/1824945‑1825237

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: