Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,441,941 G→T 11.7% E56* (GAA→TAA)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,441,9410GT11.7% 1052.7 / 72.3 319E56* (GAA→TAA) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (145/136);  new base T (21/17);  total (166/153)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 326 total aligned reads displayed.

TAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAA  >  CP009273/1441792‑1442090
                                                                                                                                                     |                                                                                                                                                     
tAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACg                                                                                                                                                       <  1:1907445/150‑1 (MQ=255)
    aGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACTAAGc                                                                                                                                                   >  1:3405441/1‑150 (MQ=255)
    aGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGc                                                                                                                                                   >  2:4350463/1‑150 (MQ=255)
    aGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGc                                                                                                                                                   >  1:2572041/1‑150 (MQ=255)
     gTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACTAAGCg                                                                                                                                                  >  2:1266005/1‑150 (MQ=255)
     gTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCg                                                                                                                                                  >  2:2474245/1‑150 (MQ=255)
      tATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCgg                                                                                                                                                 <  1:3297208/150‑1 (MQ=255)
        ttAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGTGGAt                                                                                                                                               <  2:1686434/150‑1 (MQ=255)
              cAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGTTGCCAACGAAGCGGATGTAGAt                                                                                                                                         >  2:3491089/1‑150 (MQ=255)
                ggTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATaa                                                                                                                                       >  1:3946425/1‑150 (MQ=255)
                   ccAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGc                                                                                                                                    >  2:3651604/1‑150 (MQ=255)
                    cAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCa                                                                                                                                   >  1:2007853/1‑150 (MQ=255)
                    cAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCa                                                                                                                                   >  2:4810451/1‑150 (MQ=255)
                     aaCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAg                                                                                                                                  <  2:895932/150‑1 (MQ=255)
                     aaCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAg                                                                                                                                  >  2:1519315/1‑150 (MQ=255)
                     aaCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAg                                                                                                                                  >  1:4932330/1‑150 (MQ=255)
                      aCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGt                                                                                                                                 >  1:1609858/1‑150 (MQ=255)
                                gATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGcc                                                                                                                       >  1:4122513/1‑150 (MQ=255)
                                gATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGcc                                                                                                                       >  1:314398/1‑150 (MQ=255)
                                  tCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTg                                                                                                                     <  2:1939381/150‑1 (MQ=255)
                                   cGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTgg                                                                                                                    >  1:4434034/1‑150 (MQ=255)
                                   cGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTgg                                                                                                                    <  2:4065935/150‑1 (MQ=255)
                                     tCAACACGGTCTTTATATTGCTGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                                                                                                                  <  2:4612901/150‑1 (MQ=255)
                                     tCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCTGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                                                                                                                  >  1:1298897/1‑150 (MQ=255)
                                     tCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCTGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                                                                                                                  >  1:209118/1‑150 (MQ=255)
                                     tCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                                                                                                                  >  1:4669773/1‑150 (MQ=255)
                                     tCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                                                                                                                  <  2:539974/150‑1 (MQ=255)
                                       aaCACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCggg                                                                                                                >  1:3182312/1‑150 (MQ=255)
                                       aaCACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCggg                                                                                                                >  1:3735421/1‑150 (MQ=255)
                                       aaCACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCggg                                                                                                                >  2:4055340/1‑150 (MQ=255)
                                         cacGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGcc                                                                                                              <  1:4350463/150‑1 (MQ=255)
                                          acGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCt                                                                                                             <  1:2698247/150‑1 (MQ=255)
                                            ggTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCttt                                                                                                           <  1:1078343/150‑1 (MQ=255)
                                            ggTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCttt                                                                                                           <  2:5051558/150‑1 (MQ=255)
                                            ggTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCttt                                                                                                           <  2:4715149/150‑1 (MQ=255)
                                             gTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTg                                                                                                          <  2:4145944/150‑1 (MQ=255)
                                             gTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTg                                                                                                          >  2:3151353/1‑150 (MQ=255)
                                              tCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGt                                                                                                         <  1:953961/150‑1 (MQ=255)
                                                tttATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGtct                                                                                                       <  2:2239551/150‑1 (MQ=255)
                                                   atatTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcgc                                                                                                    <  2:4465307/150‑1 (MQ=255)
                                                      ttGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACTAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTc                                                                                                 <  2:112886/150‑1 (MQ=255)
                                                          tGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTg                                                                                             <  1:4055350/150‑1 (MQ=255)
                                                              cgtcCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGcc                                                                                         <  1:2474245/150‑1 (MQ=255)
                                                                   tGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCg                                                                                    >  1:4982276/1‑150 (MQ=255)
                                                                     gCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGAt                                                                                  >  1:4606561/1‑150 (MQ=255)
                                                                       cccGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACTAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTa                                                                                >  2:722036/1‑150 (MQ=255)
                                                                       cccGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCTATTa                                                                                >  2:3032657/1‑150 (MQ=255)
                                                                        ccGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTAc                                                                               <  1:2396186/150‑1 (MQ=255)
                                                                            acaAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGc                                                                           >  1:2926337/1‑150 (MQ=255)
                                                                                aGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCTATTACCCGCagag                                                                       >  1:1358759/1‑150 (MQ=255)
                                                                                 gTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCTATTACCCGCAGAGc                                                                      <  1:2173242/150‑1 (MQ=255)
                                                                                  tGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCg                                                                     <  1:3060885/150‑1 (MQ=255)
                                                                                  tGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCg                                                                     >  2:1559405/1‑150 (MQ=255)
                                                                                   gAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGt                                                                    >  2:2728896/1‑150 (MQ=255)
                                                                                    aaaaaCGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCTATTACCCGCAGAGCGTg                                                                   <  2:2416169/150‑1 (MQ=255)
                                                                                    aaaaaCGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTg                                                                   >  1:4543566/1‑150 (MQ=255)
                                                                                      aaaCGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACTAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGaa                                                                 >  2:5303064/1‑150 (MQ=255)
                                                                                      aaaCGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGaa                                                                 >  1:5286512/1‑150 (MQ=255)
                                                                                             tGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCTATTACCCGCAGAGCGTGAACGTATTc                                                          <  2:2886150/150‑1 (MQ=255)
                                                                                             tGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTc                                                          <  2:2373614/150‑1 (MQ=255)
                                                                                             tGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTc                                                          <  1:737731/150‑1 (MQ=255)
                                                                                              ggCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACTAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCt                                                         <  1:2318363/150‑1 (MQ=255)
                                                                                              ggCGAGCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCCTGTCTGCCTGGCGGGCCTGTGTCTCGCGTCGCTGGGCCGGGCGAGTTCCCGCAGAGCGGGAACGTATTCt                                                         >  2:3316155/1‑150 (MQ=255)
                                                                                                   tCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACTAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTAc                                                    <  1:636829/150‑1 (MQ=255)
                                                                                                    cTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACg                                                   <  1:3851030/150‑1 (MQ=255)
                                                                                                     tttGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGt                                                  >  1:2989608/1‑150 (MQ=255)
                                                                                                       tGATCCGGCCACCGGGCAAGACATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCTATTACCCGCAGAGCGTGAACGTATTCTGCTACGttt                                                <  2:1727852/150‑1 (MQ=255)
                                                                                                        gATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGtttt                                               <  2:1341372/150‑1 (MQ=255)
                                                                                                        gATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGtttt                                               <  1:4390226/150‑1 (MQ=255)
                                                                                                           ccGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACTAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCt                                            >  1:4281630/1‑150 (MQ=255)
                                                                                                           ccGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACTAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCt                                            <  1:2827153/150‑1 (MQ=255)
                                                                                                           ccGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCt                                            <  1:3486407/150‑1 (MQ=255)
                                                                                                             ggCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGa                                          <  2:2815501/150‑1 (MQ=255)
                                                                                                             ggCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGa                                          >  2:3443946/1‑150 (MQ=255)
                                                                                                               ccACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATc                                        <  2:1175142/150‑1 (MQ=255)
                                                                                                                cACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCt                                       <  2:3814534/150‑1 (MQ=255)
                                                                                                                cACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCt                                       <  1:3676902/150‑1 (MQ=255)
                                                                                                                  ccGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCtgg                                     >  2:4679497/1‑150 (MQ=255)
                                                                                                                   cGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCtggt                                    >  1:509454/1‑150 (MQ=255)
                                                                                                                     ggCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCtggtgg                                  >  1:5128499/1‑150 (MQ=255)
                                                                                                                     ggCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCtggtgg                                  >  1:3454107/1‑150 (MQ=255)
                                                                                                                     ggCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCtggtgg                                  >  1:3669884/1‑150 (MQ=255)
                                                                                                                      gCAAGATATTGCGTCTACTGCTGATGCCAACGAAGCGGTTGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGa                                 <  2:1009589/150‑1 (MQ=255)
                                                                                                                             aTTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCTATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGcac                          >  1:791203/1‑150 (MQ=255)
                                                                                                                                gCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGt                       <  2:1609906/150‑1 (MQ=255)
                                                                                                                                 cGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTg                      <  1:1164943/150‑1 (MQ=255)
                                                                                                                                  gTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTga                     >  2:2604269/1‑150 (MQ=255)
                                                                                                                                   tCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTgag                    <  2:2260087/150‑1 (MQ=255)
                                                                                                                                   tCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTgag                    <  2:5184384/150‑1 (MQ=255)
                                                                                                                                    cTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTgagg                   <  1:1698629/150‑1 (MQ=255)
                                                                                                                                      aCTGCTGATGCCAACTAAGCGGATGTAGATAACGCAGTCATGACTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGATATGCTGATCTGGTGGAGCAGCACAGTgaggag                 >  1:521026/1‑150 (MQ=255)
                                                                                                                                      aCTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTgaggag                 >  2:3181002/1‑150 (MQ=255)
                                                                                                                                        tgctgATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCt               <  2:423177/150‑1 (MQ=255)
                                                                                                                                         gctgATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTg              <  2:1960619/150‑1 (MQ=255)
                                                                                                                                                ccAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCTATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAAc       >  1:3128372/1‑150 (MQ=255)
                                                                                                                                                 cAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACt      >  2:4522612/1‑150 (MQ=255)
                                                                                                                                                     gAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCTATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGaa  >  2:817631/1‑150 (MQ=255)
                                                                                                                                                     gAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGaa  >  1:3555543/1‑150 (MQ=255)
                                                                                                                                                     gAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGaa  >  1:4776104/1‑150 (MQ=255)
                                                                                                                                                     gAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGaa  >  1:2124806/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAA  >  CP009273/1441792‑1442090

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: