Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	968,252	G→A	K544K (AAG→AAA)	msbA →	ATP‑binding lipopolysaccharide transport protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	968,252	0	G	A	100.0%	141.4 / NA	45	K544K (AAG→AAA)	msbA	ATP‑binding lipopolysaccharide transport protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/21);  total (24/21)

AATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACT  >  NC_000913/968010‑968472                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                             aaTAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGa                                                                                                                                                                                                                      <  1:291233/251‑1 (MQ=255)  aTAAGCTGGCTAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATGGCCCACCGCTTGTCTACCATTGAAAAAGCAGACgc                                                                                                                                                                                                                      >  2:424312/1‑249 (MQ=255)    aaGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAAt                                                                                                                                                                                                                   >  1:219923/1‑251 (MQ=255)         ggATAACGGGCTCGAATCAGGGAGTGGTGAAAACGGCGTGCTGCTCTCTGGTGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCGTTGTTGCGTGATAGCCCGATTCTGATGCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGTGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACAGTACCTCTCTGGTTATTGCCCACCGCTTGTCTACCATTGAAAAAGCATACGAAATCGTgg                                                                                                                                                                                                              <  2:420822/251‑1 (MQ=255)                         aCAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGgtgt                                                                                                                                                                                              >  1:442588/1‑251 (MQ=255)                               aTTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATtgt                                                                                                                                                                                        <  1:65237/251‑1 (MQ=255)                                  ggTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGGGTCATTGTgg                                                                                                                                                                                      >  1:215230/1‑250 (MQ=255)                                    tGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGACCGTGCGATTCAGGCGGCACTGGATGAGTTGACGAAAAACCGTAACTCTCTGGTGATTGCCAACCGCTTGTATACCATTGGAAAAGCAGACGAAATAGTGGTCGTCGGGGGTGGTGTCATTGTGGacc                                                                                                                                                                                   >  2:333174/1‑249 (MQ=255)                                                      ctctGGCGGGCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGGTAGCCCGATTCTGATTCTGGACGGAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGCAGAAATCGTGGACGTCGAGGGTAGAGTCATTGTGGGACGCGGGACGCATAACgttt                                                                                                                                                                 >  2:130939/1‑248 (MQ=255)                                                       cctGGCGGTAAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCTGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGAttt                                                                                                                                                                <  2:230300/250‑1 (MQ=255)                                                                  gcgtcagcgtATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCAcc                                                                                                                                                     >  1:362055/1‑251 (MQ=255)                                                                           tATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGttt                                                                                                                                            >  1:395592/1‑251 (MQ=255)                                                                            aTTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTa                                                                                                                                           <  1:42067/251‑1 (MQ=255)                                                                               gCTATCGCTCGAGCCTTGGTTCGTGATAGCCCTATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTAcg                                                                                                                                         <  2:491069/250‑1 (MQ=255)                                                                                                gttgCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCCCCGCGGCGTTTACGCGCAACTTCACAAAATg                                                                                                                        <  2:356257/250‑1 (MQ=255)                                                                                                    cGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCACCTTCACCAAATGCAGtt                                                                                                                   >  2:485540/1‑251 (MQ=255)                                                                                                            cccGATTCTGATTCTGGACGACGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTAACTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGCCGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGAATACAGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGAATTTTGGACAAt                                                                                                           >  2:495223/1‑251 (MQ=255)                                                                                                               gattctgattctGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGACAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGTGCAACTTCACAACATGCAGTTTGGCCAATGAt                                                                                                        >  1:290980/1‑251 (MQ=255)                                                                                                                  tctgattctgGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCg                                                                                                      >  1:309015/1‑250 (MQ=255)                                                                                                                            gACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAAtctgg                                                                                           >  1:88952/1‑251 (MQ=255)                                                                                                                             aCGAAGCTACCTCGGCTCTGGATAACGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAAtctggt                                                                                          >  1:438329/1‑251 (MQ=255)                                                                                                                                      ccTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAAt                                                                                 <  1:207471/251‑1 (MQ=255)                                                                                                                                        ccGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATcc                                                                               <  1:266236/250‑1 (MQ=255)                                                                                                                                                   tACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGCGGATGGTGTCATTGTGGACCGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACCAAATGCAGTTTGGCCAATGATCGAACAAATCTGTTCTGGTGAATCCCCTTtgtgtcg                                                                    >  2:113427/1‑248 (MQ=255)                                                                                                                                                    aCCGAATCCGAGCGTGCTATTCATGCGGCCCTGGAAGAGTTGCAGAAAAACCGTACCTCTCTGGTGAGTGCCCGCCGCTGGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGGCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTCTACGCGCAACTGCACAAAATGCAGGTTGGCCAATGATCGAAAAAATCTGGTCTGCTTAATCCCCTTTGTggcgg                                                                   <  2:67901/251‑1 (MQ=255)                                                                                                                                                      cGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTggcgtcc                                                                 >  2:18610/1‑248 (MQ=255)                                                                                                                                                            cGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGCGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCACTGATCGAAAAAATCTGGTCTGGTGCATCCCTTTTGTGGCGGCTCTtgct                                                           >  2:224388/1‑251 (MQ=255)                                                                                                                                                             gAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgctg                                                          <  2:477606/251‑1 (MQ=255)                                                                                                                                                                cGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGcc                                                        >  1:264181/1‑250 (MQ=255)                                                                                                                                                                      aTTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCGTGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCAGTGTGGAACGCGTTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGGAGTTTGGCCAATGATCGAAAAAATCGGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTcc                                                 <  1:480347/251‑1 (MQ=255)                                                                                                                                                                          aGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGt                                             <  1:275147/251‑1 (MQ=255)                                                                                                                                                                                            ggcaaaaaaaCCGTGCCACACTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGGGTCATTGTGGAACGCGGTACGCATAACGAGTGGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTTGTCTGGTTAATCCCCTTTGTGGCGGTTATGGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTg                           <  2:33283/246‑1 (MQ=255)                                                                                                                                                                                             gCAGAAAAACCGTCCCTCTCTGGTGATTGCCCACCGCTGGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGGCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCGCCTGGTTGTATGGCCTGGTGAGTg                           <  2:442588/250‑1 (MQ=255)                                                                                                                                                                                             gCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCGTCCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTTGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGTCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCTTTTGTGGCGGCTATTTCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTg                           <  1:488245/250‑1 (MQ=255)                                                                                                                                                                                                       cGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGGGAGTGGCGCGATCCGt                >  1:140618/1‑251 (MQ=255)                                                                                                                                                                                                        gTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTCTTGCTGCCACTCTCCTGGTTGTAttgaccgggggggtgcgcgatccctc               >  2:319933/1‑226 (MQ=255)                                                                                                                                                                                                          aCCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGATTGGCGCGATCCGTCt              <  1:375521/250‑1 (MQ=255)                                                                                                                                                                                                           ccTCTCTGGGGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCAGTGTGGAACGCGGTACGCATAACGATTGGCTTGAGCACCGCTGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGGGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCttt            <  2:215230/251‑1 (MQ=255)                                                                                                                                                                                                              ctctGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGATATTTCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCt         >  1:10355/1‑251 (MQ=255)                                                                                                                                                                                                                  ggTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCtata      >  1:120131/1‑250 (MQ=255)                                                                                                                                                                                                                    tGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGAAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTTTTGCTGCCACTCTCCTGGGTGTATTGCCTGGTGAGTGGAGCGATACGTCTTTGCTATAAAc   >  1:487220/1‑251 (MQ=255)                                                                                                                                                                                                                     gATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGGCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACt  <  1:506548/251‑1 (MQ=255)                                                                                                                                                                                                                         gCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGtt                                            <  1:500195/205‑1 (MQ=255)                                                                                                                                                                                                                            cACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGaaaaa                                                                                                 >  1:174712/1‑149 (MQ=255)                                                                                                                                                                                                                            cACCGCTTGTCTACCATTGAAAAAGCAGACGAAACCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTGTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGaaaaa                                                                                                 <  2:174712/149‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                             AATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACT  >  NC_000913/968010‑968472

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