Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,677 A→G 10.3% V309A (GTC→GCC)  mutY ← A/G‑specific adenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,6770AG10.3% 178.7 / 9.7 68V309A (GTC→GCC) mutYA/G‑specific adenine glycosylase
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base G (4/3);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCT  >  NC_002947/344531‑344822
                                                                                                                                                  |                                                                                                                                                     
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                                                                                                                                              >  7:135528/1‑132 (MQ=255)
                       aTCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCACGTCGACCCGGTCGGCGAGCACgt                                                                                                                                              <  8:135528/133‑2 (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                                                                                                   >  8:84295/1‑149 (MQ=255)
                                                  agcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGt                                                                                                   >  5:58567/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                                                                                          >  5:14959/1‑148 (MQ=255)
                                                            ccGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTc                                                                                          >  1:2106/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                                                    <  6:14959/149‑1 (MQ=255)
                                                                                             cGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGc                                                    <  2:338700/149‑1 (MQ=255)
                                                                                             cGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGc                                                    <  5:102913/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)
                                                                                                                                               tCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGCTGTGGTCCAGGTCCgggggggt  >  4:266898/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGGGAGGCT  >  NC_002947/344531‑344822

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