Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,079,4680AG46.0% 67.4 / 15.1 52D55G (GAT→GGT) putPproline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (1/26);  new base G (23/0);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 2.22e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAA  >  NC_000913/1079323‑1079618
                                                                                                                                                 |                                                                                                                                                                  
aTGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGa                                                                                                                                                              <  1:543174/151‑1 (MQ=255)
  gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGCTA‑TGAgg                                                                                                                                                            >  1:976478/1‑150 (MQ=255)
  gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGc                                                                                                                                                            <  1:548244/151‑1 (MQ=255)
  gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGc                                                                                                                                                            <  1:2792/151‑1 (MQ=255)
     ggTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑gg                                                                                                                                                         <  2:472076/150‑1 (MQ=255)
          cATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTGGGGTC‑TGAggcgcagg                                                                                                                                                      >  2:233160/1‑142 (MQ=255)
           aTTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTGGGGTA‑TGAgg‑gg                                                                                                                                                         >  2:168251/1‑144 (MQ=255)
               tgtgTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGA‑TGAGG‑gggtgggcggtt                                                                                                                                               >  1:63580/1‑144 (MQ=255)
                   tCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGTCGTGTCGGGTG‑TGGGG‑Ggctgggt‑ggtta                                                                                                                                             >  1:407328/1‑143 (MQ=255)
                      atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑ggg                                                                                                                                      <  1:189711/148‑1 (MQ=255)
                             tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCggg                                                                                                                             <  2:702812/150‑1 (MQ=255)
                              ggCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGgc‑g                                                                                                                          <  1:275159/151‑1 (MQ=255)
                              ggCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGgc‑g                                                                                                                          <  2:639766/151‑1 (MQ=255)
                                 aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTA‑TGAGG‑GGGTGGCT‑G‑TT‑TAT‑GGGGGTGGCggg                                                                                                                             >  1:83644/1‑146 (MQ=255)
                                 aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTA‑TGAGG‑GGCGGGGG‑G‑TT‑TAT‑GGGGGTGGCggg                                                                                                                             >  1:597910/1‑146 (MQ=255)
                                    atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGTA‑TGAGG‑GGCTGGGT‑G‑TT‑TCT‑GGGGTTGTCggggg                                                                                                                           >  2:946835/1‑143 (MQ=255)
                                    atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTT‑TGTGCGGGCGGGCG‑T‑TT‑ATT‑GGGGTTGcgg                                                                                                                               >  2:1073555/1‑140 (MQ=255)
                                      attgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑t‑tttt                                                                                                              <  1:252573/151‑1 (MQ=255)
                                        tgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑Ttttc                                                                                                             <  1:491831/150‑1 (MQ=255)
                                                tttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTC‑TGAGC‑GGCGGGGG‑G‑CT‑GAT‑GGGGGTGGCGGGCGG‑C‑GG‑T‑TTTTTTTTCCggg                                                                                                     >  2:336864/1‑151 (MQ=255)
                                                 ttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGa                                                                                                    <  2:71300/150‑1 (MQ=255)
                                                  tATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGA‑TGAGG‑GGGTGGGT‑G‑TT‑TAT‑GGGGGTGCCGGGG‑GCC‑GG‑T‑TTTTTTTTCCgggtt                                                                                                   >  1:490088/1‑149 (MQ=255)
                                                       ccTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGa                                                                                             <  2:53237/151‑1 (MQ=255)
                                                         tGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTA‑TGAGG‑GGGTGGGT‑G‑TT‑AATGGGGGTTGCGGGGG‑G‑G‑GG‑T‑TTTTTTTTCGGGGTTTTCGGa                                                                                             >  1:862618/1‑150 (MQ=255)
                                                         tGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGA‑TTAGC‑GGGGGGGG‑G‑TTGATT‑GGGGTTTCGGGGG‑C‑C‑GG‑T‑TTTTTTTTCCGGGTTTCCCGa                                                                                             >  2:508743/1‑150 (MQ=255)
                                                           gCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGc                                                                                         <  1:250255/151‑1 (MQ=255)
                                                             gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTA‑TTGGG‑GGGTGGGT‑G‑TT‑TAT‑GGGGGTGGCGGGG‑G‑G‑GG‑G‑TTTTTTTTTCGGGGTTTTCGTAAACt                                                                                        >  2:441775/1‑150 (MQ=255)
                                                             gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATTGGGG‑GGCGGGTG‑T‑TT‑ATT‑GGGGTGGCGGGGG‑G‑G‑GT‑T‑TTTTTTTTTCGGGTTTTCCaaaa                                                                                           >  2:662894/1‑148 (MQ=255)
                                                                cAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGa                                                                                     <  1:758345/150‑1 (MQ=255)
                                                                 aaCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGA‑TGAGG‑GGGTGGCT‑G‑TT‑TAT‑GGGGGTGTCGGGG‑G‑CGGG‑T‑TTTTTTTTCCGGGTTTTCGGaaaatgggtt                                                                                    >  1:1067943/1‑144 (MQ=255)
                                                                   cGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGCGCGTCGGCTA‑TGAGC‑GGGGGGCTGG‑TT‑AAG‑GGGGTT‑CCGGGG‑G‑C‑TG‑T‑TTTTTTTTTCGGGGATTTCGaaaagag                                                                                       >  2:1027580/1‑141 (MQ=255)
                                                                     aaaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTA‑TTAGC‑GGCTGGGT‑GTTT‑ATT‑GGGGTGTCCGGGG‑C‑C‑GG‑T‑TTTTTTTTTCGGGGTTT‑CGAAAGCTGGGTC‑Gcc                                                                               >  2:989832/1‑150 (MQ=255)
                                                                           tttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTA‑TGAGC‑GGGGGGGT‑G‑TT‑TAA‑GGGGGTGGCGGGC‑G‑C‑TGTT‑TTTTTTTTCGGGGATTTCGGAAA‑CTGGGTT‑GCCATTgggc                                                                        >  1:870565/1‑149 (MQ=255)
                                                                              gacgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGa                                                                     <  1:323977/151‑1 (MQ=255)
                                                                                  acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTA‑TTAGG‑GGCCGGGT‑G‑GT‑TAT‑TGGGGTGGCGGGG‑G‑G‑TG‑T‑TTTTTTTTTCGGGGTTTTCGAAAGCTTGGGC‑GCCCCTGGCCTGGCAt                                                                  >  1:440673/1‑150 (MQ=255)
                                                                                   cTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGATA‑TGAGG‑GGCGGGGT‑G‑TT‑TAT‑GGGGGTGCGCGGC‑G‑C‑GG‑T‑TTTTTTTTTCGGGGGTTCCCCAAGCTGGATC‑GCCCTTTGCCtctca                                                                   >  1:256354/1‑144 (MQ=255)
                                                                                     atatTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTA‑TGAGG‑GGGTGGGT‑G‑TT‑TAT‑GGGGGTGGCGGGG‑G‑G‑CG‑T‑TTTTTTTTTCGGGGGTTTCGAAAACTGGGTC‑GCCCTTTGCCTGGCATTatg                                                              >  1:239305/1‑149 (MQ=255)
                                                                                                       gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGaa                                            <  2:199318/151‑1 (MQ=255)
                                                                                                       gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTA‑TTGGC‑GGGCGGGT‑GTTT‑ATT‑GGGGTTGCGGGGG‑C‑C‑GT‑T‑TTTTTTTTCGGGGGTTTCCG‑AAGGTGGGTT‑GCCACTTGGCTGGCATTTGgcggggggggttaactg                                              >  2:777027/1‑135 (MQ=255)
                                                                                                        tCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGA‑TGAGG‑GGGTGGCT‑G‑GT‑TAT‑GGGGGTGGCGGGG‑G‑C‑CG‑G‑TTTTTTTTTCGGGGTTTTCGAAAGCTGGGTC‑GCCCTTGGCCTGGCATTTGGCGCGGGGGTTAAATGGaa                                            >  1:754833/1‑150 (MQ=255)
                                                                                                           tGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCt                                         <  2:84356/150‑1 (MQ=255)
                                                                                                               ccATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtggtgg                                    <  1:978890/151‑1 (MQ=255)
                                                                                                               ccATTCGTGACGGCATTATCGGCGGGGGCGACGGGTA‑TGAGG‑GGCGGGGG‑G‑TT‑TAT‑GGGGGGGCCGGGG‑G‑G‑TG‑TGTTTTTTTTCTGGGATTT‑TGGAAGATCGGGT‑GCCAGTTTCCCGCCAATAGGCGCGGGGGTTTACTGGAAGGtgg                                       >  1:353881/1‑148 (MQ=255)
                                                                                                                  ttCGTGACGGCATTATCGGCGGGGGGGTCGGGTA‑TGAGG‑GGCTGGGT‑G‑TT‑TAT‑GGGGGTGCCGGGC‑G‑C‑TG‑T‑TTTTTTTTTCGGGGTTTTCGAAAGCTGGGGCGGCCGTTGG‑CTGTCATTATGGGCGGGGGTTTACTGGGAGCtggggg                                    >  1:539196/1‑148 (MQ=255)
                                                                                                                    cGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCggg                               <  1:689155/151‑1 (MQ=255)
                                                                                                                          ggCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTg                         <  2:63573/151‑1 (MQ=255)
                                                                                                                                aTCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtgt                    <  2:13314/150‑1 (MQ=255)
                                                                                                                                      gggTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAcc             <  2:178762/151‑1 (MQ=255)
                                                                                                                                          gCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAAt         <  2:791828/151‑1 (MQ=255)
                                                                                                                                            gTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATac       <  2:1003606/151‑1 (MQ=255)
                                                                                                                                              cGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATgca      <  2:664652/150‑1 (MQ=255)
                                                                                                                                              cGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaa     <  1:650698/151‑1 (MQ=255)
                                                                                                                                                 ata‑tGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATAcaacaa  <  2:239243/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                                  
ATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATA‑TGAGC‑GGCTGGCT‑G‑TT‑AAT‑GGGGTTGCCGGGC‑G‑C‑TG‑T‑TTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAA  >  NC_000913/1079323‑1079618

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