Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,040,0750CG42.4% 67.1 / 23.6 60L200V (CTG→GTG) yedYmembrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base C (8/26);  new base G (25/0);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 5.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGG  >  NC_000913/2039926‑2040223
                                                                                                                                                      |                                                                                                                                                     
ggTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCt                                                                                                                                                      <  2:540284/151‑1 (MQ=255)
    gATGGTGGTGCCGTGGATTGGTTTTCCGCTGCCCAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑Gaaa                                                                                                                                                 <  1:736412/151‑1 (MQ=255)
         tggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GGAata                                                                                                                                               >  2:8186/1‑148 (MQ=255)
         tggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGG‑Gaaa                                                                                                                                                 >  1:143164/1‑146 (MQ=255)
           gtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GAAATATcc                                                                                                                                            >  2:1201911/1‑149 (MQ=255)
           gtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCtt                                                                                                                                          <  1:976995/151‑1 (MQ=255)
           gtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCtt                                                                                                                                          <  1:888879/151‑1 (MQ=255)
                gTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGC‑GCAatattc                                                                                                                                            >  2:21540/1‑142 (MQ=255)
                  ggATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCg                                                                                                                                    <  2:179339/150‑1 (MQ=255)
                   gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGG‑GGGCG‑GGAAAATCCTTATGt                                                                                                                                      >  2:973542/1‑147 (MQ=255)
                         ttttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGCT‑GAATTTTCCttttttgtagggat                                                                                                                              >  1:361638/1‑138 (MQ=255)
                         ttttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GAAATATACTTTTGTCGAAgg                                                                                                                                >  1:859033/1‑147 (MQ=255)
                          tttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTg                                                                                                                            <  1:319802/150‑1 (MQ=255)
                                  gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTcgac                                                                                                                   <  1:1198459/151‑1 (MQ=255)
                                  gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTcgac                                                                                                                   <  1:319161/151‑1 (MQ=255)
                                  gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTcg                                                                                                                     <  1:430426/149‑1 (MQ=255)
                                     caAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GGAATATCCTTTTGTTGGAGGATTGCGTCTcggc                                                                                                                   >  2:67091/1‑146 (MQ=255)
                                      aaaTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCC‑GGAATATCCTTATGTCGAAGGAGTGCGTCTcgacga                                                                                                                 >  1:903312/1‑149 (MQ=255)
                                           gCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GGAAAAATCTTTTGTCCGAGGGTTGTGtgtc                                                                                                                      >  1:1135550/1‑139 (MQ=255)
                                                 gcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GGAATATCCTTTTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATcc                                                                                                     >  2:232410/1‑150 (MQ=255)
                                                  cTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GGAATATCCTTTTTTCGGAAGGTTTCGGCTCCGCGAAGCAAAGCATcc                                                                                                     >  2:876423/1‑149 (MQ=255)
                                                    tGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGc                                                                                                   <  2:678892/149‑1 (MQ=255)
                                                    tGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCt                                                                                                  <  1:819048/150‑1 (MQ=255)
                                                                  gCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GGAATATACTTTTTGCTGCGGATTTCGTCTCGCCCGAGCAATGCACCTGCTccccccgaggagc                                                                                     >  1:1148162/1‑137 (MQ=255)
                                                                      tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGG‑GGGGT‑CAAATAACCTTATCTTGTAGGGTTGCGCGCCGCCGCAACAAAGCATACACTCACCCTGCTGACCgc                                                                                   >  2:1143186/1‑146 (MQ=255)
                                                                       gCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGgt                                                                               <  1:232398/150‑1 (MQ=255)
                                                                        cGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GGAAAATACCTATGTCGGAGGATTTGGGCCCGACGAAAAAAAGCATCCCCCCCCCCTGATGAACGGAgggg                                                                              >  2:1000914/1‑148 (MQ=255)
                                                                         gAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GAAATATCCTTTTGTTGGAGGATTGCGGCTCGACGAAGCAAAGCATCCGCTCCCACTGATGACCGTAgggggt                                                                            >  2:411434/1‑147 (MQ=255)
                                                                          aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGttt                                                                           <  2:373734/151‑1 (MQ=255)
                                                                          aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGA‑GGGGTGGAAAAATCC‑TATTTTGAAGGATTGCGGCTCGACGAAACAAATCATCCGCCCACACTTGTGATCGGAgggggtt                                                                           >  2:166288/1‑146 (MQ=255)
                                                                             tATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GAGAAATTCTTTTGTCGAAGGCGTGCTTCTCCACGAAGGAAATCATCCCCTCTCACCGATGGCCGAAGGGGGTTAt                                                                         >  2:530819/1‑150 (MQ=255)
                                                                                  cGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCaa                                                                    <  2:159714/150‑1 (MQ=255)
                                                                                  cGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAg                                                                   <  2:1171525/151‑1 (MQ=255)
                                                                                      ttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCg                                                                <  2:193275/150‑1 (MQ=255)
                                                                                      ttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCg                                                                <  2:329494/150‑1 (MQ=255)
                                                                                        cGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAAAATCCTTATGTTGAAGGAGTTCGGCTCGAAGAAGCAAAGCATACCCCCACACTGATGAACGGAGGGGttttttggaaaggggt                                                              >  2:215262/1‑137 (MQ=255)
                                                                                              aaTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GAAATATCCTTTTGTCCGAGGAATGGGGTCCGACGAAGGAATACATTCCCCCCCACCGATGACCGGAGGGGGTTTTGGGAAAGGGTTACCGcc                                                        >  2:524338/1‑150 (MQ=255)
                                                                                               aTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTTCGTCTCGACGAAGCAAAGCATCCGCCCACACTGATGAACGGAGGGGTTTTTGGCAAGGCGGTACCGcccc                                                      >  2:1201122/1‑149 (MQ=255)
                                                                                                 ttATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GAAAAATCCTTTTTTCGAAGGGTTGCGTCTCGACGAAACAATGCATCCCCTCACCCTGATGACCGTAGGGGTTTTTTTCAAGGCGTTACCGcccc                                                      >  1:501655/1‑147 (MQ=255)
                                                                                                  tATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGG‑GGAATAACCTTATGTCGAAGGGTTGCGTCTCGGCGAAGCAAAGCATCCGCTCCCACTTATGACCGGAGGGGGTTTTGGCAAGGCGGTACCGcccc                                                      >  2:186909/1‑146 (MQ=255)
                                                                                                      cACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGC‑GGGGT‑GAAAGATACTTATGTCGAAGGAACGCGTCTCCATGAAGCAATGCCTCTCCTCACCCTGATAAACCTCAGGGTTTTTGGGAAGGGGTTAACGccccc                                                     >  1:785905/1‑142 (MQ=255)
                                                                                                          ggAACAGATGCCGGGCCAGCAGGCCCGTTTTATCGGCGGT‑GGGCT‑GAAATATCTTTATGTAAAAGGATTGCGTCTCGACGAAGCAATGCATCCCCTCACACGGATGCCGGTGGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGgcgc                                            <  1:468332/150‑1 (MQ=255)
                                                                                                             aCAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑TGGGT‑GAAAAATACCTTTGTCGCAGGAATGCGTCTCGGCGCAGCAAAACAACCGCTCACACCGATTACCCCCGGGGTTTTTGGGAAGGCGGTAACCCCACAAAAAGGCGCGcc                                         >  1:1001315/1‑150 (MQ=255)
                                                                                                                  tGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGGT‑GAAAAATTCTTTTGTTGAAGGGTTTGGTCTCGACGAAGCAAAGCATCCCCTCCCCCTGGTGAACGGAGGGGTTTTTGGGAAGGGGGTACCGCCCCAAAAAGGGGCGCCGGGGCg                                   >  2:83290/1‑151 (MQ=255)
                                                                                                                    ccaggccagCAGGACCGCTTTATCGGCGGT‑GGGGT‑GAAAAATCCTTATGTCGAAGGATTTCGTTTCGACGAAGCAATGCATCCGCCCCCACCGATGACCGTAGGGGGTTTTGGGAAAGCGTTACCGCCCCCAAAAGGCGCGCCGGGGCggc                                 >  2:36165/1‑149 (MQ=255)
                                                                                                                     caggccagCAGGACCGCTTTATCGGCGGC‑GGGGT‑GGAAAATCCTTTTGTCGAAGGGTTGCGTCTCGGCGAAGCAAAGCATCCCCTCCCACTGATGGCCGTAGGGGTTTTTTGGAAAGCGGTACCCCCACCAAAAGGCGGGCCGGTg                                     >  1:229036/1‑146 (MQ=255)
                                                                                                                       ggccggCAGGCCCGCTTTATCGGCGGT‑GGGCT‑GAAATTTCCTTATGTCGAGGGATTGCGTCCCGACAAAGCAATGAACCCGCTCCAACTGAGGACCGTAGGTTTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTg                               <  1:843494/150‑1 (MQ=255)
                                                                                                                        gccagCAGGACCGCTTTATCGGCGGT‑GGGGT‑TGAAAATTCTTATGTCGGAGGGGTTGGTCTCCGCGGAGCAAAGCATCCCCCCCCACTGATGGCCGGAGGGGGTTTTTGGAAGGCGGTTCCCCCCCCAAAAAGCGCGCCGGGGCGGCTg                               >  1:552970/1‑149 (MQ=255)
                                                                                                                          cagcagGACCGCTTTATCGGCGGT‑GGGCT‑GAAAAATACTTTTTTCGCAGGATTGCGTCTCGCCGGAGCAATGCCTCCGCCCACACTGATGACCGTAGGGGGTTTTTGCAAGGCGTTACCGCCCCAAAAAGGCGCGCCGGGGCGGCTGAtt                            >  1:1116240/1‑150 (MQ=255)
                                                                                                                               ggACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGcc                       <  2:214371/150‑1 (MQ=255)
                                                                                                                                gACCGCTTTATCGGCGGT‑GGGGC‑GGAATATCCTTATGTTTGAGGAATGGGTCTCGGCGAAGGAAAGGATCCCCTCACACTCATGACCGCAGGGGTTTTTTGCAAGGCGGTACCCCCACAAAAAAGCGGGCCGGGGGGACTGATtttgg                        >  1:1032947/1‑147 (MQ=255)
                                                                                                                                    gttttATCGGCGGT‑GGGCT‑AAAATATCCTTATGTCGAAGGATTGCGTCTCGAGGAAGCAATGCATCCGCTCACACTGATGACGGTGGGTGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGa                  <  1:160155/148‑1 (MQ=255)
                                                                                                                                      tttATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGaa                 <  2:514034/149‑1 (MQ=255)
                                                                                                                                       ttATCGGGGGT‑GGGCG‑AAAATACCCTTATGTCGAAGGATTGCGTCCCGACGAAGCAATGCATCCGCTCACACTGATGACCGTGGGTTTTTAGGGCAGGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAat               <  1:736/150‑1 (MQ=255)
                                                                                                                                          ccggcggT‑GGGTT‑GAAATATCTTTATGCCGAAGGATTGCGTCTGGAGGAAGAAAGGCTTCCGCGCACGCTGAGGCCCGTGGGTGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATg            <  1:423741/149‑1 (MQ=255)
                                                                                                                                           cggcggT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTTCGACTGATTGTGCCGTGGAAATATgg           >  2:272582/1‑150 (MQ=255)
                                                                                                                                             gggggggggCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACAAAGCAATGCATCCCCTCACACTGATGACCGTGGGTGTTTATGGAAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAata              <  2:676341/140‑1 (MQ=255)
                                                                                                                                             gcggT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCtt        >  1:725687/1‑151 (MQ=255)
                                                                                                                                               ggT‑GGGCT‑GAAATATCCTTATGCCGAAGGATTGCGTCTCGACGAAGCATTGCATCCGCTCACACGGAGGACCGTGGGTGTTTATGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTa      <  2:714326/151‑1 (MQ=255)
                                                                                                                                                    ggCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTAGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAgg  <  2:275644/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGT‑GGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGG  >  NC_000913/2039926‑2040223

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