Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,079,957	T→C	E560G (GAA→GGA)	tktA ←	transketolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,079,957	0	T	C	98.7%	265.6 / ‑4.4	77	E560G (GAA→GGA)	tktA	transketolase 1
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (34/42);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

AACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGG  >  NC_000913/3079821‑3080085                                                                                                                                         |                                                                                                                                 aaCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTg                                                                                                                               <  1:1099822/140‑1 (MQ=255)  aCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGa                                                                                                                              >  2:106800/1‑140 (MQ=255)  aCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGa                                                                                                                              >  2:975670/1‑140 (MQ=255)    cGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAAc                                                                                                                            >  1:411257/1‑140 (MQ=255)          cGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAg                                                                                                                      >  1:1170381/1‑140 (MQ=255)               gTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATg                                                                                                                 <  2:1194725/140‑1 (MQ=255)                    gATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAAt                                                                                                            <  1:107198/140‑1 (MQ=255)                         aCGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGtt                                                                                                       >  2:417651/1‑140 (MQ=255)                          cGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTc                                                                                                      >  1:1009279/1‑140 (MQ=255)                          cGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTc                                                                                                      >  1:794168/1‑140 (MQ=255)                          cGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTc                                                                                                      >  1:129130/1‑140 (MQ=255)                              aaGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGc                                                                                                  >  2:796740/1‑140 (MQ=255)                                 cagcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGa                                                                                               <  2:797926/140‑1 (MQ=255)                                   gcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGAcc                                                                                             >  2:716811/1‑140 (MQ=255)                                     aTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCgg                                                                                           <  2:1042016/140‑1 (MQ=255)                                     aTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCgg                                                                                           >  2:83792/1‑140 (MQ=255)                                       ccTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGgcg                                                                                         <  2:854951/140‑1 (MQ=255)                                       ccTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGgcg                                                                                         >  2:1015179/1‑140 (MQ=255)                                         tGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCATTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCa                                                                                       >  1:1143981/1‑140 (MQ=255)                                         tGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGGGCa                                                                                       <  1:4189/140‑1 (MQ=255)                                         tGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCa                                                                                       >  1:880290/1‑140 (MQ=255)                                          gCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAg                                                                                      <  1:417651/140‑1 (MQ=255)                                          gCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAg                                                                                      <  1:625731/140‑1 (MQ=255)                                          gCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAg                                                                                      <  1:1194679/140‑1 (MQ=255)                                          gCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAg                                                                                      <  2:852510/140‑1 (MQ=255)                                           cTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGt                                                                                     >  2:487424/1‑140 (MQ=255)                                           cTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGt                                                                                     >  2:679202/1‑140 (MQ=255)                                              gTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCtt                                                                                  >  2:863017/1‑140 (MQ=255)                                                cAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTc                                                                                <  2:794168/140‑1 (MQ=255)                                                cAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTc                                                                                <  1:372024/140‑1 (MQ=255)                                                  aaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAg                                                                              <  2:184357/140‑1 (MQ=255)                                                  aaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAg                                                                              <  2:1107782/140‑1 (MQ=255)                                                  aaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAg                                                                              <  1:83792/140‑1 (MQ=255)                                                  aaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAg                                                                              <  1:106800/140‑1 (MQ=255)                                                      cGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGcaca                                                                          >  1:68713/1‑140 (MQ=255)                                                        tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCa                                                                                       <  2:1015306/125‑1 (MQ=255)                                                        tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCa                                                                                       >  1:1015306/1‑125 (MQ=255)                                                        tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATa                                                                        >  2:708186/1‑140 (MQ=255)                                                         cGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATaa                                                                       >  2:729997/1‑140 (MQ=255)                                                            tAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAacca                                                                    <  1:729997/140‑1 (MQ=255)                                                                 ggCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACcgcg                                                               >  2:253730/1‑140 (MQ=255)                                                                   cATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAac                                                                      <  2:611873/131‑1 (MQ=255)                                                                   cATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAac                                                                      >  1:611873/1‑131 (MQ=255)                                                                   cATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACcgcgcg                                                             >  1:434468/1‑140 (MQ=255)                                                                     tGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACcgcgcgcg                                                           >  1:859899/1‑140 (MQ=255)                                                                      ggACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGa                                                          >  1:289551/1‑140 (MQ=255)                                                                      ggACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACcgcgcg                                                             <  2:727902/137‑1 (MQ=255)                                                                      ggACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACcgcgcg                                                             >  1:727902/1‑137 (MQ=255)                                                                        acacCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATg                                                        >  1:1109736/1‑140 (MQ=255)                                                                         caccacGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGt                                                       <  2:1009279/140‑1 (MQ=255)                                                                           ccacGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGttt                                                     >  2:829407/1‑140 (MQ=255)                                                                           ccacGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGttt                                                     <  2:880290/140‑1 (MQ=255)                                                                           ccacGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGttt                                                     <  2:411257/140‑1 (MQ=255)                                                                            cacGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTg                                                    >  1:764487/1‑140 (MQ=255)                                                                             acGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGc                                                   >  2:1181252/1‑140 (MQ=255)                                                                               gcgcgcTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCa                                                 >  1:537401/1‑140 (MQ=255)                                                                                 gcgcTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCATTTCAACTCCTGAACCGGTAGCGACGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGt                                               <  2:1143981/140‑1 (MQ=255)                                                                                   gcTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTg                                             <  1:1181252/140‑1 (MQ=255)                                                                                        cACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCtt                                        <  1:708186/140‑1 (MQ=255)                                                                                             cTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAgtt                                   <  1:391174/140‑1 (MQ=255)                                                                                                cGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAgttcgt                                <  2:129130/140‑1 (MQ=255)                                                                                                 ggCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAgttcgtt                               >  1:459500/1‑140 (MQ=255)                                                                                                       cagtTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCctgct                         >  2:1045584/1‑140 (MQ=255)                                                                                                       cagtTTTTCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCctgct                         <  2:305956/140‑1 (MQ=255)                                                                                                            tttCGTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGcc                    <  2:764487/140‑1 (MQ=255)                                                                                                                gTAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCGCCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGt                <  2:998351/140‑1 (MQ=255)                                                                                                                 tAGGCAGCAACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGtt               <  1:253730/140‑1 (MQ=255)                                                                                                                     cagcaACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGa           <  2:434468/140‑1 (MQ=255)                                                                                                                     cagcaACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGa           <  1:591783/140‑1 (MQ=255)                                                                                                                      agcaACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGAc          <  2:1060056/140‑1 (MQ=255)                                                                                                                        caACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGcc                    <  2:1036822/128‑1 (MQ=255)                                                                                                                        caACAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGcc                    >  1:1036822/1‑128 (MQ=255)                                                                                                                         aaCAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACggg       <  1:679202/140‑1 (MQ=255)                                                                                                                          aCAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGga      <  2:459500/140‑1 (MQ=255)                                                                                                                          aCAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGga      <  2:426415/140‑1 (MQ=255)                                                                                                                          aCAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGga      <  2:210350/140‑1 (MQ=255)                                                                                                                          aCAGCCAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGga      <  1:863017/140‑1 (MQ=255)                                                                                                                              ccAGTTCAACTCCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAgg  <  2:1170381/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                 AACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGG  >  NC_000913/3079821‑3080085

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