Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	2,106,741	A→T	S162S (TCT→TCA)	wbbH ←	O‑antigen polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	2,106,741	0	A	T	98.5%	218.2 / ‑3.1	67	coding (486/1167 nt)	wbbH	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (26/40);  minor base . (1/0);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTC  >  NC_000913_3_pae_tpiA/2106599‑2106885                                                                                                                                               |                                                                                                                                                  tATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaa                                                                                                                                             >  1:212312/1‑149 (MQ=255)  aTTAACAACTACGATTAATGCAGCAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAACATGAATAAAATAACATGATGTTAAGTACCTTTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaa                                                                                                                                            <  2:222415/149‑1 (MQ=255)   ttAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCaaat                                                                                                                                           <  1:407783/149‑1 (MQ=255)        acCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCGGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAaataaa‑ta                                                                                                                                     >  1:562379/1‑149 (MQ=255)         ccTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAataaa‑taa                                                                                                                                    <  1:125284/149‑1 (MQ=255)          cTACGATGAATGCAGAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAAtaaa‑taaa                                                                                                                                   <  2:258381/149‑1 (MQ=255)           tACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAg                                                                                                                                  <  1:48027/149‑1 (MQ=255)           tACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAg                                                                                                                                  >  1:757337/1‑149 (MQ=255)            aCGATGAATGCATAAGAGATGATAACCATAAAGACAATCTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TCAAGc                                                                                                                                 >  1:890463/1‑149 (MQ=255)              gATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCaa                                                                                                                               >  2:596940/1‑149 (MQ=255)                tgaatgaaTAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACATCCATTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCaaaa                                                                                                                             <  2:851901/149‑1 (MQ=255)                    tGCATAAGAGATGCTAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAagtag                                                                                                                         <  1:661385/149‑1 (MQ=255)                      cATAAGAGATGTAAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGtt                                                                                                                       <  2:920463/149‑1 (MQ=255)                      cATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGtt                                                                                                                       <  2:338101/149‑1 (MQ=255)                      cATAAGAAATGATAACCATAAAGAAAAGTTGCTTACCAGTATTCAGTATAATTGAAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGtt                                                                                                                       <  2:507438/149‑1 (MQ=255)                         atgaaatgatACCAATAAAAACAATTTCCTTACCTGTATTAAGTATAATTTCAAAGATAAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGa                                                                                                                    <  2:295791/144‑1 (MQ=255)                         aaGAGATGATAACCATAAAGACAATTTTCTTATCCGTATTCAGTATAATGGCAAACATAAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGa                                                                                                                    <  1:416308/149‑1 (MQ=255)                         aaGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCACATAAA‑TAAAGCAAAACTAGCTAGa                                                                                                                    >  1:500844/1‑149 (MQ=255)                              atgatAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAAAAAGTAAAGTAAATGTTTTACTAACCTTTGTATTAGTAAATTTTTTTGTCCAAATAAC‑TAAATCAAAAGTATTTAGaagaat                                                                                                               >  1:301635/1‑145 (MQ=255)                               tgatAACCATAAAGACAAGTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGAtt                                                                                                              <  1:170420/149‑1 (MQ=255)                                 atAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTgg                                                                                                            >  1:848863/1‑149 (MQ=255)                                   aaCCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCt                                                                                                          >  1:299685/1‑149 (MQ=255)                                     ccATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGAAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGc                                                                                                        <  2:596965/149‑1 (MQ=255)                                      cATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCa                                                                                                       <  1:368410/149‑1 (MQ=255)                                       aTAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCAt                                                                                                      >  2:922202/1‑149 (MQ=255)                                       aTAAAGACAATATGCTTACCAGTATTAAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCAt                                                                                                      <  2:228909/149‑1 (MQ=255)                                        tAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATg                                                                                                     >  2:552362/1‑149 (MQ=255)                                          aaGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTCGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTACAATGATTGGCTGCATGTa                                                                                                   >  2:1019913/1‑149 (MQ=255)                                           aGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTAt                                                                                                  <  1:858600/149‑1 (MQ=255)                                                aTTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTga                                                                                             <  1:86171/149‑1 (MQ=255)                                                    gCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAGTAAAAACAAGTCAAGAAAATGTATTACTTACCTTTGTATTAGTAAATCTTTTTGACCAAATAAA‑CAAAGCAAAAGTAGTTAGAATGCTTGGCTGCATGTATGCTGTGaaa                                                                                         >  1:889604/1‑149 (MQ=255)                                                     cTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAt                                                                                        <  2:192506/149‑1 (MQ=255)                                                       tACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAAttt                                                                                      <  1:664230/149‑1 (MQ=255)                                                         ccAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTTTTTTATTTAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTCCATGTATGCTGTAAAATTTCt                                                                                    >  2:329541/1‑149 (MQ=255)                                                         ccAGTATTCAGAATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATAATGTTCTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTGTGCTGAGAAATTTCt                                                                                    <  1:770874/149‑1 (MQ=255)                                                          cAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCtt                                                                                   <  2:896802/149‑1 (MQ=255)                                                             tATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGAt                                                                                <  1:703677/149‑1 (MQ=255)                                                                  aGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTc                                                                           >  1:968202/1‑149 (MQ=255)                                                                    tataATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCtt                                                                         >  2:329164/1‑149 (MQ=255)                                                                        aTTGCAAAGATGAAGACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAAc                                                                     <  2:871368/149‑1 (MQ=255)                                                                             aaaGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTT‑GACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGc                                                               >  2:781865/1‑149 (MQ=255)                                                                                     aTACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCtt                                                        >  1:163059/1‑149 (MQ=255)                                                                                      tACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTa                                                       >  1:547713/1‑149 (MQ=255)                                                                                         aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCa                                                    >  1:359339/1‑149 (MQ=255)                                                                                         aaTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTATTATCa                                                    >  2:392619/1‑149 (MQ=255)                                                                                          aTAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAACTTACTTAATGTGTCTTCCACATCACCATATGTTATCaa                                                   >  2:841166/1‑149 (MQ=255)                                                                                           tAAAAGGAAGTAAAGAAAATGTGGTACTTACCTTTATATTAGTAAATTTTTTTGACCATACAAA‑TAAAGCAAAAGTAGTTATAATGATTCGCTGCATGTATTCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATGaaa                                                  <  2:417651/149‑1 (MQ=255)                                                                                           tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                                  <  2:542145/149‑1 (MQ=255)                                                                                           tAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCaaa                                                  <  1:792417/149‑1 (MQ=255)                                                                                                 cAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCata                                            <  1:108968/149‑1 (MQ=255)                                                                                                   agtaaagtaaATGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCAtata                                          <  2:156495/149‑1 (MQ=255)                                                                                                      aaagtacaTGTTTTAATTACCTTTGTATTAGTAAATTTTTTTGGGCGACCAAACCGGAGC‑TAGGTTTTGAGAATAATTGGCTGCACTCATGCTGTGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCt                                       <  2:748800/149‑1 (MQ=255)                                                                                                         gtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCta                                      >  1:229527/1‑147 (MQ=255)                                                                                                         gtaaaTGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCta                                      <  2:229527/147‑1 (MQ=255)                                                                                                              tGTTTTACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGc                               <  1:670704/149‑1 (MQ=255)                                                                                                                ttttACTTACCTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATCGCTAAATAAGCTa                             >  2:290003/1‑149 (MQ=255)                                                                                                                   tacgtacCTTTGTATTAATAAAGTATTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTCGAATGATTGGCTGGACGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCACATAGCTAAGTCAGCTAGTc                          <  2:889604/149‑1 (MQ=255)                                                                                                                          cattGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGTTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt                   <  2:114179/147‑1 (MQ=255)                                                                                                                          cTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGc                                                                                                >  1:121701/1‑72 (MQ=255)                                                                                                                          cTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGc                                                                                                <  2:121701/72‑1 (MQ=255)                                                                                                                          cTTTGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACt                   <  1:695490/149‑1 (MQ=255)                                                                                                                             tGTATTAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGt                <  1:714160/149‑1 (MQ=255)                                                                                                                              gTATGACTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATATCTAAGTAAGCTAGTCCCGAACTGGTa               <  2:528112/149‑1 (MQ=255)                                                                                                                                 gtAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGCTGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGCACTGGTAttt            >  2:277807/2‑149 (MQ=255)                                                                                                                                  tAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTa           <  1:308200/149‑1 (MQ=255)                                                                                                                                  tAGTAAATTTTTTTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTa           <  2:225472/149‑1 (MQ=255)                                                                                                                                           ttttTGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTc  >  2:519575/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                  TATTAACACCTACGATGAATGCATAAGAGATGATAACCATAAAGACAATTTGCTTACCAGTATTCAGTATAATTGCAAAGATGAATACAATAAAAACAAGTAAAGTAAATGTTTTACTTACCTTTGTATTAGTAAATTTTTTAGACCAAATAAA‑TAAAGCAAAAGTAGTTAGAATGATTGGCTGCATGTATGCTGAGAAATTTCTTGATGTGTCTTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTC  >  NC_000913_3_pae_tpiA/2106599‑2106885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑