Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,662 T→A 8.8% E314V (GAG→GTG)  mutY ← A/G‑specific adenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,6620TA 8.8% 192.5 / 8.2 68E314V (GAG→GTG) mutYA/G‑specific adenine glycosylase
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base A (4/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GGCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCA  >  NC_002947/344517‑344809
                                                                                                                                                 |                                                                                                                                                       
ggCACATCACCGTCCGGGTCATCGGGCGTCTCCTGCGATCAGATCGACGGCCGCCCGATCAAGCAGCTTTATCACCGGGGCGGCCACGCCCAGGCGCGGCCGGGTCGCGAGGGTATACCAGAGCCAGTCGGCCTCGGCCACGTgcgcgc                                                                                                                                                      >  8:143239/1‑149 (MQ=255)
 gCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGcc                                                                                                                                                     >  7:274066/1‑149 (MQ=255)
      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                                                                                                                                     <  3:118294/125‑1 (MQ=255)
                                         gTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGc                                                                                                                                     >  4:118294/1‑125 (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                                                                             >  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                                       <  6:14959/149‑1 (MQ=255)
                                                                                                           cGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGc                                       <  2:338700/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                                                     >  2:164961/1‑124 (MQ=255)
                                                                                                                      cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑Gg‑ccagcca                                                     <  1:164961/124‑1 (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                                   <  3:197664/128‑1 (MQ=255)
                                                                                                                                    ctcggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGt                                   >  4:197664/1‑128 (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)
                                                                                                                                                 |                                                                                                                                                       
GGCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGC‑‑‑GG‑CCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCA  >  NC_002947/344517‑344809

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