Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,153,5730TG46.4% 78.1 / 27.1 56Y92D (TAC→GAC) pabC4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base T (0/30);  new base G (26/0);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.50e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAC  >  NC_000913/1153432‑1153720
                                                                                                                                             |                                                                                                                                                            
cGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGcaaa                                                                                                                                                      >  1:381540/1‑148 (MQ=255)
 gACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CAtt                                                                                                                                                   <  2:202796/151‑1 (MQ=255)
  aCTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTg                                                                                                                                                  <  1:294933/151‑1 (MQ=255)
     aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGca‑caat                                                                                                                                                   >  2:388867/1‑145 (MQ=255)
     aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGaa                                                                                                                                                <  1:128707/150‑1 (MQ=255)
     aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAAc                                                                                                                                               <  2:235733/151‑1 (MQ=255)
     aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAAc                                                                                                                                               <  1:230480/151‑1 (MQ=255)
      cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACa                                                                                                                                              <  2:1008906/151‑1 (MQ=255)
         gATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGca‑caattaa                                                                                                                                                >  2:289864/1‑141 (MQ=255)
          aTGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAACA‑CAATGAACa                                                                                                                                              >  1:856080/1‑147 (MQ=255)
            gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGa                                                                                                                                        <  1:983341/150‑1 (MQ=255)
            gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGAc                                                                                                                                       <  2:172237/151‑1 (MQ=255)
              ttGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCA‑CATAGaaa                                                                                                                                               >  1:595005/1‑141 (MQ=255)
               tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCg                                                                                                                                     <  2:188594/150‑1 (MQ=255)
                gTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCgg                                                                                                                                    <  1:918558/150‑1 (MQ=255)
                  cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCaa                                                                                                                                 <  2:300414/151‑1 (MQ=255)
                   aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAGGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCA‑CAATGAACAGG‑GGAGCgg                                                                                                                                    >  2:433741/1‑147 (MQ=255)
                      ggTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGGA‑CACTGAACAGC‑GGAcccg                                                                                                                                    >  2:930345/1‑142 (MQ=255)
                        ttGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAAAACA‑CATTGAACAAC‑GGACCGGGAAcgccg                                                                                                                            >  2:152605/1‑148 (MQ=255)
                           atgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATtc                                                                                                                       <  1:186462/151‑1 (MQ=255)
                                 tcatgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGACAGACGTGATCAGACGCGGCAGTGGCGGCGGCGGGGACAGGG‑CAATAAACAGCGGGTCCGGCaaaccga                                                                                                                            >  1:1047937/4‑137 (MQ=255)
                                 tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑Tcc                                                                                                                 <  1:574057/150‑1 (MQ=255)
                                 tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑Tcc                                                                                                                 <  1:417035/150‑1 (MQ=255)
                                        aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGc                                                                                                        <  1:577136/151‑1 (MQ=255)
                                                   cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGca                                                                                             <  2:387344/150‑1 (MQ=255)
                                                   cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGca                                                                                             <  1:208793/150‑1 (MQ=255)
                                                   cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGc                                                                                              <  2:147641/149‑1 (MQ=255)
                                                    aGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGcaca                                                                                           <  2:507674/151‑1 (MQ=255)
                                                       ttGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGCA‑CAATTAACAAC‑GGGCCGGCAACGC‑GGGTTTTC‑TCCC‑TTACGGCCTAT‑TCTTcac                                                                                            >  2:994891/1‑147 (MQ=255)
                                                          aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGAAGCA‑CATTGAAAAGA‑GGGGCGGGAACGC‑GGGTTCTC‑TCCG‑TTTCGGGGTAT‑TCTGcac                                                                                            >  2:245545/1‑144 (MQ=255)
                                                          aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCA‑CATTGAAAAAC‑GGACCGGCCACGC‑GGGGTCTC‑TCCC‑TTTCGGCTTTT‑CCTGcc                                                                                             >  1:614817/1‑142 (MQ=255)
                                                              agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGt                                                                                <  2:155311/151‑1 (MQ=255)
                                                                  aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGc                                                                            <  1:501735/151‑1 (MQ=255)
                                                                     aaaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑Ta                                                                        <  1:116350/151‑1 (MQ=255)
                                                                            tGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGCA‑CATTTAACAAC‑GGACCGGCGACGC‑CGGGTCTT‑TCCC‑CTACGGCTTTTCCCTGC‑CATTAAC‑ACCGCCTGCG‑TGACGAggg                                                                 >  1:848931/1‑151 (MQ=255)
                                                                                agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGCA‑CAATGAAAAAC‑GGAGCGGCAACGC‑GGGGTTTC‑TCCG‑TTTCGGGCTTT‑TCTTCCCATTATC‑ACCCCTTGCG‑GAACCAgggg                                                                >  1:143401/1‑148 (MQ=255)
                                                                                agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGAtt                                                             <  1:889561/151‑1 (MQ=255)
                                                                                     aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGATTACGtt                                                        <  2:120524/151‑1 (MQ=255)
                                                                                           aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAAGA‑CCATGAACAAC‑GGAGCGGGCACGCGGGGTTTTT‑CCCT‑TTACGCCTTTT‑CCCCCACATTCGA‑ACCCGTTG‑G‑TAAAGAGGGGGGTACGtttgggg                                                   >  1:787348/1‑145 (MQ=255)
                                                                                           aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGATTACGTTGGCGCt                                                  <  1:937472/151‑1 (MQ=255)
                                                                                            aaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGAGGGGGGGGAGCA‑AATTCAACAGG‑GGAGCGGCCACGC‑GGGGTCTC‑CCCG‑TTTCGGCCTTT‑CCTGCCCATTACA‑ACCGGTTGAGTTAACGGGGGGGTTCCTTTgcg                                                     >  1:502402/1‑146 (MQ=255)
                                                                                                  tgtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAACA‑CATTTAAAAGC‑GGACCGGGAACCC‑GGACTTTT‑TCCG‑CTATCGCCTAT‑TCTGCACAATTCT‑ACCCGTTGCG‑GTAAGAGGGGGGTATGTTGGCGAT‑CAGccc                                           >  2:696323/1‑150 (MQ=255)
                                                                                                        gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCA‑CATTTAAAAGC‑GGGCCCGCCACGC‑GGGGTCTC‑TCCC‑CTTCGGCTTAT‑CCTGCACATTACG‑ACCGTTTTCG‑TGACGAGGGGGTTACGTTTGGGCT‑AAGCCCCGGGCgg                                    >  2:658081/1‑151 (MQ=255)
                                                                                                        gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCA‑CCTTGAAAAGC‑GGGCCGGCCACGC‑GGGTTCTC‑TCCC‑TTTCGGCTTAT‑CCTCCACATTACG‑ACCCCTTGCG‑TAACGAGGGGGTTACGTTTGCGCT‑AAGCCCGGGGCgg                                    >  2:614729/1‑151 (MQ=255)
                                                                                                          aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCA‑CATTGAAAAAC‑GGAACGGGAAAGC‑GGGTTCTCTTCCG‑TTTGGGCTT‑T‑TCTGCCCACTACG‑ACCGCTTTGG‑TTACGAGGGGGTTACGTTTgcgcctagacccggggg                                      >  2:485671/1‑135 (MQ=255)
                                                                                                              cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTgggg                              <  2:442514/151‑1 (MQ=255)
                                                                                                               gTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAACACCATTGAA‑AAA‑GGAGCGGCCACAC‑GGGGTCTC‑TCCG‑TTTCGGGCTTT‑CCTGCACAATTCG‑ACCGGTTGTG‑TAAACAGGGGGATGCGTTGGCGGT‑CAACCCCGGGCGGGCgggg                              >  1:944224/1‑150 (MQ=255)
                                                                                                               gTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAACA‑CATTGAACAGA‑GGGCCGGGAACGC‑GGGTTTTC‑TCCC‑TTTCGGGTTAT‑TCCTCACATTAAGAACCGTTT‑CG‑TAACGAGGGGGGTACGTTGGCGCTAAAGCCCggggggcgggggg                              >  2:332150/1‑140 (MQ=255)
                                                                                                               gTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCA‑CATTTAACAAC‑GGACCGGCAACGC‑GGGTTCTC‑TCCC‑TTACGGCTTTT‑CCTGCCCATTACG‑ACCGCTTTCG‑GAACGAGGGGGGTACGTTTGCGCT‑CAACCCCGGGGGGCGggggg                             >  2:84231/1‑150 (MQ=255)
                                                                                                               gTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAAAGCA‑CCTTTAACAGC‑GGGCCGGCAACAC‑GGAGTCTC‑TCCC‑GTACGGCCTTT‑CCCGCACCTTAAG‑ACCGTTTTCG‑GAACGAGGGGGGTACGTTGGCGCT‑AAACCCGGCGCGGGTggggg                             >  2:955389/1‑150 (MQ=255)
                                                                                                                     aGTCGCGGTAGTGGCGGGCGGGGGGAAAGAA‑CATTGAACCAA‑GGGGCGGCCAAGC‑GGGGTGTC‑TCCGCTTAC‑GCGTAT‑TCCGCACATTTTC‑ACCGTTTGTG‑GCACGAGGGGGGTACGTTTGGGCT‑TAACCCGGGGCGGGCggggggggatc                       >  1:764486/1‑144 (MQ=255)
                                                                                                                        cgcgGTAGTGGCGGGCGGGGGGACAGCA‑CATTGAACAGC‑GGACCGGGAACGC‑GGGTTCTC‑TCCG‑TTTCGGCGTTT‑CCTTCACATTACC‑ACCGTTTTCG‑TAACGAGGGGGATAAGTTGGCGGT‑AAAGCCCGTGGGGGTGGGgggccaacc                      >  2:990391/1‑144 (MQ=255)
                                                                                                                        cgcgGTAGTGGCGGGCGAGGGTACAGCA‑CTTTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGTTGGGGCGCAATCCt                     <  1:686161/150‑1 (MQ=255)
                                                                                                                        cgcgGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTc                    <  2:345395/151‑1 (MQ=255)
                                                                                                                        cgcgGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTc                    <  1:614819/151‑1 (MQ=255)
                                                                                                                                          gggTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAc  <  1:388859/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                            
CGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCA‑CATTGAACAGC‑GGACCGGCAACGC‑GGATTCTC‑TCCG‑TTACGGCTTAT‑CCTGCACATTACG‑ACCGTTTGCG‑TAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAC  >  NC_000913/1153432‑1153720

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