Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,977,719	T→G	16.9%	K42N (AAA→AAC)	flhC ←	flagellar class II regulon transcriptional activator, with FlhD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,977,719	0	T	G	16.9%	110.7 / 2.9	58	K42N (AAA→AAC)	flhC	flagellar class II regulon transcriptional activator, with FlhD
Reads supporting (aligned to +/- strand):  ref base T (21/25);  new base G (8/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 7.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGC  >  NC_000913/1977491‑1977950                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             ggTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACTTTTTTTTCCCAGGGCCAAAACCAGTCGGTTGAGAATGCCAGCCTTCCTTTCGGCGGGGGGGGT‑CCGCCCCGT‑TTGTTATAAAAATTTGTCAAGc                                                                                                                                                                                                                                    >  2:57905/1‑239 (MQ=255)    tttctccacgatttgTCCTTCTTCTGCTTTTGGGAACTTTTCAAGGGATAACCGGTAGGCTTTGACCACCGCATCGACGCCATTCCACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTAGCAGGCGTc                                                                                                                                                                                                                                 <  2:472606/225‑1 (MQ=255)           cgccgggggccTTTTTCTGCTTTTGGGCCCTTTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATGGACGCCATTACAAAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑t                                                                                                                                                                                                                         <  2:498626/230‑1 (MQ=255)            gcagTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCCGT‑TCTTTTTAAAGTTTTTTCAGGCGTCCCCGCC‑tt                                                                                                                                                                                                                        >  1:606088/1‑239 (MQ=255)            gcagTGGTCACGCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGCCATAAACCAGTCGGTTGAGAATGGACGCATGACTTTCGGCGGGGGGGTTCCCGCCCATT‑CTTTTATAAATtttttacagatcgcccgactt                                                                                                                                                                                                                          >  2:21900/1‑222 (MQ=255)              ctgtgtcctttcttttcTTTTGGGCACTGTTCAAGGTATAAACGGGAGGCTTTGATCACCGGATCGACGCCATTACACAAACCGGTTTTCAGTAAAACCTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTG‑CCGCGCAGT‑TCGTTATAAAGTTTTATCAGGGGTCCGCGAC‑TTaa                                                                                                                                                                                                                      <  2:113340/223‑1 (MQ=255)                  gTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGATTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCAGTCCTTTCGGCGGTGGGATT‑CCGCGCAGT‑TCCTTATAAAATTTTATAAGGGGTCAGCGGC‑TTAACtg‑t                                                                                                                                                                                                                 >  2:811892/1‑239 (MQ=255)                        tttctGTTTGTTGGCCCTTTTCAAGGTTTAAACGGTAGGCTTTGATCCCCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCg                                                                                                                                                                                                           <  2:324617/238‑1 (MQ=255)                        cttctGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGATTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGAATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGACTGCACACTTGCCTTTCGGCGGTGGGGCT‑GCGCGCAGTGTGTTAATAAGGTTTT‑TCACGGGGCCGGGAC‑TAAACTGATGtttg                                                                                                                                                                                                            >  2:435020/1‑238 (MQ=255)                                           ttCAAGGAATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTTTAAAGTTTTATCAGGGGTTCGCGCC‑CTAACAG‑TGTTTTGCTTTTCATC‑ATCTTCAAAc                                                                                                                                                                                       >  1:841260/1‑239 (MQ=255)                                            tCAAGGTATAAACGGTAGGCTTTGATCCCCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACg                                                                                                                                                                                      <  1:853790/239‑1 (MQ=255)                                            tCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACg                                                                                                                                                                                      <  2:733239/239‑1 (MQ=255)                                                   ataactggTAAGCTTTGATCCCCACATAGCCGCCTTTACACAAACCGGTTTTCAGTAAAAACTGCCAGGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGTGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAgcgc                                                                                                                                                                                 <  1:525813/231‑1 (MQ=255)                                                     aaaCGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAg                                                                                                                                                                             <  1:577995/239‑1 (MQ=255)                                                     aaaCGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCCCGCAGG‑TCTTTTTAAAAGTTTTTCAGGCGGCCGCGACTTTACCTG‑TGTTT‑GCTTTCCAAC‑acattgaaaaggggaccccg                                                                                                                                                                             >  2:889756/1‑220 (MQ=255)                                                      aaacgtgggCTTTTATCCCCGCATCGACGCCATTACCAAACCCGGTTTTCAGGAAAAACTGCAATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAgg                                                                                                                                                                            <  2:600378/232‑1 (MQ=255)                                                          gTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGGTT‑CCCCGCCGT‑TCTTTTTAAAGCTTTTTCAGTGCTCCCCGGC‑TTAAATG‑CGTTTCGCTTTCCCAC‑ATTTTCAAACAAGCCCCCAgtggac                                                                                                                                                                        >  2:577995/1‑234 (MQ=255)                                                           tAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGGCT‑TCTCGCCCT‑TCTTTTTAAAGTTTTTTCAAGCGGCCCCGGC‑TTAACCC‑TGGTTCGCCTTTCAGC‑CTTTTCCAACGAACCCCCAAGGGGag                                                                                                                                                                       >  2:1529/1‑238 (MQ=255)                                                            aGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCCGGCGTCCGCGAC‑TTCACTG‑TGTTTTGCTTTCCAGC‑ATCTGCAAACGGGCGCCCAGGGGGATc                                                                                                                                                                      >  1:859882/1‑239 (MQ=255)                                                            aGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCCTTATAAAGGTTTATCAGGCGTCCGCGGC‑TTAACTG‑TGTTTCGCTTTCCAAC‑ATCTGCAAACCGACGCCCAGGGGGATc                                                                                                                                                                      >  1:514273/1‑239 (MQ=255)                                                              gCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCCT‑CCGCGCAGT‑TCTTTATAAAATTTTATCAGCCGTCCGCGAC‑TTAACTG‑TGTTTGGCTTTCCAGC‑ATCTGCCAAGGAGAGCCCCggggggac                                                                                                                                                                      >  1:698764/1‑232 (MQ=255)                                                                       cccggagcGAAGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGc                                                                                                                                                          <  1:719286/232‑1 (MQ=255)                                                                        ccGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCCGGCGTCCGCGAC‑TTTACTG‑TGTTTCGCTTTCCTGC‑ATCTGCAAAATGACGCCCCGGGGGATCA‑ATTCCCTTGc                                                                                                                                                          >  1:280904/1‑238 (MQ=255)                                                                         cGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCCCGAC‑TTAACTG‑TGTTTCGCTTTCCCGC‑CTCTTCAAACGAGGGCCCAGGGGGATCA‑ATTCCATTTcc                                                                                                                                                         >  1:420342/1‑238 (MQ=255)                                                                         cGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGCCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTAGACTT‑CCCTGCCGT‑TTCTTATAAAGGTTTACATGGCGGCCACGGC‑TTAACAG‑GGTTTTGCTTTCCACCAATTTCCAAACAAAGGCCCAggggataaaattcaattccc                                                                                                                                                          >  2:331831/1‑222 (MQ=255)                                                                                              aCCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACAAGCGCCCCGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCCCGCT‑TCCt                                                                                                                                  >  1:877074/1‑239 (MQ=255)                                                                                                             accctgccctgccTTTCAGAACATCGAAGCATTAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCtttttt                                                                                                                   <  2:420342/236‑1 (MQ=255)                                                                                                               aCTGCCCTGCATTACCGAACATCGAAACATTAACGTTTTGTTCCCAGGTCACAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCa                                                                                                                 <  1:811892/239‑1 (MQ=255)                                                                                                                    gaTGCTTTTCAGAACATTGAAGCATGAACTTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCa                                                                                                             <  2:514273/237‑1 (MQ=255)                                                                                                                           tACAGAACATCGAAGCATGGACCTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAgg                                                                                                     <  1:889756/239‑1 (MQ=255)                                                                                                                           tACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCcg                                                                                                                                                                                                                               <  2:506819/122‑1 (MQ=255)                                                                                                                           tACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGGTTTAACAGGCGTCcg                                                                                                                                                                                                                               >  1:506819/1‑122 (MQ=255)                                                                                                                              cgAACATTCAAGCATTGACGTTTTGTTCCCAGGTCATAAACCCGGCGGTTGAAATAGGCAGCTTGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGccc                                                                                                  <  2:859882/238‑1 (MQ=255)                                                                                                                                  cATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTTGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGGGATCA‑ATTCCATTGCCCGCTGAATATCCCCCCCT‑TCCTGGACAATGCTTTTTTCACTCATGACAAGGCCCtttt                                                                                              >  1:441411/1‑239 (MQ=255)                                                                                                                                  cATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCCGGGGGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTTAATCATGATCAGGCCCtttt                                                                                              >  1:576068/1‑239 (MQ=255)                                                                                                                                     cGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCCCGCAGT‑TCTTTATAAAGTTTTATCTGGCGGCCGCGGC‑TTTCCTG‑GGTTTCCCTTTTCCGC‑ATATTCAACAGCGAGACCCGGGTGTCCC‑ATTCCATTTCCAGCGGCAAATCCCGCCCT‑TCCTGAACAATCCTTTTTTATCTCCTGTTAAGACCCCtttttt                                                                                           >  2:913504/1‑239 (MQ=255)                                                                                                                                     cGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGATT‑CCGCGCAAT‑TCTTTATAAAGTTTTATTCGGGGTCCGCGGC‑CTAACTG‑TGGTTCTGTTTCCAAG‑ACCTGCAAACCAGCGCCCAGGGGGGTCA‑ATTCCACTTCCCGCTGAAAATCCCCCCCT‑TCCCGGACCAAACTTTTTTCAATCATTATCCGGGCCtttttt                                                                                            >  2:510693/1‑236 (MQ=255)                                                                                                                                     cGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGGGGGCTT‑CCGCGCCGG‑TCTTTATAAAGGTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGGTTTCCTTTCCACC‑ATCTGCAAACGAGCGCCCAGGGGGATCA‑ATTCCATTGCCAACTGAATATCCCCCGCG‑TCCTGCAAAATGCTTTTTTTACTCACTACAAGGCCCTtttttt                                                                                           >  2:714784/1‑239 (MQ=255)                                                                                                                                     cGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGCCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCCCCCACT‑TCTTTTTAAAGTTTTATTCGGTGTCCTCGAC‑CTAACTG‑TGTTTCGCTTTCTACC‑ATTTTCAAACCACGCCCCCGGGTCATCC‑ATTTCCTTTCCCGCCTACTATCCCCCCCTGTCCCG‑ACAATGCTTTTTTCCATCCTTACCAAGCCCCTTTCtt                                                                                           >  2:718272/1‑239 (MQ=255)                                                                                                                                           aTGAACTTTTGGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGGTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTgcgc                                                                                       <  1:510693/237‑1 (MQ=255)                                                                                                                                                       ttCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAgg                                                                         <  1:714784/239‑1 (MQ=255)                                                                                                                                                             ggTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGGGATCA‑ATTCCATTGCCAGCTGAATTACCCGCGCT‑TCCCGAACAATCGTTTTTTCATTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGAtt                                                                   >  1:360529/1‑239 (MQ=255)                                                                                                                                                                 aTAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCa                                                                                                                 <  2:466045/189‑1 (MQ=255)                                                                                                                                                                 aTAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCCTTTTTCa                                                                                                                 >  1:466045/1‑189 (MQ=255)                                                                                                                                                                         gTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCCTGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTTCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTTCCAGCTGAATATCCCGCGCT‑TCCTTAACCATGCTTTTTTCACTCCTGATCAGGCCCTTTTCTTGTGCCGCGCTTCTTTAGGCTGGTTAACATCATTCAg                                                       >  2:734875/1‑239 (MQ=255)                                                                                                                                                                         gTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGGTTTATCAGGCGTCCGCGAC‑TTAACTG‑GGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCCGCGCTTCTTCAGGCTGGTTAACATCATTCAg                                                       >  1:259510/1‑239 (MQ=255)                                                                                                                                                                               tGAGAATGGCAGCATGCCTTTCGGCGGTGGGCAT‑CTCCGCAGT‑TCTTTATAAAGGTTTATTAGGCGTCCGGGGC‑CTAACTG‑GGGTTCGCTTTCCAGG‑CTCTGCAAACGAGGGCCCAGGGTGAGCACATTCC‑TTGACAGCTGACTTTCCCCCGAT‑TCCCAACCAATTCTTTTTTCAACCATGTTACAGCCCTTTTCTTGGACAGGCCTTCTTTAGGATGATTAAAATCATTCacgaaacc                                                 >  2:843323/1‑232 (MQ=255)                                                                                                                                                                                    atgcccgcagtccTTTCGGCGGTGGGCTG‑TCGCGCAGG‑TCTTTATAAAGTTTTATCAGGCCTCCGCGCC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCGGCAAACGACCGTCCAGGCTGATCA‑ATTCCATTGCCAGCTGAATAACCCGCTCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTTCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGtt                                             <  1:649088/227‑1 (MQ=255)                                                                                                                                                                                        cagcaTGCCTTTCGGCGGTGGGCTT‑CCGCGCCGT‑TCTTTTTAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGCATATCCCGCGCC‑TCCTGAACAATGCTTTTTTCACTCATTATCAGGCCCTTTTTTTGCGCAGCGCCTCTTCATGTTGATTCACATAATTCAGCAAGCGTGTTgaga                                         >  2:301467/1‑238 (MQ=255)                                                                                                                                                                                        cagcaTGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTTTAAAGGTTTATCAGGCGTCCGCGAC‑TTTACTG‑TGTTTCGCTTTCCAGC‑ATCTGCCAACGAGCGCCCAGGGGGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTTGACAATGCTTTTTTTCCTCCTGTTCCAGCCCTTTTCTTTCGCACCCCTTTTTCAGGCTGTTTAACCACATTCCGCaacatggtttaggg                                        >  2:210186/1‑227 (MQ=255)                                                                                                                                                                                                ttttgcgcggTGTGCTT‑CCGCGCCGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCTAC‑TTAACTG‑TGTTTCGCTTTACAGC‑ATCTGCAAACGAGCTCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTTAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCatgatg                                 <  1:205054/232‑1 (MQ=255)                                                                                                                                                                                                cTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAgcg                                                                                   <  1:235313/188‑1 (MQ=255)                                                                                                                                                                                                 tttCGGCGGTTGGCTT‑CCGCGCAGT‑TCTTTATAAAGGTTTATCAGGCGGCCGCTGC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGGCCAGGGTGATCA‑ACTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGc                                <  1:818555/238‑1 (MQ=255)                                                                                                                                                                                                         gTGGGCTT‑CCGCGCATT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGaa                       <  2:606088/239‑1 (MQ=255)                                                                                                                                                                                                           gggCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAAt                      <  1:1529/238‑1 (MQ=255)                                                                                                                                                                                                              cTT‑CCCCTCTGT‑TCTTTATAAATTTTTTTCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGc                  <  2:841260/239‑1 (MQ=255)                                                                                                                                                                                                               tt‑CCGCGCGTT‑TTTTTATAAATTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGc                  <  2:576068/238‑1 (MQ=255)                                                                                                                                                                                                                                tataACGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGACTATCCCGCCCT‑TCCTGAACCATGCTTTTTTCACTCCTGATCAGGCCCTTTTCTTGCGCAGGGCTTCCTCAGGCTTATTCACATTCTTCAGCAAGCGAGTTGAGCGCATGCTCACGGTACGAATTAGCTGGTAGTCGCCAAcgg  >  2:701535/1‑238 (MQ=255)                                                                                                                                                                                                                                tataAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAAACTCATTCACCAAGCGGGTTGAGAGACTGATGCGGGTAGGAATTTGCTGGATAAGGTAAAcgc  >  2:95404/1‑239 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             GGTTAATGCCAGCAGTGGTCCTTCTTCTGCTTGTGGGCACTGTTCAAGGTATAAACGGTAGGCTTTGATCACCGCATCGACGCCATTACACAAACCGGTTTTCAGTAAAAACTGCCATGCATTACAGAACATCGAAGCATGAACGTTTTGTTCCCAGGTCATAAACCAGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTT‑CCGCGCAGT‑TCTTTATAAAGTTTTATCAGGCGTCCGCGAC‑TTAACTG‑TGTTTCGCTTTCCAGC‑ATCTGCAAACGAGCGCCCAGGGTGATCA‑ATTCCATTGCCAGCTGAATATCCCGCGCT‑TCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGC  >  NC_000913/1977491‑1977950

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