Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,003,639	A→G	intergenic (‑33/‑233)	fliC ← / → fliD	flagellar filament structural protein/flagellar filament capping protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,003,639	0	A	G	100.0%	238.0 / NA	67	intergenic (‑33/‑233)	fliC/fliD	flagellar filament structural protein/flagellar filament capping protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/25);  total (42/25)

ACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTAT  >  NC_000913/2003510‑2003773                                                                                                                                  |                                                                                                                                       acagacGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGt                                                                                                                              <  2:137349/140‑1 (MQ=255)    gacGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTAtt                                                                                                                                  >  1:269193/1‑133 (MQ=255)    gacGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTAtt                                                                                                                                  <  2:269193/133‑1 (MQ=255)      cGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCa                                                                                                                         <  2:940328/140‑1 (MQ=255)        cTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCAcc                                                                                                                       >  1:424734/1‑140 (MQ=255)        cTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCAcc                                                                                                                       >  2:566114/1‑140 (MQ=255)        cTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCAcc                                                                                                                       >  1:1334437/1‑140 (MQ=255)          cGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCg                                                                                                                     <  2:362947/140‑1 (MQ=255)                aaCTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCt                                                                                                               <  1:1492225/140‑1 (MQ=255)                  cTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCa                                                                                                             >  1:1487542/1‑140 (MQ=255)                       cAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGc                                                                                                        >  2:1261951/1‑140 (MQ=255)                          cgcAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGt                                                                                                     >  1:987351/1‑140 (MQ=255)                            cAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCa                                                                                                   <  2:987351/140‑1 (MQ=255)                               aCTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAAcc                                                                                                <  2:834481/140‑1 (MQ=255)                                  ggTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTg                                                                                             <  2:1487542/140‑1 (MQ=255)                                    ttCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGtt                                                                                           >  2:627407/1‑140 (MQ=255)                                         gttgATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGCCTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGt                                                                                      >  1:1066847/1‑140 (MQ=255)                                         gttgATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGt                                                                                      >  2:600834/1‑140 (MQ=255)                                          ttgATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTc                                                                                     <  2:424734/140‑1 (MQ=255)                                             atatTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTc                                                                                     >  1:310461/1‑137 (MQ=255)                                             atatTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTc                                                                                     <  2:310461/137‑1 (MQ=255)                                              tattatTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTc                                                                                 <  1:648813/140‑1 (MQ=255)                                                   ttttGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTaaa                                                                            >  1:888409/1‑140 (MQ=255)                                                     ttGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAaac                                                                          <  1:1236306/140‑1 (MQ=255)                                                      tGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAaaca                                                                         <  1:566114/140‑1 (MQ=255)                                                        aGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAaacaac                                                                       >  2:207352/1‑140 (MQ=255)                                                                      ggCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAAtttttt                                                         >  2:680087/1‑140 (MQ=255)                                                                      ggCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAAtttttt                                                         >  1:1349626/1‑140 (MQ=255)                                                                       gCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTTTATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAAttttttt                                                        >  2:38877/1‑140 (MQ=255)                                                                       gCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAAttttttt                                                        >  2:1065573/1‑140 (MQ=255)                                                                       gCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAAttttttt                                                        >  2:1423431/1‑140 (MQ=255)                                                                               tATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACa                                                <  1:680087/140‑1 (MQ=255)                                                                                   aaTGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCa                                            >  2:1363533/1‑140 (MQ=255)                                                                                   aaTGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCa                                            >  2:124909/1‑140 (MQ=255)                                                                                     tGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCAtt                                          >  2:490509/1‑140 (MQ=255)                                                                                     tGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCAtt                                          >  1:281666/1‑140 (MQ=255)                                                                                         ttGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCAtttttt                                      >  2:555340/1‑140 (MQ=255)                                                                                         ttGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCAtttttt                                      >  2:711625/1‑140 (MQ=255)                                                                                            tgCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGtt                                   >  1:815936/1‑140 (MQ=255)                                                                                            tgCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGtt                                   >  1:334280/1‑140 (MQ=255)                                                                                                 tGATTCGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCg                              <  2:1349626/140‑1 (MQ=255)                                                                                                      cGTTATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGaa                         >  1:475688/1‑140 (MQ=255)                                                                                                        ttATCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAAt                       <  1:600834/140‑1 (MQ=255)                                                                                                           tCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATAct                    >  2:336520/1‑140 (MQ=255)                                                                                                           tCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATAct                    <  1:1065573/140‑1 (MQ=255)                                                                                                           tCCTATATTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATAct                    >  2:1240409/1‑140 (MQ=255)                                                                                                              tataTTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCtt                 <  1:1363533/140‑1 (MQ=255)                                                                                                                tatTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGaaa                        >  1:227501/1‑131 (MQ=255)                                                                                                                tatTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGaaa                        <  2:227501/131‑1 (MQ=255)                                                                                                                tatTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGaa                         >  1:1351750/1‑130 (MQ=255)                                                                                                                tatTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCtttt               >  1:1430015/1‑140 (MQ=255)                                                                                                                tatTGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTCTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGaa                         <  2:1351750/130‑1 (MQ=255)                                                                                                                   tGCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTtctc            >  1:1442561/1‑140 (MQ=255)                                                                                                                    gCAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTctct           <  2:295470/140‑1 (MQ=255)                                                                                                                     cAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTGAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAAt                       <  1:1083264/127‑1 (MQ=255)                                                                                                                     cAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAAt                       >  2:1083264/1‑127 (MQ=255)                                                                                                                     cAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTg          >  2:124822/1‑140 (MQ=255)                                                                                                                     cAAGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTg          >  2:1196004/1‑140 (MQ=255)                                                                                                                      aaGTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGc         >  2:1362498/1‑140 (MQ=255)                                                                                                                        gTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGccc       <  1:124909/140‑1 (MQ=255)                                                                                                                        gTCGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGccc       <  1:336520/140‑1 (MQ=255)                                                                                                                          cGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCt     >  1:220997/1‑140 (MQ=255)                                                                                                                          cGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCt     >  2:845860/1‑140 (MQ=255)                                                                                                                          cGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTAAAAATTTTTTTCAAAACCAGCCATTTTTTGTTGGTCCCCAAAATATTCTTTT‑TCTTCCCCtt    >  2:540069/1‑140 (MQ=255)                                                                                                                          cGTTGATTGCGTATTGGGTTTCCACCCGTCGGCTCAAACGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCt     >  2:1342986/1‑140 (MQ=255)                                                                                                                            ttGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTa   >  2:1086733/1‑140 (MQ=255)                                                                                                                             tGATTGCGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTAt  <  2:1430015/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                       ACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACCTTTAGAATTTTTTTCAAAAACAGCCATTTTTTGTTAGTCGCCGAAATACTCTTTTCTCTGCCCCTTAT  >  NC_000913/2003510‑2003773

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