Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,739,087 T→G 8.4% intergenic (‑108/+49) proY ← / ← hutH proline (histidine) APC transporter/histidine ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,739,0870TG 8.4% 526.3 / 4.5 218intergenic (‑108/+49)proY/hutHproline (histidine) APC transporter/histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (103/86);  new base G (15/10);  total (121/97)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

GCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCAC  >  NC_002947/5738938‑5739236
                                                                                                                                                     |                                                                                                                                                      
  gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                                                      <  1:691691/150‑1 (MQ=255)
   ggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgtg                                                                                                                                                     <  1:6706037/150‑1 (MQ=255)
   ggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgtg                                                                                                                                                     <  1:945450/150‑1 (MQ=255)
      aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTTAATTGTTTTTGTGGTCGTGAGATTCTGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                                                                                  <  1:3558780/150‑1 (MQ=255)
      aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                                                                                  <  2:1526952/150‑1 (MQ=255)
        ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTTCACCTACTTCTGTGGGAGc                                                                                                                                                >  2:2348982/1‑150 (MQ=255)
        ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGc                                                                                                                                                >  2:5765877/1‑150 (MQ=255)
         tcgCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCg                                                                                                                                               <  2:6016864/148‑1 (MQ=255)
          aGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCGGTGGGAGCgg                                                                                                                                              >  1:6166280/1‑150 (MQ=255)
            cccGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgc                                                                                                                                            <  1:2604531/150‑1 (MQ=255)
            cccGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACGTCTGTGGGAGCGGgc                                                                                                                                            >  2:506057/1‑150 (MQ=255)
                 gcTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCAAGCCCCCTCCCACAGGGGATGGTGTTTTGCACCAACTTAGGTGGGAGCGGGCGCGcc                                                                                                                                       >  2:1817375/1‑150 (MQ=255)
                      agaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCCCGCCCACCGGGGATGGTGTTTTGCCCCTACTTCTGTGGGAGCGGGCGCGCCCCCGa                                                                                                                                  >  2:2025108/1‑150 (MQ=255)
                       gaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTACCACAGGGGAAGGTGTTTTGCACCTACATCTGTGGGAGCGGGCGCGCCCGCGaa                                                                                                                                 >  2:2244531/1‑150 (MQ=255)
                        aCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGTTGGTGGTTTGCACCTCCTTCTGTGGGTGCGGGCGCGCCCGCGAAg                                                                                                                               >  2:2549157/1‑150 (MQ=255)
                           ttGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGACGGgc                                                                                                                            >  2:5418693/1‑150 (MQ=255)
                             gAGCTTGTTGCATGTCCCGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCAGCTCCCACAGGGTATGTTTTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGgccg                                                                                                                          <  1:4722092/150‑1 (MQ=255)
                               gCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGGCCCCCCCCGCCCCCCGGGGCTGGGGGTTTGCGCCGACGTCGGGGGGGGCGGGGGCGCCCGCGGCGGggcgcc                                                                                                                        >  1:2962910/1‑150 (MQ=255)
                               gCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACGTCTGTGGGAGCGGGCGCGCCCGCGAAGGgccgcc                                                                                                                        >  2:6797297/1‑150 (MQ=255)
                               gCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCAAGCCCCCGCCCCCACGGGTTTGGGGTTTGGACCGAAATGGGTGGGGGCGGGGGCGCCCCCCCAGGgccccc                                                                                                                        >  1:2962678/1‑150 (MQ=255)
                                  tgttgCATGTCACGTCCTTGAATTTTTTTTGTGGTCGTGAGATTCTGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAg                                                                                                                     <  2:4072084/150‑1 (MQ=255)
                                            cccGTCCTTGACTTGTGTTTGTGGTCGTGAGATTCAGGGGGCTGGTGCGAAGCCCCTTCGCGGGCACGCCCGCTCCCACTGGGGATGGTGTTTTGCACCTACTTCCGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAgg                                                                                                           <  2:3837211/148‑1 (MQ=255)
                                                     taaTTTTTTTTGTGGTCGTGTGATTCTGGGGGCGGCTGCGCAGCCCCTTCGGGGGCACGCCCGCTCCCACTGGGTTTGGTGTTTTGCACCTACTTCTGTGGGAGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                                                                                  <  1:1034943/149‑1 (MQ=255)
                                                     gacttGTTTTTGGGGGCGTGAGATTCCGGGGGCTGCGGCGCCGCCCCTTCGCGGGCACGCCCGCCCCCACGGGGGAGGGTGTTTTGCACCTACTTCGGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                                                                                  <  1:1610223/147‑1 (MQ=255)
                                                        ttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGc                                                                                               <  2:5263186/150‑1 (MQ=255)
                                                        ttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGc                                                                                               <  2:6237826/150‑1 (MQ=255)
                                                        ttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCACCCACAGGGGATTGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCCCCCGCGAAGGGCCCCCCAGCGGCCCCAGGGGCGTTAAAGGc                                                                                               >  1:3403510/1‑150 (MQ=255)
                                                          gTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCCGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGTTGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCtt                                                                                             <  2:2056046/150‑1 (MQ=255)
                                                           tttttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCCCACCCACACGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCCGCGGCCCCCGGGGCGGTACACGCTTg                                                                                            >  2:6332619/1‑150 (MQ=255)
                                                           tttttGTGGTCGTGAGATGCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTg                                                                                            <  2:2723270/150‑1 (MQ=255)
                                                           tttgtGTGGTCGTGAGATTCTGGGGGCTGCTGCGCCGCCCCTTCGCGGGCACGCCCGCACCCACAGGGGAGGGTGTTTTGCACCTACTTCTGGGGGCGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTg                                                                                            <  1:4907294/150‑1 (MQ=255)
                                                              ttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                                                                         <  2:4272129/150‑1 (MQ=255)
                                                              tggtgGTCGGGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCCGGGGTTGGTGTTTTGCCCCTACTTCGGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                                                                         <  1:3889394/148‑1 (MQ=255)
                                                                 tgGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTa                                                                                      <  2:2712635/150‑1 (MQ=255)
                                                                      gTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCCCTCCCACAGGGGATGGTGTTTTGCACCTAATTCGGTGGGTGCGGGCGCGCCCCCGAAGGGCCGCCCACCGGCCCCAGGGGCGTTAAAGGCTTGGCCGTACAccc                                                                                 >  1:5005289/1‑150 (MQ=255)
                                                                          ccgcgaaggggCTGCGCAGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTGTTGCACCTCCTTCGGGGGGAGCGGGCGCGCCCGCGGCGGGCCGCCCCGCGGCCCCAGGGGCGTTACCGGCTGGGCAGTACACccgccg                                                                             >  2:4779879/8‑150 (MQ=255)
                                                                          ccgcgaaggggCTGCGCAGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccgccg                                                                             >  2:6075968/8‑150 (MQ=39)
                                                                          ccgcgaaggggCTGCGCAGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCGACGTCGGTGGGAGCGGGCGCGCCCGCGCAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCCGTACACccgccg                                                                             >  1:5029865/8‑150 (MQ=255)
                                                                           aTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGGGATGGTGTTTTGCACCTACGTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCACCGGCCCCAGGGGCGGTACAGGCTTGGCAGTACACCCGCCgg                                                                            >  2:1442999/1‑150 (MQ=255)
                                                                           aTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCCGGGGCGTTACAGGCTTGGCAGTACACCCGCCgg                                                                            >  2:3952717/1‑150 (MQ=255)
                                                                           aTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCGACTTCTGGGGGAGCGGGCGCGCCCGCCAAGGGCCGCCCAGCGGCCCCAGGGGCGGTACAGGCTTGGCAGTACACccgcccg                                                                            >  2:5891059/1‑148 (MQ=255)
                                                                                 ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGCTGGTGGTTTGCCCCGACTTCTGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCaa                                                                      >  1:4561511/1‑150 (MQ=255)
                                                                                 ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATTGTGTTTTGCCCCTACGTGTGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCCGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGgcagcac                                                                      >  2:3494753/1‑149 (MQ=39)
                                                                                  gggCTGCGCCGCCGCCCCCGGGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAg                                                                     <  2:5012743/150‑1 (MQ=39)
                                                                                    gctgctGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGGGGTGGTGGTTTGCACCGACTTCGGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCagccagcc                                                                   >  1:5574801/1‑150 (MQ=255)
                                                                                     cggctgcCCCGCCCCTGCGCGGGCACGCCCGCTCCCACCGGGGTTGGTGTTTTGCACCTACTTCTGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCg                                                                  <  1:3616967/148‑1 (MQ=255)
                                                                                       gctgcGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCGGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGt                                                                >  2:4514328/1‑150 (MQ=255)
                                                                                         tgcGCAGCCCCTTCGCGGGCACGCCCGCGCCCACTGGGGATGGTGTTTTGCACCTACGTCGGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCCGTACACCCGCCGGCCGCAAGCCGGTcc                                                              >  2:3127830/1‑149 (MQ=255)
                                                                                            gcAGCCCCTTCGCGGGCACCCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTAAACCCGCCGGCAGCACGCCGGTCAGAc                                                           >  2:3854912/1‑150 (MQ=255)
                                                                                                 cccTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCTGTGGGTGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTAccc                                                      <  2:849285/150‑1 (MQ=255)
                                                                                                   cTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCtt                                                    <  2:3793388/150‑1 (MQ=255)
                                                                                                   cTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACGTGTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCCCCCGCAGCAAGCCGGTCAGACGAcccat                                                    >  1:2496961/1‑148 (MQ=255)
                                                                                                    ttCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTg                                                   >  2:3218649/1‑150 (MQ=255)
                                                                                                     tCGCGCGCACGCCCGCTCCCACAGGGGATGGTGTTTGGCACCTACTTCGGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTgg                                                  <  2:1478417/150‑1 (MQ=255)
                                                                                                      cgggggCACGCCCGCTCCCCCCGGGGTTGGTTTTGTGCACCTACTAGGGTGGGAGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGc                                                 <  1:3699140/147‑1 (MQ=255)
                                                                                                         gggCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAgg                                                                                                <  2:3273032/100‑1 (MQ=255)
                                                                                                         gggCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAgg                                                                                                >  1:3273032/1‑100 (MQ=255)
                                                                                                          ggCACGCCCGCTCCCACGGGGGATGGTGTTTGGCACCTACTTCTGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAAc                                             <  2:1357238/150‑1 (MQ=255)
                                                                                                            cccGCCCGCTCCCCCCGGGGTGGGTTTTTTGCACCTACTTCTGTGGGTGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAg                                           <  1:674911/148‑1 (MQ=255)
                                                                                                               gCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTc                                        >  2:5601371/1‑150 (MQ=255)
                                                                                                               gCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTc                                        >  2:3613989/1‑150 (MQ=255)
                                                                                                                cccGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTcc                                       >  2:5851450/1‑150 (MQ=255)
                                                                                                                 ccGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCCGACTACCCTTGGCCAACAGGTCCa                                      >  1:2709106/1‑150 (MQ=255)
                                                                                                                 ccGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCa                                      >  1:5641356/1‑150 (MQ=255)
                                                                                                                   gcccccACTGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACg                                    <  1:994932/147‑1 (MQ=255)
                                                                                                                   gCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACg                                    >  1:5353557/1‑150 (MQ=255)
                                                                                                                    cTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACgg                                   >  1:6041233/1‑150 (MQ=255)
                                                                                                                     tCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGc                                  >  1:4544726/1‑150 (MQ=255)
                                                                                                                       ccACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCt                                >  1:3457018/1‑150 (MQ=255)
                                                                                                                        cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt                               >  1:1344526/1‑150 (MQ=255)
                                                                                                                        cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt                               >  1:1965750/1‑150 (MQ=255)
                                                                                                                         acaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCATGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCCGACTACCCTTGGCCAACAGGTCCACGGCCttt                              >  1:3142187/1‑150 (MQ=255)
                                                                                                                         acaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCttt                              >  2:6581589/1‑150 (MQ=255)
                                                                                                                          caGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtttt                             <  2:6709263/150‑1 (MQ=255)
                                                                                                                            ggggATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCg                           <  1:5083445/150‑1 (MQ=255)
                                                                                                                             gggATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGa                          <  1:1652709/150‑1 (MQ=255)
                                                                                                                               gttgggttttgGCACCTAGTTGGGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGata                        <  1:575048/143‑1 (MQ=255)
                                                                                                                               gATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGata                        <  2:6630538/150‑1 (MQ=255)
                                                                                                                                aTGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGatat                       >  2:5490785/1‑150 (MQ=255)
                                                                                                                                  ggTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCCGACTACCCTTGGCCAACAGGTCCACGGCCTGGTCGataggc                     >  2:2924045/1‑147 (MQ=255)
                                                                                                                                   gtgtTTTGCACCTACTTCTGTGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCgg                   >  2:5358597/1‑150 (MQ=255)
                                                                                                                                    tgtTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCgg                   >  2:5117940/1‑150 (MQ=255)
                                                                                                                                    tgtTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCgg                   <  1:5442213/150‑1 (MQ=255)
                                                                                                                                     gtTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGgc                  >  2:6245169/1‑150 (MQ=255)
                                                                                                                                       tttGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGgcgc                >  2:6008809/1‑150 (MQ=255)
                                                                                                                                       tttGCACCTACTGCTGGGGGAGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGgcgc                <  1:4008735/150‑1 (MQ=255)
                                                                                                                                          gCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCgaa             >  1:2116101/1‑150 (MQ=255)
                                                                                                                                          gCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCAAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCgaa             >  1:1418902/1‑150 (MQ=255)
                                                                                                                                            aCCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCgaaga           >  1:4086552/1‑150 (MQ=255)
                                                                                                                                            aCCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCgaaga           >  2:2218620/1‑150 (MQ=255)
                                                                                                                                               tagttCTGGGGGCGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAac        <  1:5988543/147‑1 (MQ=255)
                                                                                                                                                aCTTCTGTGGGCGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAacg       <  1:2348982/150‑1 (MQ=255)
                                                                                                                                                aCTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAacg       <  2:4132164/150‑1 (MQ=255)
                                                                                                                                                  ttCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAacggt     <  2:5866422/150‑1 (MQ=255)
                                                                                                                                                  ttCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAacggt     >  2:1200879/1‑150 (MQ=255)
                                                                                                                                                  ttCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAacggt     <  2:1185559/150‑1 (MQ=255)
                                                                                                                                                  ttCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAacggt     >  2:4664343/1‑150 (MQ=255)
                                                                                                                                                  ttCTGTGGGAGCGGGCGCGCACGCGAAGCGCCGCCCAGCGGCCCCATGGGCGTTACAGGCTTGGCAGTAAACACGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAacggt     >  2:1200905/1‑150 (MQ=255)
                                                                                                                                                   ccTGTGGGAGCGGGCGTGCCCGCGAAGGGGCTGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAacggt     >  1:1663685/2‑150 (MQ=255)
                                                                                                                                                    cTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAacggtca   <  1:6797297/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATATCCGGCGCGAAGAAACGGTCAC  >  NC_002947/5738938‑5739236

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: