Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,120,216 C→A 62.2% 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,2160CA62.2% 180.2 / inf 260D389Y (GAT→TAT) glcDglycolate oxidase subunit, FAD‑linked
Reads supporting (aligned to +/- strand):  ref base C (53/45);  new base A (105/57);  total (158/102)
Fisher's exact test for biased strand distribution p-value = 9.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 263 total aligned reads displayed.

TTTTTTCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAACGC  >  CP009273/3120069‑3120363
                                                                                                                                                   |                                                                                                                                                   
    ttCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGc                                                                                                                                               <  2:1774182/150‑1 (MQ=255)
     tctcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCa                                                                                                                                              <  2:4661026/150‑1 (MQ=255)
      ctcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCAt                                                                                                                                             <  1:454503/150‑1 (MQ=255)
      ctcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCAt                                                                                                                                             <  1:3157677/150‑1 (MQ=255)
      ctcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCAt                                                                                                                                             >  2:1910384/1‑150 (MQ=255)
      ctcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCAt                                                                                                                                             >  1:1495462/1‑150 (MQ=255)
      ctcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCAt                                                                                                                                             <  1:1323482/150‑1 (MQ=255)
       tcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATg                                                                                                                                            >  1:4348881/1‑150 (MQ=255)
       tcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATg                                                                                                                                            >  1:303835/1‑150 (MQ=255)
       tcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATg                                                                                                                                            >  1:2408903/1‑150 (MQ=255)
         gCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGaa                                                                                                                                          <  2:4281271/150‑1 (MQ=255)
          cccGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGaaa                                                                                                                                         <  1:2171999/150‑1 (MQ=255)
             gATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGAc                                                                                                                                      <  2:1511067/150‑1 (MQ=255)
              aTGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACg                                                                                                                                     <  2:4375286/150‑1 (MQ=255)
                gCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGtt                                                                                                                                   <  1:4603083/150‑1 (MQ=255)
                gCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGtt                                                                                                                                   >  1:1874224/1‑150 (MQ=255)
                 ccATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTg                                                                                                                                  >  2:875645/1‑150 (MQ=255)
                    tGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCa                                                                                                                               >  1:1669418/1‑150 (MQ=255)
                      ttCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAac                                                                                                                             >  2:3495863/1‑150 (MQ=255)
                           cACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTa                                                                                                                        >  2:2251897/1‑150 (MQ=255)
                           cACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTa                                                                                                                        >  2:4788781/1‑150 (MQ=255)
                              tGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAAt                                                                                                                     <  2:1536634/150‑1 (MQ=255)
                                aTGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCCTGAAAGACGTTGGCAACACGGAAATCa                                                                                                                   >  2:2505892/1‑150 (MQ=255)
                                aTGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCa                                                                                                                   >  1:1399235/1‑150 (MQ=255)
                                aTGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCa                                                                                                                   >  2:3257030/1‑150 (MQ=255)
                                   ctgcCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCatat                                                                                                                >  2:417541/1‑150 (MQ=255)
                                     gccgccAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATtg                                                                                                              <  1:3349699/150‑1 (MQ=255)
                                      ccgccAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATtgc                                                                                                             >  2:697767/1‑150 (MQ=255)
                                          cAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCg                                                                                                         <  2:4455279/150‑1 (MQ=255)
                                           aaCTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGa                                                                                                        <  2:605011/150‑1 (MQ=255)
                                                cAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAAc                                                                                                   >  2:1112864/1‑150 (MQ=255)
                                                 aaCGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAACTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACg                                                                                                  <  2:4056843/150‑1 (MQ=255)
                                                   cGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACggg                                                                                                >  1:1199822/1‑150 (MQ=255)
                                                     cAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCa                                                                                              <  1:3426457/150‑1 (MQ=255)
                                                      agagTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCaa                                                                                             >  1:408480/1‑150 (MQ=255)
                                                       gagTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAAt                                                                                            <  1:1833955/150‑1 (MQ=255)
                                                                gATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAAGTCACCGGGTTCGTTGGCATCGAATAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTTTATCATATTGCTGCGATAAACGGGCAATGCCTTCCAg                                                                                   <  2:520090/150‑1 (MQ=255)
                                                                 aTCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGt                                                                                  <  1:1108065/150‑1 (MQ=255)
                                                                    ttCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACg                                                                               >  2:2705575/1‑150 (MQ=255)
                                                                       ccGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCca                                                                            <  2:2168276/150‑1 (MQ=255)
                                                                       ccGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCca                                                                            <  1:1587320/150‑1 (MQ=255)
                                                                          cccAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcaggc                                                                         <  1:697767/150‑1 (MQ=255)
                                                                           ccAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcaggca                                                                        >  1:3214354/1‑150 (MQ=255)
                                                                           ccAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcaggca                                                                        >  2:3114986/1‑150 (MQ=255)
                                                                           ccAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcaggca                                                                        >  2:201913/1‑150 (MQ=255)
                                                                            cAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcaggcag                                                                       <  2:3760296/150‑1 (MQ=255)
                                                                             aGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcaggcagg                                                                      <  2:286605/150‑1 (MQ=255)
                                                                               ctctTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGgc                                                                    >  1:4384735/1‑150 (MQ=255)
                                                                                tctTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGgcg                                                                   <  1:3314556/150‑1 (MQ=255)
                                                                                     cgcgcgGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGAcg                                                              >  1:854978/1‑150 (MQ=255)
                                                                                        gcgGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCgg                                                           <  2:2603322/150‑1 (MQ=255)
                                                                                          gggCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGa                                                         >  2:1961553/1‑150 (MQ=255)
                                                                                            gCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATg                                                       >  2:1435587/1‑150 (MQ=255)
                                                                                              aaaTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGt                                                     >  1:2677155/1‑150 (MQ=255)
                                                                                               aaTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTg                                                    >  2:403599/1‑150 (MQ=255)
                                                                                                aTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGc                                                   >  2:3100266/1‑150 (MQ=255)
                                                                                                   cACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCCGGGCGCGACGCGGGATGGTGccat                                                >  2:3800567/1‑150 (MQ=255)
                                                                                                      cGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGccatcca                                             <  1:2343323/150‑1 (MQ=255)
                                                                                                      cGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGccatcca                                             >  1:2922588/1‑150 (MQ=255)
                                                                                                       gggTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGccatccat                                            <  2:639215/150‑1 (MQ=255)
                                                                                                         gTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGc                                          <  2:3049318/150‑1 (MQ=255)
                                                                                                          ttCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCa                                         >  1:179284/1‑150 (MQ=255)
                                                                                                            cGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCagt                                       >  1:3245485/1‑150 (MQ=255)
                                                                                                             gTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCagta                                      >  1:3498051/1‑150 (MQ=255)
                                                                                                               tGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCagtagt                                    <  2:2328445/150‑1 (MQ=255)
                                                                                                                    tCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATcc                               >  2:3980456/1‑150 (MQ=255)
                                                                                                                     cGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCg                              >  2:2833962/1‑150 (MQ=255)
                                                                                                                     cGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCg                              >  2:854437/1‑150 (MQ=255)
                                                                                                                      gAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCgg                             >  2:2915344/1‑150 (MQ=255)
                                                                                                                       aaaaGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCggg                            <  1:1628663/150‑1 (MQ=255)
                                                                                                                       aaaaGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCggg                            >  1:3716639/1‑150 (MQ=255)
                                                                                                                       aaaaGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCggg                            >  2:867727/1‑150 (MQ=255)
                                                                                                                           ggATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGgaga                        >  1:18148/1‑150 (MQ=255)
                                                                                                                            gATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGAt                       <  2:1860833/150‑1 (MQ=255)
                                                                                                                             atgaaCGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATa                      <  2:3470402/147‑1 (MQ=255)
                                                                                                                             aTTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATa                      >  2:4015158/1‑150 (MQ=255)
                                                                                                                             aTTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATa                      >  2:3621931/1‑150 (MQ=255)
                                                                                                                              ttAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATAc                     <  2:2647111/150‑1 (MQ=255)
                                                                                                                              ttAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATAGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATAc                     >  1:477606/1‑150 (MQ=255)
                                                                                                                               tAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACg                    >  1:3327195/1‑150 (MQ=255)
                                                                                                                               tAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACg                    >  1:1483148/1‑150 (MQ=255)
                                                                                                                                 aCGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTc                  >  2:4830628/1‑150 (MQ=255)
                                                                                                                                  cGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTcc                 >  2:3435687/1‑150 (MQ=255)
                                                                                                                                   gggTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCt                <  2:3260683/150‑1 (MQ=255)
                                                                                                                                       gCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccg            >  1:3881500/1‑150 (MQ=255)
                                                                                                                                         aTGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccgcc          >  1:649628/1‑150 (MQ=255)
                                                                                                                                          tGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccgccg         <  1:528997/150‑1 (MQ=255)
                                                                                                                                           gTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCgg        >  1:2549672/1‑150 (MQ=255)
                                                                                                                                            ttGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCggg       >  1:3740508/1‑150 (MQ=255)
                                                                                                                                             tGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGa      >  1:97677/1‑150 (MQ=255)
                                                                                                                                             tGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGa      >  1:3917307/1‑150 (MQ=255)
                                                                                                                                             tGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGa      >  2:4608803/1‑150 (MQ=255)
                                                                                                                                             tGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGa      >  1:1667360/1‑150 (MQ=255)
                                                                                                                                             tGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCGGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGa      >  2:1176861/1‑150 (MQ=255)
                                                                                                                                              gCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGaa     >  2:1175995/1‑150 (MQ=255)
                                                                                                                                              gCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGaa     <  2:4141941/150‑1 (MQ=255)
                                                                                                                                              gCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGaa     <  2:1687975/150‑1 (MQ=255)
                                                                                                                                              gCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGaa     >  1:1031656/1‑150 (MQ=255)
                                                                                                                                               ccATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAAc    >  2:1305815/1‑150 (MQ=255)
                                                                                                                                               ccATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAAc    >  2:1305838/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTTTTCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAACGC  >  CP009273/3120069‑3120363

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: