Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,460,345 C→T A431V (GCC→GTC)  paaH → 3‑hydroxyadipyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,460,3450CT98.3% 193.4 / NA 59A431V (GCC→GTC) paaH3‑hydroxyadipyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/27);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGA  >  NC_000913/1460099‑1460585
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
gCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACGCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTgg                                                                                                                                                                                                                                              <  2:82671/251‑1 (MQ=255)
 cAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTggg                                                                                                                                                                                                                                             >  1:43894/1‑251 (MQ=255)
                 gctgctgcAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAACGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTtggctg                                                                                                                                                                                                                             <  1:58766/251‑1 (MQ=255)
                        cAGCGGAGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCCGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACAGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTGGTATGGGGAGCGCAACTTGGCTGGCAGCg                                                                                                                                                                                                                       >  2:549804/1‑250 (MQ=255)
                            ggTGAACCCGGACTCAGCGACCCGCAAGGCCAGTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGTCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTAGCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAatat                                                                                                                                                                                                                  <  1:254401/251‑1 (MQ=255)
                               gAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTa                                                                                                                                                                                                                >  1:330864/1‑250 (MQ=255)
                                aaCCCGGACTCAGCGCCCCGCAAGGCCATTTATTAACTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATTCTGATTTGGCGAACGGTAGCTATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGAGATCGATTCCGCCATGCGTCTTGGGGTGAATTCTCCATATGGCCCACTGGTCGGGGGAGCGCAACTTGGCTGTCAGCGAATATTAAg                                                                                                                                                                                                              <  2:511909/251‑1 (MQ=255)
                                   ccGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCt                                                                                                                                                                                                           <  1:96393/251‑1 (MQ=255)
                                     ggACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGAGACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTcc                                                                                                                                                                                                         <  1:551464/251‑1 (MQ=255)
                                      gACTCAGCGACCCGCAAGACCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTCGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCGGGCAGCGAATATTAGGGCTCCt                                                                                                                                                                                                        >  1:215985/1‑251 (MQ=255)
                                         acAGCGACCCGCAGGGCCAGTTATAACATGCAACAGCAGGGCAAAACAGTGCTGCAAATTTCAGATTACCCAGGAATGCTTATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCTTTGATGCGCTTCAAAAAGGCGTGGCCTCTTAACAGGATAGCGATAGCGCCATGCGTCTTGGTGTGAATTATCCATATGGCCCACTTGGCTGGGGAGCGCAACTTGTCTGGCAGCGAATATTAAGGCTCCTTGa                                                                                                                                                                                                     <  2:239292/250‑1 (MQ=255)
                                             cGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAAt                                                                                                                                                                                                 <  1:34677/251‑1 (MQ=255)
                                                  cGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCtcc                                                                                                                                                                                             >  2:585202/1‑248 (MQ=255)
                                                   gcaaggtaaTTTATTACATGCACCCGCAGGGCAAAACAGTGCTGCATATTGCATGTTAACCCGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCACTGAAGCCCTTGATGAGCTTCACAAAGGCGTGGCCTCTGAACAGGATATCGGTACCGCAATGCGTCCTGGGGTGAATTATACACATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAAAATTAAGGCTCCCTGAGAATcaaaac                                                                                                                                                                                           >  2:530935/9‑246 (MQ=255)
                                                         ccATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCAc                                                                                                                                                                                     >  1:583731/1‑251 (MQ=255)
                                                         ccATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCAc                                                                                                                                                                                     <  1:426760/251‑1 (MQ=255)
                                                          cATTTATTACCTGCAACAGCATGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACt                                                                                                                                                                                    >  1:246915/1‑251 (MQ=255)
                                                            tttATTACCTGAAAAAGCAGGGCAAAAAAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTAt                                                                                                                                                                                  >  2:504106/1‑251 (MQ=255)
                                                               attaCCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGc                                                                                                                                                                               >  2:543983/1‑251 (MQ=255)
                                                                   ccggCAACAGCAGGACAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCAGCAAGGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCgaag                                                                                                                                                                           <  2:446789/248‑1 (MQ=255)
                                                                   ccTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCga                                                                                                                                                                             >  2:124397/1‑249 (MQ=255)
                                                                              agGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGc                                                                                                                                                                >  1:457929/1‑251 (MQ=255)
                                                                                      aCAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTc                                                                                                                                                        <  1:249300/251‑1 (MQ=255)
                                                                                                aaaTTGCAGATTACACAGGAATGCTGATTTGGCGAACGGTAGCAATGCTCATCAATGAAGCCCTTGATGCGCTTCAAACAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGAAACTTGGCTGGCAGCGAATATTACGGCTCCTTGAAAATCTACAGCATCACTCTGGCGAAGAACGATATCGCCCATGTTCATTGCTgcgc                                                                                                                                              >  2:527249/1‑251 (MQ=255)
                                                                                                 aaTTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGcc                                                                                                                                             <  2:43650/251‑1 (MQ=255)
                                                                                                  aTTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATc                                                                                                                                                                  >  1:200086/1‑229 (MQ=255)
                                                                                                  aTTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATc                                                                                                                                                                  <  2:200086/229‑1 (MQ=255)
                                                                                                      cAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACgg                                                                                                                                        >  1:5912/1‑251 (MQ=255)
                                                                                                       aGATTACCCAGGAATGCTGATTTGGCGCACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGCTATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACggg                                                                                                                                       >  2:36014/1‑251 (MQ=255)
                                                                                                       aGATTACCCAGGAATGCTGAGTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTGGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGGGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACggg                                                                                                                                       <  2:567547/251‑1 (MQ=255)
                                                                                                        gATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGgc                                                                                                                                      <  1:549804/251‑1 (MQ=255)
                                                                                                                        tGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCCTATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAACATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGt                                                                                                                      >  1:342766/1‑251 (MQ=255)
                                                                                                                          aTTTGGGGAACGGTAGCAATGATAATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCGATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTa                                                                                                                    <  2:43894/251‑1 (MQ=255)
                                                                                                                                      gTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAgg                                                                                                        <  1:182946/251‑1 (MQ=255)
                                                                                                                                         gCAATGATCACCAATGAATCCCTGGATGCGCTTCAAAAAGGCGTGTCCGCTGAACAGGATATTGATACCGCCATGCGTCTTGGGTTGGATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAAGATTAAGGCGCCGTGAAAATCTACATCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCCTAAGGCCt                                                                                                     <  2:55698/251‑1 (MQ=255)
                                                                                                                                                 catcaATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAAt                                                                                             >  1:497/1‑251 (MQ=255)
                                                                                                                                                 catcaATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAAt                                                                                             <  1:543983/251‑1 (MQ=255)
                                                                                                                                                  atcaATGAAGCCCTTGATGCGCTCCAAAAAGGCGTGGCCTTGGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGTCAAAATg                                                                                            <  2:85248/251‑1 (MQ=255)
                                                                                                                                                        gAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCAt                                                                                       <  2:146383/250‑1 (MQ=255)
                                                                                                                                                         aaGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGc                                                                                     >  1:336794/1‑251 (MQ=255)
                                                                                                                                                                     cgcTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATgaga                                                                          >  1:508255/1‑250 (MQ=255)
                                                                                                                                                                          cAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCTCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATg                                                                    <  1:45776/251‑1 (MQ=255)
                                                                                                                                                                          cAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGAt                                                                     <  1:325412/250‑1 (MQ=255)
                                                                                                                                                                           aaaaaGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGc                                                                   >  1:205086/1‑251 (MQ=255)
                                                                                                                                                                                            gAACAGGATATCGATACCGCCATGCGTCTTGGGGTGACTTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGTCGAATAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGATCGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCtt                                                   <  1:563153/250‑1 (MQ=255)
                                                                                                                                                                                            gAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTg                                                  >  2:85895/1‑251 (MQ=255)
                                                                                                                                                                                            gAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTg                                                  >  1:542047/1‑251 (MQ=255)
                                                                                                                                                                                                     atCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCCACGGGCGATTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGAc                                          >  1:16475/1‑250 (MQ=255)
                                                                                                                                                                                                        gATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACAtt                                       >  2:94289/1‑250 (MQ=255)
                                                                                                                                                                                                                   gCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGAt                           >  2:458575/1‑251 (MQ=255)
                                                                                                                                                                                                                     gTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCCTGTTCATTGCTGCGACAACGGGAGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGCGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATg                          >  1:383331/1‑250 (MQ=255)
                                                                                                                                                                                                                     gTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATg                          >  2:159766/1‑250 (MQ=255)
                                                                                                                                                                                                                         tGGGGTGAATTATCCACATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAg                       <  1:13320/249‑1 (MQ=255)
                                                                                                                                                                                                                            ggTGAATTATCCATATGTCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACCGCATCACTATGGCGAAGAACGCTATCGACCATGATCATTGCTGCGCCACCGGGCGCTTCTGGAGAGCGGTTCTGAGTCATAAGGCCTGACCACATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGACTCGCCATTATCTCAATGGCTGCAGGCttt                  >  2:469443/1‑251 (MQ=255)
                                                                                                                                                                                                                              tgcaTTATCCATAGGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCCAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCACTGGATGATGGCTTTGc                >  2:283024/4‑251 (MQ=255)
                                                                                                                                                                                                                                          atGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCAt    <  1:94289/251‑1 (MQ=255)
                                                                                                                                                                                                                                          atGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCAt    <  1:124397/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aCTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATCAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGa  >  1:99347/1‑246 (MQ=255)
                                                                                                                                                                                                                                                 aCTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGa  <  2:99347/246‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
GCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGA  >  NC_000913/1460099‑1460585

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: