Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,148,927 G→T 20.1% E72* (GAA→TAA)  yqhD → aldehyde reductase, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,148,9270GT20.1% 789.8 / 150.2 287E72* (GAA→TAA) yqhDaldehyde reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (104/125);  new base T (29/29);  total (133/154)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

TACGCGAACAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGCACATTCTGCAA  >  CP009273/3148778‑3149076
                                                                                                                                                     |                                                                                                                                                     
tACGCGAACAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATg                                                                                                                                                       >  1:2768265/1‑150 (MQ=255)
    cgAACAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAAc                                                                                                                                                   >  1:2114018/1‑150 (MQ=255)
    cgAACAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAAc                                                                                                                                                   <  2:3100622/150‑1 (MQ=255)
    cgAACAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAAc                                                                                                                                                   <  2:3100744/150‑1 (MQ=255)
    cgAACAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAAc                                                                                                                                                   >  2:1446205/1‑150 (MQ=255)
      aaCAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGc                                                                                                                                                 >  1:1935862/1‑150 (MQ=255)
       aCAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCt                                                                                                                                                <  2:2280142/150‑1 (MQ=255)
       aCAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCt                                                                                                                                                <  2:3212517/150‑1 (MQ=255)
         aaaGTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCtga                                                                                                                                              <  1:542781/150‑1 (MQ=255)
           aTTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCtgatg                                                                                                                                            >  2:839141/1‑150 (MQ=255)
               cTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACg                                                                                                                                        >  2:3981953/1‑150 (MQ=255)
               cTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACg                                                                                                                                        <  2:927168/150‑1 (MQ=255)
               cTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACg                                                                                                                                        <  2:2116892/150‑1 (MQ=255)
                    gATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTg                                                                                                                                   >  1:1628005/1‑150 (MQ=255)
                     aTGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGa                                                                                                                                  >  2:2609820/1‑150 (MQ=255)
                      tGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGaa                                                                                                                                 >  1:2507999/1‑150 (MQ=255)
                         tCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACt                                                                                                                              <  1:1958822/150‑1 (MQ=255)
                          cgcgTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGGCGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTg                                                                                                                             >  2:4461127/1‑150 (MQ=255)
                          cgcgTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTg                                                                                                                             <  2:2067612/150‑1 (MQ=255)
                           gcgTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTgg                                                                                                                            >  1:3310452/1‑150 (MQ=255)
                           gcgTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTgg                                                                                                                            <  2:2834992/150‑1 (MQ=255)
                               aTTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTcg                                                                                                                        <  2:1295339/150‑1 (MQ=255)
                               aTTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTcg                                                                                                                        <  2:659680/150‑1 (MQ=255)
                                ttGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTcgc                                                                                                                       >  1:2175817/1‑150 (MQ=255)
                                 tGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTcgcg                                                                                                                      >  1:10238/1‑150 (MQ=255)
                                    ttACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAAc                                                                                                                   <  1:3824983/150‑1 (MQ=255)
                                     tACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACa                                                                                                                  >  2:357313/1‑150 (MQ=255)
                                      aCCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAg                                                                                                                 <  1:4452821/150‑1 (MQ=255)
                                       ccTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGa                                                                                                                >  2:1671946/1‑150 (MQ=255)
                                       ccTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGa                                                                                                                >  2:1287616/1‑150 (MQ=255)
                                          aCGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAg                                                                                                             <  1:1266552/150‑1 (MQ=255)
                                            ggcggcggcAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTg                                                                                                           >  2:1069027/1‑150 (MQ=255)
                                             gcggcggcAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGa                                                                                                          <  2:2582694/150‑1 (MQ=255)
                                                  ggcATCGTTAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTc                                                                                                     <  2:1968/150‑1 (MQ=255)
                                                   gcAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGGCTTTcc                                                                                                    <  2:1843633/150‑1 (MQ=255)
                                                   gcAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTcc                                                                                                    <  2:3299139/150‑1 (MQ=255)
                                                    cAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCct                                                                                                   <  2:2101063/150‑1 (MQ=255)
                                                    cAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCct                                                                                                   <  2:3479900/150‑1 (MQ=255)
                                                            aaaaaaCCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCgg                                                                                           >  2:3171659/1‑150 (MQ=255)
                                                              aaaaCCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGtt                                                                                         >  2:381070/1‑150 (MQ=255)
                                                                 aCCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTggc                                                                                      >  2:2992086/1‑150 (MQ=255)
                                                                   cGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTggcgg                                                                                    <  2:686647/150‑1 (MQ=255)
                                                                   cGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTggcgg                                                                                    <  2:120871/150‑1 (MQ=255)
                                                                      cGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGcggcgg                                                                                 >  1:1875400/1‑150 (MQ=255)
                                                                       gTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGt                                                                                <  2:1750115/150‑1 (MQ=255)
                                                                          ctcGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCt                                                                             <  2:4274891/150‑1 (MQ=255)
                                                                            cGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGt                                                                           >  1:286044/1‑150 (MQ=255)
                                                                             gATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGTCTTATGAAAAGCTGATGAACGCCGTGAAACTGGTTCGAGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTa                                                                          >  2:3440713/1‑150 (MQ=255)
                                                                             gATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTa                                                                          >  2:2755759/1‑150 (MQ=255)
                                                                                cAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTg                                                                       >  2:371601/1‑150 (MQ=255)
                                                                                 aaGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTgg                                                                      >  2:3436604/1‑150 (MQ=255)
                                                                                  aGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGa                                                                     <  2:3890812/150‑1 (MQ=255)
                                                                                      cTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGc                                                                 >  1:880503/1‑150 (MQ=255)
                                                                                       tGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCa                                                                >  1:91992/1‑150 (MQ=255)
                                                                                        ggATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGACACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCAc                                                               >  2:2918653/1‑150 (MQ=255)
                                                                                         gATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCAcc                                                              <  1:2700425/150‑1 (MQ=255)
                                                                                          aTGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCa                                                             <  1:1451730/150‑1 (MQ=255)
                                                                                             cccTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAAt                                                          <  2:364487/150‑1 (MQ=255)
                                                                                                tGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTa                                                       >  1:2457189/1‑150 (MQ=255)
                                                                                                 gAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTAt                                                      <  2:503309/150‑1 (MQ=255)
                                                                                                  aaaGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATc                                                     <  2:774380/150‑1 (MQ=255)
                                                                                                   aaGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATcg                                                    <  1:2802315/150‑1 (MQ=255)
                                                                                                   aaGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATcg                                                    >  1:3497638/1‑150 (MQ=255)
                                                                                                     ggCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATcgcc                                                  <  1:2957050/150‑1 (MQ=255)
                                                                                                         tGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAg                                              <  1:97545/150‑1 (MQ=255)
                                                                                                            aCGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCgg                                           <  2:2186369/150‑1 (MQ=255)
                                                                                                             cGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGc                                          >  1:1289021/1‑150 (MQ=255)
                                                                                                             cGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGc                                          >  2:2333069/1‑150 (MQ=255)
                                                                                                               tGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTa                                        >  1:1896470/1‑150 (MQ=255)
                                                                                                               tGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTa                                        <  1:661740/150‑1 (MQ=255)
                                                                                                                  tGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACt                                     >  1:3723454/1‑150 (MQ=255)
                                                                                                                  tGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACt                                     >  1:2408865/1‑150 (MQ=255)
                                                                                                                   ggAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTa                                    <  1:2340244/150‑1 (MQ=255)
                                                                                                                      aTTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATcc                                 >  2:2934812/1‑150 (MQ=255)
                                                                                                                         tGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGa                              >  1:2834408/1‑150 (MQ=255)
                                                                                                                          ggCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGaa                             >  1:1523817/1‑150 (MQ=255)
                                                                                                                             ggTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAat                          <  1:445138/150‑1 (MQ=255)
                                                                                                                               tATTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAatat                        <  2:1369933/150‑1 (MQ=255)
                                                                                                                               tATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAatat                        >  1:3443830/1‑150 (MQ=255)
                                                                                                                               tATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAatat                        <  1:2637909/150‑1 (MQ=255)
                                                                                                                                aTTGAGCCAAACCCGGCTTATTAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATc                       <  1:3601675/150‑1 (MQ=255)
                                                                                                                                aTTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATc                       >  1:4101803/1‑150 (MQ=255)
                                                                                                                                  tGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGa                     >  1:2667122/1‑150 (MQ=255)
                                                                                                                                   gAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGAt                    <  2:1550404/150‑1 (MQ=255)
                                                                                                                                   gAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGAt                    <  1:810347/150‑1 (MQ=255)
                                                                                                                                    aGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATc                   >  1:2273951/1‑150 (MQ=255)
                                                                                                                                     gCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTTCTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATcc                  <  1:4488054/150‑1 (MQ=255)
                                                                                                                                        aaaCCCGGCTTATTAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTg               <  2:2506236/150‑1 (MQ=255)
                                                                                                                                          aCCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGc             <  1:3801657/150‑1 (MQ=255)
                                                                                                                                          aCCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGc             >  1:3987451/1‑150 (MQ=255)
                                                                                                                                           cccGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGca            >  2:2626093/1‑150 (MQ=255)
                                                                                                                                            ccGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGcac           >  2:4393819/1‑150 (MQ=255)
                                                                                                                                            ccGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGcac           >  1:1746333/1‑150 (MQ=255)
                                                                                                                                                cTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGCACATTc       <  1:4216156/150‑1 (MQ=255)
                                                                                                                                                 ttattaAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGCACAGTcc      <  2:3386889/150‑2 (MQ=255)
                                                                                                                                                  tATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGCACATTCTg     <  1:2566241/150‑1 (MQ=255)
                                                                                                                                                   aTGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGCACATTCTGc    >  2:913839/1‑150 (MQ=255)
                                                                                                                                                     gAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGCACATTCTGCaa  >  1:1804812/1‑150 (MQ=255)
                                                                                                                                                     gAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGCACATTCTGCaa  <  1:2235876/150‑1 (MQ=255)
                                                                                                                                                     gAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGCACATTCTGCaa  >  1:2153173/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGCGAACAAATTCCTCACGATGCTCGCGTATTGATTACCTACGGCGGCGGCAGCGTGAAAAAAACCGGCGTTCTCGATCAAGTTCTGGATGCCCTGAAAGGCATGGACGTGCTGGAATTTGGCGGTATTGAGCCAAACCCGGCTTATGAAACGCTGATGAACGCCGTGAAACTGGTTCGCGAACAGAAAGTGACTTTCCTGCTGGCGGTTGGCGGCGGTTCTGTACTGGACGGCACCAAATTTATCGCCGCAGCGGCTAACTATCCGGAAAATATCGATCCGTGGCACATTCTGCAA  >  CP009273/3148778‑3149076

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: