Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 950,317 A→T 9.4% intergenic (+59/+23) ycaK → / ← pflA putative NAD(P)H‑dependent oxidoreductase/pyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913950,3170AT 9.4% 326.8 / 16.8 123intergenic (+59/+23)ycaK/pflAputative NAD(P)H‑dependent oxidoreductase/pyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase
Reads supporting (aligned to +/- strand):  ref base A (53/58);  new base T (7/5);  total (60/63)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

TAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGCCC  >  NC_000913/950168‑950466
                                                                                                                                                     |                                                                                                                                                     
tAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGAAGAATGAAGCGCGGAATAAAAAAGCTGAAACTAAATa                                                                                                                                                       >  1:192481/1‑150 (MQ=255)
  aTACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATaaa                                                                                                                                                     >  2:207255/1‑150 (MQ=255)
  aTACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATaaa                                                                                                                                                     >  2:1719742/1‑150 (MQ=255)
    aCAGTGGTATTTGATGTTGAATAACTTCACGCTTGCTATTATCAGTCTTTATTATCTCAGGTACGGGAAATGGTAGATGCACTATAAGATTTGTTAATATCTCTGTAGCAGAATGAAGCGCGGATTAAAAAAGCGGCAACTCAATTAAGt                                                                                                                                                   <  1:125034/150‑1 (MQ=255)
    aCAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGt                                                                                                                                                   <  2:1497310/150‑1 (MQ=255)
      aGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTg                                                                                                                                                 <  2:1405371/150‑1 (MQ=255)
      aGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTg                                                                                                                                                 <  1:1036871/150‑1 (MQ=255)
       gTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGc                                                                                                                                                >  1:1515106/1‑150 (MQ=255)
       gTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGc                                                                                                                                                <  1:146073/150‑1 (MQ=255)
        tGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGcc                                                                                                                                               <  2:655610/150‑1 (MQ=255)
          ttATTCGATGGGGAAGAACTTCACGCGAGCTTTTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGc                                                                                                                                             <  2:1465584/149‑1 (MQ=255)
           tATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAATGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCt                                                                                                                                            <  1:187138/150‑1 (MQ=255)
           tATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCGCt                                                                                                                                            <  2:1515106/150‑1 (MQ=255)
           tATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCGCt                                                                                                                                            <  2:794538/150‑1 (MQ=255)
              tCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCGCtttt                                                                                                                                         <  2:103932/150‑2 (MQ=255)
               cGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTAc                                                                                                                                        >  2:1647760/1‑150 (MQ=255)
                 attgggAAGAACTTCACGCGAGCTATACTCAGTCGTTATTATCTCAGGAACGTGATATGGTAGATGCACTAAAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACgg                                                                                                                                      <  2:450428/147‑1 (MQ=255)
                 aTGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACgg                                                                                                                                      >  1:1036566/1‑150 (MQ=255)
                 aTGCGGAAGAACTGCACGCTAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGAAGCACTATAAGCTGTGTTAAAACCGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACgg                                                                                                                                      <  2:252991/150‑1 (MQ=255)
                      gaagaaCTTCACGCGAGCTATTCTCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGaaa                                                                                                                                 <  2:1569070/150‑1 (MQ=255)
                        agaaCTTCTCGCGAGCTATTATCAGTCGTTATAATATCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAAtt                                                                                                                               <  2:192481/150‑1 (MQ=255)
                        agaaCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAAtt                                                                                                                               <  2:808771/150‑1 (MQ=255)
                         gaaCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGcc                                                                                                                                               <  1:935871/133‑1 (MQ=255)
                         gaaCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTa                                                                                                                              >  1:398697/1‑150 (MQ=255)
                              tCACGCTAGCTATTATCAGTCGTTATTATCTCATGTACGGGATATGGTAGATGCACTATAAGATGTTTTAAAAACGCTGTAGCAGTATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACa                                                                                                                         <  1:1209248/150‑1 (MQ=255)
                                 cgcgAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAAT‑‑‑AAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGTGAAATTAGAACATTACCt                                                                                                                   >  2:367470/1‑150 (MQ=255)
                                      gCTATTATCAGTCGTTATTATATCAGGTACGGGATATTGTATATGCACTATAAGATGTGTTAAAAACTCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCTGCAACTAAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTa                                                                                                                 <  2:620832/150‑1 (MQ=255)
                                               aGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTAc                                                                                                        >  1:511837/1‑150 (MQ=255)
                                                   gTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCt                                                                                                    <  2:415109/150‑1 (MQ=255)
                                                     tattatCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCa                                                                                                  >  2:1494820/1‑150 (MQ=255)
                                                     tattatCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCa                                                                                                  >  2:469769/1‑150 (MQ=255)
                                                         atCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGaa                                                                                              <  1:362913/150‑1 (MQ=255)
                                                               ggTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCtt                                                                                        <  1:1087326/150‑1 (MQ=255)
                                                                   cGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCAGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGATCAAGAATGCCTTTCAc                                                                                    >  1:847884/1‑150 (MQ=255)
                                                                    gggATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCAcg                                                                                   <  1:1307553/150‑1 (MQ=255)
                                                                     ggATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCTGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCActc                                                                                  >  1:110236/1‑148 (MQ=255)
                                                                          tGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTcc                                                                             >  2:134648/1‑150 (MQ=255)
                                                                          tGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTcc                                                                             >  2:264537/1‑150 (MQ=255)
                                                                           ggTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCa                                                                            >  1:1410902/1‑150 (MQ=255)
                                                                            gTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGTGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCAt                                                                           <  1:272110/150‑1 (MQ=255)
                                                                             tAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATg                                                                          <  1:530310/150‑1 (MQ=255)
                                                                              aGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATgg                                                                         >  2:1462557/1‑150 (MQ=255)
                                                                              aGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATgg                                                                         >  1:732524/1‑150 (MQ=255)
                                                                                 tGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGtct                                                                      <  1:307030/150‑1 (MQ=255)
                                                                                  gCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGtctc                                                                     <  1:270993/150‑1 (MQ=255)
                                                                                    aCTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCtctt                                                                   <  2:732524/150‑1 (MQ=255)
                                                                                    aCTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCtctt                                                                   <  2:618995/150‑1 (MQ=255)
                                                                                     cTATAAGATGTGTTAAAAACGCTTTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCtcttt                                                                  <  1:566857/150‑1 (MQ=255)
                                                                                     cTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACTCGTTCCATGGTCtcttt                                                                  >  2:116567/1‑150 (MQ=255)
                                                                                        taAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCtctttctt                                                               >  1:582267/1‑150 (MQ=255)
                                                                                        taAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCtctttctt                                                               >  2:802027/1‑150 (MQ=255)
                                                                                         aaGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTttcttc                                                              <  2:796721/150‑1 (MQ=255)
                                                                                          aGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCg                                                             >  2:1425021/1‑150 (MQ=255)
                                                                                          aGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCg                                                             >  2:456884/1‑150 (MQ=255)
                                                                                           gATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCATGAATGACTTTCACGCGTTCCATGGTCTCTTTCTTCgg                                                            >  2:900879/1‑150 (MQ=255)
                                                                                            aTGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggt                                                           <  1:1462557/150‑1 (MQ=255)
                                                                                              gtgtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                         <  2:1372751/150‑1 (MQ=255)
                                                                                              gtgtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                         >  2:1592830/1‑150 (MQ=255)
                                                                                              gtgtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                         <  2:1221478/150‑1 (MQ=255)
                                                                                                gtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGtt                                                       >  1:967336/1‑150 (MQ=255)
                                                                                                gtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGtt                                                       >  2:864795/1‑150 (MQ=255)
                                                                                                gtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGtt                                                       <  2:61137/150‑1 (MQ=255)
                                                                                                 ttAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGttt                                                      <  1:1274238/150‑1 (MQ=255)
                                                                                                   aaaaaCGCTTTAGCACAATGAAGCGCGGAATAAAAAAGCGACAACTCAATAAAGATGCAGCTTTAAGGTGTAATTAGAACATTACCTTATGACCGAACTGCTCAATAATGCATTTCACGCGTTCGATGGTCTCTTTCTTCGGAGGTTTaa                                                    >  2:371518/1‑150 (MQ=255)
                                                                                                    aaaaCGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAac                                                   >  1:40631/1‑150 (MQ=255)
                                                                                                    aaaaCGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCATTCCGGGGAAATTAGCACATTACCTGATGACCGTACTGCTCACGAATCCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTatc                                                   >  1:253953/1‑148 (MQ=255)
                                                                                                     aaaCGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAaca                                                  >  1:1503340/1‑150 (MQ=255)
                                                                                                       aCGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACAcc                                                >  2:1323369/1‑150 (MQ=255)
                                                                                                        cGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCg                                               <  1:207255/150‑1 (MQ=255)
                                                                                                        cGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCg                                               >  2:1240571/1‑150 (MQ=255)
                                                                                                        cGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCg                                               <  1:1719742/150‑1 (MQ=255)
                                                                                                           tGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCg                                            <  1:1647760/150‑1 (MQ=255)
                                                                                                            gTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATTAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTTATGCTCAAGAATGCGTTTCACGCGGTCCATGGTCTCTTTCTTCGGTGGTTTTACACAGTCGa                                           >  2:762220/1‑150 (MQ=255)
                                                                                                              aGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGt                                         >  2:1302467/1‑150 (MQ=255)
                                                                                                               gCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGtt                                        <  2:1466418/150‑1 (MQ=255)
                                                                                                               gCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGtt                                        <  1:861979/150‑1 (MQ=255)
                                                                                                                cAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGttt                                       <  2:511837/150‑1 (MQ=255)
                                                                                                                     tGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTAct                                  <  2:1413155/150‑1 (MQ=255)
                                                                                                                       aaGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTActct                                <  2:1632970/150‑1 (MQ=255)
                                                                                                                       aaGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTActct                                <  1:1494820/150‑1 (MQ=255)
                                                                                                                          cgcgGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCa                             >  1:165868/1‑150 (MQ=255)
                                                                                                                              gAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCa                         >  2:1512658/1‑150 (MQ=255)
                                                                                                                                aTAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCAtt                       >  1:149821/1‑150 (MQ=255)
                                                                                                                                   aaaaaGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGcc                    <  1:1112177/150‑1 (MQ=255)
                                                                                                                                     aaaGCGGCAACTTTATTGAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCAc                  >  1:709603/1‑150 (MQ=255)
                                                                                                                                       aGCGGCAACTTTATTGAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTACCACCGTCGAGTTTGTACTCTTCACCCATTGCCAccc                >  2:268317/1‑150 (MQ=255)
                                                                                                                                       aGCGGCAACTTTATTGAGTTGCCGCTTTACGGGAAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTTTACTCTTCACCCATTGCCAccc                >  2:1577821/1‑150 (MQ=255)
                                                                                                                                       aGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCAccc                >  1:461496/1‑150 (MQ=255)
                                                                                                                                        gCGGCAACTTTATTGAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTAAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCa               >  1:1011479/1‑150 (MQ=255)
                                                                                                                                        gCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATTCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCa               >  2:686454/1‑150 (MQ=255)
                                                                                                                                         cGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCAt              >  2:914637/1‑150 (MQ=255)
                                                                                                                                         cGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCAt              >  2:1083659/1‑150 (MQ=255)
                                                                                                                                          ggCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAAAATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACAGTCGAGTATGTACTCTTCACACATTGCAACCCAtt             >  1:1025714/1‑150 (MQ=255)
                                                                                                                                             aaCTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtgt          <  1:609379/150‑1 (MQ=255)
                                                                                                                                              aCTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtgtg         <  2:398697/150‑1 (MQ=255)
                                                                                                                                              aCTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCGTCACCCATTGCCACCCATTtgtg         >  2:1390512/1‑150 (MQ=255)
                                                                                                                                               cTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCAATTTgtgt        >  1:487657/1‑150 (MQ=255)
                                                                                                                                                   aTAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGc    <  1:870533/150‑1 (MQ=255)
                                                                                                                                                   aTAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGc    <  1:1083659/150‑1 (MQ=255)
                                                                                                                                                     aaaGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGccc  <  1:864795/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGCCC  >  NC_000913/950168‑950466

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: