Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,665,772 G→T 86.0% E123* (GAA→TAA)  hcaD → phenylpropionate dioxygenase, ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,665,7720GT86.0% 498.0 / 67.7 223E123* (GAA→TAA) hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (21/12);  new base T (91/98);  total (113/110)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

CATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAG  >  CP009273/2665625‑2665919
                                                                                                                                                   |                                                                                                                                                   
cATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGGTGCGTGATTCACTGGgaga                                                                                                                                                    >  2:1640925/1‑149 (MQ=255)
 aTCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATaa                                                                                                                                                   <  1:2710664/149‑1 (MQ=255)
  tCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGCGCTCGCCCCCGGCCCGGGCGGGATGCACGGGGGGAAc                                                                                                                                                  >  1:2098806/1‑149 (MQ=255)
    tGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgc                                                                                                                                                <  2:1422991/149‑1 (MQ=255)
       aTTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctgc                                                                                                                                             <  2:888098/149‑1 (MQ=255)
       aTTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCCCAGCCCTTGCTTGTTGCACAGGGAGAAcccagc                                                                                                                                             >  1:3924290/1‑144 (MQ=255)
          ccGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCttt                                                                                                                                          <  1:4099099/149‑1 (MQ=255)
            ggTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTAc                                                                                                                                        <  2:1040258/149‑1 (MQ=255)
             gtgtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTAcc                                                                                                                                       <  2:34288/149‑1 (MQ=255)
                 aaCCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTAc                                                                                                                                   <  1:3801416/149‑1 (MQ=255)
                 aaCCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTAc                                                                                                                                   >  2:3539121/1‑149 (MQ=255)
                   ccATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCACCACCCCCGCCCGTGCTTGCTGCACAGGGATAACCCGGCGTTAaccaaacc                                                                                                                                 >  1:4534321/1‑144 (MQ=255)
                     aTCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCGTGATGCACTGGGATAACGCTGCGTTAACCTACGCCa                                                                                                                               >  2:621778/1‑149 (MQ=255)
                      tCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCAt                                                                                                                              <  1:163226/149‑1 (MQ=255)
                         aaaCGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCAAGcc                                                                                                                           >  2:377413/1‑149 (MQ=255)
                           aCGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCTGCCGGTGGGTGATGCACTGGGAGAACGCTGCTTTACCCAACGCCCTGCCgg                                                                                                                         >  2:3719229/1‑149 (MQ=255)
                             gCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTg                                                                                                                       <  1:1518129/149‑1 (MQ=255)
                             gCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCGACGCCATGCCGGTg                                                                                                                       >  2:3577972/1‑149 (MQ=255)
                                        acacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGa                                                                                                            <  1:4573334/149‑1 (MQ=255)
                                        acacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCAACGCCATGCCGGTGATGCCGCCAGa                                                                                                            >  2:915884/1‑149 (MQ=255)
                                        acacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCTTGCTTGATGCACGGGGAGAAAGCTGCTTTTCACTACGCCCTGCCGGTGATGCCGCCaca                                                                                                            >  2:915837/1‑147 (MQ=255)
                                        acacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGAGACCGCTGCCGTTGGGTGTTGGACAGGGCTAAAGCTGCTTTAACCGACGCCATGCCGGTGATGCCGCCAGa                                                                                                            >  2:3675000/1‑149 (MQ=255)
                                         cacacGAGAGATAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGAc                                                                                                           <  2:2373243/149‑1 (MQ=255)
                                          acacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCGTGATGCACTGGGATAACGCGGCTTTACCCTACGCCAAGCCGGTGATGCCGCCAGACt                                                                                                          >  2:708431/1‑149 (MQ=255)
                                           cacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCGTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTg                                                                                                         >  2:2712080/1‑149 (MQ=255)
                                             cGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTTCTTGATGCACTGGGATAACGCTGCTTTAACCTACGCCATGCCGGTGATGCCGCCAGACTGCg                                                                                                       >  2:4324530/1‑149 (MQ=255)
                                               agagTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCAGCCGTTGCGTGATGCACAGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCCCACAGCgag                                                                                                     >  2:4513061/1‑149 (MQ=255)
                                               agagTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGAACGCGGCCGTTGGTTGATTCAATGGGATAACGCAGCTTTACCCAAAGCCAATCCGGTGATGCCGCCACACAGCgag                                                                                                     >  1:3089769/1‑149 (MQ=255)
                                                gagTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGCGCGCGACCCCGGCCGGTGGGTGTTTGACACGGGTAAAGCTGCGTTACCCCAACCCCAGCCCGTGATGCCCCCACACAGCgaga                                                                                                    >  1:3633968/1‑149 (MQ=255)
                                                  gTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAAAGCTGGGTTACCCTACGCCCTGCCGGTGATGCCGCCCGACAGCGAGAAg                                                                                                  >  2:4257877/1‑149 (MQ=255)
                                                   ttAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGAGCGCGACCCCTTCCCTTGGTTTTTGGACTGGGATAACGCTGCGTTTCCCCAAGCCATGCCGGTGATGCCCGCCGAATGCgaggatt                                                                                                 >  1:2557772/1‑145 (MQ=255)
                                                   ttAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGt                                                                                                 <  2:1271715/149‑1 (MQ=255)
                                                    tAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGtt                                                                                                <  2:3951291/149‑1 (MQ=255)
                                                     aGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTc                                                                                               <  2:3138375/149‑1 (MQ=255)
                                                     aGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTc                                                                                               <  2:3054320/149‑1 (MQ=255)
                                                      gtgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCTCTCGACCCCTGCCGTTGCGTGATTCACGGGGGGAACGCAGCTTTTACCAACGCCCAGCCGGGGATGCCGCCAGACAGCGAGACGTTCt                                                                                              >  2:2046493/1‑149 (MQ=255)
                                                      gtgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCt                                                                                              <  1:3745197/149‑1 (MQ=255)
                                                       tgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATTACGCTGCTTTAACCTAAGCCAAGCCGGTGATGCCGCCAGACAGCGAGAAGTTCTg                                                                                             >  2:2916803/1‑149 (MQ=255)
                                                         ttAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCa                                                                                           <  1:625664/149‑1 (MQ=255)
                                                          tAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTGTTATAGCAACCGGCGCGGCAGCGCGTCCGCGGCCGGTGCGGGTTGCCCAGGGCGAACGCGGCTGTAACCAACGCCAAGCCGGTGCTGCCGCCCGACAGCGAGAACTTCTGGAg                                                                                          >  2:3226304/1‑149 (MQ=255)
                                                           aaCCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCAGCCCTTGCTTGATGGACAGGGATAACGCTGCTTTAACCCAAGCCAAGCCCGTGATGCCGCCAGACAGCGACAAATTCAGCacc                                                                                         >  2:1663885/1‑147 (MQ=255)
                                                           aaCCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGc                                                                                         <  2:145903/149‑1 (MQ=255)
                                                             ccAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGccc                                                                                       <  1:2392085/149‑1 (MQ=255)
                                                                aCGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGaa                                                                                    <  2:1986102/149‑1 (MQ=255)
                                                                 cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAAc                                                                                   <  1:1417027/149‑1 (MQ=255)
                                                                 cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAAc                                                                                   <  2:3347518/149‑1 (MQ=255)
                                                                    cGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGt                                                                                <  2:3219370/149‑1 (MQ=255)
                                                                    cGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGt                                                                                <  1:4312381/149‑1 (MQ=255)
                                                                     gAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCGGCCGGTGGTTGTTTGACAGGGGTTACGCGGCGTTTACCGACGCCCAGCCGGGGGTGCCGCCCGCCAGCCCGGAGTTGAGCAGCCCGCACGGtg                                                                               >  1:2998091/1‑148 (MQ=255)
                                                                           tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCAACGCCAAGCCGGTGATGCCGCCAGACAGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGt                                                                         >  2:310688/1‑149 (MQ=255)
                                                                               aCTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGTTGCACAGGGATAAAGCTGCGTTAACCAACCCCAAGCCGGTGATGCCCCCACACAGCGAGAAGTTATTCAGCCCCAACGGTCAGGCGGGAtt                                                                     >  1:3157815/1‑149 (MQ=255)
                                                                                 tGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCCACCGCTGCCGTTTCTTGATGCACTGGGATAACGCAGCGTTAACCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCCGCCCGAACGGTCAGTCGGGATTAt                                                                   >  2:4678143/1‑149 (MQ=255)
                                                                                  gggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATc                                                                  <  1:4480929/149‑1 (MQ=255)
                                                                                  gggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCTTGCTTGATGCACAGGGATAACGCAGCTTTACCCAACGCCCTGCCGGTGATGCCGCCAGACAGCGAGAAGTTCTGCAGCCCGAACGGTCCGTCGGGATTATc                                                                  >  2:1135666/1‑149 (MQ=255)
                                                                                   ggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCTAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCg                                                                 <  1:644803/149‑1 (MQ=255)
                                                                                   ggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACAGGGATAACGCTGCTTTAACCTACGCCAAGCCGGTGATGCCGCCAGACTGCGAGAAGTTATGCAGCCCCAAACGTCAGTCGTGATTAtag                                                                 >  2:2292534/1‑147 (MQ=255)
                                                                                   ggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCAGCCGTTTCTTGATGCACAGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATTCCGCCAGACAGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCg                                                                 >  2:1203204/1‑149 (MQ=255)
                                                                                    gATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCgg                                                                <  1:1882681/149‑1 (MQ=255)
                                                                                     aTCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGCCCGTGCGTGATTCACAGGGAGAACGCAGCGTTACCCTACGCCAAGCCCGTGATGCCGCCAGACAGCGAGAAGTTATGCAGCCCCAACGGTCAGTCGGGATTATCGGt                                                               >  2:4589106/1‑149 (MQ=255)
                                                                                        aGCTTTTTATAGCAACCGGCGCGGCAGCTCGTCCGCTGCCGTTTCTTGATTCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTATGCAGCCCCAACGGTCAGTCGTGATTATAGGTGcc                                                            >  2:2522200/1‑149 (MQ=255)
                                                                                        aGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCTTGCTTGATTCCCTGGGATAACGCGGCGTTACCCTACGCCATGCCCGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGGCAGTCGGGATTAACGGTGcc                                                            >  2:711985/1‑149 (MQ=255)
                                                                                         gCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCTGCCGTTTCTTGATGCACTGGGGTAACGCTGCTTTTACCAACGCCCAGCCGGTGATTCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATAGGTGCCg                                                           >  2:4083473/1‑149 (MQ=255)
                                                                                             tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAAc                                                       <  2:2784638/149‑1 (MQ=255)
                                                                                              ttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACt                                                      <  2:2641235/149‑1 (MQ=255)
                                                                                                ataGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                    <  1:1654692/149‑1 (MQ=255)
                                                                                                 taGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTTATGCACTGGGATAAAGCTGCTTTACCCAACGCCAAGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAtt                                                   >  2:309221/1‑149 (MQ=255)
                                                                                                 taGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAtt                                                   <  1:3257098/149‑1 (MQ=255)
                                                                                                   gCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTGACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTgg                                                 <  2:2835554/149‑1 (MQ=255)
                                                                                                       ccGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTg                                             <  1:232797/149‑1 (MQ=255)
                                                                                                          gcgcGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGaa                                          <  2:3157815/149‑1 (MQ=255)
                                                                                                          gcgcGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGaa                                          <  2:3698243/149‑1 (MQ=255)
                                                                                                           cgcgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAAc                                         <  1:2827317/149‑1 (MQ=255)
                                                                                                            gcgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACt                                        >  1:3934289/1‑149 (MQ=255)
                                                                                                              ggCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                      >  1:106160/1‑149 (MQ=255)
                                                                                                               gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGc                                     >  1:3922616/1‑149 (MQ=255)
                                                                                                                  gCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                                  >  1:1535047/1‑149 (MQ=255)
                                                                                                                   cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGcc                                 <  1:4434954/149‑1 (MQ=255)
                                                                                                                    tCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCa                                <  1:277778/149‑1 (MQ=255)
                                                                                                                    tCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCa                                >  2:1699384/1‑149 (MQ=255)
                                                                                                                      gACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgc                              <  2:2505737/149‑1 (MQ=255)
                                                                                                                       aCCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgcg                             >  1:3940008/1‑149 (MQ=255)
                                                                                                                             gCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgc                       >  1:898847/1‑149 (MQ=255)
                                                                                                                                gTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagc                    <  1:1391161/149‑1 (MQ=255)
                                                                                                                                 ttGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcg                   >  1:874154/1‑149 (MQ=255)
                                                                                                                                  tGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgc                  >  2:3881641/1‑149 (MQ=255)
                                                                                                                                   gCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGTCAGCGCCAcgcagcgca                 >  2:1782819/1‑149 (MQ=255)
                                                                                                                                   gCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgca                 >  1:4588569/1‑149 (MQ=255)
                                                                                                                                     ttGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa               <  1:4135761/149‑1 (MQ=255)
                                                                                                                                       gATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg             >  1:1473163/1‑149 (MQ=255)
                                                                                                                                       gATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg             <  2:3152259/149‑1 (MQ=255)
                                                                                                                                        aTGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGt            <  1:1926553/149‑1 (MQ=255)
                                                                                                                                         tGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTa           <  2:473045/149‑1 (MQ=255)
                                                                                                                                          gCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTaa          <  2:4080481/149‑1 (MQ=255)
                                                                                                                                           cACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAg         <  1:790494/149‑1 (MQ=255)
                                                                                                                                            aCTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAgg        >  1:479719/1‑149 (MQ=255)
                                                                                                                                            aCTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAgg        >  1:264688/1‑149 (MQ=255)
                                                                                                                                             cTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGt       >  2:4612131/1‑149 (MQ=255)
                                                                                                                                             cTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGt       <  1:1838355/149‑1 (MQ=255)
                                                                                                                                             cTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGt       <  1:21265/149‑1 (MQ=255)
                                                                                                                                                  agaACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAg  >  2:512402/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAG  >  CP009273/2665625‑2665919

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: