Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,846,605	A→G	5.5%	L119P (CTG→CCG)	topB ←	DNA topoisomerase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,846,605	0	A	G	5.5%	226.8 / 2.9	73	L119P (CTG→CCG)	topB	DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGG  >  NC_000913/1846461‑1846753                                                                                                                                                 |                                                                                                                                                     cgcgCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGa                                                                                                                                                    >  2:38451/1‑147 (MQ=255) cgcgCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATa                                                                                                                                                  <  6:377/149‑1 (MQ=255) cgcgCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATa                                                                                                                                                  <  6:48501/149‑1 (MQ=255) cgcgCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGGCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATa                                                                                                                                                  <  5:137720/149‑1 (MQ=255)    gcTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTc                                                                                                                                               <  6:197223/149‑1 (MQ=255)        gcgcCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGc                                                                                                                                           >  6:173547/1‑149 (MQ=255)                  aGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATcca                                                                                                                                 >  5:113624/1‑149 (MQ=255)                   gAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGCCGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATccac                                                                                                                                >  2:238484/1‑149 (MQ=255)                   gAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATa                                                                                                                                                  >  5:283670/1‑131 (MQ=255)                   gAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATa                                                                                                                                                  <  6:283670/131‑1 (MQ=255)                   gAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATccac                                                                                                                                <  6:56242/149‑1 (MQ=255)                      aCGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAg                                                                                                                             >  8:26246/1‑149 (MQ=255)                         cacaGCGGTACAAAGTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCaa                                                                                                                          <  1:77480/149‑1 (MQ=255)                          acaGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAAt                                                                                                                         <  8:240484/149‑1 (MQ=255)                            aGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTg                                                                                                                       <  2:340916/149‑1 (MQ=255)                              cGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGcc                                                                                                                     <  7:26246/149‑1 (MQ=255)                              cGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGcc                                                                                                                     <  6:106650/149‑1 (MQ=255)                                      aCTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTTCGGGTTCAGGTAGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTTCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGa                                                                                                             <  3:10794/149‑1 (MQ=255)                                       cTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGAt                                                                                                            >  1:264520/1‑149 (MQ=255)                                        tCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTCGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATc                                                                                                           >  1:156603/1‑149 (MQ=255)                                          aCTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATc                                                                                                           <  5:173547/147‑1 (MQ=255)                                          aCTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATc                                                                                                           <  7:55815/147‑1 (MQ=255)                                          aCTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCg                                                                                                         >  3:139573/1‑149 (MQ=255)                                            tGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCggg                                                                                                       >  1:296770/1‑149 (MQ=255)                                                ggAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTAGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTccc                                                                                                   >  2:274526/1‑149 (MQ=255)                                                 gAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTcccc                                                                                                  >  4:132781/1‑149 (MQ=255)                                                  aaCGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCg                                                                                                 <  3:132781/149‑1 (MQ=255)                                                      aaGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTCGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGt                                                                                             <  2:156603/149‑1 (MQ=255)                                                        gACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGa                                                                                                                                                    <  5:33578/92‑1 (MQ=255)                                                        gACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGa                                                                                                                                                    >  6:33578/1‑92 (MQ=255)                                                           ggTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGaa                                                                                          <  3:288950/147‑1 (MQ=255)                                                           ggTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGaa                                                                                          >  4:288950/1‑147 (MQ=255)                                                              cGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCGAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGAtt                                                                                     <  6:202426/149‑1 (MQ=255)                                                                   gCCCGCTCTACCTCCTGCGTGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCt                                                                                <  2:332230/149‑1 (MQ=255)                                                                    cccGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCCCCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTg                                                                               >  3:189340/1‑149 (MQ=255)                                                                     ccGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTgg                                                                              >  1:311559/1‑149 (MQ=255)                                                                     ccGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTgg                                                                              >  1:181889/1‑149 (MQ=255)                                                                       gCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCt                                                                            >  3:339961/1‑149 (MQ=255)                                                                           aaCCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCGGATAGTCCa                                                                                                                                             >  7:289049/1‑80 (MQ=255)                                                                           aaCCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCGGATAGTCCa                                                                                                                                             <  8:289049/80‑1 (MQ=255)                                                                           aaCCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCGGATAGTCCa                                                                                                                                             <  4:155264/80‑1 (MQ=255)                                                                           aaCCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCGGATAGTCCa                                                                                                                                             >  3:155264/1‑80 (MQ=255)                                                                           aaCCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCAGGCGTGAACGATTTCGCTGGCTTCAt                                                                        >  2:195784/1‑149 (MQ=255)                                                                              cGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCa                                                                     <  6:113624/149‑1 (MQ=255)                                                                                      ggTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGt                                                             <  4:295392/149‑1 (MQ=255)                                                                                        ttCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGttt                                                           >  4:186220/1‑149 (MQ=255)                                                                                           aGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTtgat                                                        >  2:124271/1‑149 (MQ=255)                                                                                            ggTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCGCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTttatg                                                       <  1:274526/149‑1 (MQ=255)                                                                                            ggTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTtgatg                                                       <  1:38451/149‑1 (MQ=255)                                                                                             gTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTtgatga                                                      >  1:139080/1‑149 (MQ=255)                                                                                              tCGTTTATCAAGCCACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATGCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGATCGGTTTGATGAc                                                     <  3:275555/149‑1 (MQ=255)                                                                                                gTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCATATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGt                                                   >  3:179809/1‑149 (MQ=255)                                                                                                gTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGt                                                   >  6:351320/1‑149 (MQ=255)                                                                                                    aTCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCAAGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTaa                                                <  5:65024/148‑1 (MQ=255)                                                                                                     tCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgt                                              <  1:198072/149‑1 (MQ=255)                                                                                                       aaGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgtt                                             <  2:181889/148‑1 (MQ=255)                                                                                                       aaGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgtt                                             <  2:311559/148‑1 (MQ=255)                                                                                                        aGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgttgt                                           <  4:139573/149‑1 (MQ=255)                                                                                                           aaCGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTg                                        <  3:41075/149‑1 (MQ=255)                                                                                                           aaCGCTGTACCTGCTGGAGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTg                                        >  4:41075/1‑149 (MQ=255)                                                                                                                   aCCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACg                                  >  2:105176/1‑147 (MQ=255)                                                                                                                     ctgctgGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAgg                              <  4:339961/149‑1 (MQ=255)                                                                                                                     ctgctgGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAgg                              <  1:238484/149‑1 (MQ=255)                                                                                                                     ctgctgGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAgg                              >  3:263195/1‑149 (MQ=255)                                                                                                                     ctgctgGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAgg                              >  5:325379/1‑149 (MQ=255)                                                                                                                         tgGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCg                          >  8:176402/1‑149 (MQ=255)                                                                                                                           gcgcTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCg                          <  2:296770/147‑1 (MQ=255)                                                                                                                              cTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCt                     <  2:264520/149‑1 (MQ=255)                                                                                                                                tttCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTg                    >  1:195086/1‑148 (MQ=255)                                                                                                                                 ttCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTg                    <  4:122482/147‑1 (MQ=255)                                                                                                                                             ttGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTATTTGCCACtttt      >  5:353069/1‑149 (MQ=255)                                                                                                                                             ttGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACtttt      >  6:59775/1‑149 (MQ=255)                                                                                                                                                cAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCg   <  7:176402/149‑1 (MQ=255)                                                                                                                                                 aGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCgg  >  3:327724/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGG  >  NC_000913/1846461‑1846753

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