Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,671 G→C 10.3% P311R (CCG→CGG)  mutY ← A/G‑specific adenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,6710GC10.3% 194.3 / 9.1 67P311R (CCG→CGG) mutYA/G‑specific adenine glycosylase
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base C (4/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCG  >  NC_002947/344523‑344815
                                                                                                                                                    |                                                                                                                                                    
tCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCg                                                                                                                                                      >  8:131044/1‑149 (MQ=255)
        cGGGTCATCGGTCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGAc                                                                                                                                               <  8:69857/148‑1 (MQ=255)
         gggTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCATCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGt                                                                                                                                             <  2:111270/149‑1 (MQ=255)
           gTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGCTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGc                                                                                                                                           >  8:296258/1‑149 (MQ=255)
              aTCTGGCGTCTCCTGAGATTCGGTAGTCTGCAGCGAGTACCAGGATCTTTTTAACCTTGTCGTTGAGGCCCATGTGGTGCTGGTGGGAGATGTTTTTCTAGAGCCAGTCGGCCTCGGCCACGTGCACGTCGACCCGGTCGGCGAGCAc                                                                                                                                         <  6:265692/148‑1 (MQ=255)
                       ctcCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGc                                                                                                                               <  3:72603/149‑1 (MQ=255)
                         ccTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTc                                                                                                                             >  7:10079/1‑149 (MQ=255)
                          cTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTc                                                                                                                             >  8:155352/1‑148 (MQ=255)
                           tGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGa                                                                                                                           >  4:272091/1‑149 (MQ=255)
                             cGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATc                                                                                                                         <  7:131044/149‑1 (MQ=255)
                             cGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTATCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGAt                                                                                                                          <  6:163752/148‑1 (MQ=255)
                               aTCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCACGTCGACCCGGTCGGCGAGCACgt                                                                                                                                       <  8:135528/133‑2 (MQ=255)
                               aTCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCACGTCGACCCGGTCGGCGAGCACgt                                                                                                                                       >  7:135528/1‑132 (MQ=255)
                                   gTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGc                                                                                                                                           >  4:118294/1‑125 (MQ=255)
                                   gTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGc                                                                                                                                           <  3:118294/125‑1 (MQ=255)
                                            ccGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                          >  7:161105/1‑149 (MQ=255)
                                            ccGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                          <  3:81542/149‑1 (MQ=255)
                                             cgcgcgTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGc                                                                                                         <  8:161105/149‑1 (MQ=255)
                                                 cgTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGaa                                                                                                     <  3:272091/149‑1 (MQ=255)
                                                  gTTCAAGCAGCTTTTTCACCGGGGCGGCGAGACCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAg                                                                                                    >  5:121202/1‑149 (MQ=255)
                                                   ttCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCAGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAgg                                                                                                   >  5:120260/1‑149 (MQ=255)
                                                      aaGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTAt                                                                                                >  5:232700/1‑149 (MQ=255)
                                                      aaGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTAt                                                                                                >  8:91738/1‑149 (MQ=255)
                                                       agcagcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATg                                                                                               >  3:190517/1‑149 (MQ=255)
                                                        gcagcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATgg                                                                                              >  7:243281/1‑149 (MQ=255)
                                                          agcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGt                                                                                            >  5:58567/1‑149 (MQ=255)
                                                          agcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGt                                                                                            >  8:84295/1‑149 (MQ=255)
                                                             tttttCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGCTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAg                                                                                         <  7:296258/149‑1 (MQ=255)
                                                               tttCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCAGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAgg                                                                                        <  6:120260/148‑1 (MQ=255)
                                                                   aCCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTc                                                                                   >  6:102913/1‑149 (MQ=255)
                                                                    ccGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTcc                                                                                  <  8:274066/149‑1 (MQ=255)
                                                                    ccGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTc                                                                                   >  1:2106/1‑148 (MQ=255)
                                                                    ccGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTc                                                                                   >  5:14959/1‑148 (MQ=255)
                                                                      ggggCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCAt                                                                                <  6:58567/149‑1 (MQ=255)
                                                                       gggCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCAtg                                                                               <  6:232700/149‑1 (MQ=255)
                                                                            gcgAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCAtggcgtg                                                                          >  5:305559/1‑149 (MQ=255)
                                                                                ggCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTggc‑‑‑gg                                                                   <  4:190517/149‑1 (MQ=255)
                                                                                 gCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTggc‑‑‑gg‑c                                                                 >  3:199982/1‑149 (MQ=255)
                                                                                      ggCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑Gg‑ccagc                                                             >  8:263574/1‑148 (MQ=255)
                                                                                            gcggGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGgcgc                                                      <  7:155352/149‑1 (MQ=255)
                                                                                                 gTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGcc                                                 <  8:10079/149‑1 (MQ=255)
                                                                                                    gCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCTTGGC‑‑‑GG‑CCAGCCAGGCGCATGCCggg                                              >  8:223957/1‑147 (MQ=255)
                                                                                                     cGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGc                                             <  5:102913/149‑1 (MQ=255)
                                                                                                     cGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGc                                             <  2:338700/149‑1 (MQ=255)
                                                                                                     cGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGc                                             <  6:14959/149‑1 (MQ=255)
                                                                                                      gAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑TCAGCCAGGCGCAGGCCGTGCt                                            >  1:188152/1‑149 (MQ=255)
                                                                                                      gAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCt                                            >  5:66357/1‑149 (MQ=255)
                                                                                                            ataCCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAg                                       <  7:91738/148‑1 (MQ=255)
                                                                                                             taCCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGc                                     <  6:305559/149‑1 (MQ=255)
                                                                                                               ccAGAGCCAGTCGGCCGCGGCCACGGGCTCGCCGACCGGGGCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGcc                                    <  7:143239/148‑1 (MQ=255)
                                                                                                                cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑Gg‑ccagcca                                                           <  1:164961/124‑1 (MQ=255)
                                                                                                                cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑Gg‑ccagcca                                                           >  2:164961/1‑124 (MQ=255)
                                                                                                                cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAg                                  <  4:199982/149‑1 (MQ=255)
                                                                                                                cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCGGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAg                                  >  3:223011/1‑149 (MQ=255)
                                                                                                                 agagCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAgg                                 >  3:50355/1‑149 (MQ=255)
                                                                                                                            gcctcggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTg                      <  1:197718/149‑1 (MQ=255)
                                                                                                                             cctcggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTgg                     <  7:177607/149‑1 (MQ=255)
                                                                                                                              ctcggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGt                                         >  4:197664/1‑128 (MQ=255)
                                                                                                                              ctcggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGt                                         <  3:197664/128‑1 (MQ=255)
                                                                                                                                  gccACGTGCTCGCCGACCGTGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGAt                <  1:308635/149‑1 (MQ=255)
                                                                                                                                   gctcGTGCACGTCGACCCGGTCGGCGAGCACGTGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATg               >  2:180979/4‑149 (MQ=255)
                                                                                                                                   ccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGTGGGTGCCGCCAGGCGCCCGCCGTGCTGGGAGGCCAGGTCGTCGAGCTGGGc                   >  1:283815/1‑149 (MQ=255)
                                                                                                                                    cACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATgt              >  7:181180/1‑149 (MQ=255)
                                                                                                                                     tcGTGCACGTCGACCCGGTCGGCGAGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATgtc             >  8:67743/2‑149 (MQ=255)
                                                                                                                                       gTGCACGTCGACCCGGTCGGCGAGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATgtcgt           >  8:290363/1‑149 (MQ=255)
                                                                                                                                        tGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATgtcgtc          <  7:293576/149‑1 (MQ=255)
                                                                                                                                          cTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑TCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCa        <  2:188152/149‑1 (MQ=255)
                                                                                                                                                gACCCGGTCGGCGAGCACGTGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCg  <  1:180979/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCG  >  NC_002947/344523‑344815

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