Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,134,466	C→T	A222V (GCC→GTC)	flgF →	flagellar basal‑body rod protein FlgF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,134,466	0	C	T	97.7%	126.5 / NA	43	A222V (GCC→GTC)	flgF	flagellar basal‑body rod protein FlgF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/23);  minor base G (0/1);  total (19/24)
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

AATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTA  >  NC_000913/1134227‑1134709                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    aaactcGGCGCTCACTCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGAGAAAGCCACGGGAAGCGCAGTGCCGCGCGGTGACGACGGCATTTTTCGTTTAAGAGCAGAAACCCAGGCTACGAGAGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGATGCGGCAATGAGCGACCTGATTGTCAGCGCgcggc                                                                                                                                                                                                                                          >  2:37233/4‑251 (MQ=255)      cGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCgcggcg                                                                                                                                                                                                                                         <  1:76342/247‑1 (MQ=255)      cGGCGCGCAATCCGGTCGATCCGGCAAATACGGTTGCGCCATTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGttt                                                                                                                                                                                                                                      <  1:27694/250‑1 (MQ=255)       ggCGCTCAATCCGTGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTg                                                                                                                                                                                                                                    >  1:380165/1‑251 (MQ=255)         cgcTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGa                                                                                                                                                                                                                                   >  1:555134/1‑250 (MQ=255)                           cGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATc                                                                                                                                                                                                                <  2:92746/251‑1 (MQ=255)                                   aCGGTTGCGCCAGTAGGGCGTCTGAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCg                                                                                                                                                                                                         <  2:405744/250‑1 (MQ=255)                                         gcgcCAGTAGGGCGTCTTAACCTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGCCGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGAGTTTTGAAAGGCAGATGAAGGTGATTAGCAGCGTCgatgg                                                                                                                                                                                                  <  2:73854/251‑2 (MQ=255)                                                 aGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTGTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCa                                                                                                                                                                                           <  1:517174/250‑1 (MQ=255)                                                  gggCGTCTTAACCTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCCGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGGTGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAgg                                                                                                                                                                                         >  2:562894/1‑251 (MQ=255)                                                                 gTGAAAGCCACGGGCAGCGAAGTGCAGCTCGGTGACGACGGCATTGTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAAct                                                                                                                                                                          <  2:479305/251‑1 (MQ=255)                                                                   gAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAActgc                                                                                                                                                                        <  2:233949/251‑1 (MQ=255)                                                                          aCGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCTTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATTATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGAt                                                                                                                                                                 <  1:424333/251‑1 (MQ=255)                                                                             ggCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAg                                                                                                                                                              <  1:186401/251‑1 (MQ=255)                                                                              gCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCCACTGCTGTCGATGAGt                                                                                                                                                             >  2:63371/1‑251 (MQ=255)                                                                                  cGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAAt                                                                                                                                                         <  2:503290/251‑1 (MQ=255)                                                                                    aagggCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGTGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTg                                                                                                                                                       <  2:234709/247‑1 (MQ=255)                                                                                    aaGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTCCTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTc                                                                                                                                                                                                       >  2:234568/1‑203 (MQ=255)                                                                                    aaGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTg                                                                                                                                                       <  1:324860/251‑1 (MQ=255)                                                                                               ggTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATTAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACAt                                                                                                                                            <  1:554383/251‑1 (MQ=255)                                                                                                      ccgGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCGTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGAGCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAAGTGAAAGGATACATGACAAGt                                                                                                                                     <  2:286731/250‑1 (MQ=255)                                                                                                       cgGCATTTTTCGGTTAAGCGCAGAAACCCAGGACACGCGTGGGCCGGTACTGCAGGAAGATCCAAAATTGCGTGAGATGTCGGGGGTTCTGGAAGGCAGTAAAGTAAATGCCGTTGCGGCAATGAGAGACTTGATTGTCAGCGCGAGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGACCGTGCCACACAACTGATGTCGATAAGTTAATTGAAAGAGTACCTGACAAgc                                                                                                                                     >  2:498572/1‑249 (MQ=255)                                                                                                             tgttCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGt                                                                                                                               <  2:150709/248‑1 (MQ=255)                                                                                                                   attAAGCGCAGAAACCCCGGCCACGAGCGGGCCGGTACTACAGGCAGAACCCAACTTCCGTGTGCCGTCGGGGGTTCTGGCAGGCAGTAAAGTAAATGCCGTTGCGGCAATGCGCGACATAATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGCGAACAGCAGCGTCGATGCTAACGCAGGCCGTGCCAACCAACCGCTGTCGATGAGATAATTGACAGCATACCTAACACGAATAAGTAGCCCg                                                                                                                         >  2:33485/2‑250 (MQ=255)                                                                                                                                  ccAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCAAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCAAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCg                                                                                                         >  1:10500/1‑251 (MQ=255)                                                                                                                                        cACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCt                                                                                                    >  1:501858/1‑250 (MQ=255)                                                                                                                                          agcgTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATAAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTagg                                                                                                 >  2:357906/2‑249 (MQ=255)                                                                                                                                              tGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATggggg                                                                                             >  1:525714/1‑251 (MQ=255)                                                                                                                                                          gcaggccgATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAAttt                                                                                 >  2:260938/1‑251 (MQ=255)                                                                                                                                                             ggcagATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATc                                                                              >  1:118907/1‑251 (MQ=255)                                                                                                                                                                    ccAACCTTTCGTGTGATGTCGGGGGGTCTGGAATGCAGTAACGTCAATGCCGTGGCGGCAATGATCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTTCTGTCTATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGAGGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGAtttt                                                                        <  1:328119/250‑1 (MQ=255)                                                                                                                                                                     cAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGc                                                                      >  2:315857/1‑251 (MQ=255)                                                                                                                                                                                           ggtgcTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGAGTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTTTCGATTTTGCGAGCACTGGTAggtatgataag                                                <  2:395836/247‑10 (MQ=255)                                                                                                                                                                                                      gCAGTAACGTCAATGTCGTTGCGGCAATGAGAGACAGGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGAGGATAACGCAGGCCGTGCCAACCAACTGCTGGCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTGAGCGGGTCCATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcgcc                                     <  2:380165/251‑1 (MQ=255)                                                                                                                                                                                                           aaCGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                >  1:8958/1‑251 (MQ=255)                                                                                                                                                                                                            aCGTCAATGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                               <  2:8958/251‑1 (MQ=255)                                                                                                                                                                                                                 aaTGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                          <  1:17251/251‑1 (MQ=255)                                                                                                                                                                                                                   tGCCGTTGCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGCTGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAACGGATACATGACAAGTATAAGTTGCCCGATGCGCACGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGCTTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                           >  2:258899/1‑248 (MQ=255)                                                                                                                                                                                                                          gCGGCAATGAGCGACATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACata                 <  1:562894/251‑1 (MQ=255)                                                                                                                                                                                                                                       aCATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGAGCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc    <  1:182694/251‑1 (MQ=255)                                                                                                                                                                                                                                        cATGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACt   >  1:552113/1‑251 (MQ=255)                                                                                                                                                                                                                                         aggATTGGCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGGTATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACGGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTGGCCCGATGCGGAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTAGCGATTTTACGAGCACTTGTAGGCCGGATAAGGCGTTTACTCCGCATCCGTCAATAAGACATATGCACTTTGTCACTa  <  2:613042/249‑1 (MQ=255)                                                                                                                                                                                                                                         aTGATTGTCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGACAATCGCAATTTATCGATTTTGCGGGCACTTGTAGGCCGGAAAAGGCGTTTACGCCGCATCCGGCAAGAAGCAATATGCACTTTGTCACt   >  2:367377/1‑250 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    AATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGCGTCGATGATAACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTA  >  NC_000913/1134227‑1134709

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: ‑