Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,079,4600TG46.6% 83.0 / 12.6 58G52G (GGT→GGG) putPproline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base T (4/27);  new base G (27/0);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.23e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACC  >  NC_000913/1079314‑1079607
                                                                                                                                                  |                                                                                                                                                               
aGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGt                                                                                                                                                             <  1:974394/151‑1 (MQ=255)
           gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGc                                                                                                                                                 <  1:963374/151‑1 (MQ=255)
               gTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCggc‑t                                                                                                                                            <  1:90378/151‑1 (MQ=255)
                tGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGG‑CGGGTATgggggggcg                                                                                                                                            >  1:1107497/1‑143 (MQ=255)
                   cATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCggc‑tggct                                                                                                                                        <  1:584466/151‑1 (MQ=255)
                      tttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGT‑GGGGGATGAGGggg‑tg                                                                                                                                           >  1:930331/1‑145 (MQ=255)
                          tgtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGT‑CGGGGATGAGCggg‑tggccggttat                                                                                                                                  >  1:497496/1‑143 (MQ=255)
                           gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATg                                                                                                                              <  2:234000/150‑1 (MQ=255)
                           gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGT‑CGGGTATGAGCGgg‑tgggt‑ggttat                                                                                                                                 >  2:888633/1‑144 (MQ=255)
                           gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGG‑CGGGTATGAGCgggcgggggtttat                                                                                                                                   >  2:889976/1‑137 (MQ=255)
                           gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCCTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATg                                                                                                                              <  1:793416/150‑1 (MQ=255)
                               atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGtt                                                                                                                         <  2:49755/151‑1 (MQ=255)
                               atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGt                                                                                                                          <  2:442265/150‑1 (MQ=255)
                                 atCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑g                                                                                                                       <  1:725598/150‑1 (MQ=255)
                                   cTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTACGGA‑CTGCGAGGG‑TGGGT‑G‑TTT‑ATgggg                                                                                                                           >  2:598084/1‑145 (MQ=255)
                                   cTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGT‑CGGGTATTAGGGGG‑GGGGT‑G‑TTATTTGGGGttcc                                                                                                                       >  1:650925/1‑148 (MQ=255)
                                   cTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGC‑GGGGTATGAGGggc‑tggccggtt                                                                                                                                    >  1:414041/1‑134 (MQ=255)
                                    tttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGG‑CGGGTATGTGGGGG‑GGGGT‑G‑TTA‑TTGTGGGT‑GTCgg                                                                                                                   >  2:688860/1‑151 (MQ=255)
                                     ttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCg                                                                                                                    <  1:471134/149‑1 (MQ=255)
                                     ttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCg                                                                                                                    <  1:7570/149‑1 (MQ=255)
                                      tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCggg                                                                                                                  <  2:461848/150‑1 (MQ=255)
                                      tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGT‑CGGGGATGTGCGGG‑TGGCTGG‑TTA‑TTGGGGTG‑GCgggg                                                                                                                  >  2:288397/1‑151 (MQ=255)
                                          aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGG‑CGGGTTTGAGGGGG‑CGGCG‑G‑GTT‑ATTGGGGT‑GGCggg                                                                                                                  >  1:531906/1‑146 (MQ=255)
                                            gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑tt                                                                                                      <  1:1121321/151‑1 (MQ=255)
                                             atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑ttt                                                                                                     <  1:568466/151‑1 (MQ=255)
                                             atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGG‑CGGGTATGAGGGGG‑GGGGG‑G‑GTT‑ATGGGGGT‑GGCGGGGG‑C‑‑T‑G‑ttt                                                                                                     >  2:33384/1‑151 (MQ=255)
                                               attgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑tttt                                                                                                    <  2:593113/150‑1 (MQ=255)
                                                  gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTtttctt                                                                                                <  1:684216/151‑1 (MQ=255)
                                                  gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGG‑CGGGTATTGGGGGG‑CGGGT‑G‑TTT‑ATTGGGGT‑GGCGGGGG‑C‑‑T‑G‑tttttttt                                                                                                >  2:224778/1‑151 (MQ=255)
                                                          ttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGG‑CGGCAATTAGCGGG‑CGGGC‑G‑TTG‑ATGGGGGT‑GTCGGGGG‑C‑‑C‑G‑Tttttttttttcggg                                                                                         >  1:638848/1‑144 (MQ=255)
                                                              cGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGG‑CGGGTATGAGGGGC‑TGGGG‑G‑TTT‑ATGGGGGT‑GGCGGGGG‑G‑‑TGG‑TTTTTTTTTCCGGGttt                                                                                       >  2:235561/1‑149 (MQ=255)
                                                              cGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTcc                                                                                    <  2:467758/151‑1 (MQ=255)
                                                              cGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTCGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTcc                                                                                    <  1:1014418/151‑1 (MQ=255)
                                                               gCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGG‑GGGGGATGGGGGGG‑GGGGT‑G‑TTT‑ATGGGGGG‑GGCGGGCG‑CGGT‑G‑TTTTTTTTTCggggt                                                                                         >  1:68829/1‑146 (MQ=255)
                                                                ccTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGT‑CGGGGATGAGCGGG‑TGGGT‑G‑GTT‑ATGGGGGG‑GGCGGGGG‑C‑‑G‑G‑Ttttttttttttggggtt                                                                                      >  1:367002/1‑138 (MQ=255)
                                                                ccTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGG‑CGGGTATGATCGGG‑TGGCT‑G‑TTT‑ATTGGGTTGGCGGGGGG‑C‑‑G‑G‑TTTTTTTTTCGGGGTTTTcg                                                                                    >  2:964534/1‑149 (MQ=255)
                                                                 cTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGG‑CGGGTATGGGGGGG‑GGGGT‑GTTTA‑TTGGGGTT‑G‑CGGGCG‑C‑‑T‑G‑TTTTTTTTTTCCGGGttt                                                                                      >  2:1085325/1‑146 (MQ=255)
                                                                     cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGG‑CGGGTATGAGGGGG‑TGGGTGG‑TTA‑TTGGGGTT‑GCCGGGGG‑T‑‑G‑G‑TTTTTTTTTCCGGGTTTTCCaaa                                                                                 >  2:592671/1‑148 (MQ=255)
                                                                       aTCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGG‑CGGGTATTAGCGGC‑GGGGT‑G‑TTT‑ATGGGGGT‑GGCGGGGG‑G‑‑C‑G‑TTTTTTTTTTCTGGGTTTTCGaaaaag                                                                             >  1:303964/1‑146 (MQ=255)
                                                                            cGAGAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGG‑CCTGTTTGGGCGGC‑CGTGC‑G‑TTG‑ATGGGGGT‑GCGAGGGG‑C‑‑TGG‑GTTTTTTTTCTGGGATTT‑CGGAAGCTGGTTC‑g                                                                      >  1:518805/1‑150 (MQ=255)
                                                                             gAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGG‑CGGGGATAAGGGGG‑TGGGG‑GTTTA‑AGGGGGTG‑GCGGGGGG‑G‑‑G‑G‑TTTTTTTTTCGGGGGTTTCGGAAAgggggttgcc                                                                      >  2:737393/1‑142 (MQ=255)
                                                                                aaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTTGTGACGGCATTATCGGCGGGTGCGG‑CGGGTATGAGGGGG‑TGGGT‑G‑TTT‑ATGGGGGT‑GGCGGGGG‑C‑‑G‑G‑TTTTTTTTTTCGGGGTTTTCGAAAACTGGATC‑GCCact                                                                 >  2:548893/1‑149 (MQ=255)
                                                                                 aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTg                                                                <  1:1064712/151‑1 (MQ=255)
                                                                                 aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTg                                                                <  1:1120201/151‑1 (MQ=255)
                                                                                    tttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGT‑CGGGGTTTGGGGGG‑GGGGT‑G‑TTT‑ATTGGGGT‑GCCGGGGG‑C‑‑T‑G‑TTTTTTTTTTCGGGGTTTTCCCAAGCCGGGTC‑GCCCtt                                                                 >  2:1039309/1‑147 (MQ=255)
                                                                                    tttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGG‑CGGGGATGAGCGGG‑GGGGT‑G‑TTT‑ATTGGGGT‑GGCGGGCG‑C‑‑C‑GTTTTTTTTTTCCGGGGTTTCCG‑AAGCCGGGTC‑Gccccttg                                                               >  1:225042/1‑144 (MQ=255)
                                                                                                atTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGgc                                                  <  2:50485/150‑1 (MQ=255)
                                                                                                  tCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGgcgc                                                <  1:90753/150‑1 (MQ=255)
                                                                                                    tGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGG‑CGGGTTTGGGGGGG‑TGGGT‑G‑GTT‑ATTGGGGT‑GGCGGGCGCC‑‑G‑G‑TTTTTTTTTCGGGGTTTTCCGAAAACGGGTTG‑GCCTTTGG‑CTGACCTTAgggg                                                 >  1:24785/1‑145 (MQ=255)
                                                                                                        cGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGG‑TGGGTATGAGCGGCGGGGCG‑T‑TTA‑AGGGGGTT‑G‑CGGGGG‑C‑‑T‑G‑TTTTTTTTTTCGGGGTTTTCGAAAACTCGAACTCCCATTGGCC‑GGCATTTAGCGCGTgggt                                          >  2:998620/1‑148 (MQ=255)
                                                                                                            cGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAAc                                      <  2:367016/150‑1 (MQ=255)
                                                                                                                  cTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGG‑CGGGTATGAGCGGG‑TGGGT‑G‑TTT‑ATGGGGGT‑GGCGGGCG‑C‑‑T‑G‑TTTTTTTTTTCGGGGTTTTCCAAAACTGGATC‑GCCCTTTGCCTGACATTAGGCGCGTGGGTTAAATggga                                 >  2:731094/1‑147 (MQ=255)
                                                                                                                       gCCATGCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCAGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtggtg                          <  2:1123892/151‑1 (MQ=255)
                                                                                                                        ccATTCGTGACGGCATTATCGGCGGGGGCGT‑CGGGTATGAGCGGC‑TGGGG‑G‑GTT‑TTGGGGGT‑GCCGGGGG‑C‑‑C‑G‑TTTTTTTTTTCGGGGTTTTCGAAAGCTGGGTC‑GCCATTGGGCTGACATTAGGCGCGTGGGTTAACTTGAAGCtggggg                         >  2:127428/1‑151 (MQ=255)
                                                                                                                        ccATTCGTGACGGCATTATCGGCGGGGGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGGT‑GCCGGGCG‑C‑‑T‑G‑GTTTTTTTTTCGGGGTTTTCGAAAGCTGGATC‑GCCATTTGCCTGACATTTAGCGCGGGGGTTAACTGGAAGCtggtgg                         >  2:1071707/1‑151 (MQ=255)
                                                                                                                              gTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGc                   <  2:924318/151‑1 (MQ=255)
                                                                                                                                               gggTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAcc  <  1:905186/151‑1 (MQ=255)
                                                                                                                                               gggTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAc   <  1:1007242/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                               
AGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGT‑CGGATATGAGCGGC‑TGGCT‑G‑TTA‑ATGGGGTT‑GCCGGGCG‑C‑‑T‑G‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATC‑GCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACC  >  NC_000913/1079314‑1079607

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: ‑