Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,169,382	C→G	9.9%	intergenic (+52/+30)	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,382	0	C	G	9.9%	191.5 / 11.7	61	intergenic (+52/+30)	bhsA/ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base G (3/3);  total (29/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 = 8.39e-01

GCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCA  >  NC_000913/1169235‑1169530                                                                                                                                                     |                                                                                                                                                      gcgcAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGcc                                                                                                                                                      >  1:443836/1‑149 (MQ=255)     aaaaaGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCGGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCgtgt                                                                                                                                                  >  4:24577/1‑149 (MQ=255)          gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                             <  3:649284/149‑1 (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                                                                                                                                             <  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:418547/1‑149 (MQ=255)                           aaaaTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCggg                                                                                                                           >  1:646773/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:85393/149‑1 (MQ=255)                                                      cGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtga                                                                                                <  4:49159/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                        >  2:668783/1‑149 (MQ=255)                                                                                                                              ctgA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt                        >  1:144625/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)                                                                                                                                                     ccGTGTGGCGGGGA‑TTTTTTTATCCGGGCTTTTCAGCGTTTGTGGGCTCACTTGATGTCGGCGGACATCCATCGGCATCCCGGAACGGACATCCATCACATCTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgca  <  1:338761/149‑1 (MQ=255)                                                                                                                                                     ccGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGattgca  <  2:144625/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCA  >  NC_000913/1169235‑1169530

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