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→GGTyihOputative 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

TTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATC  >  NC_000913/4064954‑4065394
                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
ttaaataatatcacgcaaTACGGGGTGTAAAGAACCTAAACCACAAATTTGACAATAAAGGCCGCCAGGGTAAACAAGTTAACAATGCATAAAACCAGCAATGGCGGATTCCGAAATGCCGTGAAAAATGATTTGAGAATCCTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATCCTCTTTGAAGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATg                                                                                                                                                                                                                       <  1:248376/211‑1 (MQ=255)
ttAATATCATCCAACACATACTGGGTGTAAAAGACCGTAATCGCCAGTTTGATATAAAATGCCGCCAGGTTACAAAAATTAGCAATGCAAAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATg                                                                                                                                                                                                                       <  2:701239/228‑1 (MQ=255)
         ttCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCTCGCAAAATGATTTGAGAATACTGGCTTTATGTCAGCTTGGCTTTGTTACACATTAATGCCTTTTGACGCTTTGTAGGCCAGCAATCTTAAAAAAAAACAGCCAATTGAGGACaaaa                                                                                                                                                                                                             >  1:26838/1‑228 (MQ=255)
         ttCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCGTTAGCACACCATCATAAAAATAAAGCCGCCAATGGAGAACATa                                                                                                                                                                                                              >  2:652447/1‑228 (MQ=255)
             cacattccTGGGTTTAAAAGACCTGAACCGCAATTTTAATATAAAATGCCGCCAGGGTAAACAAGTTAGCAAGGCATAAAACAAGAAATGGCGGATGCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTg                                                                                                                                                                                                          <  1:24096/221‑1 (MQ=255)
              cacaTACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGGTTCCAAATAATGCTCTTTGACGCCTCTTGAGCCCGGCATCTTAAAAATAAAACCACCAATGGCGAACATAAATGc                                                                                                                                                                                                         >  1:522166/1‑228 (MQ=255)
                      gggTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAgg                                                                                                                                                                                                 >  2:485269/1‑228 (MQ=255)
                          taaaaggcccTGAAACGCCAGGTTGATAGTAAATGCCGACAGGGTAAACAAGTTAGCAATGCAAAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTGGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATa                                                                                                                                                                                             <  2:193512/220‑1 (MQ=255)
                           tAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTTAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATACTCTTTGACGCCTCTGTAGCACAGCATCCTAAAAATAAAGCCGCCAATGGAGAACATAAGGACCGCACAGGCCTaa                                                                                                                                                                                            >  2:11677/1‑228 (MQ=255)
                                   cTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTAGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCCCTTTGACGCTCCTGTAGAACAGCATCATCAAAATAAAACCGCCAATGCGTACCATAAGTACCACACAGCCAAAACCAAcaga                                                                                                                                                                                    >  1:314696/1‑225 (MQ=255)
                                      aaTCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATAACGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGCCGTTGTGACCACATAATGCTCTTTGACGCCTCTTTAGCCAAGCATCAATAAAATAAACCCGCCAATGGAGAACATAAATGCCTCACAgaaacaaacaaaagatag                                                                                                                                                                                 >  1:311860/1‑211 (MQ=255)
                                               tttGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGga                                                                                                                                                                                                                     >  2:56205/1‑183 (MQ=255)
                                                     ttaaaaTGCCGCGAGGTGAACCATTTTAGCAACGAAAAAACCAAGCAAGGGATGTTCCCGAAATGTCGCGAAAATTGATTTGATCATACTGGTTTGAGGCCCGATTGGCGTTGTACCCACATAATGCTCTGTGACGCCGCTGTAGCACGGAATGAAAAAAATAAAGCCGCCAATGGAGAACAGAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGa                                                                                                                                                                  <  1:11677/225‑1 (MQ=255)
                                                     atTAAATGCCGCCAGGGTTAACAAGTTAGAAAGGCATAAAACCAGCAATGGCGGATCCCGAAATATCGCACAAAATGATTTGAGAATACTGGCTTTATGTGCGGTTGGGGTTGTGACCACATAATGCTCTTTGACGCCTCTGTAGCACAGCACCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGa                                                                                                                                                                  <  2:331562/228‑1 (MQ=255)
                                                                ccAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGGCCGGTTGGCGTTGTGTCCACATAATGCTCTTTTACGCCTCTGTAGGCCAGCATCCTAACACTAAAACCGCCAATGGGGAACATAAGTGCCGCACAGCCATAACAGAAGGTTGCACAAGAAAAAAGCGAcgccaaacaaa                                                                                                                                                       >  1:331562/1‑218 (MQ=255)
                                                                ccAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACCTGCTTTATGCCCGGTTGGCGTTGTGACCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAACCGCCAATGGAGAACATAAGGGCCGAACAGCCATAACCGACGCTTACATCGGAAAAAAGCAAcggaaaacgaa                                                                                                                                                       >  1:445207/1‑218 (MQ=255)
                                                                    ggtaacaaaGTTAGAAATTCATAAAACCATCAATGGCGGTTTCCGAAATGTCGGGCAAAATGATTTGAGAATACTGGTTTTATGTCCGGTTGGCGTTGGGTCCACATAATGCTCTTTGACGCCTCTGAAGCACAGCATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATaaa                                                                                                                                                   <  1:652447/221‑1 (MQ=255)
                                                                         acaATTTAGAAATGCATAAAACCAGCATTGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCa                                                                                                                                              <  2:522166/228‑1 (MQ=255)
                                                                         acaAGTTAGAAATGCTTTAAACCAGCAAGGGCTGTTTCCGAAGTTTCGAGGAAAATGATTTGAGAATCCGGGCCTTATGTCCGGTTGGAGTTGTGACCAAAAAATGCCCTTTGGCGCCTCTGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCa                                                                                                                                              <  1:217449/228‑1 (MQ=255)
                                                                          aataagagtaaggCATAAAACCGACAATGGCGAGTTCCGAAATGTCGCGAAAAATGATTGTAAAAGACAGGCTTTAGGTCCGGTTGGCTTTGTGTCAGCATAATGCTTTTTGACGCCCCTGAAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCa                                                                                                                                              <  2:163547/215‑1 (MQ=255)
                                                                              ttAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCCCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGAATGGAAAGGAGTAAAGGCCACggg                                                                                                                                          >  2:430861/1‑226 (MQ=255)
                                                                                aGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTAc                                                                                                                                        >  2:546369/1‑227 (MQ=255)
                                                                                     tGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGACCATAAGTGCCGCACAGGACTAACCGACGGTTCAGTCAGAAAAAAGCGACTGCAAAGGAAAAAAACCCACGGGAAAAAtaa                                                                                                                                  >  1:605594/1‑226 (MQ=255)
                                                                                     tGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATAATGGCTTTATGTCCAGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGTAGAACATAAGTGCCGCACAGGCAAAACCAACGACGCAGTAGGAAAAAAGCAATTGCACCGAAAAAAACACCACGGAAAAAACCa                                                                                                                                  >  1:302770/1‑228 (MQ=255)
                                                                                     tGCATAAAACCAGAAATGGCGGATTCCGAAATATCGCGCAAAATGATTTAAGAATACTGGCTTTATGGCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAAGCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGTAAAAAACGAACGCAACGGAATAAACGCACCGGTAAAAAACa                                                                                                                                  >  2:660740/1‑228 (MQ=255)
                                                                                      gCATAAAACCAGAAAGGGGGGATTCCGAAATAGCGCGAAAAATGATTTGAGAATACTGGTTTTATGCCCGGTTGGAGTTGTGCCCACATAATGCTCTTTGACGCCTCTGGAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCaa                                                                                                                                 <  2:314696/228‑1 (MQ=255)
                                                                                       cATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCCCAGCATCATCAAACTAAAGCCGCCAATGGAGACCATACGGACCGGACAGGCATAACCGACGGATGAGTCCGGAAAAAGCGACTGCAACGCAATAAACACCCCGGGAAAAACCAAt                                                                                                                                >  1:38197/1‑228 (MQ=255)
                                                                                       cATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAAt                                                                                                                                >  2:222383/1‑228 (MQ=255)
                                                                                                 gCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTTGCTTTATGTCCGGTTGGCGTTGTGTCCAAATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAACCGCCAATGGAGAACCTAAGTGCCCGACAGGCCTAAACGACGCCTGAGTCAGAAAAAAGAGACTGAAACGGAATAAACCCCACGGTACAAAACacgtatcttataa                                                                                                                      >  1:163547/1‑216 (MQ=255)
                                                                                                 gCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAACAATAAAGCCGCCATTGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAACGACTGCAACGAAATAAACACCAAGGTGCAAATCCATAAcccgaagg                                                                                                                      >  1:549320/1‑223 (MQ=255)
                                                                                                            ttCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCaccacc                                                                                                            >  2:421968/1‑227 (MQ=255)
                                                                                                                       cgcgCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTCGCACCGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAACAGCGACTGCAACGGAATAAACGCCACGGTACAAAACAATACCCCTATGGTAGCACCCCCCTGGCgtgagt                                                                                                 >  2:695573/1‑222 (MQ=255)
                                                                                                                       cgcgCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAATCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGCCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAAAAGCCCAATGGTTACACACACCTGCAGGTAAg                                                                                                 >  1:443082/1‑227 (MQ=255)
                                                                                                                                  aTTTGAAAATGCTGGCTTTATGTCCGGTTGGGGTTGTGTCCACAAAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTgcg                                                                                     <  2:490876/228‑1 (MQ=255)
                                                                                                                                    ttGAAAATTCTGGTTTTATCTGCGTGTGGCGTTGTACACAAAGTATCCTATTAGGCCCCTATGAAGCAACAGCATAAAAAAAAAAAGCTGCCAATGGAGAAAATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCCCCACCCTGGCGGTAAGCCGCGAGCTgcgcg                                                                                   <  1:490468/228‑1 (MQ=255)
                                                                                                                                                        cccGGTTGGCGTGGTGTACACATAATGCCCTTTGACGCCTCTGTAGCCAAGCATCATAAAAATAAAGCCGCCAATGGAGAACAAAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGt                                                               <  2:605594/227‑1 (MQ=255)
                                                                                                                                                          cGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGCCGGAAAGCCGCAAGCTCCCGGGGCCCAATCGGGGTTTTGGTaa                                                             >  1:109074/1‑228 (MQ=255)
                                                                                                                                                          cGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACAAACCTGGCGGTAAGCCGCGACCTGCCGGCCCTCATTCGGGTTTTTGGTaa                                                             >  2:598562/1‑228 (MQ=255)
                                                                                                                                                          cGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTCGAGTCAGAAAAAAGCGCCTGCACCGGAATAAAAGCCCCGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGCAAGCCGCGAGCTGACCGCACTCCTTCGGGGTTTTGGGaa                                                             >  1:419815/1‑228 (MQ=255)
                                                                                                                                                          cGGGTGGCGTTGGCACAAAATAATGCTTGTTGGAGCCTCTGGAGCACAGCAACAAAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCAAAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTaa                                                             <  1:485269/228‑1 (MQ=255)
                                                                                                                                                            gTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAAt                                                            >  2:347518/1‑227 (MQ=255)
                                                                                                                                                                          cacaTAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCAAAGGCATAACCAACGGTTGAGACAGACAAAAGCGACTGCAACGGATTAAAAGCCACAGACAAAATCAAAAACCCTATGGGTGAACAACCCAGGAGGAAAGCGCCAACCAGCACCCCCCCATCCGGGGTTTGGGTAATTGCCAGGACCATAGc                                             >  1:45562/1‑228 (MQ=255)
                                                                                                                                                                          aaaaataatgcatggagACGCCCTTGTACCACAGAAACAAAAAAATAAGGCCGCCAATGGAAAAAATAAGGGCCGAACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAAGCAATAGCCCTATGGTTGCTCCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCg                                              <  1:678992/212‑1 (MQ=255)
                                                                                                                                                                           acaTAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCTCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGatcat                                                <  2:45562/224‑1 (MQ=255)
                                                                                                                                                                                   gCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACg                                                                                                                                            >  2:284876/1‑125 (MQ=255)
                                                                                                                                                                                     ttttggcggCCGCTGAAGCAAAGCAAAAAAAAAAAAAAGCGGCCAAGGGAAAACAAAAGTGCCGGACAGGCAAAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAg                                   <  1:347518/220‑1 (MQ=255)
                                                                                                                                                                                     tctTGGACGCCTCTGTAGAACAGCATAATAAAAATAAAGCCTCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCGCCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAg                                   <  2:311860/228‑1 (MQ=255)
                                                                                                                                                                                      ctTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGAAATTGCCGGGATCATCGCACCATACGac                                   >  2:114976/1‑226 (MQ=255)
                                                                                                                                                                                       tttGGCGCCTTGTTAGCACATAATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCa                                 <  2:445207/228‑1 (MQ=255)
                                                                                                                                                                                       tttGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGGCCGCGCTCCTTCGGGCTTTTGGTAATTGCCGGGATCATCGCACCATACGAGc                                  >  2:678992/1‑227 (MQ=255)
                                                                                                                                                                                       ttggACGCCTCTGTAGCACAGCATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCa                                 <  2:109074/225‑1 (MQ=255)
                                                                                                                                                                                        aagagagagctagagcacagcacAAAAAAAAGAAAGCCGCCAAGGGAAAACAAAAGGGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCa                                 <  1:170891/214‑1 (MQ=255)
                                                                                                                                                                                         ggACGCCTCTGAAGCACAGAATAATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCa                                 <  2:594253/225‑1 (MQ=255)
                                                                                                                                                                                          gACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGAGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGACAGCCGCGACCACCACCCGCTCATTCGGGTTTTTGGTAATTGCCGGGGTCATCCCCCCCAACGAGCAGtt                              >  1:594253/1‑228 (MQ=255)
                                                                                                                                                                                                ctcagtagacaaGCAACAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTtcatca                         <  2:26838/217‑1 (MQ=255)
                                                                                                                                                                                                 tttttgacgcccaaattaaaaaaaaaaCCCCGCAAATGGAAAAGATAAGTGCCGCACAGGCATACCCGACGGTTGAGTCAGAAACAAGCGACTGCAACGGAATAAACGCCACGGTACAACTCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGTGATTACCCCACCATTCGAGCATTtcatca                         <  1:430861/209‑1 (MQ=255)
                                                                                                                                                                                                 ctGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTtcatcaa                        >  2:421347/1‑226 (MQ=255)
                                                                                                                                                                                                 ctGTAGCAAAGCAGCAAAAAAATAAGGCCGCCAATGGAGAACATAAGTGCCGCACAGGCAGAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGa                       <  2:555650/228‑1 (MQ=255)
                                                                                                                                                                                                 cgtgtgagccatgaaaaaaaaaaaaaaGGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGTTTGAGTCAGAAAAACGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTtcatca                         <  1:222383/209‑1 (MQ=255)
                                                                                                                                                                                                   tgtgtgccagaaaacaaaaaaaaaaAGCCGCCAAAGGAAAACACAAGGTCCGCACAGGCATAACCGACGGTGGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATCCGAGCAGTTCATCAGa                       <  1:452606/211‑1 (MQ=255)
                                                                                                                                                                                                          caGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGGAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGa               >  2:494195/1‑228 (MQ=255)
                                                                                                                                                                                                           ccagaacatAAAAAAAAATCCGCCAATGGAGAACATAAGTTCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGa               <  2:167684/222‑1 (MQ=255)
                                                                                                                                                                                                           aGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCACCGGAATAAACGCCACGGTACAAATCAATAGCCCTATCGTTGCACCACCCCGGCGGTAAGCCGCGAGCTGAGCCCGCTCATTCGGGTTTTGGGTATTTGCCGGGGTCATCGCCCCATACGACCAGTGCCTCAg                        >  1:555650/1‑218 (MQ=255)
                                                                                                                                                                                                             catcatAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGGGCAGTTCATCAGACGGTATGAAAg            >  2:294824/1‑228 (MQ=255)
                                                                                                                                                                                                             catcatAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGc                                  >  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)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
TTAATATCATTCAGCACATACTGGGTGTAATAGACCTGAATCGCCAGTTTGATATTAAATGCCGCCAGGGTACACAAGTTAGCAATGCATAAAACCAGCAATGGCGGATTCCGAAATATCGCGCAAAATGATTTGAGAATACTGGCTTTATGTCCGGTTGGCGTTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAACGGAATAAACGCCACGGTACAAATCAATAGCCCTATGGTTGCACCACCCTGGCGGTAAGCCGCGAGCTGCGCGCGCTCATTCGGGTTTTTGGTAATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATC  >  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: