Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,138,831	C→T	R152C (CGC→TGC)	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,831	0	C	T	97.9%	145.9 / ‑6.4	47	R152C (CGC→TGC)	flgK	flagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (22/24);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTA  >  NC_000913/1138590‑1139075                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      taCCAACCAGTTACGGGCGGAGCAGACGCAAAGTAGCGGTCTGACGGCCCGCTATGAGCAGTTGTCGAAAATCGACAATATGCTCTCCACCAGTACCCCTTCCCTGGCAACACTGATGCAGGATTTCGTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGCAAGCAGAAGGATTGGTGAAGCAGTTTAAAACCACCGACCAATATCTGTGCGACCAg                                                                                                                                                                                                                                               <  2:285766/250‑1 (MQ=255)  accaaccaGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAgg                                                                                                                                                                                                                                              <  2:15859/250‑1 (MQ=255)            tACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCCCCGATCAATATCTGTGCGACCAGGACAAACAGGTc                                                                                                                                                                                                                                   >  2:266698/1‑251 (MQ=255)             aCGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGCGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGCCCAGGACAAACAGGTc                                                                                                                                                                                                                                   >  2:18370/1‑250 (MQ=255)               gTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAat                                                                                                                                                                                                                                >  2:258084/1‑251 (MQ=255)                       cAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGAt                                                                                                                                                                                                                        <  1:455621/251‑1 (MQ=255)                             aaaaGTAGAGGTCTGACTGCCCGCTATGAGCTGATGGCGAAACGCGACAATATGCTGGCCGCCAGTGCCTGTGCGCGGGGAACAGAGGTGCGGGATTGCTGCACCAGCCTGGAAACGCTGGTGAGTAAGGGGGAAGACCCGGCAGCGGGCCTGGCGCGGAGTGGGAAAGCAGAAGGATTGGTGAATCAGTTTAAAACCACCTATCAATAGCTTTGGGACCAGGACAGACAGGTCAATATTGAGATAGGTGc                                                                                                                                                                                                                  <  2:426654/250‑1 (MQ=255)                                     cGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTg                                                                                                                                                                                                          <  1:317611/251‑1 (MQ=255)                                      ggTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAg                                                                                                                                                                                                                                               >  2:319879/1‑213 (MQ=255)                                      ggTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAg                                                                                                                                                                                                                                               <  1:319879/213‑1 (MQ=255)                                           gACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTgatcaga                                                                                                                                                                                                    <  1:229252/251‑1 (MQ=255)                                           gACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTgatca                                                                                                                                                                                                      <  1:18370/249‑1 (MQ=255)                                               gCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTgatcagatca                                                                                                                                                                                                 <  1:4801/250‑1 (MQ=255)                                               gCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCCCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGGGCCAGCGTTgatcagatca                                                                                                                                                                                                 <  1:230361/250‑1 (MQ=255)                                                 ccGCTATGAGCAGATGTCGAAAAACGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGGCCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATcaaca                                                                                                                                                                                              <  2:178991/251‑1 (MQ=255)                                                        gAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACCCAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGCAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAAAAGGTCAATATCGCGATAGGTGCCAACGTTGATAAGATCAACAACTACGc                                                                                                                                                                                       >  2:224871/1‑251 (MQ=255)                                                               tGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACa                                                                                                                                                                                >  1:276835/1‑250 (MQ=255)                                                                gTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACACACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACaa                                                                                                                                                                               >  1:170515/1‑250 (MQ=255)                                                                                           caGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTcgc                                                                                                                                                   >  2:314070/1‑251 (MQ=255)                                                                                                              acacAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGACCGATCAAATCTCGCGCCTGACAGGCGTGGGGGc                                                                                                                                >  2:262/1‑251 (MQ=255)                                                                                                                   gATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAgggg                                                                                                                           <  1:266751/251‑1 (MQ=255)                                                                                                                          gATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCAcc                                                                                                                    <  1:262/251‑1 (MQ=255)                                                                                                                              tcttcACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAAc                                                                                                                >  1:306765/1‑251 (MQ=255)                                                                                                                              tattaaCCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGATGATTGGGAAATCAGCAGGATTGGAGCATAAGTTTAAAACCACAGATCAATCTCTGTGCGACCAGGACAAACAGGACAATATCGCGATAGGTGCCAGCATTGATCAGATAAACAACTACGCACAAACAA‑TAGCAAGCCTGACCGCTCAAATCTCGCGCCTACCAGGAGTGGCGGCAGGGGCGTCACCTAAc                                                                                                                >  1:306769/6‑251 (MQ=255)                                                                                                                                          tGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATcgcg                                                                                                                                                                                                                          <  2:364459/134‑1 (MQ=255)                                                                                                                                          tGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATcgcg                                                                                                                                                                                                                          >  1:364459/1‑134 (MQ=255)                                                                                                                                           gCAAACGCTGGTGAGTAACGCGTAAGACCCGGCAGCGCGACAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGAGCGTCACCTAACAATCTGCTGGATc                                                                                                   <  2:435280/251‑1 (MQ=255)                                                                                                                                                                ggAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCg                                                                              <  1:147807/251‑1 (MQ=255)                                                                                                                                                                   aGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAATGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAAt                                                                           <  1:152535/251‑1 (MQ=255)                                                                                                                                                                        cGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAAc                                                                      <  1:266698/251‑1 (MQ=255)                                                                                                                                                                              cgcgcCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAAAtt                                                                >  2:347620/1‑251 (MQ=255)                                                                                                                                                                                 gcCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCCCCTAACAATCTGCTGGATCAACGCGATCATCTGGTGAGCGAATTAAACCAGAttgtt                                                             >  2:254581/1‑251 (MQ=255)                                                                                                                                                                                    aGGCGCTGAGTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCATCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGTCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGgt                                                          <  2:119921/251‑1 (MQ=255)                                                                                                                                                                                       cgcTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCCGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGgtgt                                                        >  2:137620/1‑250 (MQ=255)                                                                                                                                                                                       cgcTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGgtgt                                                        >  2:189258/1‑250 (MQ=255)                                                                                                                                                                                                          aGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATg                                    <  1:82286/251‑1 (MQ=255)                                                                                                                                                                                                           ggATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATgg                                   >  1:427032/1‑251 (MQ=255)                                                                                                                                                                                                           ggATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACAAGATTGTTGGAGTAGAAGTCAGCGTTCAGGATgg                                   >  1:354364/1‑251 (MQ=255)                                                                                                                                                                                                            gATTGGGGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATggc                                  <  2:276835/251‑1 (MQ=255)                                                                                                                                                                                                                ggTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCa                              <  1:209831/251‑1 (MQ=255)                                                                                                                                                                                                                 gTGAATCAGTTTAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGCCAGGCGTGGGGGCAGGGGCGTCACCTCACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGAAGGCGGCAc                             >  2:233730/1‑251 (MQ=255)                                                                                                                                                                                                                  tgaatacgTTTAAACCCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATCGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAAAAAATAGCCAGCCTGAACGATCAAATCTCGCGCCTGACACGCGAGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGCTTGTTGGTGTAGACGTCAGCGTTCAGGATGGAGGCACt                            >  2:76626/8‑251 (MQ=255)                                                                                                                                                                                                                          tttAAAACCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACAt                    >  1:454150/1‑251 (MQ=255)                                                                                                                                                                                                                             aaaaCCACCGATCAATATCTGTGCGACCAGGACAAACAGGACAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCAc                 >  1:412355/1‑251 (MQ=255)                                                                                                                                                                                                                              aaaCCACCGATCAATATCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTACCAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACg                >  2:3482/1‑251 (MQ=255)                                                                                                                                                                                                                                     cGATCAATAACGGTGCGACCAGGACAAACAGTTCAATCTCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCGCGCGCCTGACATTCGTGGGGGCAGGTGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTAGAACATCACGATGGCCa         <  1:103797/251‑1 (MQ=255)                                                                                                                                                                                                                                             atCTGTGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAAGCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTa  <  2:427032/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      TACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGCCAGGCGCTGATTGGGAAATCAGAAGGATTGGTGAATCAGTTTAAAACCACCGATCAATATCTGCGCGACCAGGACAAACAGGTCAATATCGCGATAGGTGCCAGCGTTGATCAGATCAACAACTACGC‑TAAACAAATTGCCAGCCTGAACGATCAAATCTCGCGCCTGACAGGCGTGGGGGCAGGGGCGTCACCTAACAATCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTA  >  NC_000913/1138590‑1139075

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