Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	156,068	A→G	V45A (GTA→GCA)	yadV ←	putative fimbrial chaperone YadV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	156,068	0	A	G	97.5%	108.3 / NA	40	V45A (GTA→GCA)	yadV	putative fimbrial chaperone YadV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/19);  minor base T (1/0);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATAAGTA  >  NC_000913/155835‑156311                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    gTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTg                                                                                                                                                                                                                                     >  1:71278/1‑250 (MQ=255)   tttGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACTTGCGTTGACACTTTTTTGAt                                                                                                                                                                                                                                   <  1:584218/250‑1 (MQ=255)     tggtggAACTTCCAGTACGTTAATCCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCg                                                                                                                                                                                                                                 <  1:574150/250‑1 (MQ=255)            aCTTCCAGTACGTTAAACCAGATCACGCTCTTTCTGTCTTGAGTCAGTGATGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGGTACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCGCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGGTATTCCCTTTAGTTTCCAGACGTGCGTTGACACTTTTTTGATCGTTTTtata                                                                                                                                                                                                                         <  2:338838/251‑1 (MQ=255)                  aGTACGTTAAACCATAACCCGCTCTCTCTGTCTTTAGGCAGGTAGGTGCTGGCTGGGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGTCTGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGTTTAGCGTTGTAATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCATACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTAc                                                                                                                                                                                                                   <  2:295878/251‑1 (MQ=255)                   gTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTAcg                                                                                                                                                                                                                  <  2:71278/251‑1 (MQ=255)                     aCGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGCGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGTCATCAATACGCGATACTGGCGGAGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGATGTCATCGCCAGTATCTAAACAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCATACGTTCGCTGACACTTTTTTGATCACTTTTCTATATTAcgcg                                                                                                                                                                                                                >  1:27427/1‑251 (MQ=255)                         tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGCTTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTTTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATTTATTACGCGAGt                                                                                                                                                                                                             >  2:544375/1‑250 (MQ=255)                                      gctcttagtttcgttAGGCAGTGTGGAGCTGGGTGTGGTCATTAATTTGATTTTTTGCCCGCGTTTGGAATAAGTACGCGATACTTGAGGTTTAGCATTAAAAGGGACTGTCATACGGCCCGGCATAGCTTAGTCATCGCCAGTATCTAACCAACGCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGTCGTGCGTTGACACTGTGTGTATCGCTTTTATATATTACGCGAGTACCTGACATGACaa                                                                                                                                                                                               <  2:159915/238‑1 (MQ=255)                                             ctGTCTTTAGGAAGTGAGGTGCTGGATGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGTCGGCGTAGCAGTAAAAGGGACTGTAAAACTGCCAGGCTAAGCGGTGGCTTCGCCCGTATCTAACCAACTCTGGACAAGCAACTGGTTATACACTTTATTTTCCAGACGTGCGTTGCAACGTTTTTGCTCACTTTTATATATTACGCTAGTAACCGattagacactgtccg                                                                                                                                                                                         >  2:191730/1‑236 (MQ=255)                                                    tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGCTTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTTTATATTACGCGAGTACCCGAAATGCAAATGTCCGCAATAGAc                                                                                                                                                                                 >  2:611623/1‑251 (MQ=255)                                                       gCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGAt                                                                                                                                                                              <  1:594481/251‑1 (MQ=255)                                                            gaggagCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATCCTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCCATAGACGATGa                                                                                                                                                                            >  2:383289/1‑248 (MQ=255)                                                            gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGa                                                                                                                                                                            <  1:383289/248‑1 (MQ=255)                                                                tGCTGGCTGTGTACATTAATTTGCTTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGAGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGAGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTaaa                                                                                                                                                                      >  2:612270/1‑250 (MQ=255)                                                                    ggctgTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGc                                                                                                                                                                   <  1:612270/249‑1 (MQ=255)                                                                                                    gTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGGAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATAAAGGTTACGAAGCATAAAGATGTTGTATGt                                                                                                                                 >  1:298287/1‑251 (MQ=255)                                                                                                             cAATACGCGATACTGGCGGCGTAGCCGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGtt                                                                                                                         >  1:204187/1‑250 (MQ=255)                                                                                                             cAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGtt                                                                                                                         >  1:286640/1‑250 (MQ=255)                                                                                                                 aCGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTaaaaa                                                                                                                    <  1:521546/251‑1 (MQ=255)                                                                                                                     gATACTAGCGGAGTAGCAGTAAAAGGGACTGTACTACTGCCAGGCTCAGCGTTGTCATCGCCAGTATATAACCACCTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGCCGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACACTGTCCGACATAGACGATGAACTAAAAGCCATACAGGTTACGCAGCATAAAGCTGTTGTATGTTTGGTTTTAAAAAACAt                                                                                                                >  2:429603/1‑251 (MQ=255)                                                                                                                              ggcgTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATcc                                                                                                       <  1:203668/251‑1 (MQ=255)                                                                                                                                           aaGGGACTGTAATACTTCCAGGCGCAGCGTTGTAATCGCCAGTGTCTAACCAACTCTGGACAAGCAACGTGTTATTCCCTTTAGTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTGTATATAATACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCTGAAAGCTGTTGTATGTGTGGTGTTAAAAAACATAAACCATCCCTGAGAGAGttaa                                                                                          <  2:63410/251‑1 (MQ=255)                                                                                                                                               gACTGTAATGCTGACAGGCTTAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGGTATTCCCTTTAGTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGttaattaa                                                                                      <  2:570367/251‑1 (MQ=255)                                                                                                                                                aCTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTaaa                                                                                     <  1:550472/251‑1 (MQ=255)                                                                                                                                                 cTGTAATACTGCCAGGCTAAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGCCAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCAAGCACCCGAAATGACACTGTCCGCAATAGACGAAGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTAGTGTTAAAACACATAAACCATCCCTGATAGAGATAATTaaa                                                                                     >  2:124222/1‑250 (MQ=255)                                                                                                                                                          tGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaat                                                                            >  1:496813/1‑250 (MQ=255)                                                                                                                                                          tGCCAGGCTCAGCGTTGCCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTCTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAACGCTGTTGTATGTTTGGTGTTAAAAAACATACACCATCCCTGATAGAGTTAATTAAAACAGGGaat                                                                            >  2:219458/1‑250 (MQ=255)                                                                                                                                                            ccAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaataat                                                                         >  1:299990/1‑251 (MQ=255)                                                                                                                                                                 cTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATaa                                                                     >  1:414298/1‑250 (MQ=255)                                                                                                                                                                          tgcatcttcAGTAACTAACCAACTCTGGACAAGCAACGGGTTATTCCCTGTATTTTCCAGACGTGCGTTGACTCTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACTATTAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAGACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAat                                                            <  1:672536/242‑1 (MQ=255)                                                                                                                                                                            cATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACCCTTTTTTGATCGCTTTTATATATTACGCGAGTCCCCGAAATGACAATGTCCGCAATAGACGATGAACTAACAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAAAAGTGAATAATATAAAATAAAAgttat                                                         >  2:555920/1‑251 (MQ=255)                                                                                                                                                                                    gTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCt                                                 <  1:611623/251‑1 (MQ=255)                                                                                                                                                                                             cAACTCTGGACACGCAACGGGATATGCCCTTTATTTTCCAGACGAGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCAGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAACACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGATTAATATAAAATAGAAATTATTATTCCCTGTATattc                                         >  2:242717/1‑250 (MQ=255)                                                                                                                                                                                                ctctGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAg                                                                                                                                                            <  1:575576/132‑1 (MQ=255)                                                                                                                                                                                                ctctGGACAAGCAACGGGTTATTCCCTTTATTTTCAAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCACTAGACGATGAACTAAAAGCCATACAg                                                                                                                                                            >  2:575576/1‑132 (MQ=255)                                                                                                                                                                                                     gACAAGCAACGGGTTCATCCCTTTATTTTCCAGACGTGCGTTCACCCTTTTATGATCGCTTTTATATATTACGCGAGTACCAGAAATGCCAATGTCCGCACTAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTCTGTTTGGCGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATACAATAGAAATTATTATTCCCTGTATATTCATtcaatc                                 >  2:561211/1‑250 (MQ=255)                                                                                                                                                                                                      aCAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATtcaatcaa                               <  1:544375/251‑1 (MQ=255)                                                                                                                                                                                                                ggTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGt                     <  2:204187/251‑1 (MQ=255)                                                                                                                                                                                                                                    aGACGTGCGGTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTCCCCGAAATGACAATGTCCGCAATAGCCGCTGAACTAAAAGCCATAGAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTACAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATAAGTa  >  1:38132/1‑250 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    GTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATAAGTA  >  NC_000913/155835‑156311

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