Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,290,126	A→G	D438G (GAC→GGC)	yraJ →	putative fimbrial usher protein YraJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,290,126	0	A	G	100.0%	231.6 / NA	65	D438G (GAC→GGC)	yraJ	putative fimbrial usher protein YraJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/31);  total (34/31)

ATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAAT  >  NC_000913/3289993‑3290262                                                                                                                                      |                                                                                                                                         atatatGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCa                                                                                                                                    <  2:653918/140‑1 (MQ=255)      tGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCg                                                                                                                               >  1:551214/1‑140 (MQ=255)        gTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGAt                                                                                                                             >  2:689299/1‑140 (MQ=255)         tttGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGgca                                                                                                                                        >  1:452289/1‑128 (MQ=255)         tttGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGAt                                                                                                                             >  2:508638/1‑139 (MQ=255)         tttGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGAt                                                                                                                             <  1:508638/139‑1 (MQ=255)         tttGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATa                                                                                                                            >  2:250170/1‑140 (MQ=255)         tttGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATa                                                                                                                            >  1:670506/1‑140 (MQ=255)           tGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAAc                                                                                                                          <  2:7749/140‑1 (MQ=255)            gAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTAAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTTGCGCACTGGGCGCGTTGGCGATGGATATCAACAGAGCTTGCACCCAATTCGGTAAcc                                                                                                                         >  2:128648/1‑140 (MQ=255)             aataatACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCa                                                                                                                        >  2:691139/1‑140 (MQ=255)                        tACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGTTAACCAGCACTCTTTTc                                                                                                             >  1:649392/1‑140 (MQ=255)                           gCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATg                                                                                                          <  1:689299/140‑1 (MQ=255)                            cTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATgg                                                                                                         >  2:457265/1‑140 (MQ=255)                               tATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCa                                                                                                            <  2:4492/134‑1 (MQ=255)                               tATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCa                                                                                                            >  1:4492/1‑134 (MQ=255)                               tATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTa                                                                                                      >  1:55545/1‑140 (MQ=255)                                 tGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATc                                                                                                    >  1:9824/1‑140 (MQ=255)                                  ggcggcCTGCTCGGTTCCGAAGCTTATTATGGGCTCGGGTTCGGTATCTGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCa                                                                                                   <  2:196316/140‑1 (MQ=255)                                         tGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCCCACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCGGGCACCCAGTTCGATAGCCAGCACTATTTTGATGGCTATCAATGGCgg                                                                                            >  2:212552/1‑139 (MQ=255)                                         tGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGt                                                                                            >  2:661560/1‑140 (MQ=255)                                            tCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACg                                                                                         <  1:457265/140‑1 (MQ=255)                                            tCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACg                                                                                         <  2:670506/140‑1 (MQ=255)                                              ggTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCa                                                                                       <  1:456829/140‑1 (MQ=255)                                               gTTCTGAAGATTATTCTGCGCTTGGCATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCGCTCTTTTCATGGCTATCAATGGCGTACGCAg                                                                                      <  1:212552/140‑1 (MQ=255)                                                ttCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGt                                                                                     <  1:250170/140‑1 (MQ=255)                                                ttCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGt                                                                                     <  1:169204/140‑1 (MQ=255)                                                     aaGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTAc                                                                                   <  2:624557/137‑1 (MQ=255)                                                     aaGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTAc                                                                                   >  1:624557/1‑137 (MQ=255)                                                     aaGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATc                                                                                >  1:477102/1‑140 (MQ=255)                                                          tattatGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCCCTGGGCGCGTTGTCGATGGATATCAACAGAGCAGGCACCCAATTCGATAACCAGCACTCTTTGCATGGCTATCAAAGGCGCACGCAGTACATCAAAGa                                                                           >  2:302490/1‑140 (MQ=255)                                                          tattatGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTa                                                                                           >  1:164913/1‑124 (MQ=255)                                                          tattatGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTa                                                                                           <  2:164913/124‑1 (MQ=255)                                                          tattatGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGa                                                                           >  1:473558/1‑140 (MQ=255)                                                               tGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGTCGTACGCAGTACATCAAAGATATcc                                                                      >  1:281637/1‑140 (MQ=255)                                                               tGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATcc                                                                      >  1:102427/1‑140 (MQ=255)                                                               tGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCACGGCTATCGATGACGAACGCAGCACATCAAAGATATcc                                                                      >  1:301769/1‑140 (MQ=255)                                                               tGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCAAGGCTATCAATGACGGACGCAGGACATCAAAGATATcc                                                                      >  1:190247/1‑140 (MQ=255)                                                                 cgcTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCg                                                                    >  1:580838/1‑140 (MQ=255)                                                                  gcTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTCTCAATGGCGTACGCAGTCCATCAAAGATATCCCgg                                                                   >  2:302641/1‑140 (MQ=255)                                                                  gcTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCgg                                                                   >  1:555544/1‑140 (MQ=255)                                                                   cTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTGTTTTCATGACTAACAAAGGCGTACGTAGTACATCAAAGATATCCCGGa                                                                  >  2:252493/1‑140 (MQ=255)                                                                    tGGGGATCGGTATCGGCGGCACACTTGGAGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGaa                                                                 >  1:461559/1‑140 (MQ=255)                                                                     ggggATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGaaa                                                                <  2:281637/140‑1 (MQ=255)                                                                       ggATCGGTATCGGCGGCTCACTCGGCGCACTTGGCTCGTACTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAAcc                                                              <  2:208763/140‑1 (MQ=255)                                                                         aTCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCaa                                                            <  2:651126/140‑1 (MQ=255)                                                                           cGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAaca                                                          >  2:11016/1‑140 (MQ=255)                                                                            ggTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAacac                                                         <  2:102427/140‑1 (MQ=255)                                                                            ggTATCGGCGGCACACTAGGCGGACTGGGCGGGTTATCGGTGGATATCTACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAacac                                                         <  2:301769/140‑1 (MQ=255)                                                                                tCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAat                                                     >  1:135897/1‑140 (MQ=255)                                                                                   gcggcACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCATTGGCGTACGCAGAACACCAAAGATAACCCGTAAACCAACACCAATATc                                                  >  2:225745/1‑140 (MQ=255)                                                                                    cggcACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCg                                                 <  2:477102/140‑1 (MQ=255)                                                                                          acTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCa                                           >  2:514369/1‑140 (MQ=255)                                                                                                cgcACTGGGCGCGTTATCGATGGATCTCATCAGAGCTAGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGctact                                     <  1:302490/140‑1 (MQ=255)                                                                                                    cTGGGCGCGTTGTCGCTGGATAGCACCAGCGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCg                                 <  1:225745/140‑1 (MQ=255)                                                                                                    cTGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCg                                 <  1:691139/140‑1 (MQ=255)                                                                                                     tGGGCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGc                                <  2:473558/140‑1 (MQ=255)                                                                                                      gggCGCGTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCt                               >  2:183782/1‑140 (MQ=255)                                                                                                         cgcgTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCtata                            <  2:555544/140‑1 (MQ=255)                                                                                                           cgTTGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATAcc                          <  1:183782/140‑1 (MQ=255)                                                                                                             ttGTCGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCaa                        >  1:32846/1‑140 (MQ=255)                                                                                                                 cGATGGATATCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGAt                    >  2:709013/1‑140 (MQ=255)                                                                                                                         atCAACAGAGCTGGCACCCAATTCGTTAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAttt            <  1:487524/140‑1 (MQ=255)                                                                                                                         atCAACAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTAttt            <  2:461559/140‑1 (MQ=255)                                                                                                                              cAGAGCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGAGATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGtt       <  2:43263/140‑1 (MQ=255)                                                                                                                                 agCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTa    <  1:514369/140‑1 (MQ=255)                                                                                                                                 agCTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTa    >  2:500662/1‑140 (MQ=255)                                                                                                                                   cTGGCACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAAt  <  1:128648/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         ATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGCACTCTTTTCATGGCTATCAATGGCGTACGCAGTACATCAAAGATATCCCGGAAACCAACACCAATATCGCTGTCAGCTACTATCGCTATACCAACGATGGCTATTTTAGTTTTAAT  >  NC_000913/3289993‑3290262

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