Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,169,380	T→C	9.8%	intergenic (+50/+32)	bhsA → / ← ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,169,380	0	T	C	9.8%	171.9 / 9.5	62	intergenic (+50/+32)	bhsA/ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base C (3/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

GCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTT  >  NC_000913/1169235‑1169519                                                                                                                                                   |                                                                                                                                             gcgcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGcc                                                                                                                                           >  1:443836/1‑149 (MQ=255)     aaaaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCGGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCgtgt                                                                                                                                       >  4:24577/1‑149 (MQ=255)          gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                  <  1:394299/149‑1 (MQ=255)          gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                  <  1:46079/149‑1 (MQ=255)          gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                  <  2:122637/149‑1 (MQ=255)          gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                  <  3:649284/149‑1 (MQ=255)          gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                  <  1:64716/149‑1 (MQ=255)          gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                  <  3:355345/149‑1 (MQ=255)           cGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCggg                                                                                                                                 <  2:144271/149‑1 (MQ=255)            ggATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCATGCTACTGCTGA‑TTTTTTCCTCGAGACATGCCGGGTCGCgggg                                                                                                                                >  2:612876/1‑149 (MQ=255)             gATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGa                                                                                                                               >  4:345409/1‑149 (MQ=255)                 agaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑tttt                                                                                                                          >  3:85393/1‑149 (MQ=255)                  gaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑ttttt                                                                                                                         <  2:252225/149‑1 (MQ=255)                    tGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACACGGCATGTCGCGGGGaaaaaatc                                                                                                                        <  4:581158/149‑8 (MQ=255)                           aaaaTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTTATCCgg                                                                                                                 >  2:557096/1‑149 (MQ=255)                           aaaaTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCggg                                                                                                                >  1:646773/1‑149 (MQ=255)                           aaaaTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCggg                                                                                                                >  1:418547/1‑149 (MQ=255)                               tCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCttt                                                                                                            <  3:181406/149‑1 (MQ=255)                                cTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTa                                                                                                           <  4:27207/149‑1 (MQ=255)                                  ttCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACa                                                                                                         <  4:614163/149‑1 (MQ=255)                                    ccGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGc                                                                                                       <  2:260336/149‑1 (MQ=255)                                      gTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAGTTTATAGATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTTCCCCGAGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCg                                                                                                      >  2:45625/1‑149 (MQ=255)                                      gTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGt                                                                                                     >  3:373650/1‑149 (MQ=255)                                      gTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGt                                                                                                     <  2:443836/149‑1 (MQ=255)                                           aCTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGGCATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTtgtg                                                                                                <  2:418547/149‑1 (MQ=255)                                           aCTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTtgtg                                                                                                >  1:46330/1‑149 (MQ=255)                                            cTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTgg                                                                                               >  1:571213/1‑149 (MQ=255)                                                      cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                     <  4:49159/149‑1 (MQ=255)                                                      cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                     <  4:85393/149‑1 (MQ=255)                                                       ggTCCGAATACCCTCCATGGAACAGCAGTAATTTGTAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgat                                                                                    <  4:176769/149‑1 (MQ=255)                                                        gTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatg                                                                                   >  4:592515/1‑149 (MQ=255)                                                                  ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACa                                                                         >  4:357413/1‑149 (MQ=255)                                                                    tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                  <  4:237665/91‑1 (MQ=255)                                                                    tCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                  >  3:237665/1‑91 (MQ=39)                                                                     ccATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATcc                                                                      <  2:571213/149‑1 (MQ=255)                                                                       aTGGAACTGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTTCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                     <  3:451070/132‑1 (MQ=255)                                                                       aTGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                     >  4:451070/1‑132 (MQ=255)                                                                          gAACAGTAATAATTTATAAATAATCTTTAACCCTCATCAATGCCTGCTTCTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑GTTTTTTATCAGCGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCgg                                                                 <  4:106757/149‑1 (MQ=255)                                                                           aaCAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGc                                                                >  1:34983/1‑149 (MQ=255)                                                                                cagTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATccc                                                           <  3:592515/149‑1 (MQ=255)                                                                                 attAAATTAAAAAGAAAGATTAACCCTTATTAATTGCTGCTAATGCTGA‑TTTTTTCCCCGCGACCAGCCGTGTTGCCGGGATTTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTCCGCCGTCATcaa                                                           <  1:557096/147‑3 (MQ=255)                                                                                   tAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGa                                                        >  2:599661/1‑149 (MQ=255)                                                                                       ttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACgg                                                    >  1:74074/1‑149 (MQ=255)                                                                                            aataaGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCCCACTTGATGACGGCGGACATCCACCGGCATCCCGGGACGGAcatc                                               <  2:183758/149‑1 (MQ=255)                                                                                             ataaGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc                                              <  4:310563/149‑1 (MQ=255)                                                                                               aaGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccat                                            <  1:67626/149‑1 (MQ=255)                                                                                                       aCCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                    <  3:303508/149‑1 (MQ=255)                                                                                                         ccTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGc                                  >  1:325570/1‑149 (MQ=255)                                                                                                         ccTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGc                                  <  2:507096/149‑1 (MQ=255)                                                                                                          cTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGTGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGTACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCa                                 >  1:603332/1‑149 (MQ=255)                                                                                                              ttAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCAc                             <  3:637761/149‑1 (MQ=255)                                                                                                              ttAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGATTGTGGGCTCACTTGATCACGGCGGACATCTACTGGCATCCCGGAACGGACATCCATCACATGTTGCATCAc                             >  1:502058/1‑149 (MQ=255)                                                                                                                            tgctgA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGttt               <  2:546629/149‑1 (MQ=255)                                                                                                                              ctgA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGCACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt             >  2:681691/1‑149 (MQ=255)                                                                                                                              ctgA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt             >  1:144625/1‑149 (MQ=255)                                                                                                                              ctgA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt             >  2:668783/1‑149 (MQ=255)                                                                                                                              ctgA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTCTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACCTCCACCGGCCTCCCGGAACGCACATCCATCCCATGTTGCATCACTTCAGCGTCAGTTtgt             >  2:338761/1‑149 (MQ=255)                                                                                                                               tgA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGATATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgtg            <  1:92205/149‑1 (MQ=255)                                                                                                                                 a‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGtgct          >  2:427952/1‑149 (MQ=255)                                                                                                                                          cccGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCttt  >  2:15092/1‑149 (MQ=255)                                                                                                                                          cccGCGACACGGCATGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGTGTCAGTTTGTGCTGCATCttt  >  2:136450/1‑149 (MQ=255)                                                                                                                                          cccGCGACACGGCATGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCttt  >  4:657692/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                             GCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTT  >  NC_000913/1169235‑1169519

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