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

AAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGAT  >  NC_000913/2039920‑2040392
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
aaGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCaa                                                                                                                                                                                                                                              <  1:1157343/238‑1 (MQ=255)
       gTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGAGGGCGGTAAAAATCTTATGTCTTAGGATTCTGGTTCGACGGAGCAAGGAAGCCGCCCACACCTTTGAACGTAGGTGGTTATTGCTGGGCGGTa                                                                                                                                                                                                                                       >  2:461622/1‑238 (MQ=255)
          gATGGTGGTGCCGTGGGTTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTTTGTCGAAGGAGTGCGGCTCGCCGAAGCAAAGCATCCGCTCACACTGATGACCGTTGGGGTTTTTTGCAAAGCggttccc                                                                                                                                                                                                                                    >  1:553398/1‑232 (MQ=255)
                      gTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTATTTTGGAGGAGTGCGTCTCGGCGGAGGAAAGCATCCGCTCAAACTGATTACCGCAGGTGTTTTTGGCAAAGCGTTACCGCCACaaaaaagcg                                                                                                                                                                                                                       >  1:1352871/1‑234 (MQ=255)
                         gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATGTCGAAGGAATGCGGCTCGACGCCACAAAGCCTCCGCTCCCAATGATGAACCTAGGTGTTTTTGGGCAGGCGTTACCGCCACAAAAAGgcggg                                                                                                                                                                                                                     >  1:936689/1‑236 (MQ=255)
                            tGGTTTTCCGCTGCACAATTTGCTGGCGCTTGCCGAACCTCCCAGAAAGGCAAAGTATGTCGTTTTCGAAACAATTTATGCCCCGGAACAGATGCCAGGCAAGGAGGACCGTTTTATCGGGGGGGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCCCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTGGGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCgg                                                                                                                                                                                                                 <  2:553398/239‑1 (MQ=255)
                                 ttCCGCTGCACAAATGGCTGGCGCTGCCCGAACCAACCGGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCGGCAGGCCCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGa                                                                                                                                                                                                            <  1:1007656/239‑1 (MQ=255)
                                          acaAATTGCTGGCGCTTGCCGAACCTACCAGAAAGGCGAAGTATGTCGTTTGCGAAACAATTTATGCACCGGAACAGATCCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATtgtg                                                                                                                                                                                                   <  1:1429148/239‑1 (MQ=255)
                                          aaaatttGCTGGCGCTGCCCAACCCTCCCAGCAATGCGAAGTATGTCGTTTTCAAAAAAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGTTTTAGGGGGGGTGGGCTAAAATATCCTTGTGTCAAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTAGGGCAGGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATtgtg                                                                                                                                                                                                   <  2:983481/234‑1 (MQ=255)
                                          aaaaaTTGCTGCCCCTTGCCGAACCACCCAGCAAGGCGAAGGAGGTCGCTTTCGAAACAATTTATGCACCGGAAAAGATGCCAGGCCAGCAGGACCGCTTTACGGGCGGGGGGCTGAAATTTCCTTATGTCGAAGGATGGCGTCTCGACGAAGAAATGCATCCGCTCACACTGATGACCGTGGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATtgtg                                                                                                                                                                                                   <  2:11378/237‑1 (MQ=255)
                                                   tGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATAACCTTTTGTCGAAGGGGTGCGTCTCCGCGAAACAAAGCAACCGCTCCCCCTGCTGACCGGAGGGGGTTTTTTGCAGGGGGTAACCCCCCAAAAAAGCGCCCCGGGGCGGCGGATtgtggccggggga                                                                                                                                                                                          >  2:616006/1‑230 (MQ=255)
                                                              gAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACGGGAACAGATCCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAGGGTTGGCGTCTCGACGAAGAAATGCACCCGCTCACACTGATGACGGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaa                                                                                                                                                                              <  2:570277/239‑1 (MQ=255)
                                                              gAACCTACCAGAAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTAGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaa                                                                                                                                                                              <  1:616006/239‑1 (MQ=255)
                                                              aaaCCTACAAGCATGGAGAAGTATGTCGCTTTCAAAAAAATTAATGCACGGGAACAGATGCAGGGCCAGCAGACCCTCTTTACCGGGGGTGGGCTAAAATATCCTTATGTCAAAGGATTGCGTCTCGACGAAGAAATGCATCCGCTCACACGGAGGACCGTGGGTGTTTATGGCAGGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaa                                                                                                                                                                              <  2:549915/238‑1 (MQ=255)
                                                                      cAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATCCCCCGGAAAAGATGCCAGGCCAGCAGGACCGTTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTGGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTa                                                                                                                                                                      <  1:944292/239‑1 (MQ=255)
                                                                            tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATTTCGAAGGATTTCGGCTCGCCGAAGCAAAGCATCCGCGCACACTGATGACCCCGGGGGGTTTTGGGAAAGGGGTAACCCCCCAAAAAAGGGGGCCGGGGGGACTGAGTTTGCCGGGGGAAAAATGGCTTAAAGGGGGTAAATAGa                                                                                                                                                                >  2:676972/1‑239 (MQ=255)
                                                                            tGCGAAGTATGTCGCTTTCAAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGTTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTGGGTGTTTATGGAAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGa                                                                                                                                                                <  1:1529007/239‑1 (MQ=255)
                                                                              cGAAGTATGTCGCTTTCGAAAAAATTTATGCCCCGAAACAGATGCCAGGCCAGCAGGCCCGTTTTATCGGCGGTGGGCTAAAATATCCTTAGGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGCCCGTAGGTGTTTAGGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGa                                                                                                                                                                <  1:213262/237‑1 (MQ=255)
                                                                               gAAGTATGCCGGTTCCGAAATAATTTATCCACCGAAACAGATGCCAGGGCAGCAGGACCGGTTTATGGGCGGGGGCCTGAAATTTCTTTATGCGGAAGGATGGCGTCTGGACGAAGAAATGCACCCGCTCACACTGAGGACCGTGGGTGTTTATGGAAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATa                                                                                                                                                              <  2:644045/238‑1 (MQ=255)
                                                                                aaGTGTGTGGCTTTAAAACAAATTGTTGCCCCGGAAAAGAGGCCAGGCCAGCAGGCCCCTTTTAGCGGCGGGGGGTTAAAAAATCCTTGTGTCGAGGGTTTGGGTCTCGAGGAAGAAATGCACCCGCTCACTCTGATGACCGTGGGTTTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                                            <  2:1359363/239‑1 (MQ=255)
                                                                                 aGTATGTCGTTTCCAAAACATTTTTGGCCCCGAAAAAGATGCCGGGCCAGCAGGACCGTTTTATCGGCGGTGGGTTAAAATATCCTTATGTCGAAGGAGTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTGGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                                            <  2:686675/238‑1 (MQ=255)
                                                                                     tcgcccttcggAAACAATTTAGGCACGGGAACAGACGCAAGGCCCGAAGCTCCGTTTTATCGGGGGTGGGCTGAAATACCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAGTGCTCCCGCCCACACTGATGACGGTGGGTGTTTAGGGCAGGGTGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTa                                                                                                                                                       <  2:759338/229‑1 (MQ=255)
                                                                                           tttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATACTTTAGTCGAAGGAGTGCGTCTCGACGAAGCAATGCATTCGCCCACACTGATGAGCGTAGGTGGTTATTGCAAAGGGGTACAGCCCCACAAAGGGGCGCCGGGGCGGCTGGGTGTGTCGGGGGAATATGGCTTTTAAGGGGTTAAAAAGATAGACAGAATTaa                                                                                                                                                   >  2:448727/1‑237 (MQ=255)
                                                                                           tttCAAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCGGGACCGTTTTACCGGCGGTGGGTTGAAATATCCTTATGTCGAAGGATGGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGc                                                                                                                                                 <  1:1224109/239‑1 (MQ=255)
                                                                                             tCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGCAAAATTCTTATTTCGAAAGGTTTCCTTTCGACGGAGCAATGCCTCCGCCCACACTGATGACCGTTGGGGGTTTTTGCAAGGCGTTACCACCACCAAAACGCGCGCCGGGGGGGATGATTGTGCCGTGGGAATATGGGTTTAAAAGGAGTAAATAGCTAATCAGTGTTaaattt                                                                                                                                               >  2:764446/1‑235 (MQ=255)
                                                                                                      tttATGCCCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                                            <  1:543736/217‑1 (MQ=255)
                                                                                                               cGGAACAGATGCCAGGACAGCAGGTCCGCTTTATCGGCGGTGGGCTGAAATACCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACCCTGATGACCGTAGGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                                            <  1:739785/208‑1 (MQ=255)
                                                                                                                  aaCAGATGCCAGGCCAGCAGGACCGTTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTATGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAac                                                                                                                           <  1:1263253/238‑1 (MQ=255)
                                                                                                                               caggcagGACCGCTTTTCCGGGGGTGGGCTGAAATATCTTTATGCCGAAGGTTTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTgg                                                                                                             <  1:936154/236‑1 (MQ=255)
                                                                                                                                   cagGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGAAATGCATCCCCTCCCCCTGATGACCGCAGGGGTTTTTTGCAAGGCGGTACCGCCCCCAAAAGGCGCGCCCGGGCGGCCGCTTGTGCCGTGGGAATATTGCTTTTAAGGGGTTAAATCGATAGGTAGTATTAAACTGACCCGCGCGCGTCCCCCCACCACCCGGAATCttggtgc                                                                                                         >  2:480151/1‑233 (MQ=255)
                                                                                                                                         cgttttACCGGGGGGGGTCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACAAAGCAATGCACCCGCTCACACTGATGCCCGTGGGTGTTTATGAAAAGCGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGa                                                                                                   <  1:1076157/236‑1 (MQ=255)
                                                                                                                                                       tGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGtttt                                                                                      <  1:545712/238‑1 (MQ=255)
                                                                                                                                                                  atCCTTATGTGGAAGGATTGCTTCCCGACGAAGAAATGCATCCGCTCACACTGATGACGGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                                          <  1:461622/239‑1 (MQ=255)
                                                                                                                                                                  atCCTTATGTCGAAGGATGGGGCCTCGAAGAGGCAATGCACCCGCTCACCCGGATGGCGGTGGGTGTTTTGGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                                          <  2:936689/239‑1 (MQ=255)
                                                                                                                                                                  acccTTATTGCGAAGGATTGCGTCGCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                                          <  1:480151/237‑1 (MQ=255)
                                                                                                                                                                   tCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                                                         >  1:1060553/1‑239 (MQ=255)
                                                                                                                                                                   tCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                                                         >  2:37631/1‑239 (MQ=255)
                                                                                                                                                                   tCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                                                         >  2:1358750/1‑239 (MQ=255)
                                                                                                                                                                     cTTATGTCGAGGGATTGCGTCCCGACGAAGCAATGCACCCGCTCACACTGATGACCGGAGGTTTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAAt                                                                       <  2:239051/239‑1 (MQ=255)
                                                                                                                                                                     cTTATGTCGAAGGATTGCGTCCCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAAt                                                                       <  1:764446/239‑1 (MQ=255)
                                                                                                                                                                        aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                                                    >  2:820549/1‑239 (MQ=255)
                                                                                                                                                                           tCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTAt                                                                 >  1:543419/1‑239 (MQ=255)
                                                                                                                                                                             gAAGGCTTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGAGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGAGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCGGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGTAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGt                                                               >  2:572198/1‑239 (MQ=255)
                                                                                                                                                                             gAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGttttt                                                                                     >  1:857365/1‑217 (MQ=255)
                                                                                                                                                                              aaGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                                              >  1:99378/1‑239 (MQ=255)
                                                                                                                                                                                    tGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCa                                                         >  2:1425473/1‑238 (MQ=255)
                                                                                                                                                                                         ctggaggaAGCAATGCATCCGCTCACACGGATGACGGTGGGTGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCcgc                                                    <  2:329873/232‑1 (MQ=255)
                                                                                                                                                                                                    aaTGCATCCGCTCACAGTGATGACCGTAGGTGTTTATGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTcaca                                          <  1:1351161/237‑1 (MQ=255)
                                                                                                                                                                                                    aaTGCACCCGCCCACACTGATGCCGGTGGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATc                                                                      <  2:785187/209‑1 (MQ=255)
                                                                                                                                                                                                           ccGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGa                                 >  1:635643/1‑239 (MQ=255)
                                                                                                                                                                                                              cTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGa                                                            >  1:785187/1‑209 (MQ=255)
                                                                                                                                                                                                                   acTGTTGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGATTGATTGTGTTGTGTAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGCCCGCCAACCACCTGGAATGTGTCAGCGCCTGACGAATACGGGTTTTACGCCAACGTTAATTCGTTTGCTGATCTACCGCGCTGGACACCGGCTTCAGTACGATTTa                         >  2:1178374/1‑239 (MQ=255)
                                                                                                                                                                                                                             gTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg                >  2:1056412/1‑238 (MQ=255)
                                                                                                                                                                                                                             gTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg                >  1:451451/1‑238 (MQ=255)
                                                                                                                                                                                                                                tgttttttAGGACAGGGCGTAACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAggc              <  2:101421/233‑1 (MQ=255)
                                                                                                                                                                                                                                  gtgtTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAggcggc           >  2:1349740/1‑239 (MQ=255)
                                                                                                                                                                                                                                  gtgtTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGAGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGAATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAggcggc           >  1:950138/1‑239 (MQ=255)
                                                                                                                                                                                                                                    ttttATGGAAGGAGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCAt         <  2:35808/238‑1 (MQ=255)
                                                                                                                                                                                                                                       tATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGGTAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTAACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATcc       >  2:248011/1‑238 (MQ=255)
                                                                                                                                                                                                                                         tGGAAAGGTGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCg    <  1:874395/239‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGa   >  2:1075482/1‑239 (MQ=255)
                                                                                                                                                                                                                                           gCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGGTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGAt  >  2:789088/1‑239 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
AAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGAT  >  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: