Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 123,784 G→T 6.0% D535Y (GAC→TAC)  aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273123,7840GT 6.0% 940.8 / 23.9 249D535Y (GAC→TAC) aceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (143/91);  new base T (5/10);  total (148/101)
Fisher's exact test for biased strand distribution p-value = 5.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 251 total aligned reads displayed.

CAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGG  >  CP009273/123637‑123933
                                                                                                                                                   |                                                                                                                                                     
cAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGAc                                                                                                                                                     <  1:2192937/150‑1 (MQ=255)
  gCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACgg                                                                                                                                                   <  1:1379902/150‑1 (MQ=255)
    gTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTa                                                                                                                                                 >  1:3607836/1‑150 (MQ=255)
    gTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTa                                                                                                                                                 >  1:864017/1‑150 (MQ=255)
    gTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTa                                                                                                                                                 <  2:1031746/150‑1 (MQ=255)
     tCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTaa                                                                                                                                                >  2:3332745/1‑150 (MQ=255)
     tCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTaa                                                                                                                                                >  2:3332408/1‑150 (MQ=255)
        gACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGAGACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCt                                                                                                                                             <  2:171781/150‑1 (MQ=255)
         aCCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTg                                                                                                                                            <  1:1600460/150‑1 (MQ=255)
            cttaagaaATTCATCAACCTCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCTGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAATAAAGCGCGTGACGGTAAGCTGACt                                                                                                                                         <  1:2977102/147‑1 (MQ=255)
            cTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACt                                                                                                                                         <  1:1821052/150‑1 (MQ=255)
               aagaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAATAAAGCGCGTGACGGTAAGCTGACTGCg                                                                                                                                      >  2:1625070/1‑150 (MQ=255)
                agaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCgg                                                                                                                                     <  2:111704/150‑1 (MQ=255)
                 gaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCggg                                                                                                                                    <  1:4521334/150‑1 (MQ=255)
                    aTACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGa                                                                                                                                 >  2:2338233/1‑150 (MQ=255)
                     tACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGaa                                                                                                                                >  2:775917/1‑150 (MQ=255)
                     tACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGaa                                                                                                                                <  1:327592/150‑1 (MQ=255)
                       cATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAAt                                                                                                                              >  2:2980786/1‑150 (MQ=255)
                           aaCATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAg                                                                                                                          <  2:4378529/150‑1 (MQ=255)
                           aaCATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAg                                                                                                                          <  2:563595/150‑1 (MQ=255)
                              aTCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGc                                                                                                                       <  2:2235244/150‑1 (MQ=255)
                               tCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCg                                                                                                                      <  2:1791709/150‑1 (MQ=255)
                               tCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCg                                                                                                                      >  1:1141742/1‑150 (MQ=255)
                               tCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCgg                                                                                                                     <  1:3165738/150‑1 (MQ=255)
                                cGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCgg                                                                                                                     <  1:2376889/150‑1 (MQ=255)
                                  gtgtGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGtt                                                                                                                   <  2:2179423/150‑1 (MQ=255)
                                      ggCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCtt                                                                                                               >  1:3565304/1‑150 (MQ=255)
                                        cGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCa                                                                                                             <  2:800851/150‑1 (MQ=255)
                                         ggTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTTACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCAc                                                                                                            <  2:3759700/150‑1 (MQ=255)
                                           tGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCa                                                                                                          <  1:268299/150‑1 (MQ=255)
                                            ggATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCAt                                                                                                         <  2:1238447/150‑1 (MQ=255)
                                              aTACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCAtct                                                                                                       >  1:4319488/1‑150 (MQ=255)
                                                    cGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCa                                                                                                 >  1:727786/1‑150 (MQ=255)
                                                      aaCGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATc                                                                                               >  1:2040053/1‑150 (MQ=255)
                                                      aaCGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATc                                                                                               >  2:1563083/1‑150 (MQ=255)
                                                      aaCGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATc                                                                                               >  1:542594/1‑150 (MQ=255)
                                                        cGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATcgg                                                                                             >  2:4412414/1‑150 (MQ=255)
                                                         ggTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATcggc                                                                                            >  1:2592403/1‑150 (MQ=255)
                                                         ggTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATcggc                                                                                            <  1:2455351/150‑1 (MQ=255)
                                                         ggTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATcggc                                                                                            >  2:2736913/1‑150 (MQ=255)
                                                            cTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCggcggc                                                                                         >  2:1958488/1‑150 (MQ=255)
                                                             tGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGcc                                                                                        >  1:4189072/1‑150 (MQ=255)
                                                             tGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGcc                                                                                        >  2:725841/1‑150 (MQ=255)
                                                              ggTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCt                                                                                       <  1:188830/150‑1 (MQ=255)
                                                                      cGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTTACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGtacta                                                                               >  2:2937411/1‑150 (MQ=255)
                                                                      cGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTTACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGtacta                                                                               >  2:147210/1‑150 (MQ=255)
                                                                      cGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTTACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGtacta                                                                               >  2:2447961/1‑150 (MQ=255)
                                                                      cGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGtacta                                                                               >  1:3822810/1‑150 (MQ=255)
                                                                        gTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTTACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTAcc                                                                             <  1:2591516/150‑1 (MQ=255)
                                                                         tATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTAccc                                                                            <  1:937649/150‑1 (MQ=255)
                                                                             ctaaGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCATCATCGGCGGCCTGGGTTCTACCCACtt                                                                        <  2:3243343/148‑1 (MQ=255)
                                                                             cAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACtt                                                                        >  2:3527679/1‑150 (MQ=255)
                                                                             cAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACtt                                                                        <  2:3243359/150‑1 (MQ=255)
                                                                             cAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACtt                                                                        >  2:3612320/1‑150 (MQ=255)
                                                                             cAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACtt                                                                        >  2:3527759/1‑150 (MQ=255)
                                                                              aaaGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTTACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTc                                                                       <  2:2418863/150‑1 (MQ=255)
                                                                                aGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTcgc                                                                     >  2:3064865/1‑150 (MQ=255)
                                                                                 gACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTcgcg                                                                    <  2:3391436/150‑1 (MQ=255)
                                                                                    gTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCg                                                                 <  2:890820/150‑1 (MQ=255)
                                                                                    gTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCg                                                                 <  1:3546915/150‑1 (MQ=255)
                                                                                       aacaaGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGAtt                                                              <  1:3839013/150‑1 (MQ=255)
                                                                                          aaGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATtgtg                                                           >  1:3974305/1‑150 (MQ=255)
                                                                                           aGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGa                                                          >  1:1208676/1‑150 (MQ=255)
                                                                                             aaaGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAAc                                                        >  2:185914/1‑150 (MQ=255)
                                                                                             aaaGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAAc                                                        >  2:4364840/1‑150 (MQ=255)
                                                                                              aaGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAAcg                                                       >  1:1095812/1‑150 (MQ=255)
                                                                                              aaGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAAcg                                                       <  2:746359/150‑1 (MQ=255)
                                                                                                   atcatcGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCg                                                  <  2:126299/150‑1 (MQ=255)
                                                                                                   atcatcGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCg                                                  >  1:2435626/1‑150 (MQ=255)
                                                                                                   atcatcGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCg                                                  >  2:2321893/1‑150 (MQ=255)
                                                                                                    tcatcGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCgg                                                 >  2:2880766/1‑150 (MQ=255)
                                                                                                      atcGAGCTGTCTCGCGAGCTGATGACGATTTCTAAGAATGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCAACTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGaa                                               >  1:4442070/1‑150 (MQ=255)
                                                                                                        cGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGt                                             <  2:1834738/150‑1 (MQ=255)
                                                                                                        cGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGt                                             >  2:4468981/1‑150 (MQ=255)
                                                                                                             tGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTa                                        >  2:248152/1‑150 (MQ=255)
                                                                                                             tGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTa                                        <  2:1325962/150‑1 (MQ=255)
                                                                                                             tGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTa                                        >  2:1218135/1‑150 (MQ=255)
                                                                                                                 tcGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCt                                    >  2:880335/1‑150 (MQ=255)
                                                                                                                  cgcgAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTc                                   >  1:2556789/1‑150 (MQ=255)
                                                                                                                    cgAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCgg                                 >  1:218521/1‑150 (MQ=255)
                                                                                                                      aGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCg                               >  2:4232789/1‑150 (MQ=255)
                                                                                                                        cTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGtt                             <  2:1494662/150‑1 (MQ=255)
                                                                                                                         tgatgaCTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGttt                            >  2:2028997/1‑150 (MQ=255)
                                                                                                                          gatgaCTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTc                           >  2:4078672/1‑150 (MQ=255)
                                                                                                                          gatgaCTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTc                           >  2:2671998/1‑150 (MQ=255)
                                                                                                                                    tCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCcgcg                 >  2:3022281/1‑150 (MQ=255)
                                                                                                                                      tAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGAt               >  1:594174/1‑150 (MQ=255)
                                                                                                                                       aaGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATg              >  1:3400150/1‑150 (MQ=255)
                                                                                                                                         gAAAGCGCGTTACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGa            <  2:274812/150‑1 (MQ=255)
                                                                                                                                         gAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGa            >  1:1158450/1‑150 (MQ=255)
                                                                                                                                         gAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGa            <  2:3095790/150‑1 (MQ=255)
                                                                                                                                         gAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGa            >  1:3220262/1‑150 (MQ=255)
                                                                                                                                           aaGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGc          >  1:2089712/1‑150 (MQ=255)
                                                                                                                                           aaGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGc          >  1:4226791/1‑150 (MQ=255)
                                                                                                                                             gcgcGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCg        <  2:2127244/150‑1 (MQ=255)
                                                                                                                                             gcgcGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCg        <  2:2127271/150‑1 (MQ=255)
                                                                                                                                                 gTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGgtgt    >  1:4505031/1‑150 (MQ=255)
                                                                                                                                                  tGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGtgtg   <  2:3565299/150‑1 (MQ=255)
                                                                                                                                                  tGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGtgtg   >  2:1493060/1‑150 (MQ=255)
                                                                                                                                                   gACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTgg  >  1:2853438/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGG  >  CP009273/123637‑123933

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: