Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,736,533 C→T 6.7% V305M (GTG→ATG)  hutU ← urocanate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,736,5330CT 6.7% 171.5 / 2.4 62V305M (GTG→ATG) hutUurocanate hydratase
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base T (2/2);  total (31/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 = 7.32e-01

ACAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCGGGC  >  NC_002947/5736386‑5736676
                                                                                                                                                    |                                                                                                                                                  
aCAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑tgcacgtgcac                                                                                                                                                    <  3:177927/148‑1 (MQ=255)
      ggCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATc                                                                                                                                             >  6:121715/1‑149 (MQ=255)
         ggATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGAc                                                                                                                                          >  5:139722/1‑149 (MQ=255)
           aTGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTg                                                                                                                                        <  3:105353/149‑1 (MQ=255)
                   cGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTggcgg                                                                                                                                >  5:34136/1‑149 (MQ=255)
                    gggACGAAGCCTGGGAAGTCCAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGACATCGACTGCTAggcgg                                                                                                                                >  5:393165/1‑148 (MQ=255)
                       aCGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCtt                                                                                                                            <  5:121715/149‑1 (MQ=255)
                         gatgCCTGGGAAGTCGAAGGCATTGGCGAGGCCCTCCTCCTTTGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGTGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGa                                                                                                                          <  4:101263/146‑1 (MQ=255)
                         gAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCATGGCCATCGACTGCTTGGCGGCCTTGa                                                                                                                          >  7:145781/1‑149 (MQ=255)
                         gAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCATGGCCATCGACTGCTTGGCGGCCTTGa                                                                                                                          >  3:184998/1‑149 (MQ=255)
                                   aaGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGc                                                                                                                <  5:185382/149‑1 (MQ=255)
                                        gAAGGCATTGGCCAGGCCCTCCTCGTTGGCCATCTGGCGGATGTTGTTGCCGTAATCTAAGGTCGGGATGCCCTGCTTCTGTAATTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGtt                                                                                                           <  3:365972/149‑1 (MQ=255)
                                             cATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑Tc                                                                                                     >  3:205741/1‑149 (MQ=255)
                                              aTTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCt                                                                                                    >  8:163267/1‑149 (MQ=255)
                                                   ccACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCa                                                                                               >  6:319638/1‑149 (MQ=255)
                                                    cACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCAc                                                                                              >  3:154416/1‑149 (MQ=255)
                                                      cGCCCTCGTCCTTGGCCATCTGGCGGATGTTGTTGCCTTAATCGAAGGTGGGTATGCCCTGCTTCTGGAAGTCAAGCATGGCC‑TGCACGTGCAGGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGa                                                                                            <  7:202848/149‑1 (MQ=255)
                                                         cctcctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATc‑‑g                                                                                       >  5:119665/1‑149 (MQ=255)
                                                            cctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGGCCACCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCgg                                                                                    <  4:371020/149‑1 (MQ=255)
                                                             ctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCgg                                                                                    <  5:319638/148‑1 (MQ=255)
                                                               cctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTAc                                                                                 >  3:175506/1‑149 (MQ=255)
                                                               cctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTAc                                                                                 >  4:145316/1‑149 (MQ=255)
                                                                ctTGGCCATCTGGCGGATGTTGTTGACGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTa                                                                                  <  4:205741/147‑1 (MQ=255)
                                                                     ccATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTcc                                                                           <  2:340127/149‑1 (MQ=255)
                                                                       aTCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTcc                                                                           >  2:329755/1‑147 (MQ=255)
                                                                            gCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCATGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAg                                                                        <  8:145781/145‑1 (MQ=255)
                                                                            gCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCATGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAg                                                                        <  4:184998/145‑1 (MQ=255)
                                                                               gATGTTGTTGCCGTAATCTTAGGTCTGGATGCCCTGCTTCTGGAAGTCCAGCTTGGCC‑GGATCGTGCACGTCCATCGACTGCTTGGCTGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCTGc                                                                 <  6:393165/149‑1 (MQ=255)
                                                                                    tgttgCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCt                                                            >  4:24622/1‑149 (MQ=255)
                                                                                        gCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGATCTCCAGCATGGCC‑TGCACGTGCACGGCAATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTgg                                                          >  1:255680/1‑147 (MQ=255)
                                                                                        gCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGCCCACTGCAGCCGGTACGG‑TCTGCGCACGAGC‑‑GCGGTACTGTTCCCAGGTCAACCCGGCTgg                                                          >  8:33871/1‑147 (MQ=255)
                                                                                        gCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑GCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGc                                                         <  7:163267/148‑1 (MQ=255)
                                                                                         ccGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAg                                                       <  6:139722/149‑1 (MQ=255)
                                                                                         ccGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAg                                                       >  8:114148/1‑149 (MQ=255)
                                                                                          cGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAgg                                                      >  3:76210/1‑149 (MQ=255)
                                                                                           gTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAgg                                                      >  1:275051/1‑148 (MQ=255)
                                                                                                 gAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGcc                                                 <  2:186566/147‑1 (MQ=255)
                                                                                                 gAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGcc                                                 >  1:186566/1‑147 (MQ=255)
                                                                                                  aaGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGtt                                              <  4:68175/149‑1 (MQ=255)
                                                                                                      tCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGt                                          <  1:66779/149‑1 (MQ=255)
                                                                                                              cccTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGt                                   <  3:24622/148‑1 (MQ=255)
                                                                                                               ccTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTgg                                 <  7:114148/149‑1 (MQ=255)
                                                                                                               ccTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTgg                                 >  4:357731/1‑149 (MQ=255)
                                                                                                               ccTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTgg                                 <  7:358288/149‑1 (MQ=255)
                                                                                                                      cTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctg                           <  4:76210/148‑1 (MQ=255)
                                                                                                                      cTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCAAGGTCAAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctgg                          >  5:227909/1‑149 (MQ=255)
                                                                                                                        ggAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTgg                                 >  3:5044/1‑140 (MQ=255)
                                                                                                                        ggAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTgg                                 <  4:5044/140‑1 (MQ=255)
                                                                                                                            gTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctggtctggt                    <  3:145316/149‑1 (MQ=255)
                                                                                                                             tCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCtggtctggtc                   >  4:129279/1‑149 (MQ=255)
                                                                                                                               gtcgaTGGCCGTGCACGTGCA‑GGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGATCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGAGGGCGCTGGTATGGTCgg                 >  7:206670/5‑149 (MQ=255)
                                                                                                                               cAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGt                                          <  5:330153/124‑1 (MQ=255)
                                                                                                                               cAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGt                                          >  6:330153/1‑124 (MQ=255)
                                                                                                                               cAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTc                                     >  4:225600/1‑129 (MQ=255)
                                                                                                                                 gCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGcc                                                                <  8:316272/100‑1 (MQ=255)
                                                                                                                                 gCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGcc                                                                >  7:316272/1‑100 (MQ=255)
                                                                                                                                 gCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTg               <  4:175506/149‑1 (MQ=255)
                                                                                                                                  cATGGGC‑TGCTCGTGCACGTCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGCACTGTTGCGAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGt                <  2:255680/147‑1 (MQ=255)
                                                                                                                                    tggccgtgcacgtgca‑ggCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGAc             >  4:25718/7‑148 (MQ=255)
                                                                                                                                    tGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGAcc            <  3:357731/149‑1 (MQ=255)
                                                                                                                                      ggc‑TGCACGGGCAGTGCCATCGACTGCTTGGGGGCCGTGTCCACTGCAGCCGGTTGGGTTCTGCTCACGATCCTTCAGTACTGTTCCCAGGGCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGa              <  1:329755/146‑1 (MQ=255)
                                                                                                                                      gCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCAt          >  5:141840/1‑149 (MQ=255)
                                                                                                                                      gCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCAt          <  4:154416/149‑1 (MQ=255)
                                                                                                                                       cc‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATg         >  1:269632/1‑149 (MQ=255)
                                                                                                                                               gtgcacgGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCGGGc  >  2:277652/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACAGTGGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCC‑TGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGG‑TCTGCGCACGATC‑‑GCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCGGGC  >  NC_002947/5736386‑5736676

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑