Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,198,449 T→C 7.1% pseudogene (1976/2847 nt) yehH → DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,198,4490TC 7.1% 183.9 / 4.2 56pseudogene (1976/2847 nt)yehHDUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base T (32/20);  new base C (2/2);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CTGTTGGCAAGGATGGGTTTTCAACGACGGATCGCTGGTACATATATTAAAGCTCCCGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAAC  >  NC_000913/2198175‑2198735
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                              
cTGTTGGCAAGGATGGGTTTTCAACGACGGATCGCTGGTACATATATTAAAGCTCCCGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACtt                                                                                                                                                                                                                                                                                   >  1:276040/1‑288 (MQ=255)
cTGTTGGCAAGGATGGGTTTTCAACGACGGATCGCTGGTACATATATTAAAGCTCCCGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACtt                                                                                                                                                                                                                                                                                   >  2:134362/1‑288 (MQ=255)
    tGGCAAGGATGGGTTTTCAACGACGGATCGCTGGTACATATATTAAAGCTCCCGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACg                                                                                                                                                                                                                                                                              <  1:359102/289‑1 (MQ=255)
                           cgGATCGCTGGTACATATATTAAAGCTCCCGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCa                                                                                                                                                                                                                                                       >  1:97961/1‑289 (MQ=255)
                           aggATCGCGTGTACATAGATTAAAGGTCCCGAAGCGTTAGTTGAATCATTGCTAAATGAAGATTATTCAACCGTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTAGTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCa                                                                                                                                                                                                                                                       <  2:297175/288‑1 (MQ=255)
                                        cATATATTAAAGCTCCCGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTACATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGACCATTACGAATTAcct                                                                                                                                                                                                                                          >  1:297175/1‑287 (MQ=255)
                                               tAAAGCTCCCGAAGCGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTACCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGa                                                                                                                                                                                                                                  >  2:225118/1‑289 (MQ=255)
                                                 aaGCTCCCGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAg                                                                                                                                                                                                                                 <  1:29535/289‑1 (MQ=255)
                                                      cccGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCCGCGAGTGAg                                                                                                                                                                                                                             >  2:457021/1‑288 (MQ=255)
                                                      cccGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGa                                                                                                                                                                                                                              >  1:71499/1‑287 (MQ=255)
                                                      cccGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTCAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGa                                                                                                                                                                                                                              >  1:332678/1‑287 (MQ=255)
                                                              ggTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTc                                                                                                                                                                                                                      <  2:71499/287‑1 (MQ=255)
                                                                             gCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTg                                                                                                                                                                                                       >  2:174968/1‑287 (MQ=255)
                                                                                       gATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACaaa                                                                                                                                                                                            >  2:266246/1‑288 (MQ=255)
                                                                                             tCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGaaaa                                                                                                                                                                                     <  2:276040/289‑1 (MQ=255)
                                                                                                         aGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTCAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGaaa                                                                                                                                                                         >  2:120142/1‑289 (MQ=255)
                                                                                                         aGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTCAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGaaa                                                                                                                                                                         <  2:332678/289‑1 (MQ=255)
                                                                                                              aTTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTCAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGc                                                                                                                                                                     <  1:120142/288‑1 (MQ=255)
                                                                                                                                     ccGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATAcc                                                                                                                                              >  2:413524/1‑288 (MQ=255)
                                                                                                                                                    aGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTAt                                                                                                                                                                                                                                                                          <  1:103583/149‑1 (MQ=255)
                                                                                                                                                    aGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTAt                                                                                                                                                                                                                                                                          >  2:103583/1‑149 (MQ=255)
                                                                                                                                                    aGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTa                                                                                                                               >  2:7514/1‑288 (MQ=255)
                                                                                                                                                                               tGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCaa                                                                                                   <  1:174968/289‑1 (MQ=255)
                                                                                                                                                                                      aGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGt                                                                                            >  1:141024/1‑289 (MQ=255)
                                                                                                                                                                                         aGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTg                                                                                           <  1:134362/287‑1 (MQ=255)
                                                                                                                                                                                         aGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTg                                                                                           <  2:141024/287‑1 (MQ=255)
                                                                                                                                                                                                  aTGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGa                                                                                 <  1:7514/288‑1 (MQ=255)
                                                                                                                                                                                                          cccTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGGCAAGCCCAGGATaagcccc                                                                         >  2:12327/1‑283 (MQ=255)
                                                                                                                                                                                                            cTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCcgcg                                                                      <  1:225118/289‑1 (MQ=255)
                                                                                                                                                                                                                     tctcCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGaaaa                                                                                                                                                                        <  2:311490/182‑1 (MQ=255)
                                                                                                                                                                                                                     tctcCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGaaaa                                                                                                                                                                        >  1:311490/1‑182 (MQ=255)
                                                                                                                                                                                                                      ctcCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCg                                                            <  1:266246/289‑1 (MQ=255)
                                                                                                                                                                                                                           caTACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGc                                                       <  1:170537/289‑1 (MQ=255)
                                                                                                                                                                                                                                 gataccgaATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAACCTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCCGGCGCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACACCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTCCCGGCATGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCAAAGGATACCCCCCGCGCTGCCCTGCGTACGCTTACCg                                                 >  2:380500/1‑289 (MQ=255)
                                                                                                                                                                                                                                   taccgaATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCCCGGCAAAGCCGGTGAAGcac                                                                                     >  2:272862/1‑251 (MQ=255)
                                                                                                                                                                                                                                   taccgaATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCATGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGaa                                               >  1:333629/1‑289 (MQ=255)
                                                                                                                                                                                                                                   taacgaATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTGGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGcac                                                                                     <  1:272862/248‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGaa                                                                                                                                                                          <  1:427999/161‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGaa                                                                                                                                                                          >  2:427999/1‑161 (MQ=255)
                                                                                                                                                                                                                                        aaTACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGca                                                                                      <  2:393956/245‑1 (MQ=255)
                                                                                                                                                                                                                                        aaTACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGca                                                                                      >  1:393956/1‑245 (MQ=255)
                                                                                                                                                                                                                                               cATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGtcatca                                    >  2:280243/1‑288 (MQ=255)
                                                                                                                                                                                                                                                    tACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCAtt                              >  1:212138/1‑289 (MQ=255)
                                                                                                                                                                                                                                                           tGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAgctggc                                                                                                                                                          >  1:323735/1‑158 (MQ=255)
                                                                                                                                                                                                                                                           tGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAgctggc                                                                                                                                                          <  2:323735/158‑1 (MQ=255)
                                                                                                                                                                                                                                                           tGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCaa                       >  1:486224/1‑289 (MQ=255)
                                                                                                                                                                                                                                                             gACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGa                     >  2:289424/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                tCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGaaa                   >  2:128493/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                       gTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCg            >  2:332604/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                          aCCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGt        >  1:413975/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                          aCCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCGGCCTTGCGTATGCTTCCCGAAAACGGTCATCAGCCATTACTGCAAGAAAACGCCCGCCGt        >  2:25887/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                gTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGtca                                       <  2:31862/252‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGtca                                       >  1:31862/1‑252 (MQ=255)
                                                                                                                                                                                                                                                                                gTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGccc             <  1:165535/278‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGccc             >  2:165535/1‑278 (MQ=255)
                                                                                                                                                                                                                                                                                gTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGAGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCGGCGGCGCTCGGCAAAGCCGGTGACGCACCGGATAACGCCAGCGCTGCCTTGCGTATGCTTACCGTAAACGGTCCTCAGCCCTTACTGCACGAAATCGCCCGACGTTATaa   >  2:274309/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                 ttAGTAAACGCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGAGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAAGCGCCCGACGTTATAAc  <  1:274309/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gtAGTAAACTCACTTGCACGCTAGCCACAAGAAGCCTTGCCCATCCCGAGTGACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCCAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAAc  <  2:297327/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                              
CTGTTGGCAAGGATGGGTTTTCAACGACGGATCGCTGGTACATATATTAAAGCTCCCGAAGCGGTAGTTGAGGCATGGCTAAATGAAGATTATTCAACCTTACTAAGTGAATTTAAGGTGTTTCATTCACCTACCGGGCATTATTGGCAGTTGGGGATTTTGACAACATTGCCGCTGGAGAAAGCAGTAAAAGCATGGAATGCCCTTACCCTATCTCCACATACCGATACCGAATACTCCATGTTACATTTTGGACTCAAAGGGTTACCTGGGTTAGTAAACTCACTTGCACGCTATCCACAAGAAGCCTTGCCCATCACGAATTACTTCGCAGCGAGTGAGCTGGCTCCTGCCGTCGCCCGTGCCTTCAACAAACTGAAAACGCTACGCGAAAACGCCCGTAGCTGGCTGTTGAAATACCCGGAACATGCCCTTACCGGCCTGCTGCCTGCGGCGCTCGGCAAAGCCGGTGAAGCACAGGATAACGCCCGCGCTGCCTTGCGTATGCTTACCGAAAACGGTCATCAGCCATTACTGCAAGAAATCGCCCGACGTTATAAC  >  NC_000913/2198175‑2198735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: