Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,040,129 A→C 21.3% T218P (ACA→CCA)  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,1290AC21.3% 109.9 / 3.8 59T218P (ACA→CCA) yedYmembrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base A (14/30);  new base C (7/7);  total (22/37)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTATTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAA  >  NC_000913/2039894‑2040362
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
gCGTATTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTcaca                                                                                                                                                                                                                                             <  1:93252/238‑1 (MQ=255)
 cGTATTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCGTTATCGGCGGTGGGCGGAAAAATCCTTTTGTCGAAGGGGTGTGTTTCGGCGGAGGAATAGATCCGCTCacacc                                                                                                                                                                                                                                           >  1:514380/1‑238 (MQ=255)
    atgtaacGTTTGCGCCGCGCGGCAGAGTGGGCGGTTGTGGTGCCGTGGATTGTTCATGCGGTGCACAAATTGCGGGCGCTTGCCGAACCTACGAGCAATGCGAAGTATGGCGCTTCCGAAACAAGTTAGGCACCGGAACAGATGCCGGGCCAGCAGGACCGCTTTATGGGCGGTGGGCTGAAATATCCTTATGTCGAGGGATTGCGTCTGGACGAAGCAATGCAGCCGATCACACtgca                                                                                                                                                                                                                                        <  2:831364/233‑3 (MQ=255)
    aTTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtgat                                                                                                                                                                                                                                        <  1:144214/239‑1 (MQ=255)
         tCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCg                                                                                                                                                                                                                                 <  1:163060/239‑1 (MQ=255)
         tCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCGGAAATATCCTTATTTCGGAGGATTGTGTTTCCGCGGAGGAATACAACCCCTCACACttatgaaag                                                                                                                                                                                                                                  >  1:63743/1‑236 (MQ=255)
             aTGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAGATACTTATGGCGAAGGGAGTCGTCTCGACGGAACCATGAAACCGCCAACACagatgagccttgg                                                                                                                                                                                                                              >  1:485206/1‑232 (MQ=255)
             aTGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAg                                                                                                                                                                                                                              <  1:322227/238‑1 (MQ=255)
             aTGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGGGGAATATCCTTTTGTCGAAGGAGTGCGTCTCCACGCCGGCATGCATCCCCTCCCACCGATGTGCGTAgg                                                                                                                                                                                                                             >  2:554658/1‑239 (MQ=255)
                   tGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTTTTTCTGAGGGGTGGGTCTCCGCGCAGCAAAGCATCCCCTCTCACTGATTACCgggtgggtttt                                                                                                                                                                                                                       >  2:695712/1‑229 (MQ=255)
                   tGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAAAATACTTATGTCGAAGGATTTCGTCTCGCCGAAGCAATGCATCCGCCCCCACTGTTGACCCTAGGGGGTTa                                                                                                                                                                                                                      >  1:885367/1‑239 (MQ=255)
                     cGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCCTTTTTCGCAGGATTTCGTCTCGACGAAACAAAGCATCCCCTCACACTGATGGCCCTAGGAGGTTAtt                                                                                                                                                                                                                    >  1:889083/1‑238 (MQ=255)
                        ggAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGc                                                                                                                                                                                                                  <  1:961168/238‑1 (MQ=255)
                         gAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCa                                                                                                                                                                                                                 <  2:427109/238‑1 (MQ=255)
                          aaGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCaa                                                                                                                                                                                                                <  1:1153249/238‑1 (MQ=255)
                                  tCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGCCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtgatgtgagcgg                                                                                                                                                                                                                               <  1:368676/217‑8 (MQ=255)
                                    gATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAATATCCTTATGTCGGAGGGTTGCGTCTCGACGAAGCAATGGATCCGCCCCCACTTCTTAACGTAGGGGTTTATGGCAAGGCGGTACCg                                                                                                                                                                                                      >  1:1304048/1‑238 (MQ=255)
                                            tgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAACATACTTATTTCGAAGGGTTTGGTCGCCACGAAGCAAGGCAACCCCCCACACTGCTGAGCCTAGGGGGTTTTGGGAAAGCGTTATAGCCACcaaaa                                                                                                                                                                                             >  2:826643/1‑238 (MQ=255)
                                                    gttGGTTTCCACCTGCACAAATTGCGGGCTCTGCCCAACCTTACAAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTGTGCACCGGAACAGAGCCCAGGCCAGAGGGACCGCTTTAGCGGCGGTGGGCTGAAATATCCTTATGGCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTATGGAAAGGTGTTACCGCCACAAAATGgcgcgc                                                                                                                                                                                      <  2:1351486/237‑1 (MQ=255)
                                                          tttCCGCGGCACAAATTGCTGGCGCTTGCCGACCCTACCAGAAAGGCGAAGTTGGTCGCTTCCGAAACAATTTATGCCCCGGAACAGATGCCAGGCCGGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATGGCGTCTCGACGAAGCAATGCACCCGCTCACACTGAGGACCGTGGGTGTTTAGGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCg                                                                                                                                                                               <  2:711761/239‑1 (MQ=255)
                                                                       aaTTGCTGGCGCTTGCCGAACCTACCGGAAATGCGAATTATGTCGTTTTCGAAACATTTTATGCCCGGAACCAGATGCCAGGCCAGCAGGCCCGCTTTACCGGGGGTGGGTTAAATTATCTTTATGTCGAGGGATTGCGTCTCGACGAAGCAAGCCATCCGCTCCCATGGATGCCCGTAGTTTTTTATGGCAAGCCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCg                                                                                                                                                                 <  2:1275795/239‑1 (MQ=255)
                                                                                cgctgcccGAACCTCCCAGAAATGCGAGGTATGGCGTTTTCAAAAAAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGCCCGTTTTATCGGGGGTGGGCTGAAATACCCTTATGTCGAAGGATGGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTAGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAata                                                                                                                                                         <  2:1304048/232‑1 (MQ=255)
                                                                                cgcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGGGGAAAAACCTTATGTCGAAGGAGTGTGGCTCGACGGAAGAATGCGTCCGCCCCCACTTATGACCGCAGGGGGTTTTTGCAAGGGGGTAACGCCCCCAAAAGGCGCGCCGGGGCggccgggtgtggcgcggggaaat                                                                                                                                                         >  2:772971/1‑217 (MQ=255)
                                                                                cgcTGGCCAACCCTACCAGAAAGGCAAAGTATGTCGTTTCGAAAAAATTTTAGGCCCCGGAACAGATGCCGGGCCAGAAGGCCCGCGTTATCGGGGGGGGGCTAAAATATCCTTATGTCAAAGGTTGGCGTCTCGACGAAGAAATGCATCCGCTCACCCTGATGACCGTAGGTGTTTATGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAata                                                                                                                                                         <  2:514380/238‑1 (MQ=255)
                                                                                               ccgcccaatggtatgTATCTCGATTTCGAAACAATTTAGGCACCGGAAAAGATGGCGGGCCAGCGGGACCCTTTTACCGGGGGGGGGTTGAAATATCTTTAGGGCGAAGGATGGCGTCTCGCCGAACCAATGCTTCCCCTCCCCCTGTGGCCCGTGGGGGTTTATGGCAGGTCGTGACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGAt                                                                                                                                          <  2:974946/225‑1 (MQ=255)
                                                                                               aCCAGTAATGGTAAGTAGGTCGCTTTAGAAAAAATTTAGGCACCGGAAAGGAGGCCAGGCGGGAGGTACCTCTTTATGGGTGGGGGGATGAAATATCCTTAGGTGGAAGGATGGCTCCTCGACGAAGAAAGGCATCGGCTCACATTGATGCCGGTAGGTTTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGAt                                                                                                                                          <  2:370213/239‑1 (MQ=255)
                                                                                                      cgcggAAGTAGTTCGTTTTCAAAAAAATTTATGCACCGGAACGTGTGCCGGGCAGGGAGGCCCGTTTGTCCGGCGGGGGTCTGAAAAATCTTTATGTCGAAGGTTTGCGTCTCGACGAAGAATTGCTTCCGCTCCCGCTGATGCCCGTGGGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCg                                                                                                                                   <  2:145301/235‑1 (MQ=255)
                                                                                                      agtggAAGTATGTCCCTTTCAAAACATTTTATGCACCGAAACAGATCCCAGGCGAGCAGGCCCTCTTTCGCGGCGGTGGGCTGAAATATCTTTATGTCGAAGGTTTGCGTCTCGACGAAGAAATGCATCCGCTCACACTGAGGACTGTAGGTTTTTATGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCg                                                                                                                                   <  2:615172/235‑1 (MQ=255)
                                                                                                      aTGCGAAGTATGCCGTTTTGAAAACAATTTATGCACGGGAACAGATGCCAGGCCAGCAGGACCGTTTTAGCGGGGGGGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCGCGACGAAGAAATGCATCCGCTCACACTGATGACCGTGGGTGTTTAGGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCg                                                                                                                                   <  2:1291526/239‑1 (MQ=255)
                                                                                                       tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGGAAAATCCTTTTTTTGGAGGGATGGGTTTCGGCGGAGGAATGCATCCCCTCCCCCCGGTGAGCGGAGGGGGTTTTGGGAAAGGGGTAACGCCCCCAAAAGGCGCGGCGGGGGGGCTGCTTTTTGCGGGGGAAAATTGCTTTTAAAGGGGTaaaaaga                                                                                                                                   >  2:394296/1‑236 (MQ=255)
                                                                                                       tGCGAAGGATGTCGCTTTCAAAACATTTTATGCACCGGAACAGATGCCAGGCCGGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTGGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGa                                                                                                                                  <  1:826643/239‑1 (MQ=255)
                                                                                                           aaGTATTGCGCTGTCTAAACAATTTATCCCCCGAAACAGATGCCGGCCCGGCAGGACCGTTTTATCGGTGGTGGGCTAAAATACCCTTATGGCGAAGGATTGCGTCTCGACGAGCCACTGCACCCGCTCACACGGAGGACCGTAGGTGTTTATGCAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  2:1269567/239‑1 (MQ=255)
                                                                                                            ggTATGTCGTGTTCGAAACAATTTTTGCACCGGAACAGATGCCAGGCCGGCGGGCCCGCTTTAGCGGCGGTGGGCTAAAATATCTTTATGTCGAGGGATTGCGTCTCGACAAAACAATGCACCCGCTCACACGGATGACTGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  2:917108/237‑1 (MQ=255)
                                                                                                            attATGCCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATCTCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  1:1305774/236‑1 (MQ=255)
                                                                                                            aggaTGTGGCTTGCAAAACAATGTATGCGCCGAACCAGATCCCGGGCCACGGGGACCGCTTTCCCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGTTTGCGTCTCGGCGAAGCAATGCACCCGCTCACACTGATGACGGTGGGTTTTTAGGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  2:349639/235‑1 (MQ=255)
                                                                                                            aGTATGTCGCTTGCGAAACAATTTATGCCCCGAAAAAGATGCCAGGCCAGGAGGCCCTCTTTATCGGCGGTGGGCTGAAATTTCTTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  1:502751/238‑1 (MQ=255)
                                                                                                            aGTATGGCGCTTTCAAACCAATTTATGCGCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTGTGGGCGGTGGGCTGAAATATCCTTATGGCGAAGGATGGCGTCTCGACGAAGAAATGCATCCGCGCCCACTGAGGACCTAAGGTTTTTAGGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  1:281454/238‑1 (MQ=255)
                                                                                                            aGTATGCCGTGTTGGAAACAATTTATGCCCCGAAACAGATCCCAGGCCAGCAGGCCCGCTTTATCGGCGGTGGGCTAAAATATCCTTATTTCGAAGGATTGCTTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGGAGGTTTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                              <  2:9506/238‑1 (MQ=255)
                                                                                                                tGTCGCTTTCAAAAAAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCGGAAATATCCTTATGGCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCCCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTa                                                                                                                         <  1:806743/239‑1 (MQ=255)
                                                                                                                      tttCGAAACAATTTATGCCCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGc                                                                                                                   <  1:868977/239‑1 (MQ=255)
                                                                                                                      ttagaaaaCAATTTATGCACCGAACAGGATGCCAGCCCACCAGAACCGCTTTATGGGGGGTGGGCTGAAATATCTTTATGTCGAAGGATGGGGTCTCGAGGAAGCAATGCACCCGCTCACACTGATGACCGTGGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGc                                                                                                                   <  2:29910/234‑1 (MQ=255)
                                                                                                                        tCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGCGGGAATAACCCTATGTCGAAGGACTTCGTCCCGCCGGAGGAAAGCCACCGCTCACACTGCTGACAGGTGGGGGTTTATGACAAGCGGTTACCGCACAAAAAAGCGCGGCGCGGCGCGCGGTTGTGGCGGTGGAATATTGGTTTCAAGGGGTTAAAACCAAAGTCTGCAtttaattg                                                                                                                 >  2:1127620/1‑233 (MQ=255)
                                                                                                                        acGAAACAATTTGTGGCCCGGATGAGATGCCTGGCCAGGCGGACGTGTTTGTAGGGGGTTGGCAGAAATATCTTTATGAGTAGGTATCGCGTCTCGCGAAAGCAAGCCCTCCCCGCCCAGTTATCCCCATAGGCGTGTATGCGAAGGCTTAACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTg                                                                                                                 <  2:53289/238‑1 (MQ=255)
                                                                                                                                          ccGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAATCCTTTTTTCGGAGGATTTCGTCTCCACGGAGCAAAGCATCCCCCCCCACTGTTGGACCGAAGGGTTTTTTGGAAAGGGGTACCCCCCCAAAAAAGGGCGCCGGGGGGGCTGGGTGTGGCGGGGGGAAATGGCTTTAAAGGGGATAAATAGATAAACAGGTTTAAAGTGACCcgcgcgggttccc                                                                                                 >  2:197477/1‑229 (MQ=255)
                                                                                                                                                aacgtagcccggcccagCAGGACTCTTTGAGCGGGGGTGGGCTAAAATATCCTTATGGCGAAGGATTGCGTCTCGACGAAGCAGTGCATCCGCTCACATGGATGACCGTAGGTGTTTATGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAac                                                                                             <  2:63743/225‑1 (MQ=255)
                                                                                                                                                     atgcagggccggCAGGACCGTTTTACCGGCGGTGGGCGAAAATATCCTTATGTCGAAGGATGGGTTCTCGACGAAGCAATGCACCCGCTCCCACGGTTGCCCGTAGGTTTTTATGGCAGGGCGTTCCCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAaccacc                                                                                         <  2:219254/231‑1 (MQ=255)
                                                                                                                                                                                   ggTGGGCTGGAAAATCCTTATGTCGAAGGGTTGCGTGTCGACGAAGCAATGCAACCGCTCCCCCTCCTGACCGTAGGGGTTTTTGGCCAAGCGGTACCGCCCCAAAAAAGCGCGCCGGGGCGGCTGTTTGTGCCGCGGGAATATGGCTTTATAGGGGTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGGCCCCCAACCCCCCGGGATCTCGCAGCGGCTGACGCAAACGGgtt                                                         >  2:595253/1‑239 (MQ=255)
                                                                                                                                                                                            aaaTATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGGGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCa                                                 >  2:397743/1‑238 (MQ=255)
                                                                                                                                                                                                ctCCTTATGGCGAAGGATTGCGTCTCGACGAAGCAGTGCATCCGCTCCCACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                            <  1:554658/238‑1 (MQ=255)
                                                                                                                                                                                                 tCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGAGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                           >  2:1025057/1‑239 (MQ=255)
                                                                                                                                                                                                   ctatggTCGAAGGATTGGGTCCCGACGAAGCAATCCACCCCCCCACACTGATGACGGTAGTTTTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAAt                                         <  1:394296/234‑1 (MQ=255)
                                                                                                                                                                                                      aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCCCCACAAAATGGCGCGCCGGTGCGCCTGATTGTGCCGTGGACATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCACGCGTCCGCCAACAACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                      >  2:324902/1‑239 (MQ=255)
                                                                                                                                                                                                            aaGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                >  2:844960/1‑239 (MQ=255)
                                                                                                                                                                                                            aaGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                >  1:1245871/1‑239 (MQ=255)
                                                                                                                                                                                                                  tGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCa                           >  1:292242/1‑238 (MQ=255)
                                                                                                                                                                                                                  tGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAAACACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACCTTAATCCGTATGTTGATCa                           >  2:368232/1‑238 (MQ=255)
                                                                                                                                                                                                                        gcgacgaGACAATGCACCCGCTCACCCTGAGGCCGGTGGGTTTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCgcgc                    <  2:889083/238‑1 (MQ=255)
                                                                                                                                                                                                                         cgacgaAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCcgc                      >  2:568578/1‑236 (MQ=255)
                                                                                                                                                                                                                                  aaTGCATCCGCTCACCGTGAGGACCGTAGGTTTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTcaca            <  1:814387/237‑1 (MQ=255)
                                                                                                                                                                                                                                         ccGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGa   >  2:585763/1‑239 (MQ=255)
                                                                                                                                                                                                                                          cGCTCACATGGAGGACGGTGGTGTTTTAGGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGaa  <  2:485206/239‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
GCGTATTTATCGTATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAA  >  NC_000913/2039894‑2040362

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