Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 951,030 G→A 5.5% D17D (GAC→GAT pflA ← 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_000913951,0300GA 5.5% 200.4 / 2.5 73D17D (GAC→GATpflApyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base A (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

GGTCACCACTTCCTTCATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGAT  >  NC_000913/950895‑951178
                                                                                                                                       |                                                                                                                                                    
ggTCACCACTTCCTTCATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCcac                                                                                                                                         <  1:848341/149‑1 (MQ=255)
 ggcaccacTTCCTTAATCAAATCTTCAACGGTAACATCTTAACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACACGCAGCGCATCAGGCAGCCCTTTAAAAAGGTGATAACGCGATTACCTGGGCCGTCTACGGTTCcaca                                                                                                                                        <  2:817613/147‑1 (MQ=255)
 gTCACCACTTCCTTCATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTTTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCcaca                                                                                                                                        <  2:416857/149‑1 (MQ=255)
 gTCACCACTTCCTTCATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGTTGTCGCGGTTATGACAATACAGGCAGCTCATCATGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCcaca                                                                                                                                        <  2:264551/149‑1 (MQ=255)
  tCACCACTTCCTTCATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAATCGAATACCTGGGCAGTCTACGGTTCCACAg                                                                                                                                       >  1:1001503/1‑149 (MQ=255)
       acTTCCTTCATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTc                                                                                                                                  <  1:177063/149‑1 (MQ=255)
           ccttcATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAg                                                                                                                              <  2:838127/149‑1 (MQ=255)
           ccttcATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAg                                                                                                                              <  1:871615/149‑1 (MQ=255)
           acttcATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGGGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAg                                                                                                                              <  1:83372/148‑1 (MQ=255)
              tcatcaAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAg                                                                                                                           <  1:635977/149‑1 (MQ=255)
                 tcaAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGCCCGTCTACGGTTCCACAGGATTCAAAGGAGTGa                                                                                                                        >  1:833138/1‑149 (MQ=255)
                   aaaTCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAAt                                                                                                                      >  1:1020800/1‑149 (MQ=255)
                   aaaTCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAAt                                                                                                                      >  2:262766/1‑149 (MQ=255)
                    aaTCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATAGAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATg                                                                                                                     <  2:220863/149‑1 (MQ=255)
                        atccagggTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGGGGTTAGGTCATTACAGGCAGCGCAGCAGCCAGCCCTGGAAAACGTTGATAACGCGAATACCTGGGCCGTCAACGGTACCACAGGATTCAAAGGAGTGAATGCGAc                                                                                                                 <  1:121279/143‑1 (MQ=255)
                         tCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTCGGCCGTCTACGGTTCCACAGGATTCAAAGGACTGAATGCGAcc                                                                                                                >  2:292319/1‑149 (MQ=255)
                           aaCGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCCGCCCGGGAAAAAGGTGATAAAGCGACTACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCaa                                                                                                              >  2:935525/1‑149 (MQ=255)
                             cGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATa                                                                                                            >  2:489559/1‑149 (MQ=255)
                               gTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAAc                                                                                                          <  2:493342/149‑1 (MQ=255)
                                ttaCTTCTTTAGCGCATTTCGTGTCGCTGGTGTAGCGGTTATGACAATACAGGCAGCTCATCTGGCAGCCCGGGAAAAAGGTGATAATGAGATTACCTGGGCCGTCTACGGTTCCTCAGGATTCAAAGGAGTGAATGCGACAAATAACt                                                                                                         <  2:605141/147‑1 (MQ=255)
                                aaaCTTCTATACCGCCATGCGTGTCCCCGGTGTCGCGGGTATGACAATACAGGCAGCGCAGCAGGCAGCCCTGGAAACAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCCCAGGATTCAAAGGAGTGAACGCGACCAATAACt                                                                                                         <  2:414254/148‑1 (MQ=255)
                                  aCTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGTTTATGACAATTCAGGCTGCGCATCAGGCTGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGAGCAATAACTGa                                                                                                       <  2:300864/149‑1 (MQ=255)
                                  aCTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGa                                                                                                       >  1:57200/1‑149 (MQ=255)
                                   cttcttTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCAGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATATCTgcc                                                                                                      >  2:154821/1‑147 (MQ=255)
                                    ttcttTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCATCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACa                                                                                                     >  1:99203/1‑149 (MQ=255)
                                    ttcttTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACa                                                                                                     >  2:527171/1‑149 (MQ=255)
                                    ttcttTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCAACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACa                                                                                                     >  2:995821/1‑149 (MQ=255)
                                    ttcttTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCATCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGAc                                                                                                                 <  1:391563/137‑1 (MQ=255)
                                    ttcttTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCATCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGAc                                                                                                                 >  2:391563/1‑137 (MQ=255)
                                    ttcttTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCATCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGAc                                                                                                                 >  2:770303/1‑137 (MQ=255)
                                    ttcttTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCATCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGAc                                                                                                                 <  1:770303/137‑1 (MQ=255)
                                    ttctaTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTTGAAAAAGGTGATAAATCGAATACCTTGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACa                                                                                                     <  1:910676/149‑1 (MQ=255)
                                       tttAACGCCATGCGTGACCCAGGTGTCGCGGTTATGACAAAACAGGCAGCGCATCAGGCAGCCCTGGAATAAGCTGATAAAGCGAATACCTGGTCCGTCTACGGTTCCAAAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTg                                                                                                  <  1:788342/149‑1 (MQ=255)
                                           ccGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGgt                                                                                              >  2:169723/1‑149 (MQ=255)
                                            cGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGgtg                                                                                             >  1:749500/1‑149 (MQ=255)
                                                aTGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGCATGCGACCAATAACTGACATTGCGGTGTTtc                                                                                         >  1:522907/1‑149 (MQ=255)
                                                aTGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTtc                                                                                         >  2:336084/1‑149 (MQ=255)
                                                  gCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTtctc                                                                                       >  1:342312/1‑149 (MQ=255)
                                                   cGTGTCCAAGGGGTCGCGGTTATGCCAATACTGGCAGCGCTTCAGGCAGCCATGGAAACAGGTGTTAAAGGGAATACATGGGCAGTCTACGGTTCCACAGGATTTTAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTcc                                                                                      <  2:665817/149‑1 (MQ=255)
                                                      ggcctaGGTGTCTCGGTTATGACTATACAGGCAGCGAATCAGGCTTCCCTGGAAAAAGGTCAGAAAGCGAATACCTGTGCCGTCTACGGTTCCAGAGGATTCAAAGGAGTGAATGCGCCCAATAACTGACATTGCGGTGTTTCTCCAGa                                                                                   <  2:71972/144‑1 (MQ=255)
                                                           atgtgtCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGGGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTgg                                                                              <  2:528516/147‑1 (MQ=255)
                                                           aGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGTTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGATTCTCAAGATGTgg                                                                              >  2:224399/1‑149 (MQ=255)
                                                                tCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCAt                                                                         >  1:952226/1‑149 (MQ=255)
                                                                  gcgTTTATGACTCTACAGGCAGCGCATCAGGCAGCACTGGAAAATGGTGAAAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTGAAAGGAGTGAATGCGACCAATAACTGCCATTGCGGTGTTTCTCCAGATGTGGCCCATCt                                                                       <  1:270659/149‑1 (MQ=255)
                                                                          gACAATACAGGCAGCGCATCAGGCAGCCCTGGATAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCgg                                                               >  1:210097/1‑148 (MQ=255)
                                                                          gACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCCCAGGATTCAAAGGAGTGAATGCGACCAATAACTGCCATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCgt                                                               >  1:604395/1‑149 (MQ=255)
                                                                          gACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCGACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCACCAGATGTGCACCATCTGCGGCCgt                                                               >  2:243801/1‑149 (MQ=255)
                                                                           aCAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGCTGTGGCCCATCTGAGGCCgtg                                                              >  2:595380/1‑149 (MQ=255)
                                                                            cAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCgtgt                                                             <  2:11602/149‑1 (MQ=255)
                                                                            cAATACAGGCAGCGCATCAGGCAGACCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCACTAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCgtgt                                                             <  1:224399/149‑1 (MQ=255)
                                                                                 cAGGCAGCGCATCAGGCAGCCCTGGAAAAACTGGATAAAGCGAATCCCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTg                                                        <  1:836158/149‑1 (MQ=255)
                                                                                       gcgcATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGc                                                  >  1:53612/1‑149 (MQ=255)
                                                                                           aTCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGCCATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGa                                              >  2:728552/1‑149 (MQ=255)
                                                                                           aTCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGa                                              <  2:730470/149‑1 (MQ=255)
                                                                                           aTCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGa                                              >  1:12554/1‑149 (MQ=255)
                                                                                               ggCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCAACGGTTCCACAGGAATCAAAGGAGTGAAAGCGACAAATAACAGACAATACGGAATTTCTCCAGAAGTGGCCCATCAGCAGCCGTCTTGGTGCGCAGCTCGAAgcc                                          >  1:443965/1‑147 (MQ=255)
                                                                                                gCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCCCAGGATTCAAAGGAGTGAATGCGACCCATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATATGAGGCCGGGTAGGTGCGCAGCTCGCAGGCt                                         >  1:834741/1‑149 (MQ=255)
                                                                                                     ccTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTc                                    <  1:352099/149‑1 (MQ=255)
                                                                                                     ccTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGGTGGTGCGCAGCTCGAAGGCTACGTc                                    >  2:269541/1‑149 (MQ=255)
                                                                                                       tGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGa                                  <  1:527171/149‑1 (MQ=255)
                                                                                                       tGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGa                                  <  2:69903/149‑1 (MQ=255)
                                                                                                       tGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTATCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCAACGTCGa                                  >  2:153513/1‑149 (MQ=255)
                                                                                                         gAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGt                                >  2:265656/1‑149 (MQ=255)
                                                                                                                 tGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAAAAACATACATTGCGGTGTTTCTCCAGATGTGGACCATCTGAGGCCGTGTTGGTTCGCATCTCGAAGGCAACGTCGAGACTGTTTTg                        >  2:104527/1‑149 (MQ=255)
                                                                                                                  gATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTgg                       <  1:935525/149‑1 (MQ=255)
                                                                                                                   aaaacgCGAATACCAGGGCCGTCTACGGTGCCACAGGATTCAAATGAGTGAATGCGACCAATACCTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGc                      <  2:662713/144‑1 (MQ=255)
                                                                                                                   aTAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGCCCAATAACTGACATTGCGGGGTTTCTCCAGATGTGGCCCATCTGAGGACGTGTTGGTGCGCAGCTCGATGGCTACGTCGACTCTGTTTTGGc                      >  1:786506/1‑149 (MQ=255)
                                                                                                                    tAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCa                     <  2:723831/149‑1 (MQ=255)
                                                                                                                      aaGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGt                   >  1:655424/1‑149 (MQ=255)
                                                                                                                         cGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCAc                <  2:1001503/149‑1 (MQ=255)
                                                                                                                            aTACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGAGGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGATGGCTACGTCGAGTCTGTTTTGGCAGTCACCtt             >  2:488607/1‑149 (MQ=255)
                                                                                                                                      cGTCTACGGTTCCACAGGATTCATAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTataga   <  1:294450/149‑1 (MQ=255)
                                                                                                                                      cGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGtatata   >  1:696136/1‑147 (MQ=255)
                                                                                                                                       gTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCGTAAAGTatagat  >  1:625499/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                    
GGTCACCACTTCCTTCATCAAATCTTCAACGGTAACTTCTTTACCGCCATGCGTGTCCCAGGTGTCGCGGTTATGACAATACAGGCAGCGCATCAGGCAGCCCTGGAAAAAGGTGATAAAGCGAATACCTGGGCCGTCTACGGTTCCACAGGATTCAAAGGAGTGAATGCGACCAATAACTGACATTGCGGTGTTTCTCCAGATGTGGCCCATCTGAGGCCGTGTTGGTGCGCAGCTCGAAGGCTACGTCGAGTCTGTTTTGGCAGTCACCTTAAAGTATAGAT  >  NC_000913/950895‑951178

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