Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	196,351	C→G	14.7%	T619T (ACC→ACG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	196,351	0	C	G	14.7%	120.0 / 14.5	54	T619T (ACC→ACG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (26/20);  new base G (6/2);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAC  >  AE015451/196139‑196584                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          gTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                           <  1:107142/229‑1 (MQ=17)           agaACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGc                                                                                                                                                                                                                         >  1:767523/1‑221 (MQ=39)                gCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACATACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCGAACGGCCAAACCATTACCGTTGAAGCCGGCACAACCAAGGGCa                                                                                                                                                                                                            >  1:22657/1‑229 (MQ=255)                             acttCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCGGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                           <  2:736685/199‑1 (MQ=21)                                   aGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGGCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGa                                                                                                                                                                                    <  1:831679/234‑1 (MQ=16)                                       aCTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGCGCAACCCTGCCCAGACGCCGGTGACCGTGACCCAGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                           <  2:493616/191‑1 (MQ=255)                                         tGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                           <  2:910335/189‑1 (MQ=18)                                           aCCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                           <  2:536659/187‑1 (MQ=18)                                                   cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGACCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                                                                                                                                                <  1:298036/174‑1 (MQ=38)                                                     cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCCGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                                                                                                                                                  >  2:939381/1‑234 (MQ=12)                                                        cGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                           <  1:325101/174‑1 (MQ=38)                                                           cGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                                                                                           <  1:365003/171‑1 (MQ=18)                                                              cTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCCACGGCCACAcctccac                                                                                                                                                                                                                                      >  1:603740/1‑152 (MQ=255)                                                                          tCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCAt                                                                                                                                            >  1:429796/1‑235 (MQ=14)                                                                                                   cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACAGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACaacaac                                                                                                                                                                   >  1:1124166/1‑187 (MQ=18)                                                                                                     cgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCa                                                                                                                                                                                                                 >  1:361560/1‑139 (MQ=34)                                                                                                        cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCa                                                                                                                                                      >  1:1263342/1‑195 (MQ=14)                                                                                                                      ccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGCCGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACaggcgg                                                                                                                          >  2:1137430/1‑209 (MQ=11)                                                                                                                          cGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCgc                                                                                                                                                                                                        >  1:1030860/1‑126 (MQ=255)                                                                                                                                      gcggcgTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGaa                                                                                                                                      >  1:168361/1‑181 (MQ=37)                                                                                                                                             gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGCACGCCACaggcggg                                                                                                                         >  1:380442/1‑186 (MQ=11)                                                                                                                                                     tccccGTGACCCTGAGCAATCCTGCCCAGCCGCCGGTGCCAGTGCCCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGAAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                                                                                                       <  1:211232/130‑1 (MQ=11)                                                                                                                                                              cccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTac                                                                                                                                                                         >  1:365063/1‑122 (MQ=32)                                                                                                                                                              cccTGAGCAATCCTGCCCAGACGCCGGTAACAGTGACCCTGCCCTACGACCCAACCATTACCGTTCAAGCCGGCAAAACCCAGGGCAg                                                                                                                                                                                                           >  1:417151/1‑88 (MQ=11)                                                                                                                                                                cTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGAc                                                                                                                                                                               <  2:1105472/114‑1 (MQ=12)                                                                                                                                                                cTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTAGAGCAAc                                                                                                            >  1:101458/1‑181 (MQ=37)                                                                                                                                                                cTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTCCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGGGTGTCAg                                                                                                                                                     >  1:1091229/1‑140 (MQ=255)                                                                                                                                                                   aGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTac                                                                                                                                                                         <  2:1279584/117‑1 (MQ=12)                                                                                                                                                                   aGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTac                                                                                                                                                                         >  1:1279587/1‑117 (MQ=17)                                                                                                                                                                    gCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCGGCGTCGATTTCCAGACCCCGGCGAATGACGCCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGGAATTTc                                                                                                                   >  1:308983/1‑170 (MQ=255)                                                                                                                                                                    gCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcggcg                                                                                                                           >  1:311903/1‑162 (MQ=38)                                                                                                                                                                    gCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCcaccac                                                  >  1:1096293/1‑235 (MQ=39)                                                                                                                                                                      aaTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGACTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCAc                                                                                                                                >  1:224555/1‑155 (MQ=25)                                                                                                                                                                        tCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGa                                             >  1:876372/1‑236 (MQ=38)                                                                                                                                                                           tGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGa                                                                                                                                                                                             >  1:9255/1‑89 (MQ=33)                                                                                                                                                                           tGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGa                                                                                                                                                                                             <  2:9255/89‑1 (MQ=33)                                                                                                                                                                           tGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCAAGAAAca                                                                                                           >  1:43047/1‑170 (MQ=16)                                                                                                                                                                                     ccGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcca                                                                                                                                                                                    <  2:639220/88‑3 (MQ=16)                                                                                                                                                                                     ccGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCa                                                                                                                                                                                                            >  1:906035/1‑64 (MQ=35)                                                                                                                                                                                     ccGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                                                                                  >  1:276749/1‑90 (MQ=33)                                                                                                                                                                                     ccGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                                                                                  <  2:276748/90‑1 (MQ=33)                                                                                                                                                                                     ccGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACAGTCAGTGTCACCATCGAGAACGCCACaggcg                                                                                                                           >  2:118686/1‑145 (MQ=17)                                                                                                                                                                                          gaccgTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                                                                                  >  1:1268340/5‑85 (MQ=0)                                                                                                                                                                                          gACAGTGACCCTGTCCAACGGCCAAACCATTACCGCTGTCGCCTATACAACCCAGCGCACCGACGAGTCCCAGGCCCCGCCCAATGACCTCCTCAACAACgg                                                                                                                                                                 >  1:776081/1‑102 (MQ=2)                                                                                                                                                                                             tgTGCCCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACg                                                                     <  2:1234658/190‑1 (MQ=255)                                                                                                                                                                                             tgTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCa                                                                                   <  2:1113716/176‑1 (MQ=38)                                                                                                                                                                                             ggTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCTAAACCCAGTGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCCACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGCCCCCGAATCCg                                                                                              <  1:377328/165‑1 (MQ=255)                                                                                                                                                                                                        ccAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAAt                                                                                                                                                                                  <  2:859828/71‑1 (MQ=255)                                                                                                                                                                                                          aaCGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCTTCGACATCCAGAccc                                                                                                                                                                                          >  2:1055153/1‑61 (MQ=11)                                                                                                                                                                                                          aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGa                                                                                                                                                                                    >  1:953201/1‑67 (MQ=35)                                                                                                                                                                                                          aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTa                                                                                                                        >  1:133117/1‑127 (MQ=38)                                                                                                                                                                                                          aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                              <  2:1195568/153‑1 (MQ=39)                                                                                                                                                                                                          aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCg                                                                                              >  1:1195571/1‑153 (MQ=39)                                                                                                                                                                                                           aCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc          <  1:1091524/236‑1 (MQ=21)                                                                                                                                                                                                            cGGCCAAACCATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc          <  1:982899/235‑1 (MQ=9)                                                                                                                                                                                                                    ccATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc  >  1:459348/1‑235 (MQ=39)                                                                                                                                                                                                                    ccATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGCAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACCACTCGGGCGACGCCACCACCGCGAcccc                                         >  1:278792/1‑195 (MQ=255)                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          GTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAC  >  AE015451/196139‑196584

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