Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,065,172 G→A 11.6% G198G (GGC→GGT)  yihO ← putative sulphoquinovose importer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,065,1720GA11.6% 157.4 / 2.0 75G198G (GGC→GGT) yihOputative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (7/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATC  >  NC_000913/4064954‑4065394
                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
ttaaataatatcacgcaaTACGGGGTGTAAAGAACCTAAACCACAAATTTGACAATAAAGGCCGCCAGGGTAAACAAGTTAACAATGCATAAAACCAGCAATGGCGGATTCCGAAATGCCGTGAAAAATGATTTGAGAATCCTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATCCTCTTTGAAGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATg                                                                                                                                                                                                                       <  1:248376/211‑1 (MQ=255)
ttAATATCATCCAACACATACTGGGTGTAAAAGACCGTAATCGCCAGTTTGATATAAAATGCCGCCAGGTTACAAAAATTAGCAATGCAAAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATg                                                                                                                                                                                                                       <  2:701239/228‑1 (MQ=255)
         ttCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCTCGCAAAATGATTTGAGAATACTGGCTTTATGTCAGCTTGGCTTTGT‑TACACATTAATGCCTTTTGACG‑CTTTGTAGG‑CCAGCAATCTTAAAAA‑AAAACAGCCAATTGAGGACaaaa                                                                                                                                                                                                             >  1:26838/1‑228 (MQ=255)
         ttCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCGTTAGC‑ACACC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATa                                                                                                                                                                                                              >  2:652447/1‑228 (MQ=255)
             cacattccTGGGTTTAAAAGACCTGAACCGCAATTTTAATATAAAATGCCGCCAGGGTAAACAAGTTAGCAAGGCATAAAACAAGAAATGGCGGATGCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTg                                                                                                                                                                                                          <  1:24096/221‑1 (MQ=255)
              cacaTACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGGTTCCAAA‑TAATGCTCTTTGACGCCTCTTGAGC‑CCGGC‑ATCTTAAAAATAAAACCACCAATGGCGAACATAAATGc                                                                                                                                                                                                         >  1:522166/1‑228 (MQ=255)
                      gggTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAgg                                                                                                                                                                                                 >  2:485269/1‑228 (MQ=255)
                          taaaaggcccTGAAACGCCAGGTTGATAGTAAATGCCGACAGGGTAAACAAGTTAGCAATGCAAAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTGGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATa                                                                                                                                                                                             <  2:193512/220‑1 (MQ=255)
                           tAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTTAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATACTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCCTAAAAATAAAGCCGCCAATGGAGAACATAAGGACCGCACAGGCCTaa                                                                                                                                                                                            >  2:11677/1‑228 (MQ=255)
                                   cTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTAGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCCCTTTGACGCTCCTGTAGA‑ACAGC‑ATCATCAAAATAAAACCGCCAATGCGTACCATAAGTACCACACAGCCAAAACCAAcaga                                                                                                                                                                                    >  1:314696/1‑225 (MQ=255)
                                      aaTCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATAACGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGCCGTTGTGACCACA‑TAATGCTCTTTGACGCCTCTTTAGC‑CAAGC‑ATCAATAAAATAAACCCGCCAATGGAGAACATAAATGCCTCACAgaaacaaacaaaagatag                                                                                                                                                                                 >  1:311860/1‑211 (MQ=255)
                                               tttGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGga                                                                                                                                                                                                                     >  2:56205/1‑183 (MQ=255)
                                                     ttaaaaTGCCGCGAGGTGAACCATTTTAGCAACGAAAAAACCAAGCAAGGGATGTTCCCGAAATGTCGCGAAAATTGATTTGATCATACTGGTTTGAGGCCCGATTGGCGTTGTACCCACA‑TAATGCTCTGTGACGCCGCTGTAGC‑ACGGA‑ATGAAAAAAATAAAGCCGCCAATGGAGAACAGAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGa                                                                                                                                                                  <  1:11677/225‑1 (MQ=255)
                                                     atTAAATGCCGCCAGGGTTAACAAGTTAGAAAGGCATAAAACCAGCAATGGCGGATCCCGAAATATCGCACAAAATGATTTGAGAATACTGGCTTTATGTGCGGTTGGGGTTGTGACCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ACCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGa                                                                                                                                                                  <  2:331562/228‑1 (MQ=255)
                                                                ccAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGGCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTTACGCCTCTGTAGG‑CCAGC‑ATCCTAACACTAAAACCGCCAATGGGGAACATAAGTGCCGCACAGCCATAACAGAAGGTTGCACAAGAAAAAAGCGAcgccaaacaaa                                                                                                                                                       >  1:331562/1‑218 (MQ=255)
                                                                ccAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACCTGCTTTATGCCCGGTTGGCGTTGTGACCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAACCGCCAATGGAGAACATAAGGGCCGAACAGCCATAACCGACGCTTACATCGGAAAAAAGCAAcggaaaacgaa                                                                                                                                                       >  1:445207/1‑218 (MQ=255)
                                                                    ggtaacaaaGTTAGAAATTCATAAAACCATCAATGGCGGTTTCCGAAATGTCGGGCAAAATGATTTGAGAATACTGGTTTTATGTCCGGTTGGCGTTGGGTCCACA‑TAATGCTCTTTGACGCCTCTGAAGC‑ACAGC‑ATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATaaa                                                                                                                                                   <  1:652447/221‑1 (MQ=255)
                                                                         acaATTTAGAAATGCATAAAACCAGCATTGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCa                                                                                                                                              <  2:522166/228‑1 (MQ=255)
                                                                         acaAGTTAGAAATGCTTTAAACCAGCAAGGGCTGTTTCCGAAGTTTCGAGGAAAATGATTTGAGAATCCGGGCCTTATGTCCGGTTGGAGTTGTGACCAAA‑AAATGCCCTTTGGCGCCTCTGTAGC‑ACAGC‑ATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCa                                                                                                                                              <  1:217449/228‑1 (MQ=255)
                                                                          aataagagtaaggCATAAAACCGACAATGGCGAGTTCCGAAATGTCGCGAAAAATGATTGTAAAAGACAGGCTTTAGGTCCGGTTGGCTTTGTGTCAGCA‑TAATGCTTTTTGACGCCCCTGAAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCa                                                                                                                                              <  2:163547/215‑1 (MQ=255)
                                                                              ttAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCCCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGAATGGAAAGGAGTAAAGGCCACggg                                                                                                                                          >  2:430861/1‑226 (MQ=255)
                                                                                aGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTAc                                                                                                                                        >  2:546369/1‑227 (MQ=255)
                                                                                     tGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGACCATAAGTGCCGCACAGGACTAACCGACGGTTCAGTCAGAAAAAAGCGACTGCAAAGGAAAAAAACCCACGGGAAAAAtaa                                                                                                                                  >  1:605594/1‑226 (MQ=255)
                                                                                     tGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATAATGGCTTTATGTCCAGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGTAGAACATAAGTGCCGCACAGGCAAAACCAACGACGCAGTAGGAAAAAAGCAATTGCACCGAAAAAAACACCACGGAAAAAACCa                                                                                                                                  >  1:302770/1‑228 (MQ=255)
                                                                                     tGCATAAAACCAGAAATGGCGGATTCCGAAATATCGCGCAAAATGATTTAAGAATACTGGCTTTATGGCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAAGCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGTAAAAAACGAACGCAACGGAATAAACGCACCGGTAAAAAACa                                                                                                                                  >  2:660740/1‑228 (MQ=255)
                                                                                      gCATAAAACCAGAAAGGGGGGATTCCGAAATAGCGCGAAAAATGATTTGAGAATACTGGTTTTATGCCCGGTTGGAGTTGTGCCCACA‑TAATGCTCTTTGACGCCTCTGGAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCaa                                                                                                                                 <  2:314696/228‑1 (MQ=255)
                                                                                       cATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑CCAGC‑ATCATCAAACTAAAGCCGCCAATGGAGACCATACGGACCGGACAGGCATAACCGACGGATGAGTCCGGAAAAAGCGACTGCAACGCAATAAACACCCCGGGAAAAACCAAt                                                                                                                                >  1:38197/1‑228 (MQ=255)
                                                                                       cATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAAt                                                                                                                                >  2:222383/1‑228 (MQ=255)
                                                                                                 gCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTTGCTTTATGTCCGGTTGGCGTTGTGTCCAAA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAACCGCCAATGGAGAACCTAAGTGCCCGACAGGCCTAAACGACGCCTGAGTCAGAAAAAAGAGACTGAAACGGAATAAACCCCACGGTACAAAACacgtatcttataa                                                                                                                      >  1:163547/1‑216 (MQ=255)
                                                                                                 gCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAACAATAAAGCCGCCATTGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAACGACTGCAACGAAATAAACACCAAGGTGCAAATCCATAAcccgaagg                                                                                                                      >  1:549320/1‑223 (MQ=255)
                                                                                                            ttCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCaccacc                                                                                                            >  2:421968/1‑227 (MQ=255)
                                                                                                                       cgcgCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTCGC‑ACCGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAACAGCGACTGCAACGGAATAAACGCCACGGTACAAAACAATACCCCTATGGTAGCACCCCCCTGGCgtgagt                                                                                                 >  2:695573/1‑222 (MQ=255)
                                                                                                                       cgcgCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAATCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGCCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAAAAGCCCAATGGTTACACACACCTGCAGGTAAg                                                                                                 >  1:443082/1‑227 (MQ=255)
                                                                                                                                  aTTTGAAAATGCTGGCTTTATGTCCGGTTGGGGTTGTGTCCACA‑AAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTgcg                                                                                     <  2:490876/228‑1 (MQ=255)
                                                                                                                                    ttGAAAATTCTGGTTTTATCTGCGTGTGGCGTTGTACACAAA‑GTATCCTATTAGGCCCCTATGAAGCAACAGC‑AT‑AAAAAAAAAAAGCTGCCAATGGAGAAAATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCCCCACCCTGGCGGTAAGCCGCGAGCTgcgcg                                                                                   <  1:490468/228‑1 (MQ=255)
                                                                                                                                                        cccGGTTGGCGTGGTGTACACA‑TAATGCCCTTTGACGCCTCTGTAGC‑CAAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACAAAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGt                                                               <  2:605594/227‑1 (MQ=255)
                                                                                                                                                          cGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGCCGGAAAGCCGCAAGCTCCCGGGGCCCAATCGGGGTTTTGGTaa                                                             >  1:109074/1‑228 (MQ=255)
                                                                                                                                                          cGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACAAACCTGGCGGTAAGCCGCGACCTGCCGGCCCTCATTCGGGTTTTTGGTaa                                                             >  2:598562/1‑228 (MQ=255)
                                                                                                                                                          cGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTCGAGTCAGAAAAAAGCGCCTGCACCGGAATAAAAGCCCCGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGCAAGCCGCGAGCTGACCGCACTCCTTCGGGGTTTTGGGaa                                                             >  1:419815/1‑228 (MQ=255)
                                                                                                                                                          cGGGTGGCGTTGGCACAAAA‑TAATGCTTGTTGGAGCCTCTGGAGC‑ACAGC‑AACAAAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCAAAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTaa                                                             <  1:485269/228‑1 (MQ=255)
                                                                                                                                                            gTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAAt                                                            >  2:347518/1‑227 (MQ=255)
                                                                                                                                                                          caca‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCAAAGGCATAACCAACGGTTGAGACAGACAAAAGCGACTGCAACGGATTAAAAGCCACAGACAAAATCAAAAACCCTATGGGTGAACAACCCAGGAGGAAAGCGCCAACCAGCACCCCCCCATCCGGGGTTTGGGTAATTGCCAGGACCATAGc                                             >  1:45562/1‑228 (MQ=255)
                                                                                                                                                                          aaaaataatgcatggagACGCCCTTGTACC‑ACAGA‑AACAAAAAAATAAGGCCGCCAATGGAAAAAATAAGGGCCGAACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAAGCAATAGCCCTATGGTTGCTCCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCg                                              <  1:678992/212‑1 (MQ=255)
                                                                                                                                                                           aca‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCTCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGatcat                                                <  2:45562/224‑1 (MQ=255)
                                                                                                                                                                                   gCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACg                                                                                                                                            >  2:284876/1‑125 (MQ=255)
                                                                                                                                                                                     ttttggcggCCGCTGAAGC‑AAAGC‑AAAAAAAAAAAAAAGCGGCCAAGGGAAAACAAAAGTGCCGGACAGGCAAAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAg                                   <  1:347518/220‑1 (MQ=255)
                                                                                                                                                                                     tctTGGACGCCTCTGTAGA‑ACAGC‑ATAATAAAAATAAAGCCTCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCGCCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAg                                   <  2:311860/228‑1 (MQ=255)
                                                                                                                                                                                      ctTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGAAATTGCCGGGATCATCGCACCATACGac                                   >  2:114976/1‑226 (MQ=255)
                                                                                                                                                                                       tttGGCGCCTTGTTAGC‑ACATA‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCa                                 <  2:445207/228‑1 (MQ=255)
                                                                                                                                                                                       tttGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGGCCGCGCTCCTTCGGGCTTTTGGTAATTGCCGGGATCATCGCACCATACGAGc                                  >  2:678992/1‑227 (MQ=255)
                                                                                                                                                                                       ttggACGCCTCTGTAGC‑ACAGC‑ATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCa                                 <  2:109074/225‑1 (MQ=255)
                                                                                                                                                                                        aagagagagctagagc‑acagc‑acAAAAAAAAGAAAGCCGCCAAGGGAAAACAAAAGGGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCa                                 <  1:170891/214‑1 (MQ=255)
                                                                                                                                                                                         ggACGCCTCTGAAGC‑ACAGA‑ATAATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCa                                 <  2:594253/225‑1 (MQ=255)
                                                                                                                                                                                          gACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGAGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGACAGCCGCGACCACCACCCGCTCATTCGGGTTTTTGGTAATTGCCGGGGTCATCCCCCCCAACGAGCAGtt                              >  1:594253/1‑228 (MQ=255)
                                                                                                                                                                                                ctcagtagacaaGC‑AACAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTtcatca                         <  2:26838/217‑1 (MQ=255)
                                                                                                                                                                                                 tttttgacgcccaaattaaaaaaaaaaCCCCGCAAATGGAAAAGATAAGTGCCGCACAGGCATACCCGACGGTTGAGTCAGAAACAAGCGACTGCAACGGAATAAACGCCACGGTACAACTCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGTGATTACCCCACCATTCGAGCATTtcatca                         <  1:430861/209‑1 (MQ=255)
                                                                                                                                                                                                 ctGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTtcatcaa                        >  2:421347/1‑226 (MQ=255)
                                                                                                                                                                                                 ctGTAGC‑AAAGC‑AGCAAAAAAATAAGGCCGCCAATGGAGAACATAAGTGCCGCACAGGCAGAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGa                       <  2:555650/228‑1 (MQ=255)
                                                                                                                                                                                                 cgtgtgagccatgaaaaaaaaaaaaaaGGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGTTTGAGTCAGAAAAACGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTtcatca                         <  1:222383/209‑1 (MQ=255)
                                                                                                                                                                                                   tgtgtgccagaaaacaaaaaaaaaaAGCCGCCAAAGGAAAACACAAGGTCCGCACAGGCATAACCGACGGTGGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATCCGAGCAGTTCATCAGa                       <  1:452606/211‑1 (MQ=255)
                                                                                                                                                                                                          caGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGGAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGa               >  2:494195/1‑228 (MQ=255)
                                                                                                                                                                                                           ccagaacatAAAAAAAAATCCGCCAATGGAGAACATAAGTTCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGa               <  2:167684/222‑1 (MQ=255)
                                                                                                                                                                                                           aGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCACCGGAATAAACGCCACGGTACAAATCAATAGCCCTATCGTTGCACCACCCCGGCGGTAAGCCGCGAGCTGAGCCCGCTCATTCGGGTTTTGGGTATTTGCCGGGGTCATCGCCCCATACGACCAGTGCCTCAg                        >  1:555650/1‑218 (MQ=255)
                                                                                                                                                                                                             c‑atcatAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGGGCAGTTCATCAGACGGTATGAAAg            >  2:294824/1‑228 (MQ=255)
                                                                                                                                                                                                             c‑atcatAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGc                                  >  2:170891/1‑206 (MQ=255)
                                                                                                                                                                                                                   tAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCAGCGCGCCGTCATTCGGGGATTCGGTAATTGCCGGGATCATCGCACCATCCCAGCAGTTCATCAGCCTGTCTGAAAGTCCaa       >  1:167684/1‑228 (MQ=255)
                                                                                                                                                                                                                   tAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGGTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCCCATTCGGGTTTTTGGTAATTGCCGGGATCCTCGCACCAAAAGCGCAGTTTATCAGACTGAATGAAAAGCCaa       >  1:88319/1‑228 (MQ=255)
                                                                                                                                                                                                                   tAAAAATAAAGCCGCCAATGGAGAAAATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGAGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACACAGGCGGTAAGCCGCGAGCCACGCGCGCTCCTTCGGGTTTTTGGGAATTGCCGGGATCATCGCACCATACGCACAGTTCATACAACTGGATGCAAGTCCaa       >  1:162900/1‑228 (MQ=255)
                                                                                                                                                                                                                   tAAAAAAAAGACCGCCAATGGAAAAAAGAAGGGCCGCACAGGCATAACCGCCGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATCAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCTCCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCaa       <  1:294824/228‑1 (MQ=255)
                                                                                                                                                                                                                   aaaaaaaaaAGCCGACAAAGGAGAACATAAGGGGCGCCCAGGCATAACCGACGGTTGAGTCAAAAAAAAGCGACGGCCACGGAATAAACGGCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCaa       <  1:421968/227‑1 (MQ=255)
                                                                                                                                                                                                                   aaaaaaTAAAGCCGGCAATGGAAAAAAAAAGTGCCGCACAGGCATAACCGTCGGGTGAGACAGAAAAAAGCGACTGCAACGGAATTAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCaa       <  1:598562/227‑1 (MQ=255)
                                                                                                                                                                                                                     aaaaaaaaGCCGACAATGGAAAACAAAAGGGCCGCACAGGGCTAACCGACGGTTGAGTCAGAAAAACGCGACTGCAACGGAATAAACGCCAAGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCaa       <  1:146193/226‑1 (MQ=255)
                                                                                                                                                                                                                     aaaaaaaaACCGCCAATGGAAAACAACAGGGCCACACAGGCATAACCGACGGTTGAGTCCGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTTTGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCaa       <  1:695573/226‑1 (MQ=255)
                                                                                                                                                                                                                     aaaaTAAAGCCGGCAATGGACAACATAAGTGCCCGACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCGCCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCaa       <  2:302770/226‑1 (MQ=255)
                                                                                                                                                                                                                          caaGCCGCCAAGGGTGAACATAAGTCCCGCACAGTCATAACCGACGGTTGAGTCAGAAAAAGGCGACTGCAACGGAATTAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAAcatc  <  2:666998/225‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
TTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACA‑TAATGCTCTTTGACGCCTCTGTAGC‑ACAGC‑ATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATC  >  NC_000913/4064954‑4065394

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