Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,592,468 C→A 5.6% W98C (TGG→TGT hcaD ← 3‑phenylpropionate/cinnamic acid dioxygenase ferredoxin‑‑NAD( ) reductase component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,592,4680CA 5.6% 219.8 / 2.9 72W98C (TGG→TGThcaD3‑phenylpropionate/cinnamic acid dioxygenase ferredoxin‑‑NAD( ) reductase component
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base A (2/2);  total (33/39)
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.21e-01

CCCAGCGAAGCCATGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAACCGCCGCAGACTA  >  KT2440/3592339‑3592601
                                                                                                                                 |                                                                                                                                     
cccAGCGAAGCCATGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCCCTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGAc                                                                                                                                  >  6:213948/1‑135 (MQ=255)
cccAGCGAAGCCATGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGAc                                                                                                                                  >  4:325893/1‑135 (MQ=255)
  ccgCGAAGCCATGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACtt                                                                                                                                <  1:40232/133‑1 (MQ=255)
  cAGCGACGCCATGCGCCAGCAGGGGGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACtt                                                                                                                                <  1:314929/135‑1 (MQ=255)
  cAGCGAAGCCATGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACtt                                                                                                                                <  3:206324/135‑1 (MQ=255)
   aGCGAAGCCATGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGGTCGGTTGCCGGCAGGGGCTGGTCCTCGACCCAGCCTTc                                                                                                                               >  6:118112/1‑135 (MQ=255)
       aaGCCATGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGTGCGTTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGc                                                                                                                          <  8:121521/136‑1 (MQ=255)
        aGCCATGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGc                                                                                                                          <  4:251919/135‑1 (MQ=255)
             tGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGCCCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTg                                                                                                                     >  3:346839/1‑135 (MQ=255)
              gcgcgagctgGGTGTCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGAt                                                                                                                   <  7:308313/131‑1 (MQ=255)
              gAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGa                                                                                                                    <  4:239327/135‑1 (MQ=255)
              gAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGAt                                                                                                                   >  7:220615/1‑136 (MQ=255)
               cgCCCGCTGGGGGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGAt                                                                                                                   <  5:248976/134‑1 (MQ=255)
                gCCAGCGGGGTGGCGACGGCGGCGATGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCGAGACTTCATCGCCCTTGAt                                                                                                                   <  1:279268/134‑1 (MQ=255)
                gCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGAt                                                                                                                   <  3:50067/134‑1 (MQ=255)
                gCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGAt                                                                                                                   <  6:113641/134‑1 (MQ=255)
                       gggTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGCCCCAGCCGTCATCGCCCTTGATCCACGccg                                                                                                           >  4:60913/1‑134 (MQ=255)
                       gggTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGccat                                                                                                          >  7:400555/1‑134 (MQ=255)
                             cGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCCGACTTCATCGCCCTTGATCCACGCCTTGTAgcg                                                                                                    >  7:321431/1‑136 (MQ=255)
                             cGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAgc                                                                                                     >  6:147856/1‑135 (MQ=255)
                                  gcggcgCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAgcggcg                                                                                                 >  8:135251/1‑134 (MQ=255)
                                     gcgcTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGGTGCCGGCAGGGGCTGGTCCTCGCCCCAGACGTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAgg                                                                                             >  6:165873/1‑135 (MQ=255)
                                                  cACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAg                                                                                <  1:327676/135‑1 (MQ=255)
                                                   aCCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAggc                                                                              >  8:9678/1‑136 (MQ=255)
                                                      ggCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACACAGACTTCATCGCCCTTGCTCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAggcgg                                                                            >  2:351598/1‑135 (MQ=255)
                                                      ggCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACACAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAggcgg                                                                            >  2:345594/1‑135 (MQ=255)
                                                         gcCGACCACTACGAAGCAGCGGGCATCGCTGTGACTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTc                                                                         <  1:125616/135‑1 (MQ=255)
                                                          ccGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTca                                                                        >  4:235608/1‑135 (MQ=255)
                                                           cGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCTACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTcc                                                                        >  6:423457/1‑133 (MQ=255)
                                                           cGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTcac                                                                       <  3:170768/135‑1 (MQ=255)
                                                           cGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTcac                                                                       <  1:267086/135‑1 (MQ=255)
                                                                          agcGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAggcggc                                                                           <  3:320316/116‑1 (MQ=255)
                                                                          agcGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAggcggc                                                                           >  4:320316/1‑116 (MQ=255)
                                                                          agcGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACcgtcgt                                                        >  8:313374/1‑135 (MQ=255)
                                                                              ggCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAAc                                                    <  2:271700/135‑1 (MQ=255)
                                                                              ggCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAAc                                                    <  2:148595/135‑1 (MQ=255)
                                                                               gCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACgg                                                  >  8:170944/1‑136 (MQ=255)
                                                                               gCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACg                                                   >  3:109409/1‑135 (MQ=255)
                                                                               gCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACg                                                   >  4:180948/1‑135 (MQ=255)
                                                                               gCATAGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGCGCTGGTCTTCGACCCAGACTTCATCGCGCGGGATCCACGGCGTGTAGCGGCGCAGGTCGTCGCGTGCAGGCGGCTCACAACCCCCACCGTCGTAAAccc                                                  >  8:434239/1‑134 (MQ=255)
                                                                                cATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGc                                                 <  8:220615/136‑1 (MQ=255)
                                                                                   cGCTGTGCCTGGGCGGTTCGGGTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGa                                               <  3:383252/135‑1 (MQ=255)
                                                                                     cTGTGCCGGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGa                                               <  5:277331/133‑1 (MQ=255)
                                                                                       gtgCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAgcg                                                                                                    <  1:288145/78‑1 (MQ=255)
                                                                                       gtgCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAgcg                                                                                                    >  2:288145/1‑78 (MQ=255)
                                                                                        tgCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAg                                         <  7:170944/136‑1 (MQ=255)
                                                                                           cTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCACAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCt                                       >  3:336127/1‑135 (MQ=255)
                                                                                            tGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACCAACCCCACCGTCGTCAACGGCGAACGTAGCtt                                      >  1:183998/1‑135 (MQ=255)
                                                                                            tGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGCACGTAGCtt                                      >  1:403879/1‑135 (MQ=255)
                                                                                             gggCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGCACGTAGCttt                                     >  3:389780/1‑135 (MQ=255)
                                                                                                cggttcggttGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTtgtg                                  <  5:301620/135‑1 (MQ=255)
                                                                                                     cggttGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAgg                            >  7:388870/1‑136 (MQ=255)
                                                                                                     cggttGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAgg                            <  7:358767/136‑1 (MQ=255)
                                                                                                     cggttGCCGGCAGGGGCTGGTCAGCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAgg                            <  7:266004/136‑1 (MQ=255)
                                                                                                      ggttGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAAAGGCGAACGTAGCTTTGTGCCAAgg                            >  2:277038/1‑135 (MQ=255)
                                                                                                              gCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCa                    <  3:110107/135‑1 (MQ=255)
                                                                                                                 ggggCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCa                 <  2:36062/135‑1 (MQ=255)
                                                                                                                     cTGGTCATCGACCCAGACTTCATCGCCCTTGCTCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGCCCGTCGCGTTGTGGCAAGGGCCGACCAGCAACccg            >  7:165319/1‑136 (MQ=255)
                                                                                                                      tGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAAccg            <  6:394581/135‑1 (MQ=255)
                                                                                                                      tGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAAccg            <  4:194165/135‑1 (MQ=255)
                                                                                                                      tGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAAccg            <  4:346839/135‑1 (MQ=255)
                                                                                                                      tGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGGCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAAccg            <  3:311279/135‑1 (MQ=255)
                                                                                                                      tGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGGCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAAccg            <  4:136239/135‑1 (MQ=255)
                                                                                                                      tGGTCATCGACACAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAAccg            <  1:351598/135‑1 (MQ=255)
                                                                                                                      tGGTCATCGACACAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAAccg            <  1:345594/135‑1 (MQ=255)
                                                                                                                       gggcATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAAccgcc          <  8:116693/133‑1 (MQ=255)
                                                                                                                        gTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAAccgcc          <  7:251761/135‑1 (MQ=255)
                                                                                                                           aTCGACTCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAATCCAGCAAACCGCCGCa       >  5:257920/1‑135 (MQ=255)
                                                                                                                           aTCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGACCGTAGCTTTGTGCCAAGGGCAAACCAGCAAACCGCCGCa       >  2:394850/1‑135 (MQ=255)
                                                                                                                           aTCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGCCGGCTCACAACCCACCCCGTCGTCAACGGCGCACGTAGCTTTGTGCCAAGGGCAACCCCGCAAAccgccg         >  7:272987/1‑133 (MQ=255)
                                                                                                                            tCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAACCGCCGCAGa     <  7:317866/136‑1 (MQ=255)
                                                                                                                             cGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGGCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAACCGCCGCAGa     <  3:56730/135‑1 (MQ=255)
                                                                                                                                cccAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTc                                                                                           <  4:107429/46‑1 (MQ=255)
                                                                                                                                cccAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTc                                                                                           >  3:107429/1‑46 (MQ=255)
                                                                                                                                cccAGACTTCATCGCCCTTCATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCCTACTCAACCGCGAACGTAGCTTAGTGCCAAGCGCAAACCAGCAAACCGCAGCAGACTa  >  3:189951/1‑135 (MQ=255)
                                                                                                                                 ccAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGa                                               >  4:286391/1‑89 (MQ=255)
                                                                                                                                 ccAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGa                                               <  3:286391/89‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                     
CCCAGCGAAGCCATGAGCCAGCAGGGTGGCGACGGCGGCGCTGCCGGCGGCACCGGCGCCGACCACTACGAAGCAGCGGGCATCGCTGTGCCTGGGCGGTTCGGTTGCCGGCAGGGGCTGGTCATCGACCCAGACTTCATCGCCCTTGATCCACGCCTTGTAGCGGCGCAGGTCGGCGAGTGCAGGCGGCTCACAAACCACACCGTCGTCAACGGCGAACGTAGCTTTGTGCCAAGGGCAAACCAGCAAACCGCCGCAGACTA  >  KT2440/3592339‑3592601

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