Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	167,992	C→A	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	100.0%	220.5 / NA	64	S170Y (TCC→TAC)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (30/34);  total (30/34)

CCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGAT  >  NC_000913/167846‑168137                                                                                                                                                   |                                                                                                                                                  ccaccTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGGCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACg                                                                                                                                                 <  7:279176/149‑1 (MQ=255)  caccTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGt                                                                                                                                                >  5:83138/1‑149 (MQ=255)         tCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTAc                                                                                                                                         >  5:294551/1‑149 (MQ=255)           aTTCGCGGCTTTGCGTCAGAAGGCGAAAGCCAGAATAACTCTCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACg                                                                                                                                        <  2:207331/148‑1 (MQ=255)           aTTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACg                                                                                                                                        <  1:178753/148‑1 (MQ=255)           aTTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACg                                                                                                                                        <  5:224082/148‑1 (MQ=255)           aTTCGCGGCTTTGCGGCAGAAGGACAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGAACCGAATATGCTGGGACGAGCTGAAAATATGAGTGGCCCGGTTTACGTGCTTTaagt                                                                                                                                       >  5:268360/1‑146 (MQ=255)            ttCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACgg                                                                                                                                       >  1:230674/1‑148 (MQ=255)               gcgGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACgg                                                                                                                                       <  1:316483/145‑1 (MQ=255)               gcgGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACgg                                                                                                                                       >  2:316483/1‑145 (MQ=255)                 ggCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACTGTAAAAg                                                                                                                                 <  1:243821/149‑1 (MQ=255)                 ggCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAg                                                                                                                                 <  2:142138/149‑1 (MQ=255)                  gCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAg                                                                                                                                 >  4:298626/1‑148 (MQ=255)                         ggCAGAAGGCCAAAGCCAGAATATCTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACTGTAAAAGCAGTCCTg                                                                                                                         <  7:247832/149‑1 (MQ=255)                         ggCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTg                                                                                                                         <  6:294551/149‑1 (MQ=255)                         ggCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTg                                                                                                                         <  3:156772/149‑1 (MQ=255)                         ggCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTg                                                                                                                         <  3:151490/149‑1 (MQ=255)                         ggCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTg                                                                                                                         <  2:103561/149‑1 (MQ=255)                            aGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTg                                                                                                                         <  2:131951/146‑1 (MQ=255)                            aGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATAAGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTg                                                                                                                         >  1:131951/1‑146 (MQ=255)                              aaGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTggcggc                                                                                                                    <  6:193198/149‑1 (MQ=255)                              aaGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTggcggc                                                                                                                    >  6:11874/1‑149 (MQ=255)                                 gCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTg                                                                                                                 >  3:106532/1‑149 (MQ=255)                                   cAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGtt                                                                                                               <  3:291102/149‑1 (MQ=255)                                   cAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTATATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGTACCGGTTTACGTGCTTTACGGTAAAAGCAGTCATGCCGGCCtttt                                                                                                               >  4:104353/1‑149 (MQ=255)                                     aaGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTg                                                                                                              >  6:222406/1‑148 (MQ=255)                                      aGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGaa                                                                                                            <  3:298626/149‑1 (MQ=255)                                              tAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCa                                                                                                    <  4:239557/149‑1 (MQ=255)                                                                  aGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGTCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGaa                                                                                >  3:206610/1‑149 (MQ=255)                                                                    ttGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAg                                                                                                   <  3:46927/128‑1 (MQ=255)                                                                    ttGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAg                                                                                                   >  4:46927/1‑128 (MQ=255)                                                                       cAGGGCAACTTCTATAACGATGCGGTCTTTGACCCGTATATGCTTGTACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAATGCAGTCCTTGCGGACTGTTGAATATTGTCAGCAAGCGTCCGACCACCGAACCGCt                                                                           <  6:268360/149‑1 (MQ=255)                                                                          ggCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCACTCCTGGCGGGCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGa                                                                         >  5:222959/1‑148 (MQ=255)                                                                          ggCAACTTCTATAACGATGCGGTCATT‑ACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGa                                                                         >  1:278214/1‑147 (MQ=255)                                                                                ttCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGt                                                                  >  3:325877/1‑149 (MQ=255)                                                                                  cTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGtt                                                                 <  3:183647/148‑1 (MQ=255)                                                                                    ataACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAg                                                              >  5:279757/1‑149 (MQ=255)                                                                                     taACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCa                                                               <  4:106532/147‑1 (MQ=255)                                                                                     taACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGGGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAg                                                              <  6:83138/148‑1 (MQ=255)                                                                                       aCGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGttt                                                           <  5:222406/149‑1 (MQ=255)                                                                                       aCGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGtt                                                            <  6:279757/148‑1 (MQ=255)                                                                                              ggTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGCACCGCTGAAAGAAGTTCAGTTTAAAGCCg                                                    >  8:141350/1‑149 (MQ=255)                                                                                              ggTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCTCTGAAAGAAGTTCAGTTTAAAGCCg                                                    >  3:170116/1‑149 (MQ=255)                                                                                                tCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCgg                                                   >  1:309356/1‑148 (MQ=255)                                                                                                      aCCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCATGGCGGCCTGTTGATTATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGa                                             >  7:305030/1‑148 (MQ=255)                                                                                                        ccGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAg                                          >  5:118342/1‑149 (MQ=255)                                                                                                               tGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTc                                   <  7:181205/149‑1 (MQ=255)                                                                                                                  tGGAACGCGCTGAAATTATGTGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAg                                >  2:168278/1‑149 (MQ=255)                                                                                                                  tGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTaaa                                                                                                                                   <  8:175939/50‑1 (MQ=255)                                                                                                                  tGGAACGCGCTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTaaa                                                                                                                                   >  7:175939/1‑50 (MQ=255)                                                                                                                          gcTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGc                                                      <  1:207326/119‑1 (MQ=255)                                                                                                                          gcTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGc                                                      >  2:207326/1‑119 (MQ=255)                                                                                                                           cTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGtttt                       <  4:206610/149‑1 (MQ=255)                                                                                                                           cTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGtttt                       <  4:311724/149‑1 (MQ=255)                                                                                                                           cTGAAATTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCCGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGtttt                       <  2:230674/149‑1 (MQ=255)                                                                                                                                aTTATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAACTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGCCTGGTTTTGACtt                  >  7:261238/1‑149 (MQ=255)                                                                                                                                 ttATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGCACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACtt                  >  1:316790/1‑148 (MQ=255)                                                                                                                                 ttATGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACtt                  >  7:290117/1‑148 (MQ=255)                                                                                                                                    tGCGTGGCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAg               <  2:309356/148‑1 (MQ=255)                                                                                                                                         ggCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTc         >  8:23652/1‑149 (MQ=255)                                                                                                                                          gCCCGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTAAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCg        >  3:299910/1‑149 (MQ=255)                                                                                                                                            ccGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGtt      <  4:325877/149‑1 (MQ=255)                                                                                                                                             cGGTTTACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTg     <  6:118342/149‑1 (MQ=255)                                                                                                                                                tttACGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGgat  <  4:170116/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGAT  >  NC_000913/167846‑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: ‑