Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,739,061 A→C 11.7% intergenic (‑82/+75) 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,0610AC11.7% 528.9 / 4.4 232intergenic (‑82/+75)proY/hutHproline (histidine) APC transporter/histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (86/101);  new base C (16/16);  total (104/128)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

GAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTC  >  NC_002947/5738912‑5739210
                                                                                                                                                     |                                                                                                                                                     
gAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCcaca                                                                                                                                                       <  1:228067/150‑1 (MQ=255)
 aGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCATGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                                                      <  2:4352226/150‑1 (MQ=255)
   ccGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTGTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAggg                                                                                                                                                    <  1:2839602/150‑1 (MQ=255)
     gtgggCCATGAAACGGATCTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTCTTGCATGTCACGTCCTTGATTTCTTTATGTGGTCGTGAGATTCAGGGGGCTGCAGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCTCAgggca                                                                                                                                                  <  1:556604/148‑3 (MQ=255)
     gAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGa                                                                                                                                                  <  1:3570887/150‑1 (MQ=255)
     gAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGa                                                                                                                                                  >  2:674911/1‑150 (MQ=255)
      aGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGAt                                                                                                                                                 <  1:5502078/150‑1 (MQ=255)
      aGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGAt                                                                                                                                                 <  1:3619368/150‑1 (MQ=255)
        ggCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAgaagtagg                                                                                                                                               <  1:1404448/150‑8 (MQ=255)
        ggCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATgg                                                                                                                                               <  1:2723878/150‑1 (MQ=255)
        ggCCATCAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTCTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACGGGGGATgg                                                                                                                                               >  2:5988739/1‑150 (MQ=255)
         gCCATGGAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgt                                                                                                                                              <  1:5404637/150‑1 (MQ=255)
             tGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGttt                                                                                                                                          >  2:3253905/1‑150 (MQ=255)
                aaCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCCCGCCCCCTACCACACGGGAGGGTGTTTTgg                                                                                                                                       >  1:5866422/1‑149 (MQ=255)
                 aCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCGCCCACCGGGGCTGGTGTTTTGcc                                                                                                                                      >  1:6630538/1‑149 (MQ=255)
                 aCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCCTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCa                                                                                                                                      <  2:4930734/150‑1 (MQ=255)
                  cGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCAc                                                                                                                                     <  1:1280900/150‑1 (MQ=255)
                     aTGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTa                                                                                                                                  <  2:3493724/150‑1 (MQ=255)
                            gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                           <  2:2911863/150‑1 (MQ=255)
                            gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                           <  2:3088696/150‑1 (MQ=255)
                            gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                           <  2:5451909/150‑1 (MQ=255)
                            gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                           <  2:6244065/150‑1 (MQ=255)
                            gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGCCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                           <  2:1454001/150‑1 (MQ=255)
                             ggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgtg                                                                                                                          <  1:945450/150‑1 (MQ=255)
                             ggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgtg                                                                                                                          <  1:6706037/150‑1 (MQ=255)
                              gCACTTAGCCCGCGCTTTAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTgg                                                                                                                         <  2:4247486/150‑1 (MQ=255)
                              gCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTgg                                                                                                                         <  1:5024674/150‑1 (MQ=255)
                                aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTTAATTGTTTTTGTGGTCGTGAGATTCTGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                                                       <  1:3558780/150‑1 (MQ=255)
                                aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                                                       <  2:1526952/150‑1 (MQ=255)
                                  ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGGACGCCCCCTCCCCCACGGGATGGTGTTTTGCACCTACGTCTGTGGGAGc                                                                                                                     >  1:4727743/1‑150 (MQ=255)
                                  ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTTCACCTACTTCTGTGGGAGc                                                                                                                     >  2:2348982/1‑150 (MQ=255)
                                  ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGc                                                                                                                     >  2:5765877/1‑150 (MQ=255)
                                   tcgCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCg                                                                                                                    <  2:6016864/148‑1 (MQ=255)
                                    aGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCgg                                                                                                                   <  2:583251/150‑1 (MQ=255)
                                    aGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCgg                                                                                                                   <  1:1711679/150‑1 (MQ=255)
                                      cccGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgc                                                                                                                 <  1:2604374/150‑1 (MQ=255)
                                      cccGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgc                                                                                                                 <  2:2055322/150‑1 (MQ=255)
                                      cccGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACGTCTGTGGGAGCGGgc                                                                                                                 >  2:506057/1‑150 (MQ=255)
                                         gcgcTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCCCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgcgcg                                                                                                              >  1:6086758/1‑150 (MQ=255)
                                          cgcTTGCGACCGTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTGGGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTGTGCACCTACTTCTGTGGGAGCGGgcgcgc                                                                                                             <  2:4642129/150‑1 (MQ=39)
                                           gcTCGAGACCTCGAGCTTGTGGCATGTCACGTCCTTGAAGTGTTTTTGTGGTCCTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGcc                                                                                                            <  1:4440613/150‑1 (MQ=39)
                                            gttGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTGGTGTGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGTTGGTGTTTTGCACCTACTTCTGTGGGAGGGGGCGCGccc                                                                                                           <  1:5036615/149‑1 (MQ=255)
                                             ttGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACCCCCGCTCCCACACGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCCGGCGCGCCcg                                                                                                          >  2:1357663/1‑150 (MQ=255)
                                                agaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTTCACCTACTTCTGTGGGAGCGGGCGCGCCcgcca                                                                                                       >  2:2025000/1‑148 (MQ=255)
                                                 gaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTACCACAGGGGAAGGTGTTTTGCACCTACATCTGTGGGAGCGGGCGCGCCCGCGaa                                                                                                      >  2:2244531/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)
                                                                  aTGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGGACCCCCCCTCCCACACGGGCTGGTGTTTTTCCCCTACGACGGGGGGAGCGGGCCCGCCCCCCCAAGGCCCCCCCCCGGccc                                                                                     >  1:2480938/1‑150 (MQ=255)
                                                                         gTCCTTGAATTGTTTTTGGGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCCGGGGTTGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGc                                                                              <  1:6771439/150‑1 (MQ=255)
                                                                               taaTTTTTTTTGTGGTCGTGTGATTCTGGGGGCGGCTGCGCAGCCCCTTCGGGGGCACGCCCGCTCCCACTGGGTTTGGTGTTTTGCACCTACTTCTGTGGGAGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                                                        <  1:1034943/149‑1 (MQ=255)
                                                                               gacttGTTTTTGGGGGCGTGAGATTCCGGGGGCTGCGGCGCCGCCCCTTCGCGGGCACGCCCGCCCCCACGGGGGAGGGTGTTTTGCACCTACTTCGGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                                                        <  1:1610223/147‑1 (MQ=255)
                                                                               gAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCCCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                                                        <  2:1864213/150‑1 (MQ=255)
                                                                                aaTTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCACGGGCGTTACAg                                                                       >  2:6653343/1‑150 (MQ=255)
                                                                                 aTTGTTTTTGTGGTCGTGAGTTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAgg                                                                      <  2:3585489/150‑1 (MQ=255)
                                                                                  ttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGc                                                                     <  2:5263186/150‑1 (MQ=255)
                                                                                  ttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCACCCACAGGGGATTGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCCCCCGCGAAGGGCCCCCCAGCGGCCCCAGGGGCGTTAAAGGc                                                                     >  1:3403510/1‑150 (MQ=255)
                                                                                   tGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCt                                                                    >  2:2083791/1‑150 (MQ=255)
                                                                                       ttggtgGTCGTGAGCTTCCGGGGGCGGCTGCGCAGCCCCTTCGCGGGCCCGCCCGCTCCCCCCGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGc                                                                <  2:5927414/147‑1 (MQ=255)
                                                                                        ttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                                               <  2:4791575/150‑1 (MQ=255)
                                                                                        tggtgGTCGGGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCCGGGGTTGGTGTTTTGCCCCTACTTCGGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                                               <  1:3889394/148‑1 (MQ=255)
                                                                                        ggggggtggtGAGTTTCTGGGGGGGGGGGCGCCGCCCGTGCGGGGGCCCGCCCGCGCCCCCGGGGGTTGGTGTTTGCCCCCTCCTGGGGGGGGCGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                                               <  1:6428865/146‑1 (MQ=255)
                                                                                           tgGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTa                                                            <  2:2712635/150‑1 (MQ=255)
                                                                                              tCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTAc                                                                                                                                 <  2:1753550/78‑1 (MQ=255)
                                                                                                    ccgcgaaggggCTGCGCAGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccgccg                                                   >  2:6075968/8‑150 (MQ=39)
                                                                                                     aTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCCGGGGCGTTACAGGCTTGGCAGTACACCCGCCgg                                                  >  2:3952717/1‑150 (MQ=255)
                                                                                                      ttCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCGCCCACAGGGGATGGTGTTTTGCACCAACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCCGCGGCCCCACGGGCGTTACAGGCTTGGCAGTACACCCGCCGgc                                                 >  1:4988314/1‑150 (MQ=255)
                                                                                                       tCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGgca                                                <  2:6393386/150‑1 (MQ=255)
                                                                                                        gaaggggCTGCGCAGCAGCCCCCTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCa                                                                                                                                      <  1:5566197/60‑1 (MQ=18)
                                                                                                        gaaggggCTGCGCAGCAGCCCCCTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCa                                                                                                                                      >  2:5566197/4‑63 (MQ=18)
                                                                                                           ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGCTGGTGGTTTGCCCCGACTTCTGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCaa                                            >  1:4561511/1‑150 (MQ=255)
                                                                                                           ggggCTGCGCAGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCaa                                            >  2:3817090/1‑150 (MQ=255)
                                                                                                                 gctgcGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCGGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGt                                      >  2:4514328/1‑150 (MQ=255)
                                                                                                                      gcAGCCCCTTCGCGGGCACCCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTAAACCCGCCGGCAGCACGCCGGTCAGAc                                 >  2:3854912/1‑150 (MQ=255)
                                                                                                                             cTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCtt                          <  2:3793388/150‑1 (MQ=255)
                                                                                                                              ttCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTg                         >  2:3218649/1‑150 (MQ=255)
                                                                                                                                cgggggCACGCCCGCTCCCCCCGGGGTTGGTTTTGTGCACCTACTAGGGTGGGAGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGc                       <  1:3699140/147‑1 (MQ=255)
                                                                                                                                   gggCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTCCTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACCGGCTTGGCCGTCCACCCGCCGGCAGCACGCCGGTCAGACTACCCTGGCCacc                   >  2:5376554/1‑148 (MQ=39)
                                                                                                                                   gggCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                           <  1:4078286/47‑1 (MQ=255)
                                                                                                                                       aCGCCCGCTCCCACTGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAgg                <  2:1065484/150‑1 (MQ=255)
                                                                                                                                       aCGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAgg                <  2:5693788/150‑1 (MQ=255)
                                                                                                                                         gCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTc              >  2:3613989/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCCACAGGGGATGGTGTTTTGCACCTACGTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCAGCCGGCAGCAAGCCGGTCAGACAACCCGTGGCCAACAGGTcc             >  1:5989927/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCCGACTACCCTTGGCCAACAGGTCCa            >  1:2709106/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCa            >  1:5641356/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCa            >  1:2063783/1‑150 (MQ=255)
                                                                                                                                             gcccccACTGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACg          <  1:994932/147‑1 (MQ=255)
                                                                                                                                             gCTCCCACTGGGGATGGTGTTTTGCACCTACTTCGGTGGGTGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACg          <  1:3253905/150‑1 (MQ=255)
                                                                                                                                             gCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACg          <  2:1281066/150‑1 (MQ=255)
                                                                                                                                                 ccACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCt      >  1:3457018/1‑150 (MQ=255)
                                                                                                                                                 ccACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCt      >  1:533922/1‑150 (MQ=255)
                                                                                                                                                  cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt     >  1:1352790/1‑150 (MQ=255)
                                                                                                                                                  cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt     >  2:86661/1‑150 (MQ=255)
                                                                                                                                                  cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt     >  2:3359263/1‑150 (MQ=255)
                                                                                                                                                  cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt     >  2:2569895/1‑150 (MQ=255)
                                                                                                                                                  cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt     >  2:134712/1‑150 (MQ=255)
                                                                                                                                                   acaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCATGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCCGACTACCCTTGGCCAACAGGTCCACGGCCttt    >  1:3142187/1‑150 (MQ=255)
                                                                                                                                                   acaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCttt    >  2:6581589/1‑150 (MQ=255)
                                                                                                                                                    caGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCACAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtttt   >  2:5857808/1‑150 (MQ=255)
                                                                                                                                                     aGGGGTTGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGACAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTc  <  2:5096422/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTC  >  NC_002947/5738912‑5739210

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: