Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,451,422 A→T 100% T246S (ACG→TCG)  paaE → ring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,451,4220AT98.8% 811.2 / ‑0.6 246T246S (ACG→TCG) paaEring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base T (127/116);  total (130/116)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
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 246 total aligned reads displayed.

CAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGG  >  CP009273/1451278‑1451566
                                                                                                                                                |                                                                                                                                                
          gggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAcgcg                                                                                                                                     >  1:1375559/1‑148 (MQ=255)
          gggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAcgcg                                                                                                                                     >  1:1375446/1‑148 (MQ=255)
           ggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAcgcgc                                                                                                                                    >  1:2290850/1‑148 (MQ=255)
                cGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCaa                                                                                                                               >  1:2742565/1‑148 (MQ=255)
                cGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCaa                                                                                                                               >  2:205003/1‑148 (MQ=255)
                 gCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCaaa                                                                                                                              >  1:3127035/1‑148 (MQ=255)
                  cTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAAc                                                                                                                             >  1:314116/1‑148 (MQ=255)
                  cTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAAc                                                                                                                             >  1:571603/1‑148 (MQ=255)
                   tCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACg                                                                                                                            <  2:805774/148‑1 (MQ=255)
                   tCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACg                                                                                                                            >  1:2022999/1‑148 (MQ=255)
                     aTTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTa                                                                                                                          >  2:2361818/1‑148 (MQ=255)
                     aTTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTa                                                                                                                          <  1:1002124/148‑1 (MQ=255)
                      ttAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAg                                                                                                                         >  1:2976454/1‑148 (MQ=255)
                        aaTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCg                                                                                                                       <  2:1709788/148‑1 (MQ=255)
                        aaTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCg                                                                                                                       <  1:2181409/148‑1 (MQ=255)
                          ttttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGtt                                                                                                                     >  2:2255573/1‑148 (MQ=255)
                           tttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTa                                                                                                                    <  1:954150/148‑1 (MQ=255)
                           tttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTa                                                                                                                    <  1:232215/148‑1 (MQ=255)
                            ttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTaa                                                                                                                   <  2:2824863/148‑1 (MQ=255)
                                tCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTg                                                                                                               >  2:449602/1‑148 (MQ=255)
                                 cTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGc                                                                                                              >  2:1226232/1‑148 (MQ=255)
                                   ttATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCaa                                                                                                            >  2:3516224/1‑148 (MQ=255)
                                   ttATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCaa                                                                                                            >  2:1998668/1‑148 (MQ=255)
                                      tgatgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGt                                                                                                         >  2:1563968/1‑148 (MQ=255)
                                        atgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGa                                                                                                       <  1:2797197/148‑1 (MQ=255)
                                        atgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGa                                                                                                       <  2:3456145/148‑1 (MQ=255)
                                         tgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGAc                                                                                                      >  2:2209727/1‑148 (MQ=255)
                                            gtcATTTTTTTGTGGTGCGGCGGCGATGATGGTTGACGCGGAAACCGCCTTAAAAGCACTGGGAATGGCAGATAAAACCATTCTTGTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTTACGTGCTAAGTGACGGa                                                                                                   <  2:71993/146‑1 (MQ=255)
                                            ggCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGa                                                                                                   >  1:2072537/1‑148 (MQ=255)
                                              cATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACa                                                                                                 >  2:2015156/1‑148 (MQ=255)
                                                   atttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAg                                                                                            <  1:2929026/148‑1 (MQ=255)
                                                   atttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAg                                                                                            <  2:272111/148‑1 (MQ=255)
                                                   atttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAg                                                                                            <  1:3380696/148‑1 (MQ=255)
                                                    tttGTGGTCCGGCGGCGATGATGGATGACGCTGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGt                                                                                           <  2:3538385/148‑1 (MQ=255)
                                                           tCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTa                                                                                    <  2:582530/148‑1 (MQ=255)
                                                           tCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTa                                                                                    <  1:1080152/148‑1 (MQ=255)
                                                                cggcgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCa                                                                               <  1:3053297/148‑1 (MQ=255)
                                                                cggcgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCa                                                                               <  1:3299383/148‑1 (MQ=255)
                                                                  gcgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAgg                                                                             <  1:2424653/148‑1 (MQ=255)
                                                                   cgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGa                                                                            <  2:1532423/148‑1 (MQ=255)
                                                                   cgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGa                                                                            <  2:3290463/148‑1 (MQ=255)
                                                                     atgatgGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATg                                                                          >  2:3386162/1‑148 (MQ=255)
                                                                       gatggatgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATggg                                                                        >  2:1461075/1‑148 (MQ=255)
                                                                            atgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGa                                                                   <  2:379398/148‑1 (MQ=255)
                                                                             tgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGAt                                                                  <  2:2247779/148‑1 (MQ=255)
                                                                               aCGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCg                                                                <  2:3610205/148‑1 (MQ=255)
                                                                                cgcgGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCgg                                                               >  1:1472959/1‑148 (MQ=255)
                                                                                  cgGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGa                                                             <  1:2081197/148‑1 (MQ=255)
                                                                                    gAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGaaa                                                           <  1:3535704/148‑1 (MQ=255)
                                                                                     aaaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAAt                                                          >  1:3106089/1‑148 (MQ=255)
                                                                                     aaaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAAt                                                          >  1:1678890/1‑148 (MQ=255)
                                                                                     aaaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAAt                                                          >  1:3105581/1‑148 (MQ=255)
                                                                                        ccGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGt                                                       >  1:3572547/1‑148 (MQ=255)
                                                                                         cGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTg                                                      <  1:1361974/148‑1 (MQ=255)
                                                                                         cGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTg                                                      <  1:1361926/148‑1 (MQ=255)
                                                                                            cTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCtt                                                   <  1:717276/148‑1 (MQ=255)
                                                                                             ttAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTa                                                  <  1:3686432/148‑1 (MQ=255)
                                                                                             ttAAAAGCACTGGCAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCTCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGAAAAGTGACCGACAAAAAGTGACTCTACGTCAGGATGGGCGGGATCGGGAACTCGTGCTTa                                                  >  2:3326294/1‑148 (MQ=255)
                                                                                              tAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTaa                                                 >  1:870/1‑148 (MQ=255)
                                                                                              tAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTaa                                                 >  1:1785627/1‑148 (MQ=255)
                                                                                              tAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTaa                                                 >  1:642397/1‑148 (MQ=255)
                                                                                                aaaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATg                                               <  2:1372138/148‑1 (MQ=255)
                                                                                                aaaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATg                                               <  1:2118883/148‑1 (MQ=255)
                                                                                                aaaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATg                                               <  1:3619969/148‑1 (MQ=255)
                                                                                                  aGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGcc                                             >  1:1059909/1‑148 (MQ=255)
                                                                                                   gCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcg                                            <  2:2158569/148‑1 (MQ=255)
                                                                                                       tGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcgacga                                        <  2:2679965/148‑1 (MQ=255)
                                                                                                          gAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGa                                     >  2:2443556/1‑148 (MQ=255)
                                                                                                          gAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGa                                     >  1:1764876/1‑148 (MQ=255)
                                                                                                           aaTGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaa                                    >  2:2091556/1‑148 (MQ=255)
                                                                                                            aTGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaaa                                   <  1:400044/148‑1 (MQ=255)
                                                                                                            aTGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaaa                                   <  1:1100433/148‑1 (MQ=255)
                                                                                                             tGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAg                                  >  1:1978636/1‑148 (MQ=255)
                                                                                                              gCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGc                                 <  2:2690128/148‑1 (MQ=255)
                                                                                                               ccAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCa                                >  2:1747433/1‑148 (MQ=255)
                                                                                                                 aGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAtt                              <  2:179548/148‑1 (MQ=255)
                                                                                                                 aGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAtt                              <  2:2918084/148‑1 (MQ=255)
                                                                                                                 aGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAtt                              <  1:2463395/148‑1 (MQ=255)
                                                                                                                 aGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAtt                              <  2:659729/148‑1 (MQ=255)
                                                                                                                  gATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtc                             >  2:284932/1‑148 (MQ=255)
                                                                                                                  gATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtc                             <  1:3459710/148‑1 (MQ=255)
                                                                                                                  gATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtc                             >  1:3406434/1‑148 (MQ=255)
                                                                                                                   aTAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTNACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCNGACGATGAAAGCATtct                            <  2:3127013/148‑1 (MQ=255)
                                                                                                                    tAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtctc                           >  2:1151557/1‑148 (MQ=255)
                                                                                                                    tAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtctc                           >  2:1151369/1‑148 (MQ=255)
                                                                                                                     aaaaCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCg                          <  2:3026656/148‑1 (MQ=255)
                                                                                                                        aCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                       >  2:882932/1‑148 (MQ=255)
                                                                                                                        aCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                       <  2:1101390/148‑1 (MQ=255)
                                                                                                                        aCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                       >  1:3342090/1‑148 (MQ=255)
                                                                                                                        aCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                       >  1:1132298/1‑148 (MQ=255)
                                                                                                                         ccATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGc                      >  1:3773243/1‑148 (MQ=255)
                                                                                                                         ccATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGc                      >  1:1883181/1‑148 (MQ=255)
                                                                                                                         ccATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGc                      <  2:829727/148‑1 (MQ=255)
                                                                                                                             tCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCa                  >  2:3187777/1‑148 (MQ=255)
                                                                                                                             tCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCa                  >  2:209868/1‑148 (MQ=255)
                                                                                                                                tCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTg               >  2:2516430/1‑148 (MQ=255)
                                                                                                                                tCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTg               >  2:760279/1‑148 (MQ=255)
                                                                                                                                tCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTg               <  2:953577/148‑1 (MQ=255)
                                                                                                                                tCTGGAGCGGTTTAATTAGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTg               >  2:1469938/1‑148 (MQ=255)
                                                                                                                                  tGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgcg             <  1:449602/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGG  >  CP009273/1451278‑1451566

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