Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,046,506	G→T	V523V (GTG→GTT)	yeeJ →	inverse autotransporter adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,046,506	0	G	T	95.3%	141.9 / NA	43	V523V (GTG→GTT)	yeeJ	inverse autotransporter adhesin
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (26/15);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTG‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATGGTGACGGAAGCTATACCCAGATCCTGACCACAGGAGCGATGTCTG  >  NC_000913/2046266‑2046750                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      aGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGAGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGCAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAgg                                                                                                                                                                                                                                            >  1:146673/1‑251 (MQ=255)   gagaAGTGAAATCACTGGTTTCGTCGCTCCAAACAAAATATGCCCTGAAAGGCTATAACGTCGAATCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCa                                                                                                                                                                                                                                          <  2:601852/251‑1 (MQ=255)   gagaAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCa                                                                                                                                                                                                                                          <  1:506247/251‑1 (MQ=255)     gaAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCAcc                                                                                                                                                                                                                                        >  1:368037/1‑251 (MQ=255)      aaGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCt                                                                                                                                                                                                                                       <  1:335476/251‑1 (MQ=255)        gTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTAc                                                                                                                                                                                                                                     >  1:535654/1‑251 (MQ=255)          gAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTGCCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACg                                                                                                                                                                                                                                    <  1:118616/250‑1 (MQ=255)                      ttCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGaaa                                                                                                                                                                                                                        <  1:307147/250‑1 (MQ=255)                         gtcgCTACAAACCAAATATGCCCTGAAAGGCTCTAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGAtt                                                                                                                                                                                                                    >  2:450042/1‑251 (MQ=255)                         gtcgCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGAtt                                                                                                                                                                                                                    <  2:493862/251‑1 (MQ=255)                             cTACAATCCAAATATGCCCTGAAAGGCTATAACGTCGAAGCAACCGCACTGGAAGCTGCCGGTGGCAACGTGGTCACAACGGGGAAAGATCTTCTGGATACCCTGCCGGCTTACCGGTTAACCAGTACGCCAGAAACC‑GATAACACCTGGACGATTGGAGTCACC‑GACGAAGATGTCAAAGGCAATTTGTCGACTCGTGAACAGAGCATGGTT‑GTCGTTCAGGAACCTACGTAAAGCCAGAAAGATTCCt                                                                                                                                                                                                                 >  1:680537/1‑250 (MQ=255)                                        atatGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTaa                                                                                                                                                                                                     <  1:359930/251‑1 (MQ=255)                                          atGCCCTGAAAGGCTATCACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCAAGACAGATTCCTCGGTATCGTTaatt                                                                                                                                                                                                   >  2:484200/1‑249 (MQ=255)                                                tGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAACGATATTCTGGTTACCCTGACGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTCTCGTTAAGTACCCa                                                                                                                                                                                              >  1:319470/1‑250 (MQ=255)                                                 gcaaGGCTATAACGTCGAAGCCCCCGCACGGGAAGCTGTCGGTGGCAAAGTGGTCCCAACGGGTTAAGATATTCGGGTGACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCaa                                                                                                                                                                                             <  1:379104/248‑1 (MQ=255)                                                         ataACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCt                                                                                                                                                                                                                 >  2:337853/1‑222 (MQ=255)                                                         ataACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCt                                                                                                                                                                                                                 <  1:337853/222‑1 (MQ=255)                                                                  aaGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTATCGCATCGTGACCAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTcc                                                                                                                                                                           >  1:715957/1‑251 (MQ=255)                                                                     ccACAGCACTGGAAGCTGCCGGTGGAAAAGTGGTCACAACGGGTAGAGATATTCTGGTTACCCTGCCGGCTTACCGGATCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGACAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCAt                                                                                                                                                                        >  2:34790/1‑251 (MQ=255)                                                                      cACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCAtt                                                                                                                                                                       >  2:121329/1‑251 (MQ=255)                                                                          gCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGTAACC‑GATAACACCTGGCCGCTTGAAGTCACC‑GCCGAAGGTGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTAACCAAAAATTGAACGCGGATTCCCATTCAAc                                                                                                                                                                   >  1:732495/1‑251 (MQ=255)                                                                              tGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGCTGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCAATTCAACCGcc                                                                                                                                                               >  2:308736/1‑251 (MQ=255)                                                                              gggAAGCAGCCGGTGACAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGGT‑GCCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGcc                                                                                                                                                               <  1:722667/250‑1 (MQ=255)                                                                                           tGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTAAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACtttt                                                                                                                                                   >  1:1382/1‑250 (MQ=255)                                                                                                         aaaagggTAAAGATATTCTTGTTCCCCCGCCCGCCTTCCGGTTCCCCAATAACCCAGGAACCAGAAAACCCCTGGCGGGTTGAAGTCCCC‑GCCGAAGATTTCAAAGCCAATTGGTCAAGTCGTAAAAATAGCAGTGTT‑GCGGTTCAGCCACCGCCGCTAGGCCAGAAAGATCCCGCGTTAGCGTTAAGTCCCCAAAAATTGAAGGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATgc                                                                                                                                     <  2:13620/247‑1 (MQ=255)                                                                                                                        ttCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCtgtt                                                                                                                     >  1:353293/1‑251 (MQ=255)                                                                                                                           tGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTc                                                                                                                  >  2:531073/1‑251 (MQ=255)                                                                                                                                       cGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGctc                                                                                                      >  1:266572/1‑251 (MQ=255)                                                                                                                                         gctaaccgttccaccatgcccccAGAAACC‑AATAACTCCTGGCCGATTGAAGCCCCC‑GCCGAAGATTTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCCCCTACGCTAAGCAAGAAAATTTCTTCGGTATCGTAAAGTCCCCAAACATTAAACGCGGATCCCCATCCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCtctc                                                                                                    <  1:93944/230‑1 (MQ=255)                                                                                                                                              ccGGTTCACCAGAACGCCAGACACC‑GATACCACCTGGCCGATTGAAGTCACC‑GCCGCAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGAcgc                                                                                               >  2:699641/1‑251 (MQ=255)                                                                                                                                                    caccaggccgcaaaaaaCC‑GATAACCCCTGGCCGATTGAAG‑CACCCGCCAAAGATGTAAAAGGCAATTTGTCAAATCGTAAACAGAGCGTGGTT‑GTCGTTCAGGCCCCTACGCTAAGCCAAAAAGATCCCTCGGTATCGTTAAGTCCCCAAACATTGAACGGGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACg                                                                                         <  2:353293/237‑1 (MQ=255)                                                                                                                                                    cACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCtctc                                                                                                    >  1:221882/1‑240 (MQ=255)                                                                                                                                                            cGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGtt                                                                                 >  1:174207/1‑251 (MQ=255)                                                                                                                                                                                               cccc‑gccgAAGATGTAAAAGGCAATTTTGAGAATCGTGAACAGAGCAGGGTT‑GTCGTTCAGGCACCTACGCAAAGCCAAAAAGATTCCTCGGTATCGTTAAGTCCCCAAACATGGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATg                                               <  1:450042/249‑1 (MQ=255)                                                                                                                                                                                                    gcggaagtggtaaaaGGCAATTTGTCGAATCGTAAACAGAGCATGGTT‑GTCTTACAGGCCCCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCAAAACATTGACCGCGGATCCCCATTCATCCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATGGTGa                                           <  2:547706/239‑1 (MQ=255)                                                                                                                                                                                                         aGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCCGTCTCGTTAAGCACCCCAACATTGAACGCGCATTCCCATTCAACCGCCACACTCCCTTTTATTCCGCATGATGCAGCAGGTACTCCTGTTGTCGGGCTGGTGCTCTCGACTCGCCACGAAGGTGTCCAGGACATCACCCTTTCTGACTGGAACGCTAATGGTgctgga                                       >  1:238682/1‑245 (MQ=255)                                                                                                                                                                                                           aTGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATGGTGACGGAAGc                                    >  2:172338/1‑251 (MQ=255)                                                                                                                                                                                                                   gccAATTTTTGGAATCGTGAACAGAGCATGGTG‑GTCGTCCAGGCACCTACGCTAAGCCAAAAAGTTCCCCCGTTATCTTAAAGTACCCAAACATTGAACGCGAATCCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATGGTGACGGAAGCTATACCCa                            <  1:121329/249‑1 (MQ=255)                                                                                                                                                                                                                               aaTCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAACGATAATGGTGACGGAAGCTATACCCAGATCCTGACcac                >  2:614275/1‑251 (MQ=255)                                                                                                                                                                                                                                aTCGTGAACAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATGGTGACGGAAGCTATACCCAGATCCTGACcac                >  1:275215/1‑250 (MQ=255)                                                                                                                                                                                                                                  attcaccAGAGCA‑GGTGTGTCGTCCAGCCACCACCGCTAAGCCAAAAATATCCCTGGGTATCTTTAAGTACCCAAACATGGACCGCGGATCCCCATTCACCCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATGGTGACGGAAGCTATACCCAGATCCTGACCACAgg             <  2:283467/245‑1 (MQ=255)                                                                                                                                                                                                                                       aCAGAGCATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATGGTGACGGAAGCTATACCCAGATCCTGACCACAgg             >  1:283467/1‑246 (MQ=255)                                                                                                                                                                                                                                             cATGGTT‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATGGTGACGGAAGCTATACCCAGATCCTGACCACAGGAGCGATGTCTg  >  1:603285/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      AGGAGAAGTGAAATCACTGGTTTCGTCGCTACAAACCAAATATGCCCTGAAAGGCTATAACGTCGAAGCCACCGCACTGGAAGCTGCCGGTGGCAAAGTGGTCACAACGGGTAAAGATATTCTGGTTACCCTGCCGGCTTACCGGTTCACCAGTACGCCAGAAACC‑GATAACACCTGGCCGATTGAAGTCACC‑GCCGAAGATGTCAAAGGCAATTTGTCGAATCGTGAACAGAGCATGGTG‑GTCGTTCAGGCACCTACGCTAAGCCAGAAAGATTCCTCGGTATCGTTAAGTACCCAAACATTGAACGCGGATTCCCATTCAACCGCCACACTGACTTTTATTGCGCATGATGCAGCAGGTAATCCTGTTGTCGGGCTGGTGCTCTCGACGCGTCACGAAGGTGTTCAGGACATCACCCTTTCTGACTGGAAAGATAATGGTGACGGAAGCTATACCCAGATCCTGACCACAGGAGCGATGTCTG  >  NC_000913/2046266‑2046750

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