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,7720GT99.8% 731.4 / ‑6.5 238E123* (GAA→TAA) hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (113/123);  minor base A (2/0);  total (115/123)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
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 239 total aligned reads displayed.

TTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACA  >  CP009273/2665623‑2665918
                                                                                                                                                     |                                                                                                                                                   
  cATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCGCGACCCCTGCCGTTGCTTGTTTCACAGGGATaa                                                                                                                                                   >  1:4726669/1‑150 (MQ=255)
  cATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATaa                                                                                                                                                   <  1:4384526/150‑1 (MQ=255)
   aTCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAAc                                                                                                                                                  <  1:1786147/150‑1 (MQ=255)
      tGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgct                                                                                                                                               <  1:4412484/150‑1 (MQ=255)
      tGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgct                                                                                                                                               <  2:4016486/150‑1 (MQ=255)
       gCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctg                                                                                                                                              <  2:3034301/150‑1 (MQ=255)
       gCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACGGGGGTAACgctg                                                                                                                                              >  1:2925147/1‑150 (MQ=255)
        cagtCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGGACTGGGATCAGCTTTTTAGAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctgc                                                                                                                                             <  2:4327678/147‑1 (MQ=255)
        cATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctgc                                                                                                                                             >  2:3454590/1‑150 (MQ=255)
           tCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCttt                                                                                                                                          <  2:567896/150‑1 (MQ=255)
           tCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCttt                                                                                                                                          <  2:5018222/150‑1 (MQ=255)
           tCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCttt                                                                                                                                          <  1:5015291/150‑1 (MQ=255)
            ccGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTa                                                                                                                                         <  1:4666888/150‑1 (MQ=255)
              ggTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCACGACCGCGGCCGTTGCTTGATGCACTGGGGTAAAGCGGCTTTTcc                                                                                                                                       >  1:1703416/1‑150 (MQ=255)
                tgtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCt                                                                                                                                    <  2:1026906/150‑1 (MQ=255)
                tgtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCt                                                                                                                                    <  2:3303150/150‑1 (MQ=255)
                 gtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTa                                                                                                                                   <  2:3616484/150‑1 (MQ=255)
                   aaCCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCGCGACCCCAGCCGGTGCGTGATGCACAGGGATAACGCTGCATTACCCAAcc                                                                                                                                 >  2:5131388/1‑149 (MQ=255)
                      cATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCAACGCCa                                                                                                                              >  2:2538760/1‑150 (MQ=255)
                      cATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCGGCTTTACCCTACGCCa                                                                                                                              >  2:282541/1‑150 (MQ=255)
                       aTCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCAGCCGTTGCGTGATGCACTGGGATAACGCTGCGTTACCCAACGccca                                                                                                                             >  2:1404696/1‑148 (MQ=255)
                         cAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGc                                                                                                                           >  2:3952810/1‑150 (MQ=255)
                              cGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATACCGGTg                                                                                                                      <  1:2988876/150‑1 (MQ=255)
                               gCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCTTGCTTGATGCACAGGGATAACGCTGCTTTACCCAACGCCATGCCCGTGa                                                                                                                     >  2:5227470/1‑150 (MQ=255)
                               gCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGGTTCTTGATGCACAGGGATTAAGCAGCGTTAACCAACGCCATTCCGGTGa                                                                                                                     >  2:2410646/1‑150 (MQ=255)
                                cTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGAt                                                                                                                    <  1:2898831/150‑1 (MQ=255)
                                 tGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCACGACCCCAGCCGGTGCTTGTTGCACAGGGATAAAGCTGCTTTTACCAACGCCAGGCCGGGGATg                                                                                                                   >  1:1280852/1‑150 (MQ=255)
                                   ggTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCGGCCGTTGCGTGATGCACGGGGATAACGCTGCGGTAGCCGAAGCCCTGCCGGGGCTgcc                                                                                                                 >  2:3657319/1‑150 (MQ=255)
                                    gTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATgccg                                                                                                                <  1:4496430/150‑1 (MQ=255)
                                        cgACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAg                                                                                                            <  2:3786093/150‑1 (MQ=255)
                                        cgACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCGGCCGGTGCTTGATTCACTGGGATAACGCTGGTTTAACCAACGCCATGCCGGGGAGgccgcccg                                                                                                            >  1:3723801/1‑148 (MQ=255)
                                           cacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACt                                                                                                         <  1:4899426/150‑1 (MQ=255)
                                           cacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACt                                                                                                         <  1:2623023/150‑1 (MQ=255)
                                           cacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACt                                                                                                         <  1:1621692/150‑1 (MQ=255)
                                              acGAGAGTTAGTGTTAACCAACGGCGACAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCg                                                                                                      <  1:963145/150‑1 (MQ=255)
                                              acGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCGGCCGTTGCATGATGCAATGGGATAACGCAGCGTTACCCAACGCCAAGCCGGGGATGCCCCCAAACTGCg                                                                                                      >  1:3537155/1‑150 (MQ=255)
                                                  gagTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCAt                                                                                                                             >  2:3113888/1‑123 (MQ=255)
                                                  gagTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCAt                                                                                                                             <  1:3113888/123‑1 (MQ=255)
                                                    gTTAGTGTCAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGt                                                                                                <  2:4259068/150‑1 (MQ=255)
                                                      tAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTc                                                                                              <  2:3873521/150‑1 (MQ=255)
                                                       aGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATCGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCt                                                                                             <  2:1209085/150‑1 (MQ=255)
                                                        gtgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGGTGCGTGATGCACAGGGATAACGCGGCTTTTCCCAACGCCATTCCGGGGTTGCCGCCACACTGCGAGAACTTCTg                                                                                            >  2:2696180/1‑150 (MQ=255)
                                                         tgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGc                                                                                           <  1:123995/150‑1 (MQ=255)
                                                          gtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCa                                                                                          <  1:3709558/150‑1 (MQ=255)
                                                          gtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCACGACCCCTGCCGTTGGGTGATGCACTGGGGTAAAGCAGCGTTACACAACCCCCAGCCGGGGAAGCCGCCCCAAAGCGAGAAGTTCTGCa                                                                                          >  1:613134/1‑150 (MQ=255)
                                                             aaCCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCTGCCGGTGCGTGATTCACTGGGATAACGCGGCGGTACCCGACGCCATGCCGGGGATGCCGCCCGACAGCGAGAAGTTATGCAGcc                                                                                       >  1:3025500/1‑150 (MQ=255)
                                                              aCCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCTGCCGTTGCTTGTTTCACTGGGTTAAAGCTGCTTTTCCCAACGCCATGCCGGTGGTGCCGCCAGACAGCGAGAAGTTATGCAGccc                                                                                      >  2:1137657/1‑150 (MQ=255)
                                                                cAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGTTAACGCGGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGa                                                                                    >  2:2139534/1‑150 (MQ=255)
                                                                 aaCGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGCGCGCGACCCCAGCCGTTGCTTGATTGACAGGGTTTAAGCGGCGTTAACCAACGCCCAGCCCGTGATGCCGCCCGACAGCCAGAAGTTCAGCAGCCCCaa                                                                                   >  2:5093220/1‑150 (MQ=255)
                                                                 aaCGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGTTTCGTGTTTCACAGGGATAAAGCTGCGTTTACCTACCCCCTGCCGGTGATGCCGCCAGACAGCGAGAACTTCAGCAGCCCCaa                                                                                   >  2:1068164/1‑150 (MQ=255)
                                                                   cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACg                                                                                 <  1:1725665/150‑1 (MQ=255)
                                                                      cGAAAGCAGGCACTGGGATCAGCTTTTTAGAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTc                                                                              <  2:4487817/150‑1 (MQ=255)
                                                                        aaaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAg                                                                            <  2:1618486/150‑1 (MQ=255)
                                                                        aaaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAg                                                                            <  2:222151/150‑1 (MQ=255)
                                                                        aaaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACAGCGAGAAGTTCTGCAGCCCGAACGGTCAg                                                                            >  2:3261290/1‑150 (MQ=255)
                                                                         aaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGt                                                                           <  1:2527783/150‑1 (MQ=255)
                                                                         aaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCAGCCGTTGCTTGATGCACTGGGATAAAGCTGCTTTAAACAACGCCAAGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGt                                                                           >  1:4382745/1‑150 (MQ=255)
                                                                         aaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCAGCCGTTGCGTGATTCACGGGGATAACGCAGCTTTACCCTACCCCATGCCGGTGGTGCCGCCAGACAGCGAGAAGTTCATCAGCCCGAAAGGTCAGt                                                                           >  2:2098736/1‑150 (MQ=255)
                                                                           gCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCGCGACCGCGGCCGGTGCGTTATTCACTGGGTTAAAGCGGCGTTTCCCAACGCCAAGCCGGTGATTCCGCCAGACAGCGAGAAGTTATGCAGCCCGAACGGTCAGGCg                                                                         >  1:1717872/1‑150 (MQ=255)
                                                                           gCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCg                                                                         <  1:860289/150‑1 (MQ=255)
                                                                              ggCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTAGAACGCTGCCGTTGCTTGATTCACTGGGATAACGCTGCTTTACCCTACGCCATTCCGGTGATGCCGCCAGACTGCTAGAAGTTCTGCAGCCCGAACGGTCAGTAGTGa                                                                      >  2:871613/1‑150 (MQ=255)
                                                                                  cTGGGATCATCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTAt                                                                  <  1:676210/150‑1 (MQ=255)
                                                                                  cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCCCGCGACCCCAGCCCGTTCGTGATTCACAGGGATAACCCAGCATTACCCAAACCCCAGCCCGGGGAGCCCCCACACAGCCCGAACATCAGCAGCCCCAACCGTCAGTCGTGAatat                                                                  >  2:1926559/1‑150 (MQ=255)
                                                                                  cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTAt                                                                  <  1:2290542/150‑1 (MQ=255)
                                                                                  cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTAt                                                                  <  2:2545936/150‑1 (MQ=255)
                                                                                   tGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATc                                                                 <  1:4044227/150‑1 (MQ=255)
                                                                                       aTCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCCCTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTg                                                             <  1:1060400/150‑1 (MQ=255)
                                                                                       aTCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCCCAGCCGTTGCGTGATGCCCAGGGATAACGCGGCTTTTACCGACGCCCAGCCCGTGATTCCGCCACAAAGCGAGAACTTCTGCAACCCCAACGGGCAGTCGGGATTATCgggg                                                             >  2:4578361/1‑148 (MQ=255)
                                                                                         cAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGcc                                                           <  2:4101968/150‑1 (MQ=255)
                                                                                         cAGCTTTTTATAGCAACCGGCGCGGCAGCACGAACCCAGCCCTTGAATTATGCACTGGGGTAAACCTGCTTTAAACAAACCCAAGACGGGTATTACGCCCCACAACGAGAACTTCTGCACCCCGCACGGTCAATACTGATTATAGGTggc                                                           >  2:5024555/1‑148 (MQ=255)
                                                                                          aGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTTCCCTACGCCATGCCGGTGATGCCGCCAGACTGCAACAACTTCTGCACCCCCAACCGTCAGTCGTGAATATCGGTGCCg                                                          >  2:3917843/1‑150 (MQ=255)
                                                                                               tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACt                                                     <  1:1510368/150‑1 (MQ=255)
                                                                                                 tataGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                   <  1:1617768/150‑1 (MQ=255)
                                                                                                 tataGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                   <  1:2488887/150‑1 (MQ=255)
                                                                                                     gCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGt                                               <  1:5208675/150‑1 (MQ=255)
                                                                                                       aaCCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCt                                             <  2:2999106/150‑1 (MQ=255)
                                                                                                        aCCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTg                                            <  1:1732229/150‑1 (MQ=255)
                                                                                                         ccGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTgg                                           <  2:2423128/150‑1 (MQ=255)
                                                                                                               cgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                     <  1:2668483/150‑1 (MQ=255)
                                                                                                                ggCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGc                                    <  2:1224565/150‑1 (MQ=255)
                                                                                                                 gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCt                                   >  1:5270813/1‑150 (MQ=255)
                                                                                                                          ccGCTGCCGGTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGcc                          <  2:4291797/150‑1 (MQ=255)
                                                                                                                            gCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTCCGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAc                        <  2:3379915/150‑1 (MQ=255)
                                                                                                                            gCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAc                        <  1:3493819/150‑1 (MQ=255)
                                                                                                                            gCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAc                        >  1:4612396/1‑150 (MQ=255)
                                                                                                                             cTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcg                       <  1:3789047/150‑1 (MQ=255)
                                                                                                                               gCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgca                     >  2:1694029/1‑150 (MQ=255)
                                                                                                                               gCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgca                     <  1:2242831/150‑1 (MQ=255)
                                                                                                                               gCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAATATCGGTGCCGGAACTATAGGTCTGGAACTGGCTGCAAGCGCCAcgca                     >  1:4788201/1‑150 (MQ=255)
                                                                                                                                ccGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                    >  2:4733439/1‑150 (MQ=255)
                                                                                                                                ccGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                    >  2:1502230/1‑150 (MQ=255)
                                                                                                                                 cGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagc                   <  2:1514724/150‑1 (MQ=255)
                                                                                                                                      cTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGACCTGGCTGCCAGCGCCACGCAGCGCAGa              >  1:3182225/1‑150 (MQ=255)
                                                                                                                                      cTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa              >  2:2901493/1‑150 (MQ=255)
                                                                                                                                          aTGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTa          >  1:2507184/1‑150 (MQ=255)
                                                                                                                                             cACTGGGAAAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTCTTTGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAgg       >  2:1498605/1‑150 (MQ=255)
                                                                                                                                               cTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTg     <  1:5131388/150‑1 (MQ=255)
                                                                                                                                                tGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGa    >  2:70130/1‑150 (MQ=255)
                                                                                                                                                tGGGATAACGCTGATTTACCCTACGCCATGCCGGTGATGCCGCCAGACAGCGAGAAGTTCTGCAGCCCGAACTGTCAGTCGAGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGgaga    >  2:828014/1‑150 (MQ=255)
                                                                                                                                                  ggATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACa  >  1:4065141/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACA  >  CP009273/2665623‑2665918

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: