Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 343,444 A→G 14.5% G187G (GGA→GGG yahK → broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913343,4440AG14.5% 149.7 / 2.7 60G187G (GGA→GGGyahKbroad specificity NADPH‑dependent aldehyde reductase, Zn‑containing
Reads supporting (aligned to +/- strand):  ref base A (2/48);  new base G (8/2);  total (10/50)
Fisher's exact test for biased strand distribution p-value = 7.38e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACG  >  NC_000913/343205‑343678
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
gtacccaTGCCCGGGTCCTGTAACCCGCCGACCCCGGCGGACCCGGCCATCCTCTGGGCGGCTACTCACAACAGGTCGTCGTGCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGGGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGa                                                                                                                                                                                                                                                 <  1:1288935/234‑1 (MQ=255)
      atagcggggcccTTTAGCTCCCCGCCCCCGGACGGCCAGGACCATACTCTGGGGTGTTACCCCAGACAGAACGTAGTTCATGGGGGATTGTGTCTGCGTATGCCGCACCCCTAAGGGCAGCTGGCGGCGGTGGCTCCTTGGTTGTGTGCAGGGATCACAACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACatat                                                                                                                                                                                                                                           <  1:82819/229‑1 (MQ=255)
      aTGACCGGCCCCTATACCCCGCCGACGCCGGCCGACCCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACatat                                                                                                                                                                                                                                           <  2:556724/239‑1 (MQ=255)
      aTGACCGGCCCCTATAACTCGCCGGCGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGGTCGTCGTTCATGAGCGATATGTTCTTCGTATTCGTCCCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACatat                                                                                                                                                                                                                                           <  2:407729/239‑1 (MQ=255)
      aTGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCCTCGGCGGGATGGGGCctat                                                                                                                                                                                                                                           >  2:1160197/1‑239 (MQ=255)
       tcgcaccgccccataACCCGCCGACGCGGGACGACCCGGGCCATACCCTGGGCGGGTTCTCCCAACCGACTGGCGTTCCTGAGGGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATg                                                                                                                                                                                                                                         <  1:1160197/227‑1 (MQ=255)
         aCCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCCGGGCACACCACGTATTCTCCGCTACGCGCCTGGCAGGCCGGGCCGGGCAAAAAGCAGGGGGTGGACGGCATCGGCGGGCCGGGGCagaatgg                                                                                                                                                                                                                                         >  1:1045407/1‑233 (MQ=255)
              cACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGAGGGCGTGGTCGGCATCGGCGGGCCGTGGCATTTggggggta                                                                                                                                                                                                                                  >  2:303304/1‑235 (MQ=255)
                      aCTCGCCGACGCCGGTCGAACCGGGCCATACGCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCCCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGcc                                                                                                                                                                                                                          <  2:664241/239‑1 (MQ=255)
                           ccGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCCCATGGGGAtttaacttggcccccc                                                                                                                                                                                                                     >  2:550764/1‑225 (MQ=255)
                              aCGCCGGACGAACGGGGTCATACTCTGGGCGGCTTCCCCCAACAGGGCGGCGGTCAGGGGCGATTTGTTCGGCGTTTTCGTCACCGGGAAGGGCAGCGGGCGGCGGTGGCTCCTTGGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGAt                                                                                                                                                                                                                <  1:1153440/239‑1 (MQ=255)
                                              gCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGACTGGGGCATACTGGGATTTAAGCGGCCCCCCCGCTCGGGGccccaggggtg                                                                                                                                                                                                >  2:727992/1‑229 (MQ=255)
                                                 ttACTTTGGGGGGATACTCACAACGGAGGGTGTTTCATGGGCGTTTTGTTCTGCGTATTCGTCACCCGTAAGAGCAGCTGGCGGGTGTGGCGCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCGACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCa                                                                                                                                                                                             <  1:303304/238‑1 (MQ=255)
                                                      ccggggcgttACTCACAACAGATCGGTGTTCATGAGCGTTATGTTCTGCGTATTCGTCCCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTAc                                                                                                                                                                                        <  2:711191/230‑1 (MQ=255)
                                                        gggCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTTGGGCAAATGGGGGTTAAAGTGGGCCACCGCGTTGGGGGcactgtgtgggcgttttaa                                                                                                                                                                                        >  2:457930/1‑222 (MQ=255)
                                                              cTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTTGGGACTATGGGGGTTAAGAGGGGCCCCCCGAGGGGGGGGCATGTGGTGGGGTTTATCCCTTCTg                                                                                                                                                                                >  2:1227402/1‑239 (MQ=255)
                                                                   cacaACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACatat                                                                                                                                                                                                                                           <  2:571296/178‑1 (MQ=255)
                                                                           aTCGTCGTTCATGAGCGATATGTCCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGGATATGGGGGTTAAGCTTGCCCACCCGATGGGGGGCAATATGGGGGGAGTTACCaactttgtgggaaaaaccga                                                                                                                                                                   >  2:956781/1‑220 (MQ=255)
                                                                                           gATTTGTTCGGCGTATTCGTCCCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCt                                                                                                                                                    <  2:229811/238‑1 (MQ=255)
                                                                                                 ttCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATgtggt                                                                                                                                                                                                 <  2:910213/187‑1 (MQ=255)
                                                                                                      gtttgttCGTCCCCCGCAAGAGCAGCGGGCGGGGGTGGTCCCTTTGGTGTGGGCAGGGATCACCACGTATTCGCCGCTACGTCCCTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATg                                                                                                                                         <  1:848461/234‑1 (MQ=255)
                                                                                                      gcgtttccgtgccccGCAAGGCCCGGTGGCGGCGTTGGTTCGTTTGTGGGGGCCAGGGATCACCACGTGGTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGAATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATg                                                                                                                                         <  1:550764/227‑1 (MQ=255)
                                                                                                      gCGTATTCGTCCCCCGCAAGGGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATg                                                                                                                                         <  2:661987/239‑1 (MQ=255)
                                                                                                           ttCGACCCCCGCAAGAGCAGCTGGCGTCGGTGGCTGCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACTATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttg                                                                                                                                  <  2:1037180/239‑1 (MQ=255)
                                                                                                                 taccGCAAGGGCAGGTGGCGGCGGTGGCTCCGTTGTTGTGTGCAGGGATCACCACGTATTCGCCGTGACTGCACGGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAATCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAAc                                                                                                                             <  2:136195/237‑1 (MQ=255)
                                                                                                                      cAAGGGCAGCTGGCGGGGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTGCCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACg                                                                                                                        <  2:804968/239‑1 (MQ=255)
                                                                                                                        tgagCGGCCGGCGGGGGTTGCTCCTTTGTTGTTGTCAGGGATCACCCCGGATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCGGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCa                                                                                                                      <  2:210447/238‑1 (MQ=255)
                                                                                                                        aggcccggtGGCGGCGGTGGCTCCTTAGTTGTGTGCAGGGATCACCCCGTATTCGCCGTTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCa                                                                                                                      <  2:379239/231‑1 (MQ=255)
                                                                                                                                gggcggcGGTGGCTCCTTTGGTGTGTGCAGGGATCCCCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGAt                                                                                                              <  2:353349/238‑1 (MQ=255)
                                                                                                                                 ggcggcGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGGTGGGGACAAAGGGGAGTAAACATGGCCACGCGGAGGGGGGCCCAGGGGGGGGCTTTTTCACATCTTGGGCCAAAACCGGAGCGGGAAAAACCCCGGGGGCGCAGGTAGGGTGTTAATCACGCAAAGCCggaa                                                                                                             >  1:1324274/1‑236 (MQ=255)
                                                                                                                                             tCCTTGGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCt                                                                                                  <  2:909256/238‑1 (MQ=255)
                                                                                                                                               cTTTTTTGTGTGCGGGGTTCCCCACGTATGCGGCGCTAGGTCACTGGCAGGGCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCAt                                                                                               <  1:510202/239‑1 (MQ=255)
                                                                                                                                                              ggATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGATTAGGCGGGCCCCCGCTCGGGCGGCACTTCGGGGGGCTGTTCCCACTTCTAGGGCAAACCGGGAGAGGCGAAATCCCCTGGGGCCGATCACGATCTTGACTCCACCAAAGCCCGTGAGATagcgggcgccctgaagaggagtgta                                                                                >  1:581606/1‑219 (MQ=255)
                                                                                                                                                                atccccacgtttgtgcgcgttcGTCACTGGCAGGCCGGGCCGGGTAAAAAAGGGGGCGGGGGCGGGATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGAtt                                                                               <  1:1227402/217‑1 (MQ=255)
                                                                                                                                                                 cccccgcGTTTTGGCCGGTAGGCCGCTGGCAGGACGGGTCGGGTAAAAAAGGGGGGGTGGGCGGCGTCGGCGGTCTGGGGCATATGGGGATTAAGCTGGCGCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAATGCCCTGGGTGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTc                                                                             <  1:727992/233‑1 (MQ=255)
                                                                                                                                                                  caccacGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACg                                                                                                                        <  2:188389/195‑1 (MQ=255)
                                                                                                                                                                    ccccGTTTTCGCCTGTCCGGCACTGGGGGGCGGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCGGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCAtt                                                                          <  1:954424/236‑1 (MQ=255)
                                                                                                                                                                    ccccGTATTCGCCGCTACGTCCCTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCAtt                                                                          <  2:974115/236‑1 (MQ=255)
                                                                                                                                                                    ccacGTTTTCGCCGGTAGGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCAtt                                                                          <  2:320377/239‑1 (MQ=255)
                                                                                                                                                                       cggttgccccgctccgctccTGGCGGGCCGGGCGGGGTAAAAAATGGGGCGGTGCTGGCATCGGCGGTCTGGTGCATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACCTGTGGTGGCATTTACCCCTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTg                                                                       <  1:457930/220‑1 (MQ=255)
                                                                                                                                                                          atgcgccggTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGaa                                                                     <  2:481771/235‑1 (MQ=255)
                                                                                                                                                                          atgcgccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCGGGGACTTAGTGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGaa                                                                     <  2:1278973/235‑1 (MQ=255)
                                                                                                                                                                                   tACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATg                                                                                                                                         <  2:575584/162‑1 (MQ=255)
                                                                                                                                                                                    tcGTCCCTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGGGGTCGGCATCGGCGGTCTGGGACATATGGGGATGAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTg                                                          <  2:1265255/238‑1 (MQ=255)
                                                                                                                                                                                    aCGGCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTg                                                          <  2:485615/239‑1 (MQ=255)
                                                                                                                                                                                     gggcccTGGGAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTgc                                                         <  2:1045407/234‑1 (MQ=255)
                                                                                                                                                                                     cggcACGGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTg                                                          <  2:559731/235‑1 (MQ=255)
                                                                                                                                                                                      gTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTgcg                                                        <  2:1065341/239‑1 (MQ=255)
                                                                                                                                                                                      gTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATAGGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTgcg                                                        <  2:162355/239‑1 (MQ=255)
                                                                                                                                                                                           gggcgggcgGTGCCGGGTAAAAAAGTGGGCGTGGTCGGATTCGGCGGTCTGGGAGCTATGTGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCcaca                                                   <  1:252249/238‑1 (MQ=255)
                                                                                                                                                                                              aaGGCCGGGCCGGGTAAAAAAGTGGGCGTGGGCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATa                                                 <  2:1324274/237‑1 (MQ=255)
                                                                                                                                                                                                   cgggccgggTAAAAAAGGGGGCGTGGTCGGCATCGGCGGTCTGGGACATAGGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcg                                           <  2:167692/239‑1 (MQ=255)
                                                                                                                                                                                                      ggcgggTAAAAAAGTGGGCGTGGTCGGCATCGGGGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgac                                         <  2:490784/236‑1 (MQ=255)
                                                                                                                                                                                                      gccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgac                                         <  2:1109743/238‑1 (MQ=255)
                                                                                                                                                                                                      gccgggTAAAAAAGGGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTcgac                                         <  2:1177692/238‑1 (MQ=255)
                                                                                                                                                                                                              aaaaaGGGGGCGTTGTCGGCAGCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTac                                 <  2:141958/238‑1 (MQ=255)
                                                                                                                                                                                                               aaaaGGGGGCTTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTacca                               <  2:285248/239‑1 (MQ=255)
                                                                                                                                                                                                                                     ggcggTCGGGGACATATGGGGATTAAGCTGGCCAACGCGATGGGGGCACATGTGGTGGCATTTACCCCTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTTTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCAc         <  2:834624/239‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTTTGGGAAGTTTGGGGATTAAGTGGGCCTACGCGTGGGGGGCAGATGGGGGTGCATTTACCCGTTCTGAGGCAAAACCGGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCa       <  1:1002360/238‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTCGGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCa       <  2:1133707/238‑1 (MQ=255)
                                                                                                                                                                                                                                             tGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACg  <  2:1158944/239‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
GATCACATGACCGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACG  >  NC_000913/343205‑343678

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