Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,040,157 A→G 19.4% K227R (AAG→AGG)  yedY → membrane‑anchored, periplasmic TMAO, DMSO reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,040,1570AG19.4% 117.7 / 4.2 61K227R (AAG→AGG) yedYmembrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base A (26/22);  new base G (2/11);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 2.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGAT  >  NC_000913/2039920‑2040392
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
aaGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCaa                                                                                                                                                                                                                                              <  1:1157343/238‑1 (MQ=255)
       gTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGAGGGCGGTAAAAATCTTATGTCTTAGGATTCTGGTTCGACGGAGCAAGGAAGCCGCCCACACCTTTGAACGTA‑GGTGGTTATTGCTGGGCGGTa                                                                                                                                                                                                                                       >  2:461622/1‑238 (MQ=255)
          gATGGTGGTGCCGTGGGTTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTTTGTCGAAGGAGTGCGGCTCGCCGAAGCAAAGCATCCGCTCACACTGATGACCGTT‑GGGGTTTTTTGCAAAGCggttccc                                                                                                                                                                                                                                    >  1:553398/1‑232 (MQ=255)
                      gTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTATTTTGGAGGAGTGCGTCTCGGCGGAGGAAAGCATCCGCTCAAACTGATTACCGCA‑GGTGTTTTTGGCAAAGCGTTACCGCCACaaaaaagcg                                                                                                                                                                                                                       >  1:1352871/1‑234 (MQ=255)
                         gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATGTCGAAGGAATGCGGCTCGACGCCACAAAGCCTCCGCTCCCAATGATGAACCTA‑GGTGTTTTTGGGCAGGCGTTACCGCCACAAAAAGgcggg                                                                                                                                                                                                                     >  1:936689/1‑236 (MQ=255)
                            tGGTTTTCCGCTGCACAATTTGCTGGCGCTTGCCGAACCTCCCAGAAAGGCAAAGTATGTCGTTTTCGAAACAATTTATGCCCCGGAACAGATGCCAGGCAAGGAGGACCGTTTTATCGGGGGGGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCCCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTG‑GGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCgg                                                                                                                                                                                                                 <  2:553398/239‑1 (MQ=255)
                                 ttCCGCTGCACAAATGGCTGGCGCTGCCCGAACCAACCGGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCGGCAGGCCCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGa                                                                                                                                                                                                            <  1:1007656/239‑1 (MQ=255)
                                          acaAATTGCTGGCGCTTGCCGAACCTACCAGAAAGGCGAAGTATGTCGTTTGCGAAACAATTTATGCACCGGAACAGATCCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATtgtg                                                                                                                                                                                                   <  1:1429148/239‑1 (MQ=255)
                                          aaaatttGCTGGCGCTGCCCAACCCTCCCAGCAATGCGAAGTATGTCGTTTTCAAAAAAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGTTTTAGGGGGGGTGGGCTAAAATATCCTTGTGTCAAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTAGGGCAGGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATtgtg                                                                                                                                                                                                   <  2:983481/234‑1 (MQ=255)
                                          aaaaaTTGCTGCCCCTTGCCGAACCACCCAGCAAGGCGAAGGAGGTCGCTTTCGAAACAATTTATGCACCGGAAAAGATGCCAGGCCAGCAGGACCGCTTTACGGGCGGGGGGCTGAAATTTCCTTATGTCGAAGGATGGCGTCTCGACGAAGAAATGCATCCGCTCACACTGATGACCGTG‑GGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATtgtg                                                                                                                                                                                                   <  2:11378/237‑1 (MQ=255)
                                                   tGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATAACCTTTTGTCGAAGGGGTGCGTCTCCGCGAAACAAAGCAACCGCTCCCCCTGCTGACCGGAGGGGGTTTTTTGCAGGGGGTAACC‑CCCCAAAAAAGCGCCCCGGGGCGGCGGATtgtggccggggga                                                                                                                                                                                          >  2:616006/1‑230 (MQ=255)
                                                              gAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACGGGAACAGATCCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAGGGTTGGCGTCTCGACGAAGAAATGCACCCGCTCACACTGATGACGGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTaa                                                                                                                                                                              <  2:570277/239‑1 (MQ=255)
                                                              gAACCTACCAGAAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTA‑GGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTaa                                                                                                                                                                              <  1:616006/239‑1 (MQ=255)
                                                              aaaCCTACAAGCATGGAGAAGTATGTCGCTTTCAAAAAAATTAATGCACGGGAACAGATGCAGGGCCAGCAGACCCTCTTTACCGGGGGTGGGCTAAAATATCCTTATGTCAAAGGATTGCGTCTCGACGAAGAAATGCATCCGCTCACACGGAGGACCGTG‑GGTGTTTATGGCAGGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTaa                                                                                                                                                                              <  2:549915/238‑1 (MQ=255)
                                                                      cAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATCCCCCGGAAAAGATGCCAGGCCAGCAGGACCGTTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTG‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTa                                                                                                                                                                      <  1:944292/239‑1 (MQ=255)
                                                                            tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATTTCGAAGGATTTCGGCTCGCCGAAGCAAAGCATCCGCGCACACTGATGACCCCG‑GGGGGTTTTGGGAAAGGGGTAACCCCCCAAAAAAGGGGGCCGGGGGGACTGAGTTTGCCGGGGGAAAAATGGC‑TTAAAGGGGGTAAATAGa                                                                                                                                                                >  2:676972/1‑239 (MQ=255)
                                                                            tGCGAAGTATGTCGCTTTCAAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGTTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTG‑GGTGTTTATGGAAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGa                                                                                                                                                                <  1:1529007/239‑1 (MQ=255)
                                                                              cGAAGTATGTCGCTTTCGAAAAAATTTATGCCCCGAAACAGATGCCAGGCCAGCAGGCCCGTTTTATCGGCGGTGGGCTAAAATATCCTTAGGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGCCCGTA‑GGTGTTTAGGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGa                                                                                                                                                                <  1:213262/237‑1 (MQ=255)
                                                                               gAAGTATGCCGGTTCCGAAATAATTTATCCACCGAAACAGATGCCAGGGCAGCAGGACCGGTTTATGGGCGGGGGCCTGAAATTTCTTTATGCGGAAGGATGGCGTCTGGACGAAGAAATGCACCCGCTCACACTGAGGACCGTG‑GGTGTTTATGGAAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATa                                                                                                                                                              <  2:644045/238‑1 (MQ=255)
                                                                                aaGTGTGTGGCTTTAAAACAAATTGTTGCCCCGGAAAAGAGGCCAGGCCAGCAGGCCCCTTTTAGCGGCGGGGGGTTAAAAAATCCTTGTGTCGAGGGTTTGGGTCTCGAGGAAGAAATGCACCCGCTCACTCTGATGACCGTG‑GGTTTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                                            <  2:1359363/239‑1 (MQ=255)
                                                                                 aGTATGTCGTTTCCAAAACATTTTTGGCCCCGAAAAAGATGCCGGGCCAGCAGGACCGTTTTATCGGCGGTGGGTTAAAATATCCTTATGTCGAAGGAGTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTG‑GGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                                            <  2:686675/238‑1 (MQ=255)
                                                                                     tcgcccttcggAAACAATTTAGGCACGGGAACAGACGCAAGGCCCGAAGCTCCGTTTTATCGGGGGTGGGCTGAAATACCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAGTGCTCCCGCCCACACTGATGACGGTG‑GGTGTTTAGGGCAGGGTGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTa                                                                                                                                                       <  2:759338/229‑1 (MQ=255)
                                                                                           tttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATACTTTAGTCGAAGGAGTGCGTCTCGACGAAGCAATGCATTCGCCCACACTGATGAGCGTA‑GGTGGTTATTGCAAAGGGGTACAGCCCCACAAAGGGGCGCCGGGGCGGCTGGGTGTGTC‑GGGGGAATATGGCTTTTAAGGGGTTAAAAAGATAGACAGAATTaa                                                                                                                                                   >  2:448727/1‑237 (MQ=255)
                                                                                           tttCAAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCGGGACCGTTTTACCGGCGGTGGGTTGAAATATCCTTATGTCGAAGGATGGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGc                                                                                                                                                 <  1:1224109/239‑1 (MQ=255)
                                                                                             tCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGCAAAATTCTTATTTCGAAAGGTTTCCTTTCGACGGAGCAATGCCTCCGCCCACACTGATGACCGTT‑GGGGGTTTTTGCAAGGCGTTACCACCACCAAAACGCGCGCCGGGGGGGATGATTGTGCC‑GTGGGAATATGGGTTTAAAAGGAGTAAATAGCTAATCAGTGTTaaattt                                                                                                                                               >  2:764446/1‑235 (MQ=255)
                                                                                                      tttATGCCCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                                            <  1:543736/217‑1 (MQ=255)
                                                                                                               cGGAACAGATGCCAGGACAGCAGGTCCGCTTTATCGGCGGTGGGCTGAAATACCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACCCTGATGACCGTA‑GGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                                            <  1:739785/208‑1 (MQ=255)
                                                                                                                  aaCAGATGCCAGGCCAGCAGGACCGTTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTA‑GGTGTTTATGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAac                                                                                                                           <  1:1263253/238‑1 (MQ=255)
                                                                                                                               caggcagGACCGCTTTTCCGGGGGTGGGCTGAAATATCTTTATGCCGAAGGTTTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTgg                                                                                                             <  1:936154/236‑1 (MQ=255)
                                                                                                                                   cagGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGAAATGCATCCCCTCCCCCTGATGACCGCA‑GGGGTTTTTTGCAAGGCGGTACCGCCCCCAAAAGGCGCGCCCGGGCGGCCGCTTGTGCC‑GTGGGAATATTGCTTTTAAGGGGTTAAATCGATAGGTAGTATTAAACTGACCCGCGCGCGTCCCCCCACCACCCGGAATCttggtgc                                                                                                         >  2:480151/1‑233 (MQ=255)
                                                                                                                                         cgttttACCGGGGGGGGTCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACAAAGCAATGCACCCGCTCACACTGATGCCCGTG‑GGTGTTTATGAAAAGCGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGa                                                                                                   <  1:1076157/236‑1 (MQ=255)
                                                                                                                                                       tGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGtttt                                                                                      <  1:545712/238‑1 (MQ=255)
                                                                                                                                                                  atCCTTATGTGGAAGGATTGCTTCCCGACGAAGAAATGCATCCGCTCACACTGATGACGGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                                          <  1:461622/239‑1 (MQ=255)
                                                                                                                                                                  atCCTTATGTCGAAGGATGGGGCCTCGAAGAGGCAATGCACCCGCTCACCCGGATGGCGGTG‑GGTGTTTTGGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                                          <  2:936689/239‑1 (MQ=255)
                                                                                                                                                                  acccTTATTGCGAAGGATTGCGTCGCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                                          <  1:480151/237‑1 (MQ=255)
                                                                                                                                                                   tCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                                                         >  1:1060553/1‑239 (MQ=255)
                                                                                                                                                                   tCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                                                         >  2:37631/1‑239 (MQ=255)
                                                                                                                                                                   tCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                                                         >  2:1358750/1‑239 (MQ=255)
                                                                                                                                                                     cTTATGTCGAGGGATTGCGTCCCGACGAAGCAATGCACCCGCTCACACTGATGACCGGA‑GGTTTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAAt                                                                       <  2:239051/239‑1 (MQ=255)
                                                                                                                                                                     cTTATGTCGAAGGATTGCGTCCCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAAt                                                                       <  1:764446/239‑1 (MQ=255)
                                                                                                                                                                        aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                                                    >  2:820549/1‑239 (MQ=255)
                                                                                                                                                                           tCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTAt                                                                 >  1:543419/1‑239 (MQ=255)
                                                                                                                                                                             gAAGGCTTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGAGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGAGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCGGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGTAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGt                                                               >  2:572198/1‑239 (MQ=255)
                                                                                                                                                                             gAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGttttt                                                                                     >  1:857365/1‑217 (MQ=255)
                                                                                                                                                                              aaGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                                              >  1:99378/1‑239 (MQ=255)
                                                                                                                                                                                    tGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCa                                                         >  2:1425473/1‑238 (MQ=255)
                                                                                                                                                                                         ctggaggaAGCAATGCATCCGCTCACACGGATGACGGTG‑GGTGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCcgc                                                    <  2:329873/232‑1 (MQ=255)
                                                                                                                                                                                                    aaTGCATCCGCTCACAGTGATGACCGTA‑GGTGTTTATGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTcaca                                          <  1:1351161/237‑1 (MQ=255)
                                                                                                                                                                                                    aaTGCACCCGCCCACACTGATGCCGGTG‑GGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATc                                                                      <  2:785187/209‑1 (MQ=255)
                                                                                                                                                                                                           ccGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGa                                 >  1:635643/1‑239 (MQ=255)
                                                                                                                                                                                                              cTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGa                                                            >  1:785187/1‑209 (MQ=255)
                                                                                                                                                                                                                   acTGTTGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGATTGATTGTGTT‑GTGTAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGCCCGCCAACCACCTGGAATGTGTCAGCGCCTGACGAATACGGGTTTTACGCCAACGTTAATTCGTTTGCTGATCTACCGCGCTGGACACCGGCTTCAGTACGATTTa                         >  2:1178374/1‑239 (MQ=255)
                                                                                                                                                                                                                             gTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg                >  2:1056412/1‑238 (MQ=255)
                                                                                                                                                                                                                             gTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg                >  1:451451/1‑238 (MQ=255)
                                                                                                                                                                                                                                tgttttttAGGACAGGGCGTAACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAggc              <  2:101421/233‑1 (MQ=255)
                                                                                                                                                                                                                                  gtgtTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAggcggc           >  2:1349740/1‑239 (MQ=255)
                                                                                                                                                                                                                                  gtgtTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGAGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGAATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAggcggc           >  1:950138/1‑239 (MQ=255)
                                                                                                                                                                                                                                    ttttATGGAAGGAGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCAt         <  2:35808/238‑1 (MQ=255)
                                                                                                                                                                                                                                       tATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGGTAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTAACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATcc       >  2:248011/1‑238 (MQ=255)
                                                                                                                                                                                                                                         tGGAAAGGTGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCg    <  1:874395/239‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGa   >  2:1075482/1‑239 (MQ=255)
                                                                                                                                                                                                                                           gCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGGTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGAt  >  2:789088/1‑239 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
AAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGAT  >  NC_000913/2039920‑2040392

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