Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,948,854 G→T 10.3% R55L (CGC→CTC)  ilvA → L‑threonine dehydratase, biosynthetic; also known as threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,948,8540GT10.3% 1060.6 / 62.8 291R55L (CGC→CTC) ilvAL‑threonine dehydratase, biosynthetic; also known as threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (130/131);  new base T (9/21);  total (139/152)
Fisher's exact test for biased strand distribution p-value = 5.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 310 total aligned reads displayed.

ACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGG  >  CP009273/3948705‑3949003
                                                                                                                                                     |                                                                                                                                                     
  cccTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCgcc                                                                                                                                                     >  2:3670832/1‑150 (MQ=255)
    cTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCgccag                                                                                                                                                   >  1:1426575/1‑150 (MQ=255)
      gTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCgccagcc                                                                                                                                                 >  1:2102019/1‑150 (MQ=255)
       tCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCgccagcca                                                                                                                                                <  1:309771/150‑1 (MQ=255)
        ccGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGccagccag                                                                                                                                               <  2:3919646/150‑1 (MQ=255)
          ggTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTg                                                                                                                                             <  2:1334625/150‑1 (MQ=255)
          ggTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTg                                                                                                                                             <  2:4057428/150‑1 (MQ=255)
               tCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAg                                                                                                                                        >  1:199547/1‑150 (MQ=255)
                 cGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCt                                                                                                                                      <  1:2799973/150‑1 (MQ=255)
                   gAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCttt                                                                                                                                    <  2:4129196/150‑1 (MQ=255)
                   gAAGGTGCCGAATATTGAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCAGTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCttt                                                                                                                                    <  2:3682827/150‑1 (MQ=255)
                    aaGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTa                                                                                                                                   <  2:2109494/150‑1 (MQ=255)
                    aaGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTa                                                                                                                                   >  2:3535850/1‑150 (MQ=255)
                       gTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGc                                                                                                                                <  2:3339265/150‑1 (MQ=255)
                       gTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGc                                                                                                                                <  2:1981592/150‑1 (MQ=255)
                       gTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGc                                                                                                                                <  2:3525967/150‑1 (MQ=255)
                        tGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCt                                                                                                                               <  1:2593622/150‑1 (MQ=255)
                        tGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCt                                                                                                                               <  2:2807979/150‑1 (MQ=255)
                        tGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCt                                                                                                                               >  2:1497960/1‑150 (MQ=255)
                          ccGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTgc                                                                                                                             <  1:3652/150‑1 (MQ=255)
                          ccGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTgc                                                                                                                             >  2:3889810/1‑150 (MQ=255)
                           cGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTgcg                                                                                                                            >  2:2236904/1‑150 (MQ=255)
                           cGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTgcg                                                                                                                            <  1:573520/150‑1 (MQ=255)
                             aaTATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGACGATCGCCAGCCAGTGCACAGCTTTAAGCTGcgcg                                                                                                                          >  2:1111350/1‑150 (MQ=255)
                               tatTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCgcggc                                                                                                                        <  2:1085875/150‑1 (MQ=255)
                                atTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGACGATCGCCAGCCAGTGCACAGCTTTAAGCTGCgcggcg                                                                                                                       >  2:867142/1‑150 (MQ=255)
                                atTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCgcggcg                                                                                                                       <  1:884738/150‑1 (MQ=255)
                                   tAAGAGCAGTGCTGCGCGCACCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCAt                                                                                                                    >  1:1193886/1‑150 (MQ=255)
                                        gCAGTGCAGCGCGCGCCGGTTTACGAGGCGGCGCAGGTAACGCCGCTACAAAAAATGGAAAAACTGTCCTCGCGTCTAGATAACGTCATTCTGGTGAAGCTCGAAGATCGCCAGCCAGTCCACAGCTATAAGCTGCGCGGCGCAAACGcc                                                                                                               >  1:2186887/1‑150 (MQ=255)
                                          aGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCat                                                                                                             <  2:2570003/150‑1 (MQ=255)
                                            tgctgcGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCatga                                                                                                           <  2:1162883/150‑1 (MQ=255)
                                               tgcGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATgg                                                                                                        <  1:4084405/150‑1 (MQ=255)
                                                gcgcgcgcCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGc                                                                                                       >  1:4393539/1‑150 (MQ=255)
                                                gcgcgcgcCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGc                                                                                                       >  1:2676743/1‑150 (MQ=255)
                                                 cgcgcgcCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCg                                                                                                      <  2:484085/150‑1 (MQ=255)
                                                     cgcCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATTGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGcc                                                                                                  <  2:209816/150‑1 (MQ=255)
                                                     cgcCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGcc                                                                                                  <  2:1286503/150‑1 (MQ=255)
                                                                gAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGgaagaa                                                                                       >  1:4559269/1‑150 (MQ=255)
                                                                   gcggcgCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAg                                                                                    <  1:3844013/150‑1 (MQ=255)
                                                                         cAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAgcg                                                                              >  2:2843695/1‑150 (MQ=255)
                                                                            gTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCAc                                                                           <  1:2558421/150‑1 (MQ=255)
                                                                              tACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACgg                                                                         >  1:2424774/1‑150 (MQ=255)
                                                                                cgccgcTACAAAAAATGGAAACACTGGCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCTGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCg                                                                       <  1:3491068/150‑1 (MQ=255)
                                                                                cgccgcTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCg                                                                       >  2:832339/1‑150 (MQ=255)
                                                                                      tACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCa                                                                 <  1:4247048/150‑1 (MQ=255)
                                                                                      tACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCa                                                                 <  1:3632132/150‑1 (MQ=255)
                                                                                       aCAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCAc                                                                <  1:3865882/150‑1 (MQ=255)
                                                                                         aaaaaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTg                                                              <  2:3432403/150‑1 (MQ=255)
                                                                                         aaaaaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTg                                                              >  1:4112952/1‑150 (MQ=255)
                                                                                          aaaaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGc                                                             <  2:2111636/150‑1 (MQ=255)
                                                                                          aaaaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGc                                                             >  1:7867/1‑150 (MQ=255)
                                                                                            aaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCtt                                                           >  1:2729599/1‑150 (MQ=255)
                                                                                             aaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTc                                                          <  2:2301880/150‑1 (MQ=255)
                                                                                              aTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCt                                                         >  1:3314901/1‑150 (MQ=255)
                                                                                              aTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCt                                                         <  1:4611026/150‑1 (MQ=255)
                                                                                               tGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTg                                                        >  1:4027373/1‑150 (MQ=255)
                                                                                                 gAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCg                                                      >  2:2966540/1‑150 (MQ=255)
                                                                                                   aaaaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCggg                                                    >  2:4382469/1‑150 (MQ=255)
                                                                                                   aaaaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCggg                                                    >  1:698703/1‑150 (MQ=255)
                                                                                                    aaaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGt                                                   <  2:172859/150‑1 (MQ=255)
                                                                                                    aaaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGt                                                   <  2:2748184/150‑1 (MQ=255)
                                                                                                    aaaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGt                                                   <  2:4544511/150‑1 (MQ=255)
                                                                                                     aaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTa                                                  <  2:2081510/150‑1 (MQ=255)
                                                                                                     aaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTa                                                  >  2:3653664/1‑150 (MQ=255)
                                                                                                     aaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTa                                                  >  1:1843730/1‑150 (MQ=255)
                                                                                                     aaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTa                                                  <  2:2610665/150‑1 (MQ=255)
                                                                                                      aCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTaa                                                 >  2:4497161/1‑150 (MQ=255)
                                                                                                         gtcgtcGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCa                                              >  1:1787387/1‑150 (MQ=255)
                                                                                                          tcgtcgCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCAc                                             <  2:1195358/150‑1 (MQ=255)
                                                                                                          tcgtcgCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCAc                                             <  1:2959494/150‑1 (MQ=255)
                                                                                                              cgcgTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACgcgc                                         >  2:1218332/1‑150 (MQ=255)
                                                                                                              cgcgTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACgcgc                                         >  1:1298388/1‑150 (MQ=255)
                                                                                                              cgcgTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACgcgc                                         >  1:1692/1‑150 (MQ=255)
                                                                                                                cgTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAg                                       >  2:2735097/1‑150 (MQ=255)
                                                                                                                cgTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAg                                       >  1:3359542/1‑150 (MQ=255)
                                                                                                                 gTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAgg                                      <  1:4420301/150‑1 (MQ=255)
                                                                                                                 gTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAgg                                      >  2:2662319/1‑150 (MQ=255)
                                                                                                                 gTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAgg                                      >  1:4474322/1‑150 (MQ=255)
                                                                                                                 gTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAgg                                      >  1:4315262/1‑150 (MQ=255)
                                                                                                                 gTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAgg                                      >  1:4315198/1‑150 (MQ=255)
                                                                                                                      gATAACGTCATTCTGGTGTAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGAGGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTc                                 <  1:1038300/150‑1 (MQ=255)
                                                                                                                      gATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTc                                 <  2:611853/150‑1 (MQ=255)
                                                                                                                      gATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTc                                 <  2:1031041/150‑1 (MQ=255)
                                                                                                                       aTAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTcg                                <  1:3942702/150‑1 (MQ=255)
                                                                                                                         aaCGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTcgcg                              >  2:1633758/1‑150 (MQ=255)
                                                                                                                         aaCGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACCGCTTCTGCGGGTAACCACGAGCAGGGCGTcgcg                              <  1:756956/150‑1 (MQ=255)
                                                                                                                           cGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGtt                            >  2:1404465/1‑150 (MQ=255)
                                                                                                                             tCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGtttt                          >  2:4229542/1‑150 (MQ=255)
                                                                                                                             tCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGtttt                          <  2:1517741/150‑1 (MQ=255)
                                                                                                                             tCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGtttt                          <  2:3619616/150‑1 (MQ=255)
                                                                                                                               aTTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTttct                        >  1:3009082/1‑150 (MQ=255)
                                                                                                                               aTTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTttct                        >  2:571845/1‑150 (MQ=255)
                                                                                                                               aTTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTttct                        <  2:3817797/150‑1 (MQ=255)
                                                                                                                                       gAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGcgcg                >  2:50291/1‑150 (MQ=255)
                                                                                                                                         aGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGt              <  2:2945609/150‑1 (MQ=255)
                                                                                                                                           cgcgAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTa            >  2:4603582/1‑150 (MQ=255)
                                                                                                                                            gcgAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAg           >  1:2664042/1‑150 (MQ=255)
                                                                                                                                                  aTCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGa     >  1:2740581/1‑150 (MQ=255)
                                                                                                                                                    cGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAg   >  1:3236173/1‑150 (MQ=255)
                                                                                                                                                     gccagccaGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAgg  >  1:2314823/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGG  >  CP009273/3948705‑3949003

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: