Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	854,127	G→A	6.0%	A161T (GCC→ACC)	ybiR →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	854,127	0	G	A	6.0%	213.1 / 2.7	67	A161T (GCC→ACC)	ybiR	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base A (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

TAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGAGTTGTCTGGGGCT  >  NC_000913/853981‑854272                                                                                                                                                   |                                                                                                                                                  tAAAAAGACTGTGTGAGATAACGGTTAATCGGCTGATTATTTTTTAGGCGCTGGCAGTCAACGATGGTTCGCTACTGACGCAAATTGGAAACCCTCAAAATATTCTTATCTGGGGACGATCTGGTCTTTCGTTTGCAGGATTTATTGcc                                                                                                                                                 >  6:17663/1‑149 (MQ=255)                  aTCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCgctggc                                                                                                                                >  4:211072/1‑149 (MQ=255)                  aTCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCgctggc                                                                                                                                >  2:50738/1‑149 (MQ=255)                     ccGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGctggctgg                                                                                                                             >  4:85530/1‑149 (MQ=255)                       ggTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGGCGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGgcg                                                                                                                           <  1:231428/149‑1 (MQ=255)                        gTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGgcg                                                                                                                           >  5:358587/1‑148 (MQ=255)                                gCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGc                                                                                                                  >  4:229936/1‑149 (MQ=255)                                 cTGATTATTTTTGAGTCGCTGGCAGTCAACGCTGGTTCGCTACTGACTCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCt                                                                                                                 <  8:343559/149‑1 (MQ=255)                                 cTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGGAACCAGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTGTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCt                                                                                                                 >  1:103016/1‑149 (MQ=255)                                 cTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCt                                                                                                                 >  5:100086/1‑149 (MQ=255)                                  tGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCaaa                                                                                                                                             <  3:120661/120‑1 (MQ=255)                                    attattTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGAc                                                                                                              >  1:21382/1‑149 (MQ=255)                                    attattTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTACCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGAc                                                                                                              >  7:105892/1‑149 (MQ=255)                                          tttGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACgctcct                                                                                                        >  8:68349/1‑149 (MQ=255)                                             gAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACgctcctgct                                                                                                     >  4:182434/1‑149 (MQ=255)                                                  gcTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGc                                                                                                                                         <  3:119097/108‑1 (MQ=255)                                                  gcTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGc                                                                                                                                         >  4:119097/1‑108 (MQ=255)                                                   cTGGCAGTCAACGCGGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCtgtg                                                                                               <  3:211072/149‑1 (MQ=255)                                                      gCAGTCAACGCTTGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTTGTCTTTCTTTTGCCTGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCtc                                                                                            <  5:17663/149‑2 (MQ=255)                                                      gCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTg                                                                                            <  1:184891/149‑1 (MQ=255)                                                       cAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTgg                                                                                           >  8:195693/1‑149 (MQ=255)                                                        aGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTGTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGgt                                                                                          <  1:50738/149‑1 (MQ=255)                                                        aGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTCTGCTGgt                                                                                          >  5:249721/1‑149 (MQ=255)                                                                   gTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTcc                                                                                 <  3:229936/147‑1 (MQ=255)                                                                   gTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCt                                                                               <  2:67395/149‑1 (MQ=255)                                                                       gCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGa                                                                            <  7:195693/148‑1 (MQ=255)                                                                              aCGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGa                                                                     <  6:100086/148‑1 (MQ=255)                                                                              aCGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGAt                                                                    >  8:139282/1‑149 (MQ=255)                                                                              aCGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTCGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGa                                                                     <  1:226458/148‑1 (MQ=255)                                                                               cGCCAATTGGCAACCCGGAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTACCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATg                                                                   <  8:105892/149‑1 (MQ=255)                                                                                gCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATg                                                                   >  5:222115/1‑148 (MQ=255)                                                                                 ccAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGc                                                                  >  8:182716/1‑148 (MQ=255)                                                                                  cAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCaa                                                                <  7:182716/149‑1 (MQ=255)                                                                                  cAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCaa                                                                <  8:224816/149‑1 (MQ=255)                                                                                        gCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCa                                                           >  7:265755/1‑148 (MQ=255)                                                                                            cccGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACgg                                                      <  3:182222/149‑1 (MQ=255)                                                                                              cGCAAAATCTTCTTATCTGGGGACGTGCTGGGCTGTCGTTTGCCGGCTTTCTTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCGCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACgggg                                                    <  5:358540/149‑1 (MQ=255)                                                                                              cGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACgggg                                                    <  7:139282/149‑1 (MQ=255)                                                                                              cGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTACCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACgctcctgctc                                                                                                    >  5:198111/1‑101 (MQ=255)                                                                                              cGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTACCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACgctcctgctc                                                                                                    <  6:198111/101‑1 (MQ=255)                                                                                                cAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGt                                                  >  7:240766/1‑149 (MQ=255)                                                                                                     tatTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCaaa                                             >  7:176245/1‑149 (MQ=255)                                                                                                     tatTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCaaa                                             <  1:175611/149‑1 (MQ=255)                                                                                                     tatTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCaaa                                             >  3:141152/1‑149 (MQ=255)                                                                                                      atTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCaaa                                             >  3:120515/1‑148 (MQ=255)                                                                                                          ttATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCg                                        >  3:188968/1‑149 (MQ=255)                                                                                                          ttATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCg                                        >  1:317146/1‑149 (MQ=255)                                                                                                             tCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGa                                      <  3:85530/148‑1 (MQ=255)                                                                                                                ggggACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTg                                   >  2:308521/1‑148 (MQ=255)                                                                                                                 gggACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGt                                    <  4:141152/146‑1 (MQ=255)                                                                                                                 gggACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTgg                                  <  2:320749/148‑1 (MQ=255)                                                                                                                 gggACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTgg                                  <  8:265755/148‑1 (MQ=255)                                                                                                                           ggTCATTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTCGCTGGCGCAATGATGCTGACGCTAATGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACAATACGGGGGTGCAAACACCCGAGTGGAAACCGCGGCt                       >  7:181850/1‑149 (MQ=255)                                                                                                                            gTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTg                      <  6:358587/149‑1 (MQ=255)                                                                                                                                 tCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGgt                    <  3:182434/146‑1 (MQ=255)                                                                                                                                  cGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGtgtg                 <  4:120515/148‑1 (MQ=255)                                                                                                                                   gattGCCGGATGAATTGCCGAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGa               <  2:103016/147‑1 (MQ=255)                                                                                                                                   gTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTGCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGa               <  2:21382/149‑1 (MQ=255)                                                                                                                                   gTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGa               <  8:240766/149‑1 (MQ=255)                                                                                                                                    tttGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGGGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGAg              <  6:222115/149‑1 (MQ=255)                                                                                                                                       gCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGAGtt            <  4:188968/148‑1 (MQ=255)                                                                                                                                       gCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGAGTTg           >  7:213915/1‑149 (MQ=255)                                                                                                                                          ggATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGAGTTGTCt        >  5:48978/1‑149 (MQ=255)                                                                                                                                           gATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGAGGCTGGTGTGGAGTTGTCTg       >  7:157805/1‑149 (MQ=255)                                                                                                                                            aTTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGAGTTGTCTgg      <  2:317146/149‑1 (MQ=255)                                                                                                                                              ttATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGAGTTGTCTggg     <  6:48978/148‑1 (MQ=255)                                                                                                                                                 ttGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGAGTTGTCTGGGGCt  <  1:308521/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCGTTTGCCGGATTTATTGCCCAAATGGCACCGCTGGCTGGCGCAATGATGCTGACGCTCCTGCTCCTGTGCTGGTGTTGTTTCCCTGGAAAGGCGATGCAATACCATACGGGGGTGCAAACACCGGAGTGGAAACCGCGGCTGGTGTGGAGTTGTCTGGGGCT  >  NC_000913/853981‑854272

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