Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	5,664,432	G→C	6.8%	intergenic (+27/+24)	PP_4971 → / ← PP_4972	putative Outer membrane‑bound lytic murein transglycolase A/conserved exported protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	5,664,432	0	G	C	6.8%	185.3 / 4.4	56	intergenic (+27/+24)	PP_4971/PP_4972	putative Outer membrane‑bound lytic murein transglycolase A/conserved exported protein of unknown function
Reads supporting (aligned to +/- strand):  ref base G (27/25);  new base C (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

CGGACGGCACGCCGGTGGTACGACCGGTTGGCGCGCAGGACACCGGCGGGGCGATTACCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGCCCACGATATCCGCCACAAGCACCGCCGCAGCAAATGCCGGAAACG  >  AE015451/5664209‑5664661                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      cGGACGGCACGCCGGTGGTACGACCGGTTGGCGCGCAGGACACCGGCGGGGCGATTACCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATg                                                                                                                                                                                                                           >  1:557072/1‑236 (MQ=255)    aCGGCACGCCGGTGGTACGACCGGTTGGCGCGCAGGACACCGGCGGGGCGATTACCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATAGcc                                                                                                                                                                                                                        <  1:576437/236‑1 (MQ=255)                   tACGACCGGTTGGCGCGCAGGACACCGGCGGGGCGATTACCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCcgc                                                                                                                                                                                                                                  >  1:575691/1‑211 (MQ=255)                                    cAGGACACCGGCGGGGCGATTACCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACg                                                                                                                                                                                         <  1:4213/235‑1 (MQ=255)                                     aGGACACCGGCGGGGCGATTACCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGGCA‑GC‑CCCCGCGAAGATGGCCAGCCGGTTTCAGATGGaaa                                                                                                                                                                                                     <  2:1071840/222‑1 (MQ=255)                                     aGGACACCGGCGGGGCGATTACCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑CAGC‑CCCCGCGAAGATGGCCAGCCGGTTTCAGATGGaaa                                                                                                                                                                                                     >  1:1071843/1‑222 (MQ=255)                                               cggGGCGATTACCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGt                                                                                                                                                                              >  1:723044/1‑235 (MQ=255)                                                         aCCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGc                                                                                                                                                                   <  1:920021/236‑1 (MQ=255)                                                                        gcCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCCGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACga                                                                                                                                                                                               >  1:962003/1‑193 (MQ=255)                                                                         ccGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCcg                                                                                                                                                                                                                                   <  2:482453/156‑1 (MQ=255)                                                                                              ggCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACgaaga                                                                                                                                                                                            >  2:806415/1‑174 (MQ=255)                                                                                              ggCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACg                                                                                                                                                                                                >  1:231006/1‑170 (MQ=255)                                                                                              ggCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCg                                                                                                                                 >  2:370490/1‑233 (MQ=255)                                                                                              ggCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGCACCTTGATCGTTGTGGCGATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGa                                                                                                                                                                                                           >  1:255949/1‑159 (MQ=255)                                                                                              ggCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAAGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTGGCGGG‑G‑GCGGCCcgc                                                                                                                                                                                                                                  >  2:759898/1‑136 (MQ=255)                                                                                               gCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCa                                                                                                                             <  1:1005262/236‑1 (MQ=255)                                                                                                     cgCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCa                                                                                                                                          <  1:741487/216‑1 (MQ=255)                                                                                                            gAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCCGG‑G‑GCTGCCCGCGGAGATGGCCAGCCGGTTTCCGATGGAAACCACGAAGAACGACCCg                                                                                                                                                                                    >  2:149616/1‑168 (MQ=255)                                                                                                                gggCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCg                                                                                                                                 <  1:618112/214‑1 (MQ=255)                                                                                                                                 cggcaAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGc                                                                                                                                                                       <  2:537508/158‑1 (MQ=255)                                                                                                                                 cagcaAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCgg                                                                                                   >  2:559532/1‑228 (MQ=255)                                                                                                                                       gggCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGgcagca                                                                                                                      >  2:353172/1‑203 (MQ=255)                                                                                                                                             aTCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCg                                                                                                                                 <  2:378733/186‑1 (MQ=255)                                                                                                                                               cTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCAccccc                                                                              <  1:472206/235‑1 (MQ=255)                                                                                                                                                   aTGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGCGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGcac                                                                         <  1:430543/236‑1 (MQ=255)                                                                                                                                                       tgtgGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGcac                                                                         <  2:723041/232‑1 (MQ=255)                                                                                                                                                            cccAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTg                                                                                                                                                                     >  2:477816/1‑135 (MQ=255)                                                                                                                                                            cccAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTg                                                                                                                                                                     <  1:477819/135‑1 (MQ=255)                                                                                                                                                            cccAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTg                                                                                                >  1:873470/1‑204 (MQ=255)                                                                                                                                                                  ggCCAGCCACGGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAg                                                                                                                                             <  2:185198/153‑1 (MQ=255)                                                                                                                                                                   gccagccaCTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGGCAGCCGGGTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCCgg                                                                                                                           >  2:183133/1‑170 (MQ=255)                                                                                                                                                                        ccaCTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGAcc                                                                                                                                                  >  1:9393/1‑142 (MQ=255)                                                                                                                                                                           gtGCATGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCa                                                                                                                             <  2:346324/159‑1 (MQ=255)                                                                                                                                                                             gCCTGAGGTACCCAAGGTACCTTGATCGTGGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGt                                                                                              <  2:99474/189‑1 (MQ=255)                                                                                                                                                                                tGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTggg                                                                 >  1:88789/1‑215 (MQ=255)                                                                                                                                                                                tGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGCCCa                                            >  1:993176/1‑236 (MQ=255)                                                                                                                                                                                    gTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCCGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGTAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGGAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAaagcaagc                                                                                                     <  2:568173/175‑1 (MQ=255)                                                                                                                                                                                     tACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGCCCACg                                          >  2:645870/1‑233 (MQ=255)                                                                                                                                                                                       cccAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCg                                                      <  1:380827/219‑1 (MQ=255)                                                                                                                                                                                          aaGGTACCTTGATCGTTGTGGCTATCTTCGCGGGCA‑GC‑CCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGCCCACGATATc                                     >  2:931924/1‑233 (MQ=255)                                                                                                                                                                                           aGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCAtt                                                 >  2:450715/1‑220 (MQ=255)                                                                                                                                                                                               aCCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGCCCACGATATCCGCcaca                              >  1:676677/1‑235 (MQ=255)                                                                                                                                                                                                ccTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGt                                                                                                                          >  1:989079/1‑142 (MQ=255)                                                                                                                                                                                                ccTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGt                                                                                                                          <  2:989076/142‑1 (MQ=255)                                                                                                                                                                                                    gATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAaagcaag                                                                                                      >  1:632242/1‑158 (MQ=255)                                                                                                                                                                                                                     tCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGCCCACGATATCCGCCACAAGCACCGCCGCAGCAAATGcc         >  2:457013/1‑234 (MQ=255)                                                                                                                                                                                                                       gcgGG‑G‑GCTGCCCGCGAAGATGGCCGGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGt                                                                <  1:166226/177‑1 (MQ=255)                                                                                                                                                                                                                       gcgGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTGTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCCACCGCAAgg                                                                                                                                            >  2:76445/1‑101 (MQ=255)                                                                                                                                                                                                                       gcgGG‑G‑GCTGCCCGCGAAGAGGGCCAGCCGGTTGCAGATGGAAGCCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAA‑CACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATAc                                                                                       <  2:48019/153‑1 (MQ=255)                                                                                                                                                                                                                         ggg‑g‑gCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAATAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCg                                                      >  2:440900/1‑185 (MQ=255)                                                                                                                                                                                                                         ggg‑g‑gCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGccc                                             <  2:769680/194‑1 (MQ=255)                                                                                                                                                                                                                          gg‑g‑gCTCCCCGCGAAGATGCCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGa                                                                                                                                                    <  2:272424/90‑1 (MQ=255)                                                                                                                                                                                                                           g‑g‑gCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCg                                                      <  2:231003/183‑1 (MQ=255)                                                                                                                                                                                                                           g‑g‑gCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGCCCACGATATCCGCCACAAGCAccgc                      >  2:1055607/1‑215 (MQ=255)                                                                                                                                                                                                                             g‑gCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGCCCACGATATCCGCCACAAGCACCGCCGCAGCAAATGCCGGAAACg  <  2:580518/234‑1 (MQ=255)                                                                                                                                                                                                                                cTGCCCGCGACGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGACAAACCACACCAGCGACGGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCgg                                                                                                   <  2:955373/135‑1 (MQ=255)                                                                                                                                                                                                                                  gCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACg                                                                                                                                                                                                <  2:531225/40‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      CGGACGGCACGCCGGTGGTACGACCGGTTGGCGCGCAGGACACCGGCGGGGCGATTACCGGCGAAGTGCGTGCCGACCTGTTCTGGGGCACCGGGCCGGAAGCAGGCGAACTGGCCGGGAACATGAAGCAGCAAGGGCAGATCTGGATGCTGTGGCCCAAGGGCCAGCCACTGCCTGAGGTACCCAAGGTACCTTGATCGTTGTGGCTATCTTCGCGGG‑G‑GCTGCCCGCGAAGATGGCCAGCCGGTTTCAGATGGAAACCACGAAGAACGACACGATCAGTGCCAGGCTGGCAAACCACACCAGCGACCGCAAGGCAGCCCAGTCGGCCAGGTAGCAGATGATGTAAAGCAAGCGGCTGGTGATATACATCACCCCCAGCACATCCTGGGTCACCTGCTCGGCATTGCCCACGATATCCGCCACAAGCACCGCCGCAGCAAATGCCGGAAACG  >  AE015451/5664209‑5664661

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