Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,360 A→C I58S (ATC→AGC)  ymfD ← e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,3600AC100.0% 369.3 / NA 94I58S (ATC→AGC) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (46/48);  total (46/48)

TATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTT  >  NC_000913/1197134‑1197564
                                                                                                                                                                                                                                  |                                                                                                                                                                                                            
tATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAg                                                                                                                                                                                                  >  2:1127669/1‑239 (MQ=255)
  ttttGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGaa                                                                                                                                                                                                >  1:555082/1‑239 (MQ=255)
      gTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATc                                                                                                                                                                                            <  1:86263/239‑1 (MQ=255)
       tCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCt                                                                                                                                                                                           >  1:673792/1‑239 (MQ=255)
            tGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGtt                                                                                                                                                                                      >  2:346350/1‑239 (MQ=255)
                        aGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGACaa                                                                                                                                                                          >  2:289204/1‑239 (MQ=255)
                                        aCATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGaa                                                                                                                                                                                                >  2:1094875/1‑201 (MQ=255)
                                        aCATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGaa                                                                                                                                                                                                <  1:1094875/201‑1 (MQ=255)
                                        aCATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTAAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTgg                                                                                                                                                            <  1:1124467/237‑1 (MQ=255)
                                               agCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtttt                                                                                                                                                     >  2:546062/1‑237 (MQ=255)
                                                   aaGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCaa                                                                                                                                               >  2:938854/1‑239 (MQ=255)
                                                   aaGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCa                                                                                                                                                >  2:563206/1‑238 (MQ=255)
                                                   aaGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCa                                                                                                                                                <  1:563206/238‑1 (MQ=255)
                                                    aGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAAt                                                                                                                                              <  2:879007/239‑1 (MQ=255)
                                                       aTGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCaa                                                                                                                                               <  1:1127669/235‑1 (MQ=255)
                                                       aTGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGc                                                                                                                                            <  1:901381/238‑1 (MQ=255)
                                                        tGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCtt                                                                                                                                          <  1:289204/239‑1 (MQ=255)
                                                                tAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAatat                                                                                                                                   >  1:1026676/1‑238 (MQ=255)
                                                                tAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAatat                                                                                                                                   >  2:979771/1‑238 (MQ=255)
                                                                    ccccAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCg                                                                                                                              <  1:36368/239‑1 (MQ=255)
                                                                             ttttGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGatat                                                                                                                     <  2:748509/239‑1 (MQ=255)
                                                                              tttGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCa                                                                                                                                                <  1:48011/211‑1 (MQ=255)
                                                                              tttGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCa                                                                                                                                                >  2:48011/1‑211 (MQ=255)
                                                                                 gTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTag                                                                                                                 <  1:346350/239‑1 (MQ=255)
                                                                                       aTCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCa                                                                                                           <  2:555082/239‑1 (MQ=255)
                                                                                          cTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAggg                                                                                                             <  1:712796/234‑1 (MQ=255)
                                                                                            cGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTa                                                                                                      <  1:849921/239‑1 (MQ=255)
                                                                                            cGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTa                                                                                                      >  1:269536/1‑239 (MQ=255)
                                                                                            cGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATACCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtt                                                                                                                                                       >  2:274028/1‑190 (MQ=255)
                                                                                            cGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATACCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCAtt                                                                                                                                                       <  1:274028/190‑1 (MQ=255)
                                                                                                  cAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGcc                                                                                                                                                                            <  1:948298/163‑1 (MQ=255)
                                                                                                  cAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGcc                                                                                                                                                                            >  2:948298/1‑163 (MQ=255)
                                                                                                     cTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGtttt                                                                                             >  1:680063/1‑239 (MQ=255)
                                                                                                       gTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGtttt                                                                                             >  1:522755/1‑237 (MQ=255)
                                                                                                            tGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCAt                                                                                      >  2:347582/1‑239 (MQ=255)
                                                                                                                tttttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCAt                                                                                      >  2:986373/1‑235 (MQ=255)
                                                                                                                  tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCt                                                                                                                                                                                           >  2:1059752/1‑132 (MQ=255)
                                                                                                                  tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCt                                                                                                                                                                                           <  1:1059752/132‑1 (MQ=255)
                                                                                                                  tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATtc                                                                                    >  2:1168083/1‑235 (MQ=255)
                                                                                                                  tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATtc                                                                                    >  2:1106531/1‑235 (MQ=255)
                                                                                                                  tttATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATtc                                                                                    >  2:327805/1‑235 (MQ=255)
                                                                                                                    tataATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAg                                                                              <  1:327805/239‑1 (MQ=255)
                                                                                                                             gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCg                                                                                                                              <  2:50686/182‑1 (MQ=255)
                                                                                                                             gTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCg                                                                                                                              >  1:50686/1‑182 (MQ=255)
                                                                                                                                    aTAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGt                                                              <  1:979771/239‑1 (MQ=255)
                                                                                                                                      aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAAt                                                                                                                                              <  1:495675/157‑1 (MQ=255)
                                                                                                                                      aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAAt                                                                                                                                              >  2:495675/1‑157 (MQ=255)
                                                                                                                                      aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCaaa                                                                                                                          >  1:544040/1‑177 (MQ=255)
                                                                                                                                      aGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCaaa                                                                                                                          <  2:544040/177‑1 (MQ=255)
                                                                                                                                         cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTg                                                             <  1:347582/235‑1 (MQ=255)
                                                                                                                                         cAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCa                                                           <  1:1168083/237‑1 (MQ=255)
                                                                                                                                            taatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCa                                                                                                                            >  2:696544/1‑169 (MQ=255)
                                                                                                                                            taatTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCa                                                                                                                            <  1:696544/169‑1 (MQ=255)
                                                                                                                                                  tAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGtac                                                <  1:546062/239‑3 (MQ=255)
                                                                                                                                                       ttAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGc                                           <  1:1106531/239‑1 (MQ=255)
                                                                                                                                                         aaTTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCaa                                                                                                                                               >  2:568017/1‑137 (MQ=255)
                                                                                                                                                         aaTTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCaa                                                                                                                                               <  1:568017/137‑1 (MQ=255)
                                                                                                                                                          aTTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAg                                        <  1:986373/239‑1 (MQ=255)
                                                                                                                                                              ctctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTAcc                                    <  2:522755/239‑1 (MQ=255)
                                                                                                                                                               tctAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCa                                   <  2:1026676/239‑1 (MQ=255)
                                                                                                                                                                   ataatTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCa                               >  1:410594/1‑239 (MQ=255)
                                                                                                                                                                             ctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTcc                                                                <  2:884184/196‑1 (MQ=255)
                                                                                                                                                                             ctctTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTcc                                                                >  1:884184/1‑196 (MQ=255)
                                                                                                                                                                                tttCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCa                                                                                                           <  1:1032570/150‑1 (MQ=255)
                                                                                                                                                                                tttCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCa                                                                                                           >  2:1032570/1‑150 (MQ=255)
                                                                                                                                                                                     aGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGt             >  2:584527/1‑239 (MQ=255)
                                                                                                                                                                                        tGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGcc                                                                                                                                                                            >  2:328920/1‑77 (MQ=255)
                                                                                                                                                                                        tGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGcc                                                                                                                                                                            <  1:328920/77‑1 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGa                                                                                                                                                                                                 <  2:176218/54‑1 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGa                                                                                                                                                                                                 >  1:176218/1‑54 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCa                                                                                                                            <  1:524830/123‑1 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCa                                                                                                                            >  2:524830/1‑123 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGc                                                                                                            >  2:813039/1‑139 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGc                                                                                                            <  1:813039/139‑1 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCa                                                                                                           >  2:352040/1‑140 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCa                                                                                                           <  1:352040/140‑1 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAg                                                                              <  2:725133/169‑1 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAg                                                                              >  1:725133/1‑169 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTAc                                     <  2:822342/210‑1 (MQ=255)
                                                                                                                                                                                          ccTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTAc                                     >  1:822342/1‑210 (MQ=255)
                                                                                                                                                                                                      tAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCaa                                                                                                                                               >  2:834282/1‑92 (MQ=255)
                                                                                                                                                                                                      tAGAAAAGTAACTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCaa                                                                                                                                               <  1:834282/92‑1 (MQ=255)
                                                                                                                                                                                                                aCTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtt  >  2:713665/1‑223 (MQ=255)
                                                                                                                                                                                                                aCTTCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtt  <  1:713665/223‑1 (MQ=255)
                                                                                                                                                                                                                  ttCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTcc                                                                <  1:985584/159‑1 (MQ=255)
                                                                                                                                                                                                                  ttCATCAAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTcc                                                                >  2:985584/1‑159 (MQ=255)
                                                                                                                                                                                                                     atcaAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCa                                                                                                                            >  1:16112/1‑96 (MQ=255)
                                                                                                                                                                                                                     atcaAATTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCa                                                                                                                            <  2:16112/96‑1 (MQ=255)
                                                                                                                                                                                                                        aaaTTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGc                                                                                                                                                                             <  2:795061/44‑1 (MQ=255)
                                                                                                                                                                                                                        aaaTTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGc                                                                                                                                                                             >  1:795061/1‑44 (MQ=255)
                                                                                                                                                                                                                        aaaTTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCaaa                                                                                                                          <  2:1056396/95‑1 (MQ=255)
                                                                                                                                                                                                                        aaaTTTACTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCaaa                                                                                                                          >  1:1056396/1‑95 (MQ=255)
                                                                                                                                                                                                                              aCTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCa                                                           <  1:300826/152‑1 (MQ=255)
                                                                                                                                                                                                                              aCTGCTTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCa                                                           >  2:300826/1‑152 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                            
TATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAATTCCTCTAATAATTTGCTCTCTTTCAAGTTGCCTTTTCGAGTCTAGAAAAGTAACTTCATCAAATTTACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTT  >  NC_000913/1197134‑1197564

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