Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,079,4680AG47.5% 74.3 / 19.5 60D55G (GAT→GGT) putPproline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (3/28);  new base G (28/0);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 9.70e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079325‑1079616
                                                                                                                                                |                                                                                                                                                               
gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGc                                                                                                                                                          <  1:963374/151‑1 (MQ=255)
    gTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCggct                                                                                                                                                     <  1:90378/151‑1 (MQ=255)
     tGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATgggggggcg                                                                                                                                                     >  1:1107497/1‑143 (MQ=255)
        cATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCggctggct                                                                                                                                                 <  1:584466/151‑1 (MQ=255)
           tttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTGGGGGATGAGGgggtg                                                                                                                                                    >  1:930331/1‑145 (MQ=255)
               tgtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGGATGAGCgggtggccggttat                                                                                                                                           >  1:497496/1‑143 (MQ=255)
                gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATg                                                                                                                                       <  2:234000/150‑1 (MQ=255)
                gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGCGggtgggtggttat                                                                                                                                          >  2:888633/1‑144 (MQ=255)
                gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGCgggcgggggtttat                                                                                                                                            >  2:889976/1‑137 (MQ=255)
                gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCCTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATg                                                                                                                                       <  1:793416/150‑1 (MQ=255)
                    atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGtt                                                                                                                                  <  2:49755/151‑1 (MQ=255)
                    atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGt                                                                                                                                   <  2:442265/150‑1 (MQ=255)
                      atCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTg                                                                                                                                <  1:725598/150‑1 (MQ=255)
                        cTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTACGGACTGCGAGGGTGGGTGTTTATgggg                                                                                                                                    >  2:598084/1‑145 (MQ=255)
                        cTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATTAGGGGGGGGGTGTTATTTGGGGttcc                                                                                                                                >  1:650925/1‑148 (MQ=255)
                        cTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGCGGGGTATGAGGggctggccggtt                                                                                                                                             >  1:414041/1‑134 (MQ=255)
                         tttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGTGGGGGGGGGTGTTATTGTGGGTGTCgg                                                                                                                            >  2:688860/1‑151 (MQ=255)
                          ttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCg                                                                                                                             <  1:471134/149‑1 (MQ=255)
                          ttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCg                                                                                                                             <  1:7570/149‑1 (MQ=255)
                           tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCggg                                                                                                                           <  2:461848/150‑1 (MQ=255)
                           tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGATGTGCGGGTGGCTGGTTATTGGGGTGGCgggg                                                                                                                           >  2:288397/1‑151 (MQ=255)
                               aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTTTGAGGGGGCGGCGGGTTATTGGGGTGGCggg                                                                                                                           >  1:531906/1‑146 (MQ=255)
                                 gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGtt                                                                                                               <  1:1121321/151‑1 (MQ=255)
                                  atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGttt                                                                                                              <  1:568466/151‑1 (MQ=255)
                                  atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGGGGGGGGGTTATGGGGGTGGCGGGGGC‑‑TGttt                                                                                                              >  2:33384/1‑151 (MQ=255)
                                    attgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGtttt                                                                                                             <  2:593113/150‑1 (MQ=255)
                                       gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTtttctt                                                                                                         <  1:684216/151‑1 (MQ=255)
                                       gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATTGGGGGGCGGGTGTTTATTGGGGTGGCGGGGGC‑‑TGtttttttt                                                                                                         >  2:224778/1‑151 (MQ=255)
                                               ttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGCAATTAGCGGGCGGGCGTTGATGGGGGTGTCGGGGGC‑‑CGTttttttttttcggg                                                                                                  >  1:638848/1‑144 (MQ=255)
                                                   cGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGAGGGGCTGGGGGTTTATGGGGGTGGCGGGGGG‑‑TGGTTTTTTTTTCCGGGttt                                                                                                >  2:235561/1‑149 (MQ=255)
                                                   cGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTcc                                                                                             <  2:467758/151‑1 (MQ=255)
                                                   cGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTCGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTcc                                                                                             <  1:1014418/151‑1 (MQ=255)
                                                    gCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGGGGGGATGGGGGGGGGGGTGTTTATGGGGGGGGCGGGCGCGGTGTTTTTTTTTCggggt                                                                                                  >  1:68829/1‑146 (MQ=255)
                                                     ccTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGATGAGCGGGTGGGTGGTTATGGGGGGGGCGGGGGC‑‑GGTtttttttttttggggtt                                                                                               >  1:367002/1‑138 (MQ=255)
                                                     ccTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGATCGGGTGGCTGTTTATTGGGTTGGCGGGGGGC‑‑GGTTTTTTTTTCGGGGTTTTcg                                                                                             >  2:964534/1‑149 (MQ=255)
                                                      cTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGGGGGGGGGGGTGTTTATTGGGGTTGCGGGCGC‑‑TGTTTTTTTTTTCCGGGttt                                                                                               >  2:1085325/1‑146 (MQ=255)
                                                          cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGGTGGGTGGTTATTGGGGTTGCCGGGGGT‑‑GGTTTTTTTTTCCGGGTTTTCCaaa                                                                                          >  2:592671/1‑148 (MQ=255)
                                                            aTCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATTAGCGGCGGGGTGTTTATGGGGGTGGCGGGGGG‑‑CGTTTTTTTTTTCTGGGTTTTCGaaaaag                                                                                      >  1:303964/1‑146 (MQ=255)
                                                                 cGAGAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCCTGTTTGGGCGGCCGTGCGTTGATGGGGGTGCGAGGGGC‑‑TGGGTTTTTTTTCTGGGATTTCGGAAGCTGGTTCg                                                                               >  1:518805/1‑150 (MQ=255)
                                                                  gAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATAAGGGGGTGGGGGTTTAAGGGGGTGGCGGGGGGG‑‑GGTTTTTTTTTCGGGGGTTTCGGAAAgggggttgcc                                                                               >  2:737393/1‑142 (MQ=255)
                                                                     aaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTTGTGACGGCATTATCGGCGGGTGCGGCGGGTATGAGGGGGTGGGTGTTTATGGGGGTGGCGGGGGC‑‑GGTTTTTTTTTTCGGGGTTTTCGAAAACTGGATCGCCact                                                                          >  2:548893/1‑149 (MQ=255)
                                                                      aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTg                                                                         <  1:1064712/151‑1 (MQ=255)
                                                                      aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTg                                                                         <  1:1120201/151‑1 (MQ=255)
                                                                         tttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGTTTGGGGGGGGGGTGTTTATTGGGGTGCCGGGGGC‑‑TGTTTTTTTTTTCGGGGTTTTCCCAAGCCGGGTCGCCCtt                                                                          >  2:1039309/1‑147 (MQ=255)
                                                                         tttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATGAGCGGGGGGGTGTTTATTGGGGTGGCGGGCGC‑‑CGTTTTTTTTTTCCGGGGTTTCCGAAGCCGGGTCGccccttg                                                                        >  1:225042/1‑144 (MQ=255)
                                                                                     atTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGgc                                                           <  2:50485/150‑1 (MQ=255)
                                                                                       tCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGgcgc                                                         <  1:90753/150‑1 (MQ=255)
                                                                                         tGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTTTGGGGGGGTGGGTGGTTATTGGGGTGGCGGGCGCC‑‑GGTTTTTTTTTCGGGGTTTTCCGAAAACGGGTTGGCCTTTGGCTGACCTTAgggg                                                          >  1:24785/1‑145 (MQ=255)
                                                                                             cGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGTGGGTATGAGCGGCGGGGCGTTTAAGGGGGTTGCGGGGGC‑‑TGTTTTTTTTTTCGGGGTTTTCGAAAACTCGAACTCCCATTGGCCGGCATTTAGCGCGTgggt                                                   >  2:998620/1‑148 (MQ=255)
                                                                                                 cGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAAc                                               <  2:367016/150‑1 (MQ=255)
                                                                                                       cTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGGTGGGTGTTTATGGGGGTGGCGGGCGC‑‑TGTTTTTTTTTTCGGGGTTTTCCAAAACTGGATCGCCCTTTGCCTGACATTAGGCGCGTGGGTTAAATggga                                          >  2:731094/1‑147 (MQ=255)
                                                                                                            gCCATGCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCAGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtggtg                                   <  2:1123892/151‑1 (MQ=255)
                                                                                                             ccATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGCGGCTGGGGGGTTTTGGGGGTGCCGGGGGC‑‑CGTTTTTTTTTTCGGGGTTTTCGAAAGCTGGGTCGCCATTGGGCTGACATTAGGCGCGTGGGTTAACTTGAAGCtggggg                                  >  2:127428/1‑151 (MQ=255)
                                                                                                             ccATTCGTGACGGCATTATCGGCGGGGGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGGTGCCGGGCGC‑‑TGGTTTTTTTTTCGGGGTTTTCGAAAGCTGGATCGCCATTTGCCTGACATTTAGCGCGGGGGTTAACTGGAAGCtggtgg                                  >  2:1071707/1‑151 (MQ=255)
                                                                                                                   gTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGc                            <  2:924318/151‑1 (MQ=255)
                                                                                                                                    gggTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAcc           <  1:905186/151‑1 (MQ=255)
                                                                                                                                    gggTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAc            <  1:1007242/150‑1 (MQ=255)
                                                                                                                                        gCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGaa        <  1:453354/150‑1 (MQ=255)
                                                                                                                                         cGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAAt       >  1:1080312/1‑150 (MQ=255)
                                                                                                                                              ggATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaac  <  2:199744/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                               
GTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGC‑‑TGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079325‑1079616

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: