Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,750,574 T→G 5.6% L207R (CTG→CGG)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,750,5740TG 5.6% 566.3 / 24.1 196L207R (CTG→CGG) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (84/101);  new base G (6/5);  total (90/106)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 199 total aligned reads displayed.

TGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGT  >  CP009273/1750427‑1750721
                                                                                                                                                   |                                                                                                                                                   
   aCCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGaa                                                                                                                                                  <  2:2238616/148‑1 (MQ=255)
    cctgcctgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGaaa                                                                                                                                                 <  1:3391524/148‑1 (MQ=255)
    cctgcctgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGaaa                                                                                                                                                 >  1:2175284/1‑148 (MQ=255)
     ctgcctgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAg                                                                                                                                                <  1:1465746/148‑1 (MQ=255)
      tgcctgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTTCTGCTTCCTTTATTAGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAgc                                                                                                                                               >  1:2708429/1‑148 (MQ=255)
       gcctgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAgcg                                                                                                                                              >  2:479189/1‑148 (MQ=255)
         ctgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCa                                                                                                                                            <  2:1491964/148‑1 (MQ=255)
          tgGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCAc                                                                                                                                           >  2:2453459/1‑148 (MQ=255)
           ggCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCAGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATACGTGAGCACCTGAAAGCGCAcg                                                                                                                                          >  1:631526/1‑148 (MQ=255)
            gCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCAcgg                                                                                                                                         >  1:1583986/1‑148 (MQ=255)
               tttCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCAcggcgg                                                                                                                                      <  2:1320023/148‑1 (MQ=255)
                    aTTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGaaa                                                                                                                                 >  2:3237305/1‑148 (MQ=255)
                    aTTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGaaa                                                                                                                                 >  1:2461310/1‑148 (MQ=255)
                     ttGCTCTTCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGAGTTCTGACGTTATCGAAATCCGTGAGCAACTGAAAGCGCACGGCGGCGaaaa                                                                                                                                >  2:292165/1‑148 (MQ=255)
                     ttGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCGGAAAGCGCACGGCGGCGaaaa                                                                                                                                >  1:38783/1‑148 (MQ=255)
                      tGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAAc                                                                                                                               <  1:2769872/148‑1 (MQ=255)
                          ctGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATcc                                                                                                                           <  1:722437/148‑1 (MQ=255)
                          ctGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATcc                                                                                                                           <  1:2504653/148‑1 (MQ=255)
                                 cACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCAcatcat                                                                                                                    >  2:2683978/1‑148 (MQ=255)
                                 cACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCAcatcat                                                                                                                    >  2:3229317/1‑148 (MQ=255)
                                 cACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCAcatcat                                                                                                                    >  2:1175948/1‑148 (MQ=255)
                                      gCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCa                                                                                                               <  1:109798/148‑1 (MQ=255)
                                      gCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCa                                                                                                               >  1:592045/1‑148 (MQ=255)
                                         gAAAAAGACAAACGGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCaaaa                                                                                                            <  2:1807527/148‑1 (MQ=255)
                                         gAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCaaaa                                                                                                            >  1:909097/1‑148 (MQ=255)
                                           aaaaGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATc                                                                                                          >  2:3366094/1‑148 (MQ=255)
                                           aaaaGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATc                                                                                                          >  1:1415285/1‑148 (MQ=255)
                                            aaaGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCg                                                                                                         >  1:14558/1‑148 (MQ=255)
                                             aaGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGa                                                                                                        <  1:3090765/148‑1 (MQ=255)
                                             aaGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGa                                                                                                        <  1:1316988/148‑1 (MQ=255)
                                              aGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTTCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGaa                                                                                                       <  1:2712398/148‑1 (MQ=255)
                                              aGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGaa                                                                                                       <  1:2296232/148‑1 (MQ=255)
                                              aGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGaa                                                                                                       >  2:2416290/1‑148 (MQ=255)
                                              aGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGaa                                                                                                       <  1:415567/148‑1 (MQ=255)
                                               gACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCTCATCATCTCCAAAATCGaaa                                                                                                      <  2:390359/148‑1 (MQ=255)
                                               gACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCGGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGaaa                                                                                                      <  1:659042/148‑1 (MQ=255)
                                                aCAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGaaaa                                                                                                     >  2:1837669/1‑148 (MQ=255)
                                                  aaaCAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAAcc                                                                                                   >  1:1088146/1‑148 (MQ=255)
                                                     cAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAgg                                                                                                <  2:379513/148‑1 (MQ=255)
                                                     cAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAgg                                                                                                <  2:2170791/148‑1 (MQ=255)
                                                     cAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAgg                                                                                                <  1:976875/148‑1 (MQ=255)
                                                       ggACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGaa                                                                                              <  1:231791/148‑1 (MQ=255)
                                                       ggACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGaa                                                                                              >  1:2244725/1‑148 (MQ=255)
                                                         aCCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAgg                                                                                            <  2:982996/148‑1 (MQ=255)
                                                         aCCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAgg                                                                                            >  2:431631/1‑148 (MQ=255)
                                                          ccTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGc                                                                                           <  1:933774/148‑1 (MQ=255)
                                                            tGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCt                                                                                         <  1:550493/148‑1 (MQ=255)
                                                             gATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTc                                                                                        >  1:524073/1‑148 (MQ=255)
                                                                cTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTcaac                                                                                     >  1:1723531/1‑148 (MQ=255)
                                                                  ttGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTcaacaa                                                                                   >  2:1458763/1‑148 (MQ=255)
                                                                   tGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCaacaac                                                                                  >  2:493742/1‑148 (MQ=255)
                                                                    ggTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACt                                                                                 <  1:2562923/148‑1 (MQ=255)
                                                                       tGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTcg                                                                              <  2:2512526/148‑1 (MQ=255)
                                                                       tGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTcg                                                                              <  2:2434597/148‑1 (MQ=255)
                                                                       tGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTcg                                                                              <  1:994466/148‑1 (MQ=255)
                                                                           aaCAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTcgacga                                                                          <  1:2366861/148‑1 (MQ=255)
                                                                             cAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGaaa                                                                        >  2:2347281/1‑148 (MQ=255)
                                                                             cAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGaaa                                                                        <  2:675461/148‑1 (MQ=255)
                                                                              aaGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAAt                                                                       <  1:1928479/148‑1 (MQ=255)
                                                                                ggCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATcc                                                                     >  2:1611886/1‑148 (MQ=255)
                                                                                 gCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCGGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCt                                                                    >  2:2312721/1‑148 (MQ=255)
                                                                                   gTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCg                                                                  >  2:2995545/1‑148 (MQ=255)
                                                                                    tAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGa                                                                 <  2:1105783/148‑1 (MQ=255)
                                                                                     aGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGaa                                                                <  1:1669081/148‑1 (MQ=255)
                                                                                      gACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAg                                                               >  2:614878/1‑148 (MQ=255)
                                                                                         tttGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCct                                                            <  1:1269646/148‑1 (MQ=255)
                                                                                          ttgttgCTGCTTCCTTTTTTCTTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACTTCATCTCCAAAATCGAAAACCAGTAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCctc                                                           <  1:611925/148‑1 (MQ=255)
                                                                                          ttgttgCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCctc                                                           <  1:2433236/148‑1 (MQ=255)
                                                                                                ctgctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACgg                                                     <  1:2713223/148‑1 (MQ=255)
                                                                                                  gctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGca                                                   >  2:101961/1‑148 (MQ=255)
                                                                                                    ttCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcatc                                                 <  2:697653/148‑1 (MQ=255)
                                                                                                       cTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACCTCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATg                                              <  1:3445053/148‑1 (MQ=255)
                                                                                                       cTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATg                                              >  2:709395/1‑148 (MQ=255)
                                                                                                        tttATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATgg                                             >  1:789107/1‑148 (MQ=255)
                                                                                                          tATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCGGAAAGCGCACTGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGtt                                           >  1:482334/1‑148 (MQ=255)
                                                                                                             tCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTgcg                                        <  2:1144049/148‑1 (MQ=255)
                                                                                                               gTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGcgcg                                      <  1:3237305/148‑1 (MQ=255)
                                                                                                                 aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTg                                    >  2:2614778/1‑148 (MQ=255)
                                                                                                                      ttCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCCTCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGAc                               <  2:400285/148‑1 (MQ=255)
                                                                                                                      ttCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGAc                               <  2:1286636/148‑1 (MQ=255)
                                                                                                                      ttCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGAc                               <  1:1765423/148‑1 (MQ=255)
                                                                                                                        cTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCt                             <  1:866251/148‑1 (MQ=255)
                                                                                                                          gACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTgg                           >  1:1971262/1‑148 (MQ=255)
                                                                                                                            cGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgt                         <  2:1455880/148‑1 (MQ=255)
                                                                                                                             gTTATCGAAATCCGTGAGCACCGGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgtg                        <  2:1097347/148‑1 (MQ=255)
                                                                                                                             gTTATCGAAATCCGTGAGCACCGGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgtg                        <  2:2871815/148‑1 (MQ=255)
                                                                                                                              ttATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgtgt                       >  2:3112979/1‑148 (MQ=255)
                                                                                                                                    aaaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAAt                 <  2:1723531/148‑1 (MQ=255)
                                                                                                                                       tCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATccc              >  2:2830980/1‑148 (MQ=255)
                                                                                                                                        ccGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCg             >  2:1229954/1‑148 (MQ=255)
                                                                                                                                         cGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCgg            <  1:3513295/148‑1 (MQ=255)
                                                                                                                                             aGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTaga        <  2:3196705/148‑1 (MQ=255)
                                                                                                                                             aGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTaga        >  2:862037/1‑148 (MQ=255)
                                                                                                                                              gCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAACTCGAAAACCAGTAAGGCCTAAACAACTTCGACGAAATACTCGAAGCCTCTGACGGCATCATGGTTGCGAGTGGCGACCTGGGTGTAGAAATCCCGGTagaa       <  1:288792/148‑1 (MQ=255)
                                                                                                                                              gCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTagaa       <  1:3199996/148‑1 (MQ=255)
                                                                                                                                               cACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGTCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTagaag      <  1:1611886/148‑1 (MQ=255)
                                                                                                                                                 ccTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAgaagaa    <  1:1141075/148‑1 (MQ=255)
                                                                                                                                                 ccTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAgaagaa    <  2:2708429/148‑1 (MQ=255)
                                                                                                                                                  cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGaagaag   <  1:138399/148‑1 (MQ=255)
                                                                                                                                                   tGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGaagaagc  <  2:2461310/148‑2 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGAACCTGCCTGGCGTTTCCATTGCTCTGCCAGCACTGGCTGAAAAAGACAAACAGGACCTGATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGT  >  CP009273/1750427‑1750721

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: ‑