Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,412,258 G→A 6.3% E68K (GAG→AAG)  tusA‑I → sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,412,2580GA 6.3% 186.1 / 2.1 63E68K (GAG→AAG) tusA‑Isulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

CCTTTGCCACTGCTCAAGGCCAAGATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCGGTCCTGTTGTTCCTGGCTTTTACC  >  NC_002947/1412111‑1412406
                                                                                                                                                   |                                                                                                                                                    
ccTTTGCCACTGCTCAAGGCCAAGATGGAGCTCAATGGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGa                                                                                                                                                     <  5:210995/149‑1 (MQ=255)
 cTTTGCCACTGCTCAAGGCCAAGATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAg                                                                                                                                                    >  5:119119/1‑149 (MQ=255)
       cACTGCTCAAGGCCAAGATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCgg                                                                                                                                               >  1:138328/1‑148 (MQ=255)
             tCAAGGCCAAGATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGt                                                                                                                                              <  1:43828/143‑1 (MQ=255)
             tCAAGGCCAAGATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGt                                                                                                                                              >  2:43828/1‑143 (MQ=255)
                   ccAAGATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTAtacttac                                                                                                                                  >  4:268032/1‑149 (MQ=255)
                    cAAGATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGTTACCTAtacttact                                                                                                                                 >  3:288154/1‑149 (MQ=255)
                       gATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCg                                                                                                                                                      <  4:341236/125‑1 (MQ=255)
                       gATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCg                                                                                                                                                      >  3:341236/1‑125 (MQ=255)
                           gAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTgc                                                                                                                           <  3:268032/148‑1 (MQ=255)
                              cTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTgcgc                                                                                                                         >  5:331783/1‑147 (MQ=255)
                              cTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCaa                                                                                                                       <  6:119119/149‑1 (MQ=255)
                                 aaTCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGTTACCTATACTTACTGGCTGCGCAAGGc                                                                                                                    <  4:288154/149‑1 (MQ=255)
                                   tCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCt                                                                                                                  <  4:94030/149‑1 (MQ=255)
                                         ggCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCt                                                                                                            >  5:198211/1‑149 (MQ=255)
                                          gCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCtt                                                                                                           <  6:8714/149‑1 (MQ=255)
                                                 gcgcGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCaa                                                                                                     <  6:331783/148‑1 (MQ=255)
                                                 gcgcGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAg                                                                                                    >  4:25021/1‑149 (MQ=255)
                                                       tGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGt                                                                                              <  6:198211/149‑1 (MQ=255)
                                                            aaGGTTATCGCCACGGACGCAGGCTCCAAGCGCGACTTCCTCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCaa                                                                                           <  3:25021/147‑1 (MQ=255)
                                                             aGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGt                                                                                        >  7:128179/1‑149 (MQ=255)
                                                                      ccACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGc                                                                               <  8:128179/149‑1 (MQ=255)
                                                                        aCGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCtt                                                                             >  5:12837/1‑149 (MQ=255)
                                                                          ggACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCaaa                                                                                   >  1:334481/1‑141 (MQ=255)
                                                                          ggACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCaaa                                                                                   <  2:334481/141‑1 (MQ=255)
                                                                          ggACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTcg                                                                           <  6:12837/149‑1 (MQ=255)
                                                                             cGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCt                                                                                                                  >  4:65350/1‑107 (MQ=255)
                                                                             cGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCt                                                                                                                  <  3:65350/107‑1 (MQ=255)
                                                                                aGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTg                                                                     >  4:7694/1‑149 (MQ=255)
                                                                                 ggCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTgg                                                                    <  4:134738/149‑1 (MQ=255)
                                                                                 ggCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTgg                                                                    <  4:129625/149‑1 (MQ=255)
                                                                                      ccAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTgg                                                                                                                               <  7:264738/85‑1 (MQ=255)
                                                                                      ccAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTgg                                                                                                                               >  8:264738/1‑85 (MQ=255)
                                                                                             gACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTAc                                                        >  8:241816/1‑149 (MQ=255)
                                                                                                ttCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCAAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTtc                                                     >  3:154165/1‑149 (MQ=255)
                                                                                                ttCCGCACTTTCGCCCAGTTGGCCGGTAATACCCTGCTGCAAGAAACGGCCAATGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTtc                                                     >  3:159458/1‑149 (MQ=255)
                                                                                                        tttCTCCGAGTTGTCCGGTCATACCCTGCTGCAAGTAACGGCCGAGGCCGGTACCTATACTTAGTGTCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCa                                                               <  8:143689/131‑1 (MQ=255)
                                                                                                        tttCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGAt                                               <  7:41887/147‑1 (MQ=255)
                                                                                                        tttCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGAt                                               <  7:47210/147‑1 (MQ=255)
                                                                                                        tttCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGAt                                               >  8:41887/1‑147 (MQ=255)
                                                                                                        tttCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGAt                                               >  8:47210/1‑147 (MQ=255)
                                                                                                        tttCGCCCAGTTGGCAGGTCATACCCTGATGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGCCTGCGCAAGGCCTGAGTCATTTCCAACGATCGTCAATCTTCAAAGTGCTTCGCGACTGGCTGCa                                                               >  7:143689/1‑131 (MQ=255)
                                                                                                          tCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGa                                             <  8:18620/147‑1 (MQ=255)
                                                                                                             cccAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCt                                                          >  3:117333/1‑131 (MQ=255)
                                                                                                             cccAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCt                                                          <  4:117333/131‑1 (MQ=255)
                                                                                                                  ttGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCggtg                                   <  7:241816/149‑1 (MQ=255)
                                                                                                                         gTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTAc                                                        >  6:19615/1‑121 (MQ=255)
                                                                                                                         gTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTAc                                                        <  5:19615/121‑1 (MQ=255)
                                                                                                                          tCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGc                           >  5:183556/1‑149 (MQ=255)
                                                                                                                           cATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCg                          >  1:128937/1‑149 (MQ=255)
                                                                                                                                ccTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCTCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCggtg                                   <  6:263130/135‑1 (MQ=255)
                                                                                                                                ccTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCggtg                                   >  5:263130/1‑135 (MQ=255)
                                                                                                                                ccTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCgg                                     >  4:65583/1‑133 (MQ=255)
                                                                                                                                ccTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCgg                                     <  3:65583/133‑1 (MQ=255)
                                                                                                                                ccTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGc                               <  7:29859/139‑1 (MQ=255)
                                                                                                                                ccTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGc                               >  8:29859/1‑139 (MQ=255)
                                                                                                                                ccTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCGGTCCt                     >  6:255203/1‑149 (MQ=255)
                                                                                                                                    cttcAAGAAACGGCCAATGACGGTACCTATAATTACTGGATGCGAAAGTCATGAGTCTTTTCCAATGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGTTTGTTCTGGCGTTCCtgttg                 <  4:159458/146‑1 (MQ=255)
                                                                                                                                    ctgcAAGAAACGGCCAAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCGGTCCtgttg                 <  4:154165/149‑1 (MQ=255)
                                                                                                                                        aaGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTg                             <  6:88725/133‑1 (MQ=255)
                                                                                                                                        aaGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTg                             >  5:88725/1‑133 (MQ=255)
                                                                                                                                        aaGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCGGTCCTGTTGTTcc             >  1:320758/1‑149 (MQ=255)
                                                                                                                                            aaCGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCGGTCCTGTTGTTCCTgg          <  6:183556/148‑1 (MQ=255)
                                                                                                                                                 ccGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCGGTCCTGTTGTTCCTGGCtttt     >  8:29190/1‑148 (MQ=255)
                                                                                                                                                 ccGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCGGTCCTGTTGTTCCTGGCTTTTa    >  8:107478/1‑149 (MQ=255)
                                                                                                                                                   gAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCGGTCCTGTTGTTCCTGGCTTTTAcc  <  2:128937/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTTTGCCACTGCTCAAGGCCAAGATGGAGCTCAATCGCCTGGCCAGTGGCGCGGTGCTCAAGGTTATCGCCACGGACGCAGGCTCCCAGCGCGACTTCCGCACTTTCGCCCAGTTGGCCGGTCATACCCTGCTGCAAGAAACGGCCGAGGCCGGTACCTATACTTACTGGCTGCGCAAGGCCTGAGTCTTTTCCAAGGATCGTCAATGTTCAAAGTGCTTCGCGACTGGGTGCAGCGCTACTTCTCCGATGAGGAAGCGGTGGTGCTGGCGGTCCTGTTGTTCCTGGCTTTTACC  >  NC_002947/1412111‑1412406

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