Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,002,222	G→A	5.5%	A462V (GCG→GTG)	fliC ←	flagellar filament structural protein (flagellin)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,002,222	0	G	A	5.5%	231.9 / 2.2	72	A462V (GCG→GTG)	fliC	flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (31/37);  new base A (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

TTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAG  >  NC_000913/2002076‑2002357                                                                                                                                                   |                                                                                                                                        ttGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCa                                                                                                                                       <  2:118446/149‑1 (MQ=255)      ggcAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGCGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTc                                                                                                                                  <  6:221786/147‑1 (MQ=255)         aGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGc                                                                                                                               <  6:231828/149‑1 (MQ=255)         aGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGc                                                                                                                               <  2:339324/149‑1 (MQ=255)         aGTCTCAGATAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGc                                                                                                                               >  8:299318/1‑149 (MQ=255)             tcAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTcc                                                                                                                           >  2:273523/1‑149 (MQ=255)             tcAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTAGCATAGTCGGCGTcc                                                                                                                           >  8:268539/1‑149 (MQ=255)               aGTTAATCCGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCt                                                                                                                          <  8:237890/148‑1 (MQ=255)               aGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCt                                                                                                                          <  6:363497/148‑1 (MQ=255)                    aTCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATa                                                                                                                     <  1:273523/148‑1 (MQ=255)                         gTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTc                                                                                                                                  <  2:179864/130‑1 (MQ=255)                         gTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTc                                                                                                                                  >  1:179864/1‑130 (MQ=255)                            aCAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTgc                                                                                                            >  1:266132/1‑149 (MQ=255)                             cAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTgcg                                                                                                           <  2:140986/149‑1 (MQ=255)                              aaCGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTgcgc                                                                                                          >  5:173770/1‑149 (MQ=255)                              aaCGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTgcgc                                                                                                          >  7:277262/1‑149 (MQ=255)                              aaCGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTgcgc                                                                                                          >  6:342983/1‑149 (MQ=255)                              aaCGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTgcgc                                                                                                          >  3:56456/1‑149 (MQ=255)                              aaCGATTAACCCTGCAGCAGAGACAGAAACTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTgcgc                                                                                                          >  6:60963/1‑149 (MQ=255)                                 gATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTc                                                                                                       <  3:116386/149‑1 (MQ=255)                                   ttAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAg                                                                                                     >  3:226793/1‑149 (MQ=255)                                      aCCCTGCAGCAGAGACAGAACGTGCTGCGGTACGTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTTAATACGGGACTGCGCTTCAGACa                                                                                                  <  5:375181/149‑1 (MQ=255)                                      aCCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGa                                                                                                    >  5:262720/1‑147 (MQ=255)                                            cagcagAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTgg                                                                                            <  7:299318/149‑1 (MQ=255)                                              gcagAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTgg                                                                                            >  4:300121/1‑147 (MQ=255)                                                    aCAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAgtggtg                                                                                    >  3:236558/1‑149 (MQ=255)                                                      aGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGtgtt                                                                                  <  8:277262/149‑1 (MQ=255)                                                        aaCCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTa                                                                                          <  1:84300/139‑1 (MQ=255)                                                        aaCCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTa                                                                                          >  2:84300/1‑139 (MQ=255)                                                           ctgctgCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTgttgtt                                                                               >  6:159782/1‑147 (MQ=255)                                                           ctgctgCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTgttgtt                                                                               >  7:149672/1‑147 (MQ=255)                                                           ctgctgCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCa                                                                             <  2:266132/149‑1 (MQ=255)                                                           ctgctgCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCa                                                                             <  5:60963/149‑1 (MQ=255)                                                            tgctgcGGTACCTGGTTAGCTCTTTCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTACGACATATTGGACACTTCGGTCGCATAGTCGGCGTACTGAATACTGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAg                                                                            >  6:173754/1‑149 (MQ=255)                                                              ctgcGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGt                                                                          <  8:149672/149‑1 (MQ=255)                                                              ctgcGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGATCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGt                                                                          <  6:308027/149‑1 (MQ=255)                                                                 cGGTACCTGGTTAGGTTTTTCCTACACGGTGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACTTATTGGTCTCTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCTGGtt                                                                         <  6:208445/147‑1 (MQ=255)                                                                 cGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTg                                                                        <  6:173770/148‑1 (MQ=255)                                                                    tACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCa                                                                                                      >  5:138103/1‑115 (MQ=255)                                                                    tACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCa                                                                                                      <  6:138103/115‑1 (MQ=255)                                                                              gCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGGGGTGTTGTTCAGGTTGGTAACCGCGGAATc                                                          <  5:342983/149‑1 (MQ=255)                                                                                  ttGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTGCAGACAGGTTTGTAGTGGGGTTGTTCAGGTTGGTAACCGCGGAATCCa                                                        <  5:173754/147‑1 (MQ=255)                                                                                   tGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGa                                                      <  3:300121/148‑1 (MQ=255)                                                                                     ccAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACgg                                                   <  5:159782/149‑1 (MQ=255)                                                                                      cAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACGTATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGAc                                                                                                   <  4:349259/100‑1 (MQ=255)                                                                                      cAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGt                                                  >  5:358008/1‑149 (MQ=255)                                                                                      cAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGt                                                  >  6:35829/1‑149 (MQ=255)                                                                                       aaCACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGtt                                                 <  7:268539/149‑1 (MQ=255)                                                                                       aaCACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCACATAGTCGGCGTCCTGAATACGGGACTGCGCTTc                                                                                                       >  7:331762/1‑95 (MQ=255)                                                                                       aaCACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCACATAGTCGGCGTCCTGAATACGGGACTGCGCTTc                                                                                                       <  8:331762/95‑1 (MQ=255)                                                                                               gTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTgcaccg                                         <  6:262720/149‑1 (MQ=255)                                                                                                     ggCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTgcaccgcacc                                     >  1:173692/1‑147 (MQ=255)                                                                                                      gCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAg                                  <  7:189115/149‑1 (MQ=255)                                                                                                        ctgctgGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGa                                                                                                    <  6:91427/81‑1 (MQ=255)                                                                                                        ctgctgGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGa                                                                                                    >  5:91427/1‑81 (MQ=255)                                                                                                                  atCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGTGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGa                       <  4:226793/148‑1 (MQ=255)                                                                                                                   tCTGCGCTTTAGACATATTGGTCACTTCGGTCGCATAGTCGGCGTCCTGAATACGGTACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCAACGCACCGAGGGAAGATCGGAAt                     >  6:275712/1‑149 (MQ=255)                                                                                                                     tGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAAttt                   >  7:326921/1‑149 (MQ=255)                                                                                                                      gcgcTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTg                  >  3:133081/1‑149 (MQ=255)                                                                                                                       cgcTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTg                  >  8:323400/1‑148 (MQ=255)                                                                                                                       cgcTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGt                 <  4:284528/149‑1 (MQ=255)                                                                                                                            tCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTa              <  6:358008/147‑1 (MQ=255)                                                                                                                               aCATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGgaagacc                          >  2:11774/1‑130 (MQ=255)                                                                                                                               aCATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAAc                          <  1:11774/132‑1 (MQ=255)                                                                                                                                   cttGGACAGTTCGGTCGCATAGTCGGCGTCCTGAATACGGGCCTGCGCTTCAGACAGGATGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACCGAATTTGTCTACAGATGCGa     <  4:56456/148‑1 (MQ=255)                                                                                                                                   atTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGa     >  4:268979/1‑149 (MQ=255)                                                                                                                                      ggACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATa   <  7:323400/148‑1 (MQ=255)                                                                                                                                       gACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACAGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAg  >  5:173408/1‑148 (MQ=255)                                                                                                                                          acTTCGGTCACATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGtt                                                 >  5:73361/1‑98 (MQ=255)                                                                                                                                          acTTCGGTCACATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGtt                                                 <  6:73361/98‑1 (MQ=255)                                                                                                                                            ttCGGTCGCATAGTCGGCGTGCTGAATACGGGACTGAGCTTCAGAGAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATa   <  8:137292/142‑1 (MQ=255)                                                                                                                                            ttCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATa   >  7:137292/1‑142 (MQ=255)                                                                                                                                                   |                                                                                                                                        TTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAACACGGAGTTACCGGCCTGCTGGATGATCTGCGCTTTCGACATATTGGACACTTCGGTCGCATAGTCGGCGTCCTGAATACGGGACTGCGCTTCAGACAGGTTGGTAGTGGTGTTGTTCAGGTTGGTAACCGCGGAATCCAGACGGTTTTGCACCGCACCGAGGGAAGAACGGAATTTGTCTACAGATGCGATAG  >  NC_000913/2002076‑2002357

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