Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,665,772 G→T 100% 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,7720GT98.2% 367.7 / ‑1.2 127E123* (GAA→TAA) hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (44/78);  total (47/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 127 total aligned reads displayed.

CTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTG  >  CP009273/2665628‑2665915
                                                                                                                                                 |                                                                                                                                               
cTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCCCGGCCGTTGCTTGATGCACTGGGATaaa                                                                                                                                              >  1:2311372/1‑147 (MQ=255)
 tGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACg                                                                                                                                             <  1:570304/148‑1 (MQ=255)
 tGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACg                                                                                                                                             <  2:162832/148‑1 (MQ=255)
 tGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACg                                                                                                                                             <  1:570870/148‑1 (MQ=255)
 tGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACg                                                                                                                                             <  2:3068015/148‑1 (MQ=255)
 tGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACg                                                                                                                                             <  2:163186/148‑1 (MQ=255)
  gCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGCCCGCTGCCGGTGCTTGTTGCACGGGGATAAAgc                                                                                                                                            >  2:1806479/1‑148 (MQ=255)
   cATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCCCTGGGATAACgct                                                                                                                                           <  2:1253780/148‑1 (MQ=255)
    aTTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctg                                                                                                                                          <  1:2355635/148‑1 (MQ=255)
       ccGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACGGGGATAACGCGGCtt                                                                                                                                       >  2:2905508/1‑148 (MQ=255)
          gtgtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTAc                                                                                                                                    <  2:2342869/148‑1 (MQ=255)
              aaCCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTa                                                                                                                                <  1:2116027/148‑1 (MQ=255)
                 cATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGc                                                                                                                             <  2:1870116/148‑1 (MQ=255)
                 cATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCACCTACACCGCTGCCGGTGATTGGTGCAAAGGGATAACGCTGCATTAACCTAcca                                                                                                                             >  2:1614145/1‑146 (MQ=255)
                  aTCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATTCACTGGGATAACGCTGCTTTACCCTAAGcc                                                                                                                            >  2:2597864/1‑148 (MQ=255)
                     aaaaCGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATg                                                                                                                         <  2:216126/148‑1 (MQ=255)
                       aaCGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGcc                                                                                                                       <  1:2716301/148‑1 (MQ=255)
                       aaCGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGGTGCTTGATGCCCTGGGATAACGCTGCTTTACCCTACGCCATGcc                                                                                                                       >  2:3337553/1‑148 (MQ=255)
                         cGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCgg                                                                                                                     <  2:2553317/148‑1 (MQ=255)
                           cTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTg                                                                                                                   <  2:3330599/148‑1 (MQ=255)
                           cTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCGTGATGCACTGGGATAACGCGGCGTTACCCTACGCCATGCCGGTg                                                                                                                   >  2:805892/1‑148 (MQ=255)
                              ggTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCGGCGCGACCGCGGCCCGTGCTTTATGCACTGGGATAAAGCTGCGTTTCCCTACGCCAAGCCCGTGATg                                                                                                                >  1:1282139/1‑148 (MQ=255)
                              ggTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATg                                                                                                                <  1:626589/148‑1 (MQ=255)
                              ggTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATg                                                                                                                <  2:2032689/148‑1 (MQ=255)
                               gTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATgc                                                                                                               <  1:2507408/148‑1 (MQ=255)
                                     acacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAg                                                                                                         <  1:1549716/148‑1 (MQ=255)
                                     acacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCCCGGCCGGTGGGTGTTGCACAGGGGTAAACCCGCTGTTACCAACGCCAACCCGGGGATGGCGCCAg                                                                                                         >  1:2489329/1‑148 (MQ=255)
                                      cacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGa                                                                                                        <  2:136583/148‑1 (MQ=255)
                                                  aGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGtt                                                                                            <  1:2335071/148‑1 (MQ=255)
                                                   gtgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTc                                                                                           <  1:2521739/148‑1 (MQ=255)
                                                   gtgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTc                                                                                           <  1:1695978/148‑1 (MQ=255)
                                                   gtgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTc                                                                                           <  1:734426/148‑1 (MQ=255)
                                                     gtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAAAGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAACTTCTg                                                                                         >  2:835048/1‑148 (MQ=255)
                                                      ttAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGc                                                                                        <  1:105033/148‑1 (MQ=255)
                                                      ttAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGc                                                                                        <  1:105475/148‑1 (MQ=255)
                                                      ttAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGc                                                                                        <  1:105575/148‑1 (MQ=255)
                                                      ttAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGc                                                                                        <  1:105030/148‑1 (MQ=255)
                                                             aCGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGa                                                                                 <  1:1707409/148‑1 (MQ=255)
                                                              cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGaa                                                                                >  2:1477602/1‑148 (MQ=255)
                                                               ggCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAAc                                                                               <  1:3190366/148‑1 (MQ=255)
                                                               ggCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAAc                                                                               <  1:733710/148‑1 (MQ=255)
                                                                 cGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACgg                                                                             <  1:2714654/148‑1 (MQ=255)
                                                                 cGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACgg                                                                             <  1:2714454/148‑1 (MQ=255)
                                                                   aaaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTc                                                                           <  1:3530880/148‑1 (MQ=255)
                                                                   aaaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCTGCCCTTGGGTGATTCCCATGGGTAAACCAGGGTTAACCAACGCCAAGCCCGTGTTGGCGCCCCACAGCGAGAAGTTATTCACCCCCAAccggc                                                                           >  2:1206985/1‑144 (MQ=255)
                                                                      gCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGt                                                                        <  1:2684942/148‑1 (MQ=255)
                                                                       cTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTc                                                                       <  2:3159691/148‑1 (MQ=255)
                                                                        tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCg                                                                      <  1:2729137/148‑1 (MQ=255)
                                                                             cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAtt                                                                 <  1:2711054/148‑1 (MQ=255)
                                                                              tGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGCCCGCGGCCCGTGCTTGATTGACTGGGATAAAGCTGCTTTAACCAACGCCATGCCGGTGATGCCGCCCGACTGCGAGAAGGTCTGCAGCCCGAACGGTGCGTCGTGAttt                                                                >  2:401287/1‑147 (MQ=255)
                                                                               gggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTAt                                                               <  2:848128/148‑1 (MQ=255)
                                                                               gggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTAt                                                               <  1:432730/148‑1 (MQ=255)
                                                                                 gATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCg                                                             <  2:2311372/148‑1 (MQ=255)
                                                                                   tCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGTTGGTTGATTCACTGGGATAAAGCAGCTTTACCCTACGCCAAGCCCGTGATGGCGCCAGACAGCGAGAAATTATGCAGCCCGAACGGTCACTGGTGATTATAGGt                                                           >  2:2426207/1‑148 (MQ=255)
                                                                                    cAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTg                                                          <  2:1164746/148‑1 (MQ=255)
                                                                                      tcTTATTATAGCAACCGCCGCGACAACTCTAACGATCACGTTTCTATATGCCATGGGAAAAACCAACTATAACCTACGCCATGCCGGTGATGCCGCCAGACTGCGCGAAGTTCTGCAGCCCGCACGGTCAGTCGTGGTTATCgggggc                                                        >  2:2309902/2‑144 (MQ=255)
                                                                                       cTTTTTATAGCAACCGGCGCGGGAGCTCGACCGCGGGCGGTGGGTTATGGACAGGGGTAAACCGGCGTTTCCCGAAGCCATGCCCGTGATTCCGCCCGACTGCGCGAAGTTCTGCCGCCCGAACGGGCAGTCGTGGTTTTCGGTGCCg                                                       >  2:646199/1‑148 (MQ=255)
                                                                                        tttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTTGTTGATGCACTGGGATAACGCAGCTTTTCCCAACGCCATGCCGGTGATGCCGCCCGACTGCGAGAAGTTCTGCACCCCGCACCGTCAGTCGTGATTATCGGGGCCgg                                                      >  2:845572/1‑148 (MQ=255)
                                                                                         ttttATAGCAACCGGCGCGGCCGCGCGACCGCTGCCGTTGCTTGATGCACTGGGGTAACGCAGCTTTACCCTACCCCCTGCCGGTGATGCCGCCACACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGa                                                     >  2:2170985/1‑148 (MQ=255)
                                                                                         ttttATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGTTGCTTGTTGCACTGGGATAACGCTGCTTTAACCAACGCCAAGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGa                                                     >  2:1460395/1‑148 (MQ=255)
                                                                                             ataGCAACCGGCGCGGCCGCTAGACCGCTGCCGTTGCTTTATGCACTGGGATAACGCTGCTTTACCCTTCGCCATTCCGGTGATGCCGCCACACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGGTTATCGGTGGCGGAACTa                                                 >  2:2617817/1‑148 (MQ=255)
                                                                                               aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGGTTGATGCACTGGGATAACGCAGCTTTACCCAACGCCATGCCGGTGATGCCGCCACACAGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGCACTAtt                                               >  2:1094583/1‑148 (MQ=255)
                                                                                               aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGGTGCACGGGGCTAAAGCGGCTTTAACCTACCCCCAGCCGGGGATTCCCCCCGACTGCGAGAACTTCTGCAGCCCCAACCGGCACGCGGGCTTATCGGTGCCGGAACGAtt                                               >  2:372578/1‑148 (MQ=255)
                                                                                               aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAtt                                               <  2:1650375/148‑1 (MQ=255)
                                                                                               aGCAACCGGCGCGGCAGCTCGACCGCGGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCAACGCCAAGCCCGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCCAACGGTCAGTCGGGATTATCGGTGCCGGAACTAtt                                               >  2:1090517/1‑148 (MQ=255)
                                                                                                gCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCCCGGGGGTAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCCGACTGCGAGAAGGTCTGCCGCCCGCACGGTCCGGCGTGATTATAGGGGCCGGAACGATTg                                              >  2:3232825/1‑148 (MQ=255)
                                                                                                  aaCCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGt                                            <  2:1985104/148‑1 (MQ=255)
                                                                                                   aCCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTc                                           <  2:2046826/148‑1 (MQ=255)
                                                                                                        cgcgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGaa                                      <  2:1260680/148‑1 (MQ=255)
                                                                                                        cgcgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGaa                                      >  2:204459/1‑148 (MQ=255)
                                                                                                        cgcgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGaa                                      <  2:721220/148‑1 (MQ=255)
                                                                                                            gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                  <  1:625118/148‑1 (MQ=255)
                                                                                                            gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                  >  1:2352514/1‑148 (MQ=255)
                                                                                                            gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                  <  2:2520945/148‑1 (MQ=255)
                                                                                                               gCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTg                               >  1:1143511/1‑148 (MQ=255)
                                                                                                               gCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTg                               >  2:1501970/1‑148 (MQ=255)
                                                                                                                cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                              <  1:535393/148‑1 (MQ=255)
                                                                                                                cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                              >  1:2495392/1‑148 (MQ=255)
                                                                                                                cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                              <  1:1617240/148‑1 (MQ=255)
                                                                                                                cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                              <  1:3177568/148‑1 (MQ=255)
                                                                                                                 tCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGcc                             <  2:619665/148‑1 (MQ=255)
                                                                                                                    aCCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgc                          >  2:477210/1‑148 (MQ=255)
                                                                                                                    aCCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgc                          >  2:206700/1‑148 (MQ=255)
                                                                                                                    aCCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgc                          >  2:477187/1‑148 (MQ=255)
                                                                                                                    aCCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgc                          >  2:1878588/1‑148 (MQ=255)
                                                                                                                      cGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgcgc                        >  1:335539/1‑148 (MQ=255)
                                                                                                                       gCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGcc                       <  2:3196826/148‑1 (MQ=255)
                                                                                                                          gCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcg                    >  2:181377/1‑148 (MQ=255)
                                                                                                                             gTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAATTTATGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                 <  1:1614145/148‑1 (MQ=255)
                                                                                                                             gTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                 >  2:3287314/1‑148 (MQ=255)
                                                                                                                               tGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcg               >  2:330855/1‑148 (MQ=255)
                                                                                                                                  ttGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgcag            <  1:1581241/148‑1 (MQ=255)
                                                                                                                                   tGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa           <  1:2662521/148‑1 (MQ=255)
                                                                                                                                   tGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa           <  1:3529974/148‑1 (MQ=255)
                                                                                                                                   tGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa           <  1:463216/148‑1 (MQ=255)
                                                                                                                                     aTGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg         >  1:3232487/1‑148 (MQ=255)
                                                                                                                                     aTGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg         <  1:1244540/148‑1 (MQ=255)
                                                                                                                                       gCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTa       <  2:1497125/148‑1 (MQ=255)
                                                                                                                                        cACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTaa      <  2:2024295/148‑1 (MQ=255)
                                                                                                                                          cTGTGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGNTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAgg    <  1:1806479/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTG  >  CP009273/2665628‑2665915

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: