Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,736,522:1 +G 6.2% coding (924/1674 nt) hutU ← urocanate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,736,5221.G 6.2% 260.1 / 10.2 65coding (924/1674 nt)hutUurocanate hydratase
Reads supporting (aligned to +/- strand):  ref base . (33/28);  new base G (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

CGCCGCGGCAGAACAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACC  >  NC_002947/5736374‑5736666
                                                                                                                                                     |                                                                                                                                                
cgccgcGGCAGAACAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGcc                                                                                                                                                   <  4:43368/149‑1 (MQ=255)
      ggCAGAACAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCtgcacg                                                                                                                                            >  7:260484/1‑149 (MQ=255)
          gAACAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCtgcacgtgca                                                                                                                                        <  5:46159/149‑1 (MQ=255)
           aaCAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCtgcacgtgcac                                                                                                                                       >  5:148232/1‑149 (MQ=255)
             cAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACgg                                                                                                                                     <  7:68335/149‑1 (MQ=255)
                              ccGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTggcg                                                                                                                    <  5:366619/149‑1 (MQ=255)
                               cGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTg                                                                                                                           >  1:65052/1‑141 (MQ=255)
                               cGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTg                                                                                                                           <  2:65052/141‑1 (MQ=255)
                                    cGAAGCCTGGGAATTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCtt                                                                                                               >  8:167189/1‑148 (MQ=255)
                                    cGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTg                                                                                                              <  2:157651/149‑1 (MQ=255)
                                     gAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGa                                                                                                             >  6:358186/1‑149 (MQ=255)
                                       aGCCTGGTAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGTATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGAc                                                                                                            <  5:334906/148‑1 (MQ=255)
                                        gCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCa                                                                                                          <  1:240899/149‑1 (MQ=255)
                                            gggAATTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGc                                                                                                      <  7:167189/149‑1 (MQ=255)
                                            gggAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTg                                                                                                       >  6:274042/1‑148 (MQ=255)
                                             ggAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCa                                                                                                     <  4:323649/149‑1 (MQ=255)
                                                aGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGcc                                                                                                  >  8:25401/1‑149 (MQ=255)
                                                     aaGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTc                                                                                             <  6:148232/149‑1 (MQ=255)
                                                        gCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGt                                                                                          >  7:36162/1‑149 (MQ=255)
                                                            tGGCCACACCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTgc                                                                                      >  4:323284/1‑149 (MQ=255)
                                                              gCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTgcgc                                                                                    >  6:39265/1‑149 (MQ=255)
                                                              gCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTgcg                                                                                     >  5:40949/1‑148 (MQ=255)
                                                               ccACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCAAGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCa                                                                                   >  2:318442/1‑149 (MQ=255)
                                                                         ctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATcgcg                                                                           <  3:329163/147‑1 (MQ=255)
                                                                         ctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCgg                                                                          <  3:289687/148‑1 (MQ=255)
                                                                                 ccATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTcc                                                                 <  6:86024/149‑1 (MQ=255)
                                                                                   aTCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCGTGCACGTGCAGGCCAtgctg                                                                                                                             <  1:54323/87‑5 (MQ=38)
                                                                                   aTCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCGTGCACGTGCAGGCCAtgctg                                                                                                                             >  2:54323/1‑83 (MQ=38)
                                                                                       ggCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTc                                                           >  2:192108/1‑149 (MQ=255)
                                                                                        gCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCTGCCTTGACAACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTc                                                           >  7:7985/1‑148 (MQ=255)
                                                                                        gCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAg                                                              <  5:274042/145‑1 (MQ=255)
                                                                                        gCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTc                                                           >  5:2398/1‑148 (MQ=255)
                                                                                          ggATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAg                                                        >  7:274191/1‑149 (MQ=255)
                                                                                                     ccGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAg                                             <  1:192108/149‑1 (MQ=255)
                                                                                                       gTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGt                                           <  3:323284/149‑1 (MQ=255)
                                                                                                         aaTCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAg                                         >  5:251216/1‑149 (MQ=255)
                                                                                                           tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTc                                                           <  6:230396/129‑1 (MQ=255)
                                                                                                           tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTc                                                           >  5:230396/1‑129 (MQ=255)
                                                                                                             gAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGt                                     >  5:343901/1‑149 (MQ=255)
                                                                                                             gAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGt                                     >  3:300175/1‑149 (MQ=255)
                                                                                                                 gTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTcc                                                                 <  5:320138/117‑1 (MQ=255)
                                                                                                                 gTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTcc                                                                 >  6:320138/1‑117 (MQ=255)
                                                                                                                    gggATGCCCTGCTTCTGGAAGTCTAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGACGGTTCGGTCTGCGCACGATCGAGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCAGTTCAGTgg                              >  5:297870/1‑149 (MQ=255)
                                                                                                                    gggATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTgg                              >  3:249372/1‑149 (MQ=255)
                                                                                                                    gggATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTgg                              <  6:2398/149‑1 (MQ=255)
                                                                                                                           ccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTgg                       >  5:68944/1‑149 (MQ=255)
                                                                                                                           ccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTgg                       >  1:295448/1‑149 (MQ=255)
                                                                                                                           ccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTgg                       >  3:341887/1‑149 (MQ=255)
                                                                                                                                ttCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCgg                                                    <  6:316717/115‑1 (MQ=255)
                                                                                                                                ttCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCgg                                                    >  5:316717/1‑115 (MQ=255)
                                                                                                                                ttCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGt                         >  5:46499/1‑142 (MQ=255)
                                                                                                                                ttCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGt                         <  6:46499/142‑1 (MQ=255)
                                                                                                                                  cTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctg                 <  5:358186/148‑1 (MQ=255)
                                                                                                                                      aaGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctggtctg            <  7:80762/149‑1 (MQ=255)
                                                                                                                                        gTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctggtctggt          >  2:102168/1‑149 (MQ=255)
                                                                                                                                         tCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCtggtctggtc         >  5:188777/1‑149 (MQ=255)
                                                                                                                                         tCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTGCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCtggtctggtc         <  6:273340/149‑1 (MQ=255)
                                                                                                                                          ccAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCg        >  8:243296/1‑149 (MQ=255)
                                                                                                                                           cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCa                                                                                                          >  8:134663/1‑50 (MQ=255)
                                                                                                                                           cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCa                                                                                                          <  7:134663/50‑1 (MQ=255)
                                                                                                                                           cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTTGGTCGTGGGCGCTGGTCTGGTCgg       >  8:161428/1‑149 (MQ=255)
                                                                                                                                           cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCgg       >  3:174111/1‑149 (MQ=255)
                                                                                                                                              cATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTg                                                 >  2:147325/1‑104 (MQ=255)
                                                                                                                                              cATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTg                                                 <  1:147325/104‑1 (MQ=255)
                                                                                                                                               aTGGCCGTGCACGTGCAGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGgc                     >  5:230870/1‑131 (MQ=255)
                                                                                                                                               aTGGCCGTGCACGTGCAGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGgc                     <  6:230870/131‑1 (MQ=255)
                                                                                                                                                tGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGAcc  <  5:39265/149‑1 (MQ=255)
                                                                                                                                                tGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGAcc  <  8:260484/149‑1 (MQ=255)
                                                                                                                                                tGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGAcc  <  8:36162/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CGCCGCGGCAGAACAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACC  >  NC_002947/5736374‑5736666

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: