Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,169,384	G→A	9.7%	intergenic (+54/+28)	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,384	0	G	A	9.7%	174.5 / 7.4	62	intergenic (+54/+28)	bhsA/ldtC	biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base A (3/3);  total (28/34)
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.95e-01

AAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCA  >  NC_000913/1169239‑1169530                                                                                                                                                   |                                                                                                                                                    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)                                                                                                                                                   |                                                                                                                                                    AAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGA‑TTTTTTCCCCGCGACATGCCGTGTCGCGGGGA‑TTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATTGCA  >  NC_000913/1169239‑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: ‑