Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,703,763	A→C	L182L (CTT→CTG)	ydgJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,703,763	0	A	C	98.7%	310.7 / ‑4.0	78	L182L (CTT→CTG)	ydgJ	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (35/42);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

AGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTT  >  NC_000913/1703620‑1703910                                                                                                                                                |                                                                                                                                                    aGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCcaga                                                                                                                                                  >  2:22110/1‑147 (MQ=255)   cATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGAtgt                                                                                                                                            <  1:235633/151‑1 (MQ=255)      ggTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATgtggt                                                                                                                                         <  1:706372/151‑1 (MQ=255)       gTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATgtggtg                                                                                                                                        <  2:857227/151‑1 (MQ=255)        tACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGc                                                                                                                                       >  2:579929/1‑151 (MQ=255)         aCCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCt                                                                                                                                      >  2:783027/1‑151 (MQ=255)           cGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTa                                                                                                                                     >  2:779804/1‑150 (MQ=255)               tAAAATGACTCCCCGCTGTGGATATTACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGAGCCAGAAGATGTGGTGCTAAAt                                                                                                                                  <  1:427729/149‑1 (MQ=255)                aaaaTGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGt                                                                                                                               <  2:336619/151‑1 (MQ=255)                aaaaGGCCTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGt                                                                                                                               <  1:482233/151‑1 (MQ=255)                      aCTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTAcc                                                                                                                            >  2:1105243/1‑148 (MQ=255)                          gccgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATAcc                                                                                                                     >  2:913230/1‑151 (MQ=255)                          gccgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATAcc                                                                                                                     >  1:498040/1‑151 (MQ=255)                           ccgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCg                                                                                                                    >  2:1051977/1‑151 (MQ=255)                              cTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTg                                                                                                                 >  1:341897/1‑151 (MQ=255)                               tgtgGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGc                                                                                                                >  2:799849/1‑151 (MQ=255)                                  ggATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTg                                                                                                             >  2:997952/1‑151 (MQ=255)                                     tAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGgc                                                                                                          >  2:985631/1‑151 (MQ=255)                                       ggACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTGCCAGATACCGCTGCCTGGgccg                                                                                                        <  2:263280/151‑1 (MQ=255)                                          cAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGccc                                                                                                     >  2:1087753/1‑151 (MQ=255)                                           aaGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCt                                                                                                    <  1:130195/151‑1 (MQ=255)                                           aaGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCt                                                                                                    <  2:1132947/151‑1 (MQ=255)                                            aGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTg                                                                                                   <  1:317216/151‑1 (MQ=255)                                            aGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTg                                                                                                   <  1:969748/151‑1 (MQ=255)                                             gATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGt                                                                                                  <  1:869552/151‑1 (MQ=255)                                              aTGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGtt                                                                                                 >  2:601772/1‑151 (MQ=255)                                               tGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTc                                                                                                >  2:1301198/1‑151 (MQ=255)                                                 gCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCAc                                                                                              >  1:366170/1‑151 (MQ=255)                                                 gCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCAc                                                                                              <  1:347498/151‑1 (MQ=255)                                                      gAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCa                                                                                          >  2:65474/1‑150 (MQ=255)                                                       aaaTAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAAc                                                                                        <  2:439045/151‑1 (MQ=255)                                                        aataatCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACg                                                                                       <  2:1268630/151‑1 (MQ=255)                                                        aataatCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACg                                                                                       >  1:392729/1‑151 (MQ=255)                                                         ataatCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGa                                                                                      <  2:479004/151‑1 (MQ=255)                                                              cGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATcgcg                                                                                 >  1:731289/1‑151 (MQ=255)                                                               ggTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCgcgc                                                                                >  1:239476/1‑151 (MQ=255)                                                               ggTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCgcgc                                                                                <  2:990944/151‑1 (MQ=255)                                                                  cGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCAcc                                                                             >  1:94494/1‑151 (MQ=255)                                                                   gACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCt                                                                            <  2:1085738/151‑1 (MQ=255)                                                                   gACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCt                                                                            <  1:647545/151‑1 (MQ=255)                                                                            ccGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAACTCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGcggacgga                                                                   <  2:949829/151‑1 (MQ=255)                                                                            ccGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGcggacgg                                                                    <  2:1072520/150‑1 (MQ=255)                                                                            ccGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGcggacg                                                                     <  1:311205/149‑1 (MQ=255)                                                                            ccGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGcggacg                                                                     <  2:531007/149‑1 (MQ=255)                                                                                gCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCg                                                               <  2:1295324/151‑1 (MQ=255)                                                                                    tAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTc                                                            >  2:32943/1‑150 (MQ=255)                                                                                    tAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTc                                                            >  2:1226869/1‑150 (MQ=255)                                                                                    tAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTc                                                            >  2:136855/1‑150 (MQ=255)                                                                                     aaCTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCa                                                           <  2:112594/150‑1 (MQ=255)                                                                                         gtgCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCaaaa                                                        <  2:963313/149‑1 (MQ=255)                                                                                         gtgCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCaaaa                                                        <  1:779775/149‑1 (MQ=255)                                                                                           gCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATgag                                                    <  2:33436/151‑1 (MQ=255)                                                                                           gCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATga                                                     <  2:1194282/150‑1 (MQ=255)                                                                                            ccAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATg                                                      >  1:503512/1‑148 (MQ=255)                                                                                              aaaTCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACt                                                 >  2:173803/1‑151 (MQ=255)                                                                                                 tCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCa                                               >  2:1016514/1‑150 (MQ=255)                                                                                                  cTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCaa                                              >  1:777349/1‑150 (MQ=255)                                                                                                     ccGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAg                                            <  1:236653/149‑1 (MQ=255)                                                                                                     ccGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGt                                           >  2:20813/1‑150 (MQ=255)                                                                                                         cATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAATATGTGGGGCGAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATTAGTCTCAAAGTAAGc                                       <  1:534809/150‑1 (MQ=255)                                                                                                         cATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGc                                       <  2:169675/150‑1 (MQ=255)                                                                                                          aTGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCaa                                     <  1:603805/151‑1 (MQ=255)                                                                                                               gACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTc                                 >  2:559611/1‑150 (MQ=255)                                                                                                               gACCGGTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCa                                <  1:216262/151‑1 (MQ=255)                                                                                                                    gTAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAg                           <  1:1161522/151‑1 (MQ=255)                                                                                                                     tAAACCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGc                          >  2:420133/1‑151 (MQ=255)                                                                                                                      aaaCCAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGc                          >  1:1080279/1‑150 (MQ=255)                                                                                                                          cAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGcc                     >  1:872222/1‑151 (MQ=255)                                                                                                                          cAAAAAGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGcc                     >  1:1102169/1‑151 (MQ=255)                                                                                                                               aGCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCcg                <  1:144249/151‑1 (MQ=255)                                                                                                                                gCGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCcgc               <  2:267580/151‑1 (MQ=255)                                                                                                                                 cGTAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCcgcg              >  2:967555/1‑151 (MQ=255)                                                                                                                                   tAATGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCcgcg              >  1:327154/1‑149 (MQ=255)                                                                                                                                    aaTGGCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGt           <  1:1061581/151‑1 (MQ=255)                                                                                                                                       ggCCTGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTaaa        <  2:913683/151‑1 (MQ=255)                                                                                                                                           tGATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCtt    >  2:1057068/1‑151 (MQ=255)                                                                                                                                            gATCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCttt   <  2:907939/151‑1 (MQ=255)                                                                                                                                             aTCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCCCCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCt     <  2:779905/148‑1 (MQ=255)                                                                                                                                             aTCCAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCtttt  <  2:127916/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    AGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCAGATACCGCTGCCTGGGCCGCCCTGTTCACGCCAACGATCGCGCACCTGCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTT  >  NC_000913/1703620‑1703910

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