Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	396,540	G→A	A208T (GCA→ACA)	USA300HOU_0364 →	possible lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	396,540	0	G	A	100.0%	123.5 / NA	38	A208T (GCA→ACA)	USA300HOU_0364	possible lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/19);  total (19/19)

GATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCA  >  CP000730/396408‑396677                                                                                                                                     |                                                                                                                                          gATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACaaaaa                                                                                                                                   <  1:799991/141‑1 (MQ=255)    atCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAtt                                                                                                                                >  1:452531/1‑141 (MQ=255)             aaaTTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTa                                                                                                                       >  1:604194/1‑141 (MQ=255)             aaaTTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTa                                                                                                                       >  1:30465/1‑141 (MQ=255)             aaaTTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTa                                                                                                                       <  1:274094/141‑1 (MQ=255)               aTTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAACCTCGCAGGTGAAGACGTGATTGATTCACATACAGAGTTATATGATTGTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAtt                                                                                                                     >  1:915357/1‑141 (MQ=255)                ttAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAttt                                                                                                                    >  1:320987/1‑141 (MQ=255)                  aaTGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAACAAATTTATGACTTATTTaa                                                                                                                  <  1:692600/141‑1 (MQ=255)                            ggTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAACAATTTATGACTTATTTAAACCTATTTTa                                                                                                        <  1:597766/141‑1 (MQ=255)                            ggTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTa                                                                                                        <  1:275773/141‑1 (MQ=255)                            ggTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTa                                                                                                        <  1:606535/141‑1 (MQ=255)                                 tGTTGACCTATTAAAAGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaa                                                                                                   >  1:686426/1‑141 (MQ=255)                                                     ttGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtg                                                                               <  1:77689/141‑1 (MQ=255)                                                       gCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATGTTAGAGAAAAAAGATAAAAAATAAAGtgat                                                                             >  1:404102/1‑141 (MQ=255)                                                           cTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGata                                                                         <  1:412222/141‑1 (MQ=255)                                                              cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATcc                                                                      >  1:534550/1‑141 (MQ=255)                                                                aaaaTTACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaa                                                                    <  1:716300/141‑1 (MQ=255)                                                                aaaaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaa                                                                    >  1:241746/1‑141 (MQ=255)                                                                 aaaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaaa                                                                   >  1:646155/1‑141 (MQ=255)                                                                   aTCACAGGTGAAGAACAAATTTATTCACATACAGCTTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATATAAACCTCTTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATg                                                                 <  1:564328/141‑1 (MQ=255)                                                                                    aaTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGaa                                                <  1:205693/141‑1 (MQ=255)                                                                                     aTTTATTCACATACAGATTGATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAAt                                               <  1:652656/141‑1 (MQ=255)                                                                                        tATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCaa                                            <  1:95618/141‑1 (MQ=255)                                                                                          ttCACATACAGATTTATATCATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAatt                                          <  1:249604/141‑1 (MQ=255)                                                                                                tACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGa                                    <  1:324577/141‑1 (MQ=255)                                                                                                  cAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATa                                  <  1:589753/141‑1 (MQ=255)                                                                                                  cAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATa                                  >  1:497929/1‑141 (MQ=255)                                                                                                           atGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAg                         >  1:798046/1‑141 (MQ=255)                                                                                                                   aaaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGCATAAATATAAAGATaacaac                 <  1:565233/141‑1 (MQ=255)                                                                                                                    aaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcg                >  1:743326/1‑141 (MQ=255)                                                                                                                     aGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATGAATTATTGGATAAATATAAAGATAACAAcgg               >  1:251196/1‑141 (MQ=255)                                                                                                                     aGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcgg               >  1:115252/1‑141 (MQ=255)                                                                                                                      gCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcggc              >  1:278486/1‑141 (MQ=255)                                                                                                                            gTTGAAGGCACACAAAACATTTATCAGTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTAt        >  1:449127/1‑141 (MQ=255)                                                                                                                             ttGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATg       >  1:315918/1‑141 (MQ=255)                                                                                                                             ttGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATg       >  1:316120/1‑141 (MQ=255)                                                                                                                                  ggCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCa  <  1:890028/141‑1 (MQ=255)                                                                                                                                  ggCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTCGATAAATATAAAGATAACAACGGCGGTTATGAGTCa  <  1:10903/141‑1 (MQ=255)                                                                                                                                     |                                                                                                                                          GATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATGAGTCA  >  CP000730/396408‑396677

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