Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 950,317 A→T 5.8% 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 5.8% 819.7 / 23.6 229intergenic (+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 (122/93);  new base T (4/9);  total (127/102)
Fisher's exact test for biased strand distribution p-value = 8.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 234 total aligned reads displayed.

TAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGCC  >  NC_000913/950168‑950465
                                                                                                                                                     |                                                                                                                                                    
tAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATa                                                                                                                                                      >  1:102654/1‑150 (MQ=255)
  aTACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATaaa                                                                                                                                                    <  2:1843856/150‑1 (MQ=255)
     ctgTGGTATTCGTGGGGTAATAACTTCACGCGAGTTATTATCAGTCTTTTTTATCTCTGGTACGGTATATGTTTGATGAACTTTAATATGTTTTAAAAACGCTGTAGATTAATGAAGCGCGTAATAAAAAAGCGGCAACTCAATAAAGtt                                                                                                                                                 <  1:4004727/148‑1 (MQ=255)
      aGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTg                                                                                                                                                >  1:3790885/1‑150 (MQ=255)
        tGGTATTCGATGGGGAAGAACTTCACGCGAGATATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGcc                                                                                                                                              <  1:2382348/150‑1 (MQ=255)
         ggTATTGGTGGGGGAAGAACTTCACGCGCGCTATTATCAGTCGTTATTATCTCCGGTACGGGATATGGTAGATGCACTATAAGATTTGTTAAAAACGCTTTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCg                                                                                                                                             <  1:3687593/150‑1 (MQ=255)
          gTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAACCGCTGTAGCAGAATGACGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCGc                                                                                                                                            <  1:231803/150‑1 (MQ=255)
          gTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCGc                                                                                                                                            <  2:627385/150‑1 (MQ=255)
          gTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCGc                                                                                                                                            <  2:3256681/150‑1 (MQ=255)
           tATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATTTGTTAAAAAAGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTTCCGCt                                                                                                                                           >  1:748661/1‑150 (MQ=255)
           tATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCt                                                                                                                                           >  2:2724274/1‑150 (MQ=255)
             ttCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGTTGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAATTGCGGCAACTCAATAAAGTTTCCGCttt                                                                                                                                         >  1:674507/1‑150 (MQ=255)
             ttCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTTTATTGAGTTGCCGCttt                                                                                                                                         <  2:2219988/150‑1 (MQ=255)
               cGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTAc                                                                                                                                       <  2:2144897/150‑1 (MQ=255)
               cGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTAc                                                                                                                                       >  2:888132/1‑150 (MQ=255)
               cGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCGGTAGCAGAATGAAGCGCGGAATAAAAAATCGGCAACTCAATAAAGTTGCCGCTTTAc                                                                                                                                       >  1:1440785/1‑150 (MQ=255)
                gATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACg                                                                                                                                      <  2:1528284/150‑1 (MQ=255)
                 aTGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACgg                                                                                                                                     <  2:1501247/150‑1 (MQ=255)
                      gaagaaCTTCACGCGCGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCCCGATAACATGTGTTAAAAACGCGGTAGCAGAATGAAGCGCGGCATAAAAAAGCGGCAACGCAAAAAAGTTGCCGCGTTACGGGGaaa                                                                                                                                >  1:2630700/1‑150 (MQ=255)
                         gaaCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGCGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCAGTAGCAGAATGAAGAGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTa                                                                                                                             >  2:77065/1‑150 (MQ=255)
                           aCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGa                                                                                                                           >  2:3869448/1‑150 (MQ=255)
                           aCTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGa                                                                                                                           <  2:4158861/150‑1 (MQ=255)
                            cTTCACTCTAGCTATTATCTGTCGTTATTATCTCAGGTACGGGATATGGTAGATGAACTATAAGATGTGTTATAAACGCTGTATCTTATTGAAGCGCGGAATAATAAAGAGGCTACTCAATAAAGTTGCCGCTTTACGGGTAAATTTGaa                                                                                                                          <  2:3598737/150‑1 (MQ=255)
                            cTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGaa                                                                                                                          <  2:2297361/150‑1 (MQ=255)
                             ttCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGACGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAAc                                                                                                                         >  2:3955397/1‑150 (MQ=255)
                             ttCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAAc                                                                                                                         >  2:4230016/1‑150 (MQ=255)
                             ttCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAAc                                                                                                                         >  2:48471/1‑150 (MQ=255)
                             ttCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGAACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAAc                                                                                                                         >  1:3456230/1‑150 (MQ=255)
                             ttCACGCGAGCTATTATCAGTCGTTATTATCTAAGGTACCGGATATGGTAGATGCACGATAAGATTTTTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACACAAAAAAGTTTCCGCATTACGGGGAAAATAGAAc                                                                                                                         >  1:1038944/1‑150 (MQ=255)
                                aCGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGGTATGGTAGGTGCACAATAATATGTGTTAAAAAAGCTGTAGCAGAAAGAAGCGCGGAATAAAAAATCGGCAACACAATAAAGTTGCCGCTTTAAGGGGAAAATATaaaaat                                                                                                                      >  1:841934/1‑146 (MQ=255)
                                aCGCGAGCTATTATCAGTCGTGATTATCTCAGGGACGGGATATGGGAGATGCACAATAAGATTTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCCACACAATAAACTTTCCGCTTTAAGGGGAACTTAGAACaat                                                                                                                      >  1:1253654/1‑148 (MQ=255)
                                 cgcgAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTa                                                                                                                     <  2:2042291/150‑1 (MQ=255)
                                   cgAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTTGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTAcc                                                                                                                   <  1:1429282/150‑1 (MQ=255)
                                     aGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCtt                                                                                                                 >  2:2006189/1‑150 (MQ=255)
                                       cTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTAt                                                                                                               <  1:3834207/150‑1 (MQ=255)
                                        tattatCAGTCGTTATTATCTCAGGTACGGGATATTGTAGATGCACTATAAGAAGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTAACGGGGAAATTAGAACAATACCTTATg                                                                                                              >  2:3250462/1‑150 (MQ=255)
                                             tCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGt                                                                                                         >  2:92639/1‑150 (MQ=255)
                                             tCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGt                                                                                                         <  1:2498145/150‑1 (MQ=255)
                                                 tCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACtt                                                                                                     <  2:3950825/150‑2 (MQ=255)
                                                    ttattaTCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTc                                                                                                  >  1:2700636/1‑150 (MQ=255)
                                                    ttattaTCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTc                                                                                                  >  1:3054610/1‑150 (MQ=255)
                                                    ttattaTCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTc                                                                                                  >  1:552947/1‑150 (MQ=255)
                                                      attatCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCaa                                                                                                >  2:3982869/1‑150 (MQ=255)
                                                        tatCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGa                                                                                              >  1:4011915/1‑150 (MQ=255)
                                                            tcAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGc                                                                                          >  2:3295420/1‑150 (MQ=255)
                                                            tcAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGc                                                                                          <  1:706698/150‑1 (MQ=255)
                                                            tcAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGc                                                                                          >  2:3016387/1‑150 (MQ=255)
                                                            tcAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGc                                                                                          >  1:2242450/1‑150 (MQ=255)
                                                              aGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCt                                                                                        <  2:3541582/150‑1 (MQ=255)
                                                              aGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCt                                                                                        >  2:3418528/1‑150 (MQ=255)
                                                                  aCGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCa                                                                                    >  2:22489/1‑150 (MQ=255)
                                                                  aCGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCa                                                                                    >  2:2182540/1‑150 (MQ=255)
                                                                    gggATATGGTAGATGCACTATAAGATGTGATAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCAcg                                                                                  <  2:3689148/150‑1 (MQ=255)
                                                                        tatGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGtt                                                                              <  1:3955385/150‑1 (MQ=255)
                                                                           ggTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCa                                                                           <  1:2845101/150‑1 (MQ=255)
                                                                            gTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCAt                                                                          >  1:3341128/1‑150 (MQ=255)
                                                                            gTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCAt                                                                          >  1:2227475/1‑150 (MQ=255)
                                                                                  gCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGtctc                                                                    >  1:3590863/1‑150 (MQ=255)
                                                                                   cACTATAAGATGTGTTAAAAACGCTGTAGAAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTctct                                                                   <  2:726883/150‑1 (MQ=255)
                                                                                   cACTATAAGATGTGATAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCATTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTctct                                                                   <  1:130481/150‑1 (MQ=255)
                                                                                     cTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCtcttt                                                                 <  2:4060941/150‑1 (MQ=255)
                                                                                          aGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCg                                                            >  2:532687/1‑150 (MQ=255)
                                                                                           gATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCgg                                                           >  2:298840/1‑150 (MQ=255)
                                                                                            aTGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggt                                                          >  1:975428/1‑150 (MQ=255)
                                                                                              gtgtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATTCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                        <  2:3458336/150‑1 (MQ=255)
                                                                                              gtgtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCggtgg                                                        >  1:2884749/1‑150 (MQ=255)
                                                                                                gtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGtt                                                      >  1:1925680/1‑150 (MQ=255)
                                                                                                gtTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACAGTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGtt                                                      <  1:4167430/150‑1 (MQ=255)
                                                                                                 ttAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGttt                                                     >  2:1181073/1‑150 (MQ=255)
                                                                                                 ttAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGttt                                                     >  1:2006738/1‑150 (MQ=255)
                                                                                                 ttAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGttt                                                     >  2:2344974/1‑150 (MQ=255)
                                                                                                 ttAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGttt                                                     >  1:917384/1‑150 (MQ=255)
                                                                                                  tAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTa                                                    <  2:248834/150‑1 (MQ=255)
                                                                                                   aataaCGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGATCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTaa                                                   <  1:2532585/147‑1 (MQ=255)
                                                                                                          cTGTAGCAGAATGAAGCTCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTc                                            <  1:1674674/150‑1 (MQ=255)
                                                                                                          cTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTc                                            <  2:3831727/150‑1 (MQ=255)
                                                                                                               gCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAACTTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGtt                                       >  2:1331277/1‑150 (MQ=255)
                                                                                                               gCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGtt                                       <  2:2827607/150‑1 (MQ=255)
                                                                                                               gCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGtt                                       >  1:3008504/1‑150 (MQ=255)
                                                                                                                  gAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGt                                    <  1:3213174/150‑1 (MQ=255)
                                                                                                                  gAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAATTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGt                                    <  2:952872/150‑1 (MQ=255)
                                                                                                                   aaTGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTa                                   <  1:1813318/150‑1 (MQ=255)
                                                                                                                        aGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCtt                              <  2:3501172/150‑1 (MQ=255)
                                                                                                                         gcgcGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTc                             <  1:3869448/150‑1 (MQ=255)
                                                                                                                          cgcgGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCa                            <  2:3194310/150‑1 (MQ=255)
                                                                                                                                aTAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCAtt                      <  1:2236902/150‑1 (MQ=255)
                                                                                                                                aTAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCAtt                      <  1:1198158/150‑1 (MQ=255)
                                                                                                                                 ccgtaaaGCGGCAACTTTATTGAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTg                     >  2:12451/5‑150 (MQ=255)
                                                                                                                                  aaaaaaGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACAGTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGc                    <  2:2768644/150‑1 (MQ=255)
                                                                                                                                    aaaaGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACAATACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTGCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGccc                  >  1:1596393/1‑149 (MQ=255)
                                                                                                                                       aGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCAccc               >  2:3202296/1‑150 (MQ=255)
                                                                                                                                        gCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTGTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGACACCCa              >  2:3135267/1‑150 (MQ=255)
                                                                                                                                           gCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCAttt           >  1:117045/1‑150 (MQ=255)
                                                                                                                                             aaCTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtgt         >  1:2771211/1‑150 (MQ=255)
                                                                                                                                             aaCTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtgt         >  2:1653308/1‑150 (MQ=255)
                                                                                                                                              aCTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTtgtg        >  1:2352387/1‑150 (MQ=255)
                                                                                                                                                tCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGtt      >  2:3420747/1‑150 (MQ=255)
                                                                                                                                                tCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGtt      <  1:1217146/150‑1 (MQ=255)
                                                                                                                                                   aTAAAGTTGCCGCATTACGGGGAAATTACAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGc   <  1:1540397/150‑1 (MQ=255)
                                                                                                                                                    tAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGcc  <  1:3599933/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTATTATCTCAGGTACGGGATATGGTAGATGCACTATAAGATGTGTTAAAAACGCTGTAGCAGAATGAAGCGCGGAATAAAAAAGCGGCAACTCAATAAAGTTGCCGCTTTACGGGGAAATTAGAACATTACCTTATGACCGTACTGCTCAAGAATGCCTTTCACGCGTTCCATGGTCTCTTTCTTCGGTGGTTTAACACCGTCGAGTTTGTACTCTTCACCCATTGCCACCCATTTGTGTTTGCC  >  NC_000913/950168‑950465

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: