Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,825,091 G→C 7.5% 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.5% 201.5 / 2.8 70R55G (CGC→GGC) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base C (4/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGT  >  NC_002947/1824945‑1825236
                                                                                                                                                   |                                                                                                                                                 
cACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgg                                                                                                                                                  <  5:182008/148‑1 (MQ=255)
cACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACAGACTCGCCGAGCTTGACCAGCTAGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcggg                                                                                                                                                 >  7:297532/1‑149 (MQ=255)
   cGGGCCCTTGCTGGTGGAGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgcgcgg                                                                                                                                              >  5:127478/1‑149 (MQ=255)
      gCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCa                                                                                                                                           <  6:345424/149‑1 (MQ=255)
       cacTTGCTGGTCGCGGGGTAGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCCGCCGGTCAg                                                                                                                                          >  6:399502/3‑149 (MQ=255)
         cTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAggg                                                                                                                                        >  6:380151/1‑149 (MQ=255)
         cTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCCCGCGGTCAgtc                                                                                                                                        <  2:192704/149‑3 (MQ=255)
          ttGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACAAGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGt                                                                                                                                       >  7:85870/1‑149 (MQ=255)
           tGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTg                                                                                                                                      >  5:205602/1‑149 (MQ=255)
              tggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCAAATGGGGATGGACTGGCCGGCCGACTGCACGCCGGCGGTCACGGTGat                                                                                                                                    >  8:343835/1‑148 (MQ=255)
                gtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATc                                                                                                                                 >  8:173525/1‑149 (MQ=255)
                gtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATc                                                                                                                                 <  2:205537/149‑1 (MQ=255)
                     ggggTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGgtataccg                                                                                                                              >  2:332670/1‑143 (MQ=255)
                            ccGATCTGCTCACAGATCACCACCGACTCGCCGATCTTGACCAGCTTGGCCAGGTATCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTTATATCCAGCAGTTTcgc                                                                                                                     <  8:297532/149‑1 (MQ=255)
                                   gCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCttt                                                                                                              >  3:322413/1‑149 (MQ=255)
                                    cTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgcg                                                                                                                    <  6:205602/142‑1 (MQ=255)
                                    cTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTcgcg                                                                                                                    <  7:337465/142‑1 (MQ=255)
                                    cTCACAGATCACCACCGACTAGACGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGTGGATGGACTGGCCGGCCGACTGACCGCTGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCttt                                                                                                              <  5:366580/148‑1 (MQ=255)
                                     tCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCttttt                                                                                                            >  5:293934/1‑149 (MQ=255)
                                     tCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCtttt                                                                                                             >  4:252665/1‑148 (MQ=255)
                                      agcaGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                           <  6:182231/147‑1 (MQ=255)
                                       acaGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc                                                                                                           <  2:79900/148‑1 (MQ=255)
                                       acaGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCg                                                                                                          >  5:91032/1‑149 (MQ=255)
                                       acaGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGCGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCg                                                                                                          >  3:333198/1‑149 (MQ=255)
                                              accaccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTc                                                                                                   <  1:286439/149‑1 (MQ=255)
                                               ccaccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATTTCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCg                                                                                                  >  6:159088/1‑149 (MQ=255)
                                                caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCg                                                                                                  >  1:57956/1‑149 (MQ=255)
                                                caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGt                                                                                                 >  1:317053/1‑149 (MQ=255)
                                                 accGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa                                                                                                <  5:159088/149‑1 (MQ=255)
                                                    gACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTagaa                                                                                             <  3:107687/149‑1 (MQ=255)
                                                           cGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCTACGAAGAtgcga                                                                                      <  4:333198/149‑5 (MQ=255)
                                                             aGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATc                                                                                                      >  6:93828/1‑131 (MQ=255)
                                                             aGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATc                                                                                                      <  5:93828/131‑1 (MQ=255)
                                                               cTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCTGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAg                                                                                   >  7:365525/1‑148 (MQ=255)
                                                                  gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGaa                                                                                 <  7:173525/147‑1 (MQ=255)
                                                                  gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGt                                                                               >  2:40500/1‑149 (MQ=255)
                                                                  gACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGt                                                                               >  6:99568/1‑149 (MQ=255)
                                                                  gACCAGCGTGGCCAGGTAGCCTTCCAGCGAGTGGATCGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGt                                                                               <  8:305464/148‑1 (MQ=255)
                                                                     caccTTGCCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGCCTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGc                                                                            <  2:284041/146‑1 (MQ=255)
                                                                     cAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGc                                                                            <  4:384197/149‑1 (MQ=255)
                                                                     cAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGc                                                                            <  1:40500/149‑1 (MQ=255)
                                                                           ggCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATg                                                                       <  5:380151/148‑1 (MQ=255)
                                                                             ccAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTTACCGCGGGCGGTCAGGGTGATATCAAGCAGTTTACCGGCCATTTTCGCATAATATTAGAACCTCTCGTAGAAGTCGCCCAttccg                                                                    >  8:316259/1‑145 (MQ=255)
                                                                             caaGGTAGCCTTCCAGCGAATGTAACGGTATCCCGCACATGGGGATGTACTGGCCGGGCGACTGACAGCGGGCTGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCgg                                                                    <  5:399502/147‑1 (MQ=255)
                                                                                     gCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAAc                                                              <  3:252665/148‑1 (MQ=255)
                                                                                        ttCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGCCCGCGGGCGGGCAGGGTGATCTCCAGCAGTTTCGCGGCCTGTTTCGCATCTTCGTAGACGATCGCGTAGAAGTCGCCCATGCGGTAGCACAt                                                            >  2:48789/1‑147 (MQ=255)
                                                                                          ccAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCa                                                          <  2:317053/147‑1 (MQ=255)
                                                                                             gCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTgg                                                     >  7:112171/1‑149 (MQ=255)
                                                                                             gCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTgg                                                     <  7:316259/149‑1 (MQ=255)
                                                                                                  tGGAACGTAATCCAGCACATGGGGATGGACTGGCCGGCCGACTGACCGCTGGCGGTCAGGGTTATATCCAGCAGGTTCGCGGCCTTTTTCGCATCTTCGTAGAATATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCgg                                                <  7:343835/149‑1 (MQ=255)
                                                                                                    gcaCGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGAAATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCggg                                               <  6:206512/146‑1 (MQ=255)
                                                                                                    gAACGGAATCCCGCACATGCGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCCTGCGGTAGAACATCAGCTGGTCCGGGt                                              <  1:212686/149‑1 (MQ=255)
                                                                                                     aaCGGAATCCTGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCt                                                                                                     >  6:261436/1‑93 (MQ=255)
                                                                                                     aaCGGAATCCTGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCt                                                                                                     <  5:261436/93‑1 (MQ=255)
                                                                                                        ggAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTg                                          <  4:322413/149‑1 (MQ=255)
                                                                                                            tCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGttt                                      >  1:12919/1‑149 (MQ=255)
                                                                                                              ccGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGttttt                                    <  5:99568/149‑1 (MQ=255)
                                                                                                               cGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCgg                                                >  6:140901/1‑136 (MQ=255)
                                                                                                               cGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCgg                                                <  5:140901/136‑1 (MQ=255)
                                                                                                               cGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTc                                   <  8:112171/149‑1 (MQ=255)
                                                                                                                            ggACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTc                             >  4:340315/1‑142 (MQ=255)
                                                                                                                            ggACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTc                             <  3:340315/142‑1 (MQ=255)
                                                                                                                              aCTGTCCGGCCGACTGACCGCGGGCTGTCAGGGGGATATCCAGCAGTTTCGAGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCTCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActgc                    <  6:127478/149‑1 (MQ=255)
                                                                                                                               cTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActgct                   <  2:57956/149‑1 (MQ=255)
                                                                                                                                tGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTActgctg                  >  3:147924/1‑149 (MQ=255)
                                                                                                                                  gccggccgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGca                >  6:150655/1‑149 (MQ=255)
                                                                                                                                       ccgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCatcatc           >  2:89147/1‑149 (MQ=255)
                                                                                                                                          ccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGCAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTc                                   <  6:346838/121‑1 (MQ=255)
                                                                                                                                          ccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGCAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTc                                   >  5:346838/2‑122 (MQ=255)
                                                                                                                                          aCTGCCCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGAATATACGTAGAAAATCTCGAAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCggg        >  5:271661/1‑149 (MQ=255)
                                                                                                                                           cTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTGTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCgggg       >  1:280351/1‑149 (MQ=255)
                                                                                                                                              aCCGCGGGCGGTCAGGGTGATATCCAGCAGATTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgt      <  8:85870/147‑1 (MQ=255)
                                                                                                                                              aCCGCGGGCGGTAAGGGTGATATCCAGCAGTGTAGCGGCCTTTTTCGCATCTTCGTATAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtg     <  8:158303/148‑1 (MQ=255)
                                                                                                                                                cgcgGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtgtgt  >  7:327184/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGT  >  NC_002947/1824945‑1825236

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