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	100.0%	153.6 / NA	43	V45A (GTA→GCA)	yadV	putative fimbrial chaperone YadV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/21);  total (22/21)

TTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATA  >  NC_000913/155823‑156307                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                tttCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGtt                                                                                                                                                                                                                                             <  1:180767/250‑1 (MQ=255)     tgcatctagttttccagcagcTTCCAGAACGATAAACCCGAACACGCTCTCTCTGTCTATACGCAGTGACGTGCTGGCTGTGTACATCAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGCCGGCGTAGCAGTAACAGGGACTGTAATCCTCCCAGGCTCAGCGTTGTCATCGCCCGTATCTAACCAAATCTGGACAAACAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGacac                                                                                                                                                                                                                                        >  2:244739/21‑251 (MQ=255)            ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACtttt                                                                                                                                                                                                                                    <  2:197355/248‑1 (MQ=255)               tttGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGAt                                                                                                                                                                                                                               <  1:527659/250‑1 (MQ=255)                     ggAACTTCCAGGACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGGAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACCGGCGGCGTAGCAGTAAAAGGGACTGTAATCCTGCCAGGCTCAGCTTTGTCATCGCCAGTATCTCACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGaacgattt                                                                                                                                                                                                                         <  2:321721/250‑8 (MQ=255)                       aaCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTtat                                                                                                                                                                                                                      >  2:414518/1‑251 (MQ=255)                           tCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTtatata                                                                                                                                                                                                                   <  1:385720/250‑1 (MQ=255)                            ccAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATAtt                                                                                                                                                                                                                 >  2:133576/1‑251 (MQ=255)                                 aCGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTAcgcg                                                                                                                                                                                                            >  2:349009/1‑251 (MQ=255)                                                            tCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCaa                                                                                                                                                                                  <  2:478717/250‑1 (MQ=255)                                                              tttAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATa                                                                                                                                                                                >  2:300669/1‑250 (MQ=255)                                                                  ggCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACg                                                                                                                                                                            >  2:324393/1‑250 (MQ=255)                                                                   gCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGAt                                                                                                                                                                          <  2:196803/251‑1 (MQ=255)                                                                            tGCTGGCTGTGTACATTACTTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCCGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTaaa                                                                                                                                                                  >  2:367205/1‑250 (MQ=255)                                                                             gctggctgTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTaaa                                                                                                                                                                  >  1:471417/1‑249 (MQ=255)                                                                             gctggctgTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAg                                                                                                                                                                <  2:575167/251‑1 (MQ=255)                                                                                   tgtgTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATAc                                                                                                                                                          >  1:246198/1‑251 (MQ=255)                                                                                        aCATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGt                                                                                                                                                      <  2:263162/250‑1 (MQ=255)                                                                                            tAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATAATGGCGGAGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCCGAGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCACAGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATCTATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACg                                                                                                                                                  >  1:548135/1‑250 (MQ=255)                                                                                                  gATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCAt                                                                                                                                            >  2:487905/1‑250 (MQ=255)                                                                                                         tGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCtgt                                                                                                                                    >  1:228246/1‑251 (MQ=255)                                                                                                          gCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCtgtt                                                                                                                                   <  2:123719/251‑1 (MQ=255)                                                                                                              aCGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTAt                                                                                                                               <  1:321062/251‑1 (MQ=255)                                                                                                               cGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATg                                                                                                                              <  2:574118/251‑1 (MQ=255)                                                                                                                    ggCATCAATAAGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACTAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTgg                                                                                                                         >  1:500421/1‑251 (MQ=255)                                                                                                                           aTACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTaaa                                                                                                                  <  1:578490/251‑1 (MQ=255)                                                                                                                               gcgATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAAc                                                                                                              <  2:449854/251‑1 (MQ=255)                                                                                                                                    aCTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATaa                                                                                                          <  2:77061/250‑1 (MQ=255)                                                                                                                                          ggcgTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATcc                                                                                                   >  1:525333/1‑251 (MQ=255)                                                                                                                                                     aaaaGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGtt                                                                                        <  1:280624/251‑1 (MQ=255)                                                                                                                                                        aGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTTCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGttaat                                                                                     >  2:559080/1‑251 (MQ=255)                                                                                                                                                               gTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTaaaa                                                                                <  1:319489/249‑1 (MQ=255)                                                                                                                                                                  aTACTGCCACGCTCAGCGCTGTCATCGCCAGTATCTAACCACCTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACTTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCTCAATAGACGATGAACTAACAGCCATACAGGTTACGAAGCATGAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAggg                                                                           >  2:371563/1‑251 (MQ=255)                                                                                                                                                                       gCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaata                                                                       <  1:154086/250‑1 (MQ=255)                                                                                                                                                                        ccAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTTTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaataat                                                                     <  1:549066/251‑1 (MQ=255)                                                                                                                                                                        ccAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGaataat                                                                     >  1:312960/1‑251 (MQ=255)                                                                                                                                                                                         aTCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAattatt                                                    >  2:331215/1‑251 (MQ=255)                                                                                                                                                                                             ccAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTc                                                >  2:582041/1‑251 (MQ=255)                                                                                                                                                                                                 tatatAACCAACTCTGGACAAGCAACGGGTTATTACCTTTATTTTCCAGACAAGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGATAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTAGGTTTGGTGTTAAAAACCATAAACCATCCCTGATAGCGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTg                                            >  2:135510/5‑251 (MQ=255)                                                                                                                                                                                                       aCCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATatt                                      >  1:19988/1‑251 (MQ=255)                                                                                                                                                                                                                      gCAACGGGTTATTCCCTTTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTa                       <  2:150198/251‑1 (MQ=255)                                                                                                                                                                                                                          cGGGTTATTCCCTGTATTTTCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTg                   <  2:111019/251‑1 (MQ=255)                                                                                                                                                                                                                                             tCCAGACGTGCGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATa  >  2:313651/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                TTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGACGTACGTTGACACTTTTTTGATCGCTTTTATATATTACGCGAGTACCCGAAATGACAATGTCCGCAATAGACGATGAACTAAAAGCCATACAGGTTACGAAGCATAAAGCTGTTGTATGTTTGGTGTTAAAAAACATAAACCATCCCTGATAGAGTTAATTAAAACAGGGAATAATATAAAATAGAAATTATTATTCCCTGTATATTCATTCAATCAATTTAACTGGTGAATACTTACTGATA  >  NC_000913/155823‑156307

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