Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,754 G→A 5.4% intergenic (+34/‑205) serS → / → dmsA seryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913940,7540GA 5.4% 298.2 / 5.9 92intergenic (+34/‑205)serS/dmsAseryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base G (41/46);  new base A (2/3);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-03

ACCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGAGCAAGCTTGAGTGAGCCATTATGAAAACGAAAATCCCTGATGCGGTATTGGCTGCTG  >  NC_000913/940485‑940995
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
aCCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgc                                                                                                                                                                                                                                                >  2:267295/1‑273 (MQ=255)
   gTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGAAGGCTGATGGTCGTATTGAAGTACAAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATATAAAAAAGGCCTACGGGCGCttt                                                                                                                                                                                                                                           >  4:185689/1‑275 (MQ=255)
         tcctgtcatGCTCCAACGTTTGGGATGTCCAGGGAGGTCGGAGGCAGGAACGTTGCCGCAGCAAGTCGGACAAGCAAACCCGTCTGGTTCATACCCTGAACGGGTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGt                                                                                                                                                                                                                                     <  2:394716/267‑1 (MQ=255)
           tctTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                                                                                                                                                                       >  2:234136/1‑271 (MQ=255)
               ccTGCTCCAACGTTTGTGATTTTCAGGAACGTCGTATGCAGGTACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtct                                                                                                                                                                                                                                   <  4:234132/271‑1 (MQ=255)
                      caactgttGGGGTTTCCAGGATCGTCGTATGCGGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATa                                                                                                                                                                                                                         <  4:267291/268‑1 (MQ=255)
                      cAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATa                                                                                                                                                                                                                         <  2:251252/274‑1 (MQ=255)
                                     ccAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAct                                                                                                                                                                                                         >  2:144743/1‑275 (MQ=255)
                                      cAGGCACGTCGTATGCAGGCACGTGTCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                                                                                                                        <  2:73979/275‑1 (MQ=255)
                                         gatcgtcgtATGGAGGCACGGTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTTGTCGTACGCTGGTTGCAGTAAGGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                                                                                                                                                                      <  1:19409/271‑1 (MQ=255)
                                                gtATGCAGGCACGTGGCCGCAGCAAGTCGTACAAGCAAACCCGTTTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                                              <  2:407869/275‑1 (MQ=255)
                                                cagtGCAGGCACGTTGCCGCAGCGAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                                              <  2:379317/272‑1 (MQ=255)
                                                  gtgtaGGCCCTTTGCCGCCGCAATCCGGACAAGAAAACCCTTCTGGTTCATCCCCGGAACGTTTCTGTTATGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCATGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACg                                                                                                                                                                                            <  4:226907/271‑1 (MQ=255)
                                                              tGCCGCAGCGAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGATAATACCCACTTTTTCTGAATCTAAAAAGCGCCTGAGGGCGCTTTTTTTGTCTCACTTTGCTACCGAACAATAATTACTCCCCACTTACAAGTAAtacttact                                                                                                                                                                                 >  4:379312/1‑270 (MQ=255)
                                                              tGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTa                                                                                                                                                                                                >  2:226911/1‑259 (MQ=255)
                                                                                                         tGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTCATCCTCACTTACACGTAATCCCTCTTTCGAGTGAAAATTCACCTATTCCTTTTATTTTCAAATTATTcca                                                                                                                                     >  4:398394/1‑273 (MQ=255)
                                                                                                         tGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGa                                                                                                                                     >  2:211371/1‑275 (MQ=255)
                                                                                                          gAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATAcc                                                                                                                                                                                                                       >  4:405010/1‑192 (MQ=255)
                                                                                                              ggTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                                                                                                                                                                   <  2:405015/192‑1 (MQ=255)
                                                                                                                ttCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTAACAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTTACTTACACGTCATACTACTTTCGAGTCAAAATACCCCTATCACTTTTATTTTCAAATTTTTCGATGtatcc                                                                                                                              >  4:456644/1‑273 (MQ=255)
                                                                                                                 tctggtctgtCTGTTTGTCGTACGCTGGTGGGAGGAATGGAAAACGATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGGACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATGTTCAAATTATTCGATGTATACa                                                                                                                             <  2:280680/275‑1 (MQ=255)
                                                                                                                 tctggtctggCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACt                                                                                                                                                                                                  >  1:113593/1‑206 (MQ=255)
                                                                                                                  gtggtatggCCGTTGGTCGTACGCGGGATGCAGTAATGGGAAACTATCAGCAGGCTGCTGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTATCCCTGTGATACCGAACAATAATTACTCCACACTTACACGTAATACTACTTTCGAGTGAAAATCTTCCTATCTCTTTGATTTTCAAATTATTCGATGTATACa                                                                                                                             <  2:204235/273‑1 (MQ=255)
                                                                                                                    ggtctggCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGTACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCg                                                                                                                                      <  2:456649/263‑1 (MQ=255)
                                                                                                                     gtctggCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                                              <  3:113593/206‑1 (MQ=255)
                                                                                                                          gCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                                                                                                                                                                            >  2:366891/1‑219 (MQ=255)
                                                                                                                              ttGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaa                                                                                                                                                                        <  4:366886/219‑1 (MQ=255)
                                                                                                                                ggTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGAAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa                                                                                                              <  2:185691/275‑1 (MQ=255)
                                                                                                                                    gTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAAtt                                                                                                                                                                                                            >  2:19319/1‑177 (MQ=255)
                                                                                                                                     tACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttt                                                                                                                                                                                                                                         >  4:193988/1‑147 (MQ=255)
                                                                                                                                      aCGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttaga                                                                                                                                                                                                                                      <  7:458559/149‑4 (MQ=255)
                                                                                                                                        tcTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                                                                                                                                                                        <  4:19319/176‑1 (MQ=255)
                                                                                                                                        gCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcc                                                                                                                                                                                                       <  2:300066/178‑1 (MQ=255)
                                                                                                                                         cTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGt                                                                                                                                                                                                                                     <  2:193990/147‑1 (MQ=255)
                                                                                                                                          tGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGa                                                                                                    <  2:398399/275‑1 (MQ=255)
                                                                                                                                                     aTGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATAcc                                                                                                                                                                                                                       >  7:73669/1‑149 (MQ=255)
                                                                                                                                                      cggAAAAGTTTCAGCAGGCTGAGTGTTGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATNTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCt                                                                                                                        <  3:104808/242‑1 (MQ=255)
                                                                                                                                                       ggAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGa                                                                                                                                                                                                                     >  5:48487/1‑149 (MQ=255)
                                                                                                                                                        gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttt                                                                                                                                                                                                                                           <  7:432901/126‑1 (MQ=255)
                                                                                                                                                        gAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttt                                                                                                                                                                                                                                           >  8:432901/1‑126 (MQ=255)
                                                                                                                                                              tATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                                                                                                                                                                              <  8:73669/149‑1 (MQ=255)
                                                                                                                                                               ttCAGCAGGCGGATGGTCGTATTGAAGTACCAGAAGGTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGtt                                                                                <  4:211368/273‑1 (MQ=255)
                                                                                                                                                                tCAGCAGGCTGATGGTCGTATTGAGGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACt                                                                                                                                                                                                  >  4:525122/1‑159 (MQ=255)
                                                                                                                                                                tCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGTACAATAATTAc                                                                                                                                                                                                          >  4:545172/1‑151 (MQ=255)
                                                                                                                                                                   gcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAct                                                                                                                                                                                                         >  7:462926/1‑149 (MQ=255)
                                                                                                                                                                   gcagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACaa                                                                                                                            >  1:104808/1‑226 (MQ=255)
                                                                                                                                                                    cagGCTGATGGTCGTATTGAGGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                                                                                                                                                                              <  2:525127/159‑1 (MQ=255)
                                                                                                                                                                    cagGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGTACAATAATTActcct                                                                                                                                                                                                      <  2:545177/151‑1 (MQ=255)
                                                                                                                                                                       gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                                                                                                                                                                       >  4:225217/1‑115 (MQ=255)
                                                                                                                                                                       gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcctc                                                                                                                                                                                                     >  7:213523/1‑149 (MQ=255)
                                                                                                                                                                       gCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcctc                                                                                                                                                                                                     >  3:10127/1‑149 (MQ=255)
                                                                                                                                                                           aTGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtct                                                                                                                                                                                                                                   <  2:225221/115‑1 (MQ=255)
                                                                                                                                                                           aTGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                                                                                                                                                                 <  1:10127/149‑1 (MQ=255)
                                                                                                                                                                              gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                                                                                                                                                                       >  2:2569/1‑108 (MQ=255)
                                                                                                                                                                                 gTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcctc                                                                                                                                                                                                     >  1:130065/1‑139 (MQ=255)
                                                                                                                                                                                  tATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtct                                                                                                                                                                                                                                   <  4:2569/108‑1 (MQ=255)
                                                                                                                                                                                     tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                                                                                                                                                                 <  3:130065/139‑1 (MQ=255)
                                                                                                                                                                                     tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAta                                                                                                                                                                                       <  8:462926/149‑1 (MQ=255)
                                                                                                                                                                                     tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAta                                                                                                                                                                                       <  8:213523/149‑1 (MQ=255)
                                                                                                                                                                                                 aGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGt                                                                                                                                                                           <  6:48487/149‑1 (MQ=255)
                                                                                                                                                                                                          tCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTa                                                                                                                                                                  >  6:382618/1‑149 (MQ=255)
                                                                                                                                                                                                                        gACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                                                                                                                                                                       >  3:141151/1‑66 (MQ=255)
                                                                                                                                                                                                                          cTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                                                                                                                                                                             >  2:546043/1‑122 (MQ=255)
                                                                                                                                                                                                                           tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                                    <  6:84235/82‑1 (MQ=255)
                                                                                                                                                                                                                           tGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                                    >  5:84235/1‑82 (MQ=255)
                                                                                                                                                                                                                            ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtct                                                                                                                                                                                                                                   <  1:141151/66‑1 (MQ=255)
                                                                                                                                                                                                                             gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                                                                                                               >  6:175916/1‑117 (MQ=255)
                                                                                                                                                                                                                             gAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                                                                                                                                                                               <  5:175916/117‑1 (MQ=255)
                                                                                                                                                                                                                              aaTATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGa                                                                                                                                                                         <  4:546038/122‑1 (MQ=255)
                                                                                                                                                                                                                               atatatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtacta                                                                                                                                                                                    >  5:318507/1‑110 (MQ=255)
                                                                                                                                                                                                                               atatatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtacta                                                                                                                                                                                    <  6:318507/110‑1 (MQ=255)
                                                                                                                                                                                                                                   atTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGaa                                                                                                                                                                                                                    >  2:102024/1‑74 (MQ=255)
                                                                                                                                                                                                                                       gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaat                                                                                                                                                                                                                <  4:102023/74‑1 (MQ=255)
                                                                                                                                                                                                                                       gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCta                                                                                                        >  2:386653/1‑178 (MQ=255)
                                                                                                                                                                                                                                       gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTACCTTTGATACCg                                                                                                                                                                                                                      >  4:522441/1‑68 (MQ=255)
                                                                                                                                                                                                                                           aTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACa                                                                                                                                                                                                                  <  2:522446/68‑1 (MQ=255)
                                                                                                                                                                                                                                           aTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGa                                                                                                    <  4:386648/178‑1 (MQ=255)
                                                                                                                                                                                                                                            tACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcctc                                                                                                                                                                                                     >  4:173554/1‑80 (MQ=255)
                                                                                                                                                                                                                                            tACACAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCACTTTATTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTCACACGTAATACTACTTTAGAGTGAAAACCTACCAATCTCTTTGATTTTCAAATTATTAGATGTTTACAAGCCTAGATACAGAAGTCATATAGAAATAATTACACAATACGGTTTGTTACTGGAACAAATCGTGAGCAAGCTTTAGGTAGCTATTCTGAAAACAAAAATCCCTGATACggagttgtcccagg  >  4:25141/1‑263 (MQ=255)
                                                                                                                                                                                                                                               ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACg                                                                                                                                                                                            <  6:136988/86‑1 (MQ=255)
                                                                                                                                                                                                                                               ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACg                                                                                                                                                                                            >  5:136988/1‑86 (MQ=255)
                                                                                                                                                                                                                                                cAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                                                                                                                                                                 <  2:173556/80‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGt                                                                                                                                                                           >  2:522795/1‑101 (MQ=255)
                                                                                                                                                                                                                                                     tttCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaa                                                                                                                                                                       <  4:522790/101‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAt                                                                                                                                                    <  6:263955/112‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAt                                                                                                                                                    >  5:263955/1‑112 (MQ=255)
                                                                                                                                                                                                                                                                aaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                                                                                                                 <  5:131421/96‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                                                                                                                                                                 >  6:131421/1‑96 (MQ=255)
                                                                                                                                                                                                                                                                 aaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATAc                                                                                                                              >  2:233928/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                     gcgcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGc                                                                                                                          <  4:233924/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGAGCAAGCTtgagtg                                              >  4:118115/1‑204 (MQ=255)
                                                                                                                                                                                                                                                                           gCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGAGCAAGCTTGAGTGAGcc                                          <  2:118116/204‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
ACCGTGAGATCTCTTCCTGCTCCAACGTTTGGGATTTCCAGGCACGTCGTATGCAGGCACGTTGCCGCAGCAAGTCGGACAAGAAAACCCGTCTGGTTCATACCCTGAACGGTTCTGGTCTGGCTGTTGGTCGTACGCTGGTTGCAGTAATGGAAAACTATCAGCAGGCTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTACTGGAATCAATCGTGAGCAAGCTTGAGTGAGCCATTATGAAAACGAAAATCCCTGATGCGGTATTGGCTGCTG  >  NC_000913/940485‑940995

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: