Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,120,216 C→A 92.4% D389Y (GAT→TAT)  glcD ← glycolate oxidase subunit, FAD‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,120,2160CA92.4% 761.2 / 48.8 260D389Y (GAT→TAT) glcDglycolate oxidase subunit, FAD‑linked
Reads supporting (aligned to +/- strand):  ref base C (10/10);  new base A (100/140);  total (110/150)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

GATTTTTTCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAAC  >  CP009273/3120067‑3120361
                                                                                                                                                     |                                                                                                                                                 
   tttttCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATcc                                                                                                                                                <  1:4456157/150‑1 (MQ=255)
   tttttCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATcc                                                                                                                                                <  2:3863087/150‑1 (MQ=255)
     tttCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCgg                                                                                                                                              <  2:3609009/150‑1 (MQ=255)
        ctcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCAt                                                                                                                                           >  1:20858/1‑150 (MQ=255)
           gCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGaa                                                                                                                                        <  1:4222531/150‑1 (MQ=255)
             ccGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAg                                                                                                                                      <  1:2363854/150‑1 (MQ=255)
               gATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGAc                                                                                                                                    >  2:1323101/1‑150 (MQ=255)
               gATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGAc                                                                                                                                    <  1:2144222/150‑1 (MQ=255)
                  gCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGtt                                                                                                                                 <  2:2211296/150‑1 (MQ=255)
                  gCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGtt                                                                                                                                 <  1:2442437/150‑1 (MQ=255)
                  gCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGtt                                                                                                                                 <  2:4473692/150‑1 (MQ=255)
                   ccATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTg                                                                                                                                >  1:4377607/1‑150 (MQ=255)
                     aTGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGTAAATTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACATGGTTCGTTGTCATCGAAAAGGATAAACGGGTGCATGAATCCATATCCTGCATGAAAGACTATGGc                                                                                                                              >  2:3976536/1‑150 (MQ=255)
                      tGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCa                                                                                                                             >  2:3371226/1‑150 (MQ=255)
                       gTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCaa                                                                                                                            >  2:2842306/1‑150 (MQ=255)
                       gTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCaa                                                                                                                            <  2:2372657/150‑1 (MQ=255)
                         tCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAaca                                                                                                                          <  2:163618/150‑1 (MQ=255)
                         tCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAaca                                                                                                                          <  2:2462419/150‑1 (MQ=255)
                           gCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACg                                                                                                                        <  1:614912/150‑1 (MQ=255)
                           gCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACg                                                                                                                        <  1:3131934/150‑1 (MQ=255)
                            ccACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGt                                                                                                                       >  2:1840912/1‑150 (MQ=255)
                               cTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTaaa                                                                                                                    >  1:828780/1‑150 (MQ=255)
                                tGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAAt                                                                                                                   <  2:3620037/150‑1 (MQ=255)
                                  aTGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCa                                                                                                                 <  2:4163659/150‑1 (MQ=255)
                                   tgctgcCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCat                                                                                                                <  2:1076527/150‑1 (MQ=255)
                                     ctgcCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCatat                                                                                                              <  2:3651901/150‑1 (MQ=255)
                                      tgcCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATAtt                                                                                                             <  2:50544/150‑1 (MQ=255)
                                      tgcCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATAtt                                                                                                             <  1:4233864/150‑1 (MQ=255)
                                         caccAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATtgct                                                                                                          <  1:2464946/148‑1 (MQ=255)
                                             aaCTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGa                                                                                                      <  1:1729098/150‑1 (MQ=255)
                                               cTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATa                                                                                                    <  2:3035082/150‑1 (MQ=255)
                                               cTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATa                                                                                                    <  1:1005510/150‑1 (MQ=255)
                                                   aaCGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACg                                                                                                >  1:122262/1‑150 (MQ=255)
                                                        agagTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCaa                                                                                           >  2:1495172/1‑150 (MQ=255)
                                                          agTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATg                                                                                         <  2:2312117/150‑1 (MQ=255)
                                                          agTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATg                                                                                         <  2:1023496/150‑1 (MQ=255)
                                                          agTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATg                                                                                         <  1:3957608/150‑1 (MQ=255)
                                                           gTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGc                                                                                        <  1:3047409/150‑1 (MQ=255)
                                                            ttCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGcc                                                                                       <  1:3767239/150‑1 (MQ=255)
                                                            ttCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGcc                                                                                       <  2:1592848/150‑1 (MQ=255)
                                                                 ggATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCa                                                                                  <  1:101363/150‑1 (MQ=255)
                                                                 ggATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCa                                                                                  >  2:2994299/1‑150 (MQ=255)
                                                                    tCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTa                                                                               >  2:2615667/1‑150 (MQ=255)
                                                                     cTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGAGTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTAc                                                                              >  2:3071153/1‑150 (MQ=255)
                                                                        cccGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGcc                                                                           <  1:3904010/150‑1 (MQ=255)
                                                                        cccGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGcc                                                                           <  1:1646211/150‑1 (MQ=255)
                                                                        cccGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCGTCCAGTACGcc                                                                           >  1:1439439/1‑150 (MQ=255)
                                                                          cGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcag                                                                         >  2:3305030/1‑150 (MQ=255)
                                                                                    ttCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGcgcg                                                               >  2:3649186/1‑150 (MQ=255)
                                                                                      ccGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGAc                                                             <  2:2406834/150‑1 (MQ=255)
                                                                                      ccGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGAc                                                             <  1:3643409/150‑1 (MQ=255)
                                                                                       cgcgcgGGCAACTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGAGCGAcg                                                            <  1:1751624/150‑1 (MQ=255)
                                                                                       cgcgcgGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGAcg                                                            >  1:3509650/1‑150 (MQ=255)
                                                                                       cgcgcgGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGAcg                                                            >  1:4418167/1‑150 (MQ=255)
                                                                                       cgcgcgGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGAGCGAcg                                                            <  1:1822191/150‑1 (MQ=255)
                                                                                          gcgGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCgg                                                         <  2:884828/150‑1 (MQ=255)
                                                                                           cgGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCggg                                                        <  2:3180147/150‑1 (MQ=255)
                                                                                           cgGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCggg                                                        <  1:2388176/150‑1 (MQ=255)
                                                                                           cgGGCAAATTCACCGGGTACGTGGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCggg                                                        <  2:640408/150‑1 (MQ=255)
                                                                                            gggCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGTCGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGa                                                       <  2:763644/150‑1 (MQ=255)
                                                                                            gggCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGa                                                       <  2:2564835/150‑1 (MQ=255)
                                                                                             ggCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCCGGCAGGGCGCGACGCGGGAt                                                      >  2:4074219/1‑150 (MQ=255)
                                                                                             ggCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGAt                                                      <  2:1921071/150‑1 (MQ=255)
                                                                                             ggCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGAt                                                      <  1:2892006/150‑1 (MQ=255)
                                                                                                 aaTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTg                                                  >  2:1056022/1‑150 (MQ=255)
                                                                                                 aaTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTg                                                  <  2:2893943/150‑1 (MQ=255)
                                                                                                    tCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGcca                                               >  1:4034883/1‑150 (MQ=255)
                                                                                                     cACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGccat                                              <  2:2793003/150‑1 (MQ=255)
                                                                                                     cACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGccat                                              <  1:502057/150‑1 (MQ=255)
                                                                                                      aCCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGccatc                                             >  2:517475/1‑150 (MQ=255)
                                                                                                      aCCGGGTTCGTTGGCATCGAAAAGGATTAACGGGAGCATGTTGCCATAACCGGCATGACAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCACTGCCTTCCAGTCCGCCAGGCAGGGCGCGACGCGGGATGGTGccatc                                             >  2:4286395/1‑150 (MQ=255)
                                                                                                        ctgggtCGTGGGCATCGAAAAGGTTAAACGCGTGCACGGTCCCATATCCGGCATGTAAGACGTGGGCAACACGTAAATCATTTTGCTGCGATTACCGGGCAAGGCCTTCCAGCACGCCAGGCAGGGCGAGACGCGGGATGGTGccatcca                                           <  2:3990169/145‑1 (MQ=255)
                                                                                                          ggTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATg                                         >  1:2484473/1‑150 (MQ=255)
                                                                                                            ttCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCa                                       >  2:2118144/1‑150 (MQ=255)
                                                                                                             tCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCag                                      <  1:2273945/150‑1 (MQ=255)
                                                                                                                 tGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGTCATGATAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCagtagt                                  <  1:97264/150‑1 (MQ=255)
                                                                                                                 tGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCagtagt                                  <  1:2397819/150‑1 (MQ=255)
                                                                                                                  ggCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAgtagta                                 <  2:4449530/150‑1 (MQ=255)
                                                                                                                   gCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTaa                                >  1:3090155/1‑150 (MQ=255)
                                                                                                                     aTCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATc                              <  2:2637837/150‑1 (MQ=255)
                                                                                                                       cGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCg                            >  1:3269230/1‑150 (MQ=255)
                                                                                                                       cGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCg                            >  1:3269285/1‑150 (MQ=255)
                                                                                                                       cGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCg                            <  1:3631495/150‑1 (MQ=255)
                                                                                                                        gAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCgg                           >  2:3457774/1‑150 (MQ=255)
                                                                                                                        gAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCgg                           >  1:439189/1‑150 (MQ=255)
                                                                                                                         aaaaGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCggg                          >  1:2371319/1‑150 (MQ=255)
                                                                                                                          aaaGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCgggg                         <  1:224123/150‑1 (MQ=255)
                                                                                                                           aaGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGga                        <  1:2162648/150‑1 (MQ=255)
                                                                                                                            aGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGgag                       >  1:2598754/1‑150 (MQ=255)
                                                                                                                              gATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGAt                     <  1:2894964/150‑1 (MQ=255)
                                                                                                                              gATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGAt                     >  2:672907/1‑150 (MQ=255)
                                                                                                                                 tAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACg                  <  2:4362804/150‑1 (MQ=255)
                                                                                                                                    cGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTcc               <  1:3988671/150‑1 (MQ=255)
                                                                                                                                      ggTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTa             <  1:2402325/150‑1 (MQ=255)
                                                                                                                                       gTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAc            >  2:3493328/1‑150 (MQ=255)
                                                                                                                                       gTGCATGATGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAc            >  2:1998609/1‑150 (MQ=255)
                                                                                                                                        tGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAcc           >  2:1735353/1‑150 (MQ=255)
                                                                                                                                          cATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccgc         >  1:1346150/1‑150 (MQ=255)
                                                                                                                                           aTGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccgcc        >  2:3775826/1‑150 (MQ=255)
                                                                                                                                                 ccATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAAc  >  1:967917/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GATTTTTTCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAAC  >  CP009273/3120067‑3120361

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: