Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	167,992	C→A	79.5%	S170Y (TCC→TAC)	fhuA →	ferrichrome outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	167,992	0	C	A	79.5%	174.5 / 47.8	83	S170Y (TCC→TAC)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (8/9);  new base A (32/34);  total (40/43)
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.94e-01

ACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGAT  >  NC_000913/167848‑168137                                                                                                                                                 |                                                                                                                                                  accTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTg                                                                                                                                               >  8:6349/1‑149 (MQ=255) accTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTg                                                                                                                                               <  4:340292/149‑1 (MQ=255)  ccTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGc                                                                                                                                              >  8:35200/1‑149 (MQ=255)     gATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCtt                                                                                                                                            >  8:385104/1‑148 (MQ=255)      aTCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTa                                                                                                                                          <  8:205857/149‑1 (MQ=255)      aTCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGACAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTa                                                                                                                                          >  7:225744/1‑149 (MQ=255)         aTTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACgg                                                                                                                                       >  7:319593/1‑149 (MQ=255)         aTTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACg                                                                                                                                        <  6:364357/148‑1 (MQ=255)             gcgGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTaa                                                                                                                                    >  5:13023/1‑148 (MQ=255)              cgGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTaaaa                                                                                                                                  <  3:172560/149‑1 (MQ=255)                         cAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGCAGTTGCAGGGCAACTTCTATATCGATGCGGTCATTGACCCGTATATGCTGGACCGCGCTGAAATTATGCGTGTCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTggc                                                                                                                       <  3:210421/149‑1 (MQ=255)                                ccAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGt                                                                                                                <  7:35200/149‑1 (MQ=255)                                 cAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGtt                                                                                                               <  7:6349/149‑1 (MQ=255)                                 cAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGACAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGtt                                                                                                               <  8:225744/149‑1 (MQ=255)                                  aaaGCCAGAATAACTATCTGATTGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTg                                                                                                              <  6:13023/149‑1 (MQ=255)                                     gCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAat                                                                                                           <  5:309082/149‑1 (MQ=255)                                     gCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAat                                                                                                           <  8:319593/149‑1 (MQ=255)                                         gAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGtt                                                                                                               <  2:225750/141‑1 (MQ=255)                                         gAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGtt                                                                                                               >  1:225750/1‑141 (MQ=255)                                             aaCTATCTGAATGGACTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGAAACGTATATGCTGGAACGAGCTGAAATTATGCGTTGACCGGTTTAAGTTCTTTACGGTAAAAGCAGTCCTGGAGGCCTGTTGAATATGGt                                                                                                      >  1:160150/1‑146 (MQ=255)                                                tATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGc                                                                                                  >  1:307174/1‑147 (MQ=255)                                                 aTCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAg                                                                                                                                 >  8:183271/1‑115 (MQ=255)                                                 aTCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAg                                                                                                                                 <  7:183271/115‑1 (MQ=255)                                                  tCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGACCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGc                                                                                              >  1:61536/1‑149 (MQ=255)                                                     gAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTc                                                                                           >  6:242271/1‑149 (MQ=255)                                                      aaTGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGCACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTcc                                                                                          >  2:227931/1‑149 (MQ=255)                                                      aaTGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGt                                                                                            <  7:385104/147‑1 (MQ=255)                                                      aaTGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACACGTATATGCTGGAACGCGCTGAAATTATGCCTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGACCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTcc                                                                                          >  7:330720/1‑149 (MQ=255)                                                          gCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGac                                                                                       >  6:187617/1‑148 (MQ=255)                                                          gCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCATGTTGAATATGGTCAGCAAGCGTCCGac                                                                                       <  2:307174/148‑1 (MQ=255)                                                              gAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGaccacc                                                                                   >  6:364576/1‑148 (MQ=255)                                                               aaGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGtt                                                                                                               >  6:351001/1‑119 (MQ=255)                                                               aaGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGtt                                                                                                               <  5:351001/119‑1 (MQ=255)                                                                        ggCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGAAAGCGTCCGACCACCGAACCGCTGa                                                                         >  1:14606/1‑148 (MQ=255)                                                                         gCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGa                                                                         >  5:256319/1‑147 (MQ=255)                                                                           aaCTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAaag                                                                      >  7:118284/1‑148 (MQ=255)                                                                            aCTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAaagaa                                                                    <  6:256319/149‑1 (MQ=255)                                                                            aCTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGGCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAaagaa                                                                    <  5:97073/149‑1 (MQ=255)                                                                             cTTCTATAACGATTCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACTGTAAAAGCAGTCCTGGCGGCTTGTTGAATATTGTCATCAATCTTCCGACCACCGAACCGCTTAaagaat                                                                   <  2:160150/149‑2 (MQ=255)                                                                             cTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAaagaag                                                                   <  6:298480/149‑1 (MQ=255)                                                                                cTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGtt                                                                 >  2:371982/1‑148 (MQ=255)                                                                                 tataACGAGGCGGGCATTGTCCCGTATATGCTGGAACTCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTTCGGTAAAATCAGTCCTGTCGGCCTGTTGAATATGGTCa                                                                                                    <  1:61396/112‑1 (MQ=255)                                                                                  ataACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAg                                                              >  7:322547/1‑149 (MQ=255)                                                                                   taACGATGCGGTCATTGCCCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTTCTTTACGGTAAAAGCAGTCATGGCGGCCTGTTGAATATGGTCAGAAAGCGTCCGACCACCTAACAGCAGAAAGAAGTTCAGt                                                             >  5:175002/1‑149 (MQ=255)                                                                                   taACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCt                                                                                                                  >  5:57383/1‑96 (MQ=255)                                                                                   taACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTc                                                                                                     >  8:336457/1‑109 (MQ=255)                                                                                   taACGATGCGGGCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCt                                                                                                                  <  6:57383/96‑1 (MQ=255)                                                                                      cGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGttt                                                           >  2:45347/1‑148 (MQ=255)                                                                                       gATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTa                                                          >  2:88736/1‑148 (MQ=255)                                                                                        aTGCGGTCATTGCCCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGCACCGCTGAAAGAAGTTCAGTTTa                                                          >  3:87014/1‑147 (MQ=255)                                                                                            ggTCATTAACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCa                                                               >  2:387029/1‑138 (MQ=255)                                                                                            ggTCATTAACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCa                                                               <  1:387029/138‑1 (MQ=255)                                                                                              tCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCgg                                                   >  6:302287/1‑148 (MQ=255)                                                                                                aTTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTa                                                 <  5:302287/148‑1 (MQ=255)                                                                                                  tGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACt                                               <  4:87014/148‑1 (MQ=255)                                                                                                    aCCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGa                                             <  1:88736/148‑1 (MQ=255)                                                                                                    aCCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGa                                             >  8:251684/1‑148 (MQ=255)                                                                                                     cccGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACa                                           <  8:118284/149‑1 (MQ=255)                                                                                                          atatGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGtt                                                                                                               <  4:279341/76‑1 (MQ=255)                                                                                                          atatGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGtt                                                                                                               >  3:279341/1‑76 (MQ=255)                                                                                                           tatGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTg                                      <  8:330720/148‑1 (MQ=255)                                                                                                           tatGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGt                                     >  7:386043/1‑149 (MQ=255)                                                                                                              gCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAaagaag                                                                   >  6:335087/1‑116 (MQ=255)                                                                                                              gCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAaagaag                                                                   <  5:335087/116‑1 (MQ=255)                                                                                                                 ggAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAg                                >  5:296652/1‑148 (MQ=255)                                                                                                                        gcTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCt                                                                                                                          >  7:78607/1‑51 (MQ=255)                                                                                                                        gcTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCt                                                                                                                          <  8:78607/51‑1 (MQ=255)                                                                                                                         cTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTcc                                  >  3:65918/1‑138 (MQ=255)                                                                                                                         cTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTcc                                  <  4:65918/138‑1 (MQ=255)                                                                                                                         cTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGtttt                       <  2:61536/149‑1 (MQ=255)                                                                                                                         cTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGtttt                       <  5:242271/149‑1 (MQ=255)                                                                                                                         cTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGtttt                       <  5:364576/149‑1 (MQ=255)                                                                                                                            aaaTTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAATGAAGTTCAGTTTAAAGCCTGTACTGACAGCCTGTTCCAGACTGGTTTTGa                     <  1:227931/148‑1 (MQ=255)                                                                                                                             aaTTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTTGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGCACCGCTGAAAGAAGTTCAGATTAAAGCCGGTACTGACCTCCTGTTCCAGCCTGGTTTTGACt                   >  5:89580/1‑149 (MQ=255)                                                                                                                                 aTGCGTGGCCCGGTTTACGTGGTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTa                <  6:296652/148‑1 (MQ=255)                                                                                                                                          ccGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGtt      <  1:371982/149‑1 (MQ=255)                                                                                                                                          ccGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGt       >  8:204942/1‑148 (MQ=255)                                                                                                                                          ccGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGt       <  7:204942/148‑1 (MQ=255)                                                                                                                                           cGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTcc                                  <  2:370650/120‑1 (MQ=255)                                                                                                                                           cGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTcc                                  >  1:370650/1‑120 (MQ=255)                                                                                                                                            ggTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTgg    <  1:45347/149‑1 (MQ=255)                                                                                                                                              tttACGTTCTTTACGGTAAAAGCTGTCCTTGAGGACTGTTGAATATGTTCAGCAATCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGgat  <  6:175002/149‑1 (MQ=255)                                                                                                                                              tttACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGgat  <  8:322547/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  ACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGAT  >  NC_000913/167848‑168137

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: ‑