Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,138,533	C→T	G52G (GGC→GGT)	flgK →	flagellar hook‑filament junction protein 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,138,533	0	C	T	97.5%	102.8 / NA	40	G52G (GGC→GGT)	flgK	flagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/17);  minor base G (1/0);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTG  >  NC_000913/1138290‑1138776                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     tGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTTGGTTg                                                                                                                                                                                                                                               >  2:178823/1‑251 (MQ=255)                tGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGGTGGGTTGGCACTGGTGGCTACg                                                                                                                                                                                                                                >  2:76381/1‑251 (MQ=255)                tGTTCTGAATAACTCAAGTCCGGCGGGTCGCTACCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCCGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCttcg                                                                                                                                                                                                                                >  1:234342/1‑248 (MQ=255)                    cTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTc                                                                                                                                                                                                                            >  1:492815/1‑251 (MQ=255)                       aaTAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGGGTCTACGTTTCTgg                                                                                                                                                                                                                         >  1:24519/1‑251 (MQ=255)                         tAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGACGGCTTATAAGGAACCTCCATGTCCAGCTTGATTACTAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCCGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTTGGAATTGTGTCTCCGTTTCTGgtt                                                                                                                                                                                                                       >  1:261770/1‑250 (MQ=255)                                 gTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCgg                                                                                                                                                                                                               >  1:437220/1‑250 (MQ=255)                                    cggcggGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTCAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAg                                                                                                                                                                                                            <  1:109208/251‑1 (MQ=255)                                         ggTCGCTGCCGATAATACTCGTTATTTAAAGTGTTATAAGGAACCTCCATGCCCAGCTTGATTAATAACGCCTTGAGCGGACTGATCGCGGCCAAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACTTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTatga                                                                                                                                                                                                       <  2:174433/251‑1 (MQ=255)                                            cGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCCGGCGGCGTTAACTACGGCAAGTAATAATATCTCCTGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTTGCAATGGGGTCTACGTTTCTGGTGTGccgggggaggatgaagg                                                                                                                                                                                                    >  2:234016/1‑235 (MQ=255)                                               tGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCCATGGTGTCTACGTTTCTGGGGTGCCGCGTGAGTATGATGCGtt                                                                                                                                                                                                 >  2:191392/1‑251 (MQ=255)                                                   gATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGGGTGCAGCGTGAGTATGATGCGTttatt                                                                                                                                                                                             >  1:227790/1‑251 (MQ=255)                                                         aCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGATGGACGCTGGCGGTTGGGGTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGGGGGGATGGTGAGTTTCTTaccaac                                                                                                                                                                                       >  2:531502/1‑251 (MQ=255)                                                              gTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTaccaac‑ca                                                                                                                                                                                    <  1:156919/248‑1 (MQ=255)                                                                         cTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGGGCAGCGTGAGTATGATGCGTTTATTTACAAC‑CAGTTACGTGCGgcgc                                                                                                                                                                      >  1:91650/1‑251 (MQ=255)                                                                                                    gATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGcc                                                                                                                                            <  1:30108/250‑1 (MQ=255)                                                                                                    gATTAATAACGCCATGAGCGGACGGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGccc                                                                                                                                           <  1:204270/251‑1 (MQ=255)                                                                                                          taaCGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTCTAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATGGGGTCTACGTTTCTGGTGTGCAGCGGGAGTAAGGTGCGTTTATTACGAAC‑CATTTTCGTGCGGCCCAGACCCAAGAAAGCGGTCTGCCCGACCGCTAtt                                                                                                                                     >  2:369073/1‑250 (MQ=255)                                                                                                           atgacatgcaGCGGACTGACGGCGCCACAGTCGTCGTTAATTACGGCATGTAATAATCACTCCCGCAACAACGTGGCCGGATATCCCCGCCAAACCACTATTATGGCGCAGTCCAAAAGCACGTTGGGCGCTGGCGGTTGGGTTGGCATTGGTGTCTACGTTTCTTGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGAGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATgt                                                                                                                                    <  2:17439/242‑2 (MQ=255)                                                                                                               ccATGAGCGGACTGAACGCGACCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTTGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGACTGCCCGcctttagc                                                                                                                                  >  2:471356/1‑242 (MQ=255)                                                                                                                       ggACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGGATGATGCGGTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGa                                                                                                                         >  2:406705/1‑250 (MQ=255)                                                                                                                                      cAGGCGGCGTTAAACATGGCAAGTAATACTATCTCCAGCTATAACGTTCCCGGATATACCCGCCAAACCACTATTATGGCGAAGGCACATAGCACGTTGGGCGCTGGCGGTTTGGTTCGCCATGGAGTCTCAGTTTCTGGTGTGCCGCGCGAGGATGATTCGATTATT‑CCAACACAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTATGACAGACCGCTATAAGCAGAAGGCGAAAATCGACAATATGct                                                                                                         >  2:603709/1‑251 (MQ=255)                                                                                                                                          cggcgTTAAATACGGCAAGTAATATTATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTTGGTTGGCAATGGTGTCTACCTTTCTTGTGTGCAGCGTGAGTATTATGCGTTTATTACCAAC‑CCGTTTCGTTCGGCGCAGACGCAAAGTAGCGGTTTTACTTCCCTCTTTGTGCTGATGTCGGATATCTACAATATGttcttc                                                                                                     >  2:340908/1‑249 (MQ=255)                                                                                                                                             cgTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCCGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCCGAAGCAAAGTAGCGGTCTGCCTGCCCGCTATGAGCAGATGTCGAAAATCGACCATATGCTTTccacc                                                                                                  >  2:541007/1‑251 (MQ=255)                                                                                                                                               atAAATACGGCAATTAATAACATCGCCAGCTAAAACATGGCCGGAAATACCCGCCAAACATCTATGATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCCACGGTTCTGGGGTGCAGCGTGAGCATGATGCGTTTATTACCAAG‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGACGGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCGGCACCAg                                                                                                <  2:261770/250‑1 (MQ=255)                                                                                                                                                tAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAg                                                                                                <  1:234016/250‑1 (MQ=255)                                                                                                                                                     aCGGCAAGTAATAATATCTCCCGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAAAGGTGTCTACGTTTCTGGGGTGCAGCGTGGGTATGATGCGTTTATTACCAAC‑CCGTTACGTGCGGCGCCGGCGCAAAGTAGCGGTAGGACTTCCCGCCATGAGCAGATGTCGAAAAACGAAAATATTCACTCCACCAGTAcccc                                                                                          >  2:529759/1‑249 (MQ=255)                                                                                                                                                                   tatCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGGCGAAAATCGACAATATGCTCTCCACCAATACCTCTTCGCTGGCAacac                                                                            >  1:351686/1‑251 (MQ=255)                                                                                                                                                                                  cGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATttct                                                             <  1:346428/251‑1 (MQ=255)                                                                                                                                                                                                 cccccAAACCACTATTATGGCGCAGGCCAATAGCCCGTGGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGGTACGTGCGGCGCAGACGCGAAGTAGCTGTCTGACTGCCCGCTATGAGCAGATGTCGAAAAGCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCaaa                                              <  2:10397/248‑1 (MQ=255)                                                                                                                                                                                                 ccGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCaaa                                              <  1:620685/251‑1 (MQ=255)                                                                                                                                                                                                  cGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAAc                                             <  1:132428/251‑1 (MQ=255)                                                                                                                                                                                                      aaaCCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTg                                         <  1:261577/251‑1 (MQ=255)                                                                                                                                                                                                          cACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCGCCACCAGTACCTCGTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGa                                     <  1:541007/251‑1 (MQ=255)                                                                                                                                                                                                            cTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCCATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTAACAAC‑CAGTTACGTGTGGCGCAGACGCCAAGGAGCGGGCTGCCTGCCCGCCATGAGCAGATGTCGAAAATCGACCATATGCTCTCCACCCGTACCTCTTCGCTGGCCACACAGATGCAGGATTTCTTCACCAGCCTGCCAAAGCTGGTgggt                                   >  2:187686/1‑248 (MQ=255)                                                                                                                                                                                                                      gcAGGCCAATAGCGCGTTGGGCGCGGGCGGTGGGGTTGGCTATGGGGTCGACTTTTCTGGGGTGCTGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGGGGCTCAGACGCAAAGTAGCGTTCTGACTGCCCGCTAGTAGCAGATGTCGAAAATCGACAATATTGTCTCCGCCAGGCCCTCTGTTCTGGCAACACAGATGGATGATTTCTTGGCCAGCCTGCAGTCGCTGCTGAGGAACGCGGAAg                         <  2:558034/251‑1 (MQ=255)                                                                                                                                                                                                                        aGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGAc                       <  1:178823/251‑1 (MQ=255)                                                                                                                                                                                                                                  cACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAgcgc             <  1:549566/251‑1 (MQ=255)                                                                                                                                                                                                                                     gTTGGGCGCTGGCGGTTGGGTTGGCAATGGGGTCTACGTTTCTGGTGGGCAGCGTGAGTATGATGCGTTTATTAACAAC‑CAGTTTCGTGCGGCGCCGACGCAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTTTTCACCCGCCTGCAAACGCTGGGGAGGAACGCGGAAGACCCGGCAGCGCGcc          >  1:186522/1‑251 (MQ=255)                                                                                                                                                                                                                                         agcgcTGGCGGTTGGGATGGACATGGTGTCTACGCTTCTTGTGTGGAGCGGGAGGAAGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGGCTGAATGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTTCACCAGTACCTCTTCGCTTGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGAACCGGCAGCGCGCCAggg      >  2:105428/2‑250 (MQ=255)                                                                                                                                                                                                                                               ggcggtTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTg  <  2:227790/249‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     TGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAAC‑CAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTG  >  NC_000913/1138290‑1138776

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