Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,274,674 G→T 7.2% S196I (AGC→ATC)  glcD → glycolate oxidase FAD‑linked subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,274,6740GT 7.2% 303.1 / 8.2 96S196I (AGC→ATC) glcDglycolate oxidase FAD‑linked subunit
Reads supporting (aligned to +/- strand):  ref base G (41/48);  new base T (2/5);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

ACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGATACTGGCCAGTTTCGACAGCGTCGAG  >  NC_002947/4274527‑4274822
                                                                                                                                                   |                                                                                                                                                    
aCCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCg                                                                                                                                                    >  1:1659334/1‑150 (MQ=255)
     tCCTCGCAAATTGCCTGCTCGAGCGGTGGCAACGGCGCCAAAAACGCCCGTACCGTCCTCTGCCTCCAGACCCCCCAGACACTAGACAATGTTCTCCGGGTGGATTTCGCTGTGGTCGATGGGAAACGCCTGTTGCTGGGCAGCGATGcc                                                                                                                                               <  2:276380/150‑1 (MQ=255)
        tCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCATCGCTGcccag                                                                                                                                            >  1:996760/1‑148 (MQ=255)
        tCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCATCGCTGcccag                                                                                                                                            <  2:996760/150‑3 (MQ=255)
        tCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCATCGCTGcccag                                                                                                                                            <  2:1026212/150‑3 (MQ=255)
        tCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTg                                                                                                                                            <  1:250430/150‑1 (MQ=255)
        tCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTg                                                                                                                                            <  2:1659334/150‑1 (MQ=255)
        tCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGCGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTg                                                                                                                                            >  1:1728159/1‑150 (MQ=255)
        tCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCAAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTg                                                                                                                                            >  1:467180/1‑150 (MQ=255)
        tCGCAAATTGCCTGCTCGATCGGTGGCAACGGCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCGCTg                                                                                                                                            >  1:2253463/1‑150 (MQ=255)
         cGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCATCGCTGCCCAgg                                                                                                                                           <  2:511038/150‑1 (MQ=255)
          gCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCATCGCTGCCCAggc                                                                                                                                          <  1:2150420/150‑2 (MQ=255)
          gCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGa                                                                                                                                          <  1:1409301/150‑1 (MQ=255)
          gCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGa                                                                                                                                          <  1:405265/150‑1 (MQ=255)
             aaTTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTCACCCATGTGCTCAAGGTAGCCAGGGCGGTGGCGAAATGCGAACTCCTGAGCCTGGGCAGCGATGCCCTGGACAg                                                                                                                                       <  1:1022615/150‑1 (MQ=255)
                 gCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGcccc                                                                                                                                   >  1:144817/1‑150 (MQ=255)
                 gCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGcccc                                                                                                                                   >  2:2058642/1‑150 (MQ=255)
                   caacTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCgg                                                                                                                                 >  1:84327/4‑150 (MQ=255)
                      cTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCtt                                                                                                                              >  1:1585875/1‑150 (MQ=255)
                         gATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGa                                                                                                                           >  1:544842/1‑150 (MQ=255)
                         gATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGa                                                                                                                           >  2:117712/1‑150 (MQ=255)
                         gATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGa                                                                                                                           <  2:84327/150‑1 (MQ=255)
                          aTCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCTAc                                                                                                                          <  2:1903842/150‑1 (MQ=255)
                          aTCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGAc                                                                                                                          <  2:2138702/150‑1 (MQ=255)
                                ggCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACctgctg                                                                                                                    <  2:1339951/150‑1 (MQ=255)
                                 gCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTgg                                                                                                                   >  2:1733366/1‑150 (MQ=255)
                                     cGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCAtt                                                                                                               <  1:70780/150‑1 (MQ=255)
                                         gCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTc                                                                                                           >  1:895527/1‑150 (MQ=255)
                                         gCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTc                                                                                                           >  1:979251/1‑150 (MQ=255)
                                         gCCGAAAACGCCGGCGGCGTGCATTGCCTCAAGTAGGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACCCCTGAGGCTGGGCAGCGATGCTCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTc                                                                                                           <  2:336233/150‑1 (MQ=255)
                                          ccGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACctgc                                                                                                                      >  2:1001043/1‑138 (MQ=255)
                                          ccGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACctgc                                                                                                                      <  1:1001043/138‑1 (MQ=255)
                                           cGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCAt                                                                                                         <  1:1733366/150‑1 (MQ=255)
                                            gAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCATCGCTGcccagg                                                                                                                                           >  2:2126084/1‑112 (MQ=255)
                                            gAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCATCGCTGCCCAgg                                                                                                                                           <  1:2126084/115‑1 (MQ=255)
                                            gAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATc                                                                                                        >  2:2065921/1‑150 (MQ=255)
                                            gAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATc                                                                                                        <  1:1840862/150‑1 (MQ=255)
                                               aaCGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGc                                                                                                     <  1:1651110/150‑1 (MQ=255)
                                                       tGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAgg                                                                                             <  1:994979/150‑1 (MQ=255)
                                                       tGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAgg                                                                                             <  2:544842/150‑1 (MQ=255)
                                                         gCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCa                                                                                           >  1:1294656/1‑150 (MQ=255)
                                                           gTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATg                                                                                         <  2:895527/150‑1 (MQ=255)
                                                               aCTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCg                                                                                     >  2:2376263/1‑150 (MQ=255)
                                                                cTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCgg                                                                                    >  1:1774436/1‑150 (MQ=255)
                                                                 tGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGt                                                                                   >  2:1704679/1‑150 (MQ=255)
                                                                  gCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTa                                                                                  <  2:1288110/150‑1 (MQ=255)
                                                                  gCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTa                                                                                  <  1:1540319/150‑1 (MQ=255)
                                                                    cTCAAGTACGGCGTGAACGTACACAATGTGCTCAAGGTGGACATCTTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATCGTTCATCGGCTTCGAACGCATGATCGGTATc                                                                                <  2:144817/150‑1 (MQ=255)
                                                                    cTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATc                                                                                <  1:703474/150‑1 (MQ=255)
                                                                    cTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATc                                                                                >  2:2115608/1‑150 (MQ=255)
                                                                     tCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCg                                                                               >  1:2156539/1‑150 (MQ=255)
                                                                     tCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCCGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCg                                                                               >  2:343006/1‑150 (MQ=255)
                                                                      cAAGGACCGCCTCACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGAGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGt                                                                              <  2:467180/150‑1 (MQ=255)
                                                                       aaGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCCTc                                                                             >  2:668471/1‑150 (MQ=255)
                                                                         gTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCAc                                                                           <  2:1774436/150‑1 (MQ=255)
                                                                         gTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATACTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCGCTGGTTGGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATACTCGGTATCGTCAc                                                                           <  2:1417236/150‑1 (MQ=255)
                                                                             ggCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGaa                                                                       <  2:2156539/150‑1 (MQ=255)
                                                                                 tGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCt                                                                                                    <  1:1523517/117‑1 (MQ=255)
                                                                                 tGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCt                                                                                                    >  2:1523517/1‑117 (MQ=255)
                                                                                 tGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCa                                                                   <  1:818869/150‑1 (MQ=255)
                                                                                      gTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTc                                                              >  2:822301/1‑150 (MQ=255)
                                                                                          acaATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCCTTGTTCCTAGGCCGCGAGGGCCATCTCGGATTTAGCATCGAAGTCACCGTCAAgc                                                          >  2:1231218/1‑150 (MQ=255)
                                                                                          acaATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCAtt                                                                                                               <  2:1038280/97‑1 (MQ=255)
                                                                                          acaATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCAtt                                                                                                               >  1:1038280/1‑97 (MQ=255)
                                                                                          acaATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAgc                                                          >  1:515134/1‑150 (MQ=255)
                                                                                            aaTGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCt                                                                                                                                             <  1:2250666/65‑1 (MQ=255)
                                                                                            aaTGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCt                                                                                                                                             >  2:2250666/1‑65 (MQ=255)
                                                                                              tgtgCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAgctgct                                                      >  1:210931/1‑150 (MQ=255)
                                                                                                   tCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCa                                                 >  1:2233289/1‑150 (MQ=255)
                                                                                                     aaGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATTGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAg                                               >  2:2000591/1‑150 (MQ=255)
                                                                                                        gTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGccc                                            <  2:1294656/150‑1 (MQ=255)
                                                                                                             cATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGCCCTGCTGGCATTGTTCCTCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCACGCCCCAGGt                                       <  1:1324377/150‑1 (MQ=255)
                                                                                                              aTCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATTGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTg                                      <  1:2000591/150‑1 (MQ=255)
                                                                                                              aTCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTg                                      <  1:117712/150‑1 (MQ=255)
                                                                                                                 cTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTc                                                                                      <  2:1566203/99‑1 (MQ=255)
                                                                                                                 cTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTc                                                                                      >  1:1566203/1‑99 (MQ=255)
                                                                                                                  ttACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCg                                                                                     >  2:1934948/1‑99 (MQ=255)
                                                                                                                  ttACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCg                                                                                     <  1:1934948/99‑1 (MQ=255)
                                                                                                                  ttACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGgcgc                                  >  1:513361/1‑150 (MQ=255)
                                                                                                                    aCGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCgg                                <  2:1422269/150‑1 (MQ=255)
                                                                                                                    aCGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCgg                                <  2:513361/150‑1 (MQ=255)
                                                                                                                    aCGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCgg                                <  2:979251/150‑1 (MQ=255)
                                                                                                                     cGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCggg                               >  2:550435/1‑150 (MQ=255)
                                                                                                                       gTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTg                             <  1:116524/150‑1 (MQ=255)
                                                                                                                        tCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGa                            >  1:91350/1‑150 (MQ=255)
                                                                                                                            gggCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGATACt                        <  2:2233289/150‑1 (MQ=255)
                                                                                                                                gAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCa                                                                            <  1:890839/94‑1 (MQ=255)
                                                                                                                                gAAAGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCa                                                                            >  2:890839/1‑94 (MQ=255)
                                                                                                                                   cGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGATACTGGCCAGt                 <  1:822301/150‑1 (MQ=255)
                                                                                                                                      cTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGATAc                         <  1:1328042/139‑1 (MQ=255)
                                                                                                                                      cTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGATAc                         >  2:1328042/1‑139 (MQ=255)
                                                                                                                                           ccTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAg                                                                                              <  1:1521874/65‑1 (MQ=255)
                                                                                                                                           ccTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAg                                                                                              >  2:1521874/1‑65 (MQ=255)
                                                                                                                                            cTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCGCCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGATACTGGCCAGTTTCGACAGc        >  2:1665262/1‑150 (MQ=255)
                                                                                                                                                  aGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCGCCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGATACTGGCCAGTTTCGACAGCGTCGAg  <  1:1665262/150‑1 (MQ=255)
                                                                                                                                                  aGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGATACTGGCCAGTTTCGACAGCGTCGAg  <  1:1704679/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCCTTCCTCGCAAATTGCCTGCTCGATCGGTGGCAACGTCGCCGAAAACGCCGGTGGCGTGCACTGCCTCAAGTACGGCCTGACCGTACACAATGTGCTCAAGGTGGACATCCTTACGGTCGAGGGCGAACGCCTGAGCCTGGGCAGCGATGCCCTGGACAGCCCCGGCTTCGACCTGCTGGCATTGTTCATCGGCTCCGAAGGCATGCTCGGTATCGTCACCGAAGTCACCGTCAAGCTGCTGCCCAAGCCCCAGGTGGCGCGGGTGATACTGGCCAGTTTCGACAGCGTCGAG  >  NC_002947/4274527‑4274822

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