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

GGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGG  >  NC_000913/2039926‑2040223
                                                                                                                                                      |                                                                                                                                                     
ggTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCt                                                                                                                                                      <  2:540284/151‑1 (MQ=255)
    gATGGTGGTGCCGTGGATTGGTTTTCCGCTGCCCAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGaaa                                                                                                                                                 <  1:736412/151‑1 (MQ=255)
         tggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAata                                                                                                                                               >  2:8186/1‑148 (MQ=255)
         tggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGaaa                                                                                                                                                 >  1:143164/1‑146 (MQ=255)
           gtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATcc                                                                                                                                            >  2:1201911/1‑149 (MQ=255)
           gtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCtt                                                                                                                                          <  1:976995/151‑1 (MQ=255)
           gtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCtt                                                                                                                                          <  1:888879/151‑1 (MQ=255)
                gTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGCGCAatattc                                                                                                                                            >  2:21540/1‑142 (MQ=255)
                  ggATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCg                                                                                                                                    <  2:179339/150‑1 (MQ=255)
                   gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGCGGGAAAATCCTTATGt                                                                                                                                      >  2:973542/1‑147 (MQ=255)
                         ttttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGCTGAATTTTCCttttttgtagggat                                                                                                                              >  1:361638/1‑138 (MQ=255)
                         ttttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATACTTTTGTCGAAgg                                                                                                                                >  1:859033/1‑147 (MQ=255)
                          tttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTg                                                                                                                            <  1:319802/150‑1 (MQ=255)
                                  gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTcgac                                                                                                                   <  1:1198459/151‑1 (MQ=255)
                                  gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTcgac                                                                                                                   <  1:319161/151‑1 (MQ=255)
                                  gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTcg                                                                                                                     <  1:430426/149‑1 (MQ=255)
                                     caAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTTTGTTGGAGGATTGCGTCTcggc                                                                                                                   >  2:67091/1‑146 (MQ=255)
                                      aaaTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCCGGAATATCCTTATGTCGAAGGAGTGCGTCTcgacga                                                                                                                 >  1:903312/1‑149 (MQ=255)
                                           gCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAAATCTTTTGTCCGAGGGTTGTGtgtc                                                                                                                      >  1:1135550/1‑139 (MQ=255)
                                                 gcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTTTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATcc                                                                                                     >  2:232410/1‑150 (MQ=255)
                                                  cTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTTTTTCGGAAGGTTTCGGCTCCGCGAAGCAAAGCATcc                                                                                                     >  2:876423/1‑149 (MQ=255)
                                                    tGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGc                                                                                                   <  2:678892/149‑1 (MQ=255)
                                                    tGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCt                                                                                                  <  1:819048/150‑1 (MQ=255)
                                                                  gCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATACTTTTTGCTGCGGATTTCGTCTCGCCCGAGCAATGCACCTGCTccccccgaggagc                                                                                     >  1:1148162/1‑137 (MQ=255)
                                                                      tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTCAAATAACCTTATCTTGTAGGGTTGCGCGCCGCCGCAACAAAGCATACACTCACCCTGCTGACCgc                                                                                   >  2:1143186/1‑146 (MQ=255)
                                                                       gCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGgt                                                                               <  1:232398/150‑1 (MQ=255)
                                                                        cGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATACCTATGTCGGAGGATTTGGGCCCGACGAAAAAAAGCATCCCCCCCCCCTGATGAACGGAgggg                                                                              >  2:1000914/1‑148 (MQ=255)
                                                                         gAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTTTGTTGGAGGATTGCGGCTCGACGAAGCAAAGCATCCGCTCCCACTGATGACCGTAgggggt                                                                            >  2:411434/1‑147 (MQ=255)
                                                                          aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGttt                                                                           <  2:373734/151‑1 (MQ=255)
                                                                          aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGAGGGGTGGAAAAATCCTATTTTGAAGGATTGCGGCTCGACGAAACAAATCATCCGCCCACACTTGTGATCGGAgggggtt                                                                           >  2:166288/1‑146 (MQ=255)
                                                                             tATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAGAAATTCTTTTGTCGAAGGCGTGCTTCTCCACGAAGGAAATCATCCCCTCTCACCGATGGCCGAAGGGGGTTAt                                                                         >  2:530819/1‑150 (MQ=255)
                                                                                  cGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCaa                                                                    <  2:159714/150‑1 (MQ=255)
                                                                                  cGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAg                                                                   <  2:1171525/151‑1 (MQ=255)
                                                                                      ttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCg                                                                <  2:193275/150‑1 (MQ=255)
                                                                                      ttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCg                                                                <  2:329494/150‑1 (MQ=255)
                                                                                        cGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAAAATCCTTATGTTGAAGGAGTTCGGCTCGAAGAAGCAAAGCATACCCCCACACTGATGAACGGAGGGGttttttggaaaggggt                                                              >  2:215262/1‑137 (MQ=255)
                                                                                              aaTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTTTGTCCGAGGAATGGGGTCCGACGAAGGAATACATTCCCCCCCACCGATGACCGGAGGGGGTTTTGGGAAAGGGTTACCGcc                                                        >  2:524338/1‑150 (MQ=255)
                                                                                               aTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTTCGTCTCGACGAAGCAAAGCATCCGCCCACACTGATGAACGGAGGGGTTTTTGGCAAGGCGGTACCGcccc                                                      >  2:1201122/1‑149 (MQ=255)
                                                                                                 ttATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTTTCGAAGGGTTGCGTCTCGACGAAACAATGCATCCCCTCACCCTGATGACCGTAGGGGTTTTTTTCAAGGCGTTACCGcccc                                                      >  1:501655/1‑147 (MQ=255)
                                                                                                  tATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGGAATAACCTTATGTCGAAGGGTTGCGTCTCGGCGAAGCAAAGCATCCGCTCCCACTTATGACCGGAGGGGGTTTTGGCAAGGCGGTACCGcccc                                                      >  2:186909/1‑146 (MQ=255)
                                                                                                      cACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTGAAAGATACTTATGTCGAAGGAACGCGTCTCCATGAAGCAATGCCTCTCCTCACCCTGATAAACCTCAGGGTTTTTGGGAAGGGGTTAACGccccc                                                     >  1:785905/1‑142 (MQ=255)
                                                                                                          ggAACAGATGCCGGGCCAGCAGGCCCGTTTTATCGGCGGTGGGCTGAAATATCTTTATGTAAAAGGATTGCGTCTCGACGAAGCAATGCATCCCCTCACACGGATGCCGGTGGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGgcgc                                            <  1:468332/150‑1 (MQ=255)
                                                                                                             aCAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTTGGGTGAAAAATACCTTTGTCGCAGGAATGCGTCTCGGCGCAGCAAAACAACCGCTCACACCGATTACCCCCGGGGTTTTTGGGAAGGCGGTAACCCCACAAAAAGGCGCGcc                                         >  1:1001315/1‑150 (MQ=255)
                                                                                                                  tGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATTCTTTTGTTGAAGGGTTTGGTCTCGACGAAGCAAAGCATCCCCTCCCCCTGGTGAACGGAGGGGTTTTTGGGAAGGGGGTACCGCCCCAAAAAGGGGCGCCGGGGCg                                   >  2:83290/1‑151 (MQ=255)
                                                                                                                    ccaggccagCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTCGAAGGATTTCGTTTCGACGAAGCAATGCATCCGCCCCCACCGATGACCGTAGGGGGTTTTGGGAAAGCGTTACCGCCCCCAAAAGGCGCGCCGGGGCggc                                 >  2:36165/1‑149 (MQ=255)
                                                                                                                     caggccagCAGGACCGCTTTATCGGCGGCGGGGTGGAAAATCCTTTTGTCGAAGGGTTGCGTCTCGGCGAAGCAAAGCATCCCCTCCCACTGATGGCCGTAGGGGTTTTTTGGAAAGCGGTACCCCCACCAAAAGGCGGGCCGGTg                                     >  1:229036/1‑146 (MQ=255)
                                                                                                                       ggccggCAGGCCCGCTTTATCGGCGGTGGGCTGAAATTTCCTTATGTCGAGGGATTGCGTCCCGACAAAGCAATGAACCCGCTCCAACTGAGGACCGTAGGTTTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTg                               <  1:843494/150‑1 (MQ=255)
                                                                                                                        gccagCAGGACCGCTTTATCGGCGGTGGGGTTGAAAATTCTTATGTCGGAGGGGTTGGTCTCCGCGGAGCAAAGCATCCCCCCCCACTGATGGCCGGAGGGGGTTTTTGGAAGGCGGTTCCCCCCCCAAAAAGCGCGCCGGGGCGGCTg                               >  1:552970/1‑149 (MQ=255)
                                                                                                                          cagcagGACCGCTTTATCGGCGGTGGGCTGAAAAATACTTTTTTCGCAGGATTGCGTCTCGCCGGAGCAATGCCTCCGCCCACACTGATGACCGTAGGGGGTTTTTGCAAGGCGTTACCGCCCCAAAAAGGCGCGCCGGGGCGGCTGAtt                            >  1:1116240/1‑150 (MQ=255)
                                                                                                                               ggACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGcc                       <  2:214371/150‑1 (MQ=255)
                                                                                                                                gACCGCTTTATCGGCGGTGGGGCGGAATATCCTTATGTTTGAGGAATGGGTCTCGGCGAAGGAAAGGATCCCCTCACACTCATGACCGCAGGGGTTTTTTGCAAGGCGGTACCCCCACAAAAAAGCGGGCCGGGGGGACTGATtttgg                        >  1:1032947/1‑147 (MQ=255)
                                                                                                                                    gttttATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGAGGAAGCAATGCATCCGCTCACACTGATGACGGTGGGTGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGa                  <  1:160155/148‑1 (MQ=255)
                                                                                                                                      tttATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGaa                 <  2:514034/149‑1 (MQ=255)
                                                                                                                                       ttATCGGGGGTGGGCGAAAATACCCTTATGTCGAAGGATTGCGTCCCGACGAAGCAATGCATCCGCTCACACTGATGACCGTGGGTTTTTAGGGCAGGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAat               <  1:736/150‑1 (MQ=255)
                                                                                                                                          ccggcggTGGGTTGAAATATCTTTATGCCGAAGGATTGCGTCTGGAGGAAGAAAGGCTTCCGCGCACGCTGAGGCCCGTGGGTGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATg            <  1:423741/149‑1 (MQ=255)
                                                                                                                                           cggcggTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTTCGACTGATTGTGCCGTGGAAATATgg           >  2:272582/1‑150 (MQ=255)
                                                                                                                                             gggggggggCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACAAAGCAATGCATCCCCTCACACTGATGACCGTGGGTGTTTATGGAAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAata              <  2:676341/140‑1 (MQ=255)
                                                                                                                                             gcggTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCtt        >  1:725687/1‑151 (MQ=255)
                                                                                                                                               ggTGGGCTGAAATATCCTTATGCCGAAGGATTGCGTCTCGACGAAGCATTGCATCCGCTCACACGGAGGACCGTGGGTGTTTATGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTa      <  2:714326/151‑1 (MQ=255)
                                                                                                                                                    ggCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTAGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAgg  <  2:275644/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGG  >  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: