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

CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  NC_000913/1153423‑1153722
                                                                                                                                                      |                                                                                                                                                        
cATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                                                                                                                                                          <  3:380187/151‑1 (MQ=255)
   atCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTAc                                                                                                                                                        <  4:1514535/150‑1 (MQ=255)
         cGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑At                                                                                                                                                <  4:409444/151‑1 (MQ=255)
          gACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGcac‑aat                                                                                                                                               >  2:85970/1‑149 (MQ=255)
              aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGa                                                                                                                                             <  4:22918/149‑1 (MQ=255)
              aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAAc                                                                                                                                           <  2:1123127/151‑1 (MQ=255)
                aGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGCAGCAC‑AATGaaa                                                                                                                                           >  1:308517/1‑148 (MQ=255)
                   aTGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGGAC‑ATTTaaaacggg                                                                                                                                      >  4:679745/1‑145 (MQ=255)
                           cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCAC‑ACTTAAccaccgaaccg                                                                                                                                 >  4:1473275/1‑138 (MQ=255)
                           cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCa                                                                                                                              <  1:436172/150‑1 (MQ=255)
                           cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCa                                                                                                                              <  4:175185/150‑1 (MQ=255)
                              cGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGAAC‑ATTTAACAAC‑GGAcc                                                                                                                                  >  4:1649382/1‑143 (MQ=255)
                                  tgatgaTTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGa                                                                                                                       <  2:1070197/150‑1 (MQ=255)
                                    atgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGAtt                                                                                                                     <  1:189518/150‑1 (MQ=255)
                                       atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGAGGCCC‑ATTTTACAGC‑GGGCCGGCAACGCgggat                                                                                                                     >  3:207733/1‑144 (MQ=255)
                                       atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAAAGCAC‑ATTTGACAGA‑GGACCGGCacag                                                                                                                           >  4:1475295/1‑138 (MQ=255)
                                           ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTTAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCg                                                                                                             <  1:97757/150‑1 (MQ=255)
                                           ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCg                                                                                                             <  4:71243/150‑1 (MQ=255)
                                           ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCg                                                                                                             <  3:125063/150‑1 (MQ=255)
                                             tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAAAAC‑AGTGAACAGC‑GGACCGGCAACGAGGGGtctc                                                                                                                  >  3:7628/1‑144 (MQ=255)
                                             tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTGAACAAC‑GGACCGGCAACGCGGGTTCTCT‑CCGt                                                                                                            >  4:386257/1‑149 (MQ=255)
                                              gtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTa                                                                                                          <  2:192831/150‑1 (MQ=255)
                                               tgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGGAA‑ATTTAACAAG‑GGACCGGCAACGCGGCTTTTTT‑CCCTTAc                                                                                                         >  4:1050746/1‑150 (MQ=255)
                                                gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACgg                                                                                                       <  3:48750/151‑1 (MQ=255)
                                                 aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGCGGGGAGAACGC‑ACTGAAAAAC‑GGCCCGGCAACGGGGGTTCTCT‑CCGTTAggg                                                                                                       >  2:222377/1‑150 (MQ=255)
                                                    tCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCtt                                                                                                    <  4:597450/150‑1 (MQ=255)
                                                    tCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCt                                                                                                     <  4:142923/149‑1 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCGAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCAC‑ATTTAAAAGC‑CGACCGGCAACGCGGAGTTTCT‑CCGTTACGGCtttttc                                                                                                >  4:784692/1‑147 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGTTCAGCAC‑ACTTAAAAAC‑GGGGCGGGCACGCGGGGTTTCT‑CCCTTTCGGCtttt                                                                                                  >  4:97682/1‑147 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑TGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTAt                                                                                                  <  2:1071064/149‑1 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTa                                                                                                   <  4:989248/148‑1 (MQ=255)
                                                                ttGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCAC‑AATTAAAAAC‑GGAACGGCAACGCGGGGTCTCT‑CCGTTTCGGCCTTTCCttccccat                                                                                        >  3:216465/1‑143 (MQ=255)
                                                                 tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCAC‑ACTGAAAAGC‑GGAGCGGCAACACGGGGTCTTT‑CCGCTTCGGCTTTTCCTTcaccatac                                                                                      >  2:135345/1‑146 (MQ=255)
                                                                 tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAGAGCAA‑AATGAACAGC‑GGACCGGCAACGCGGGGTCTCT‑CCCTTACCGCCTTTTCTGCACATTa                                                                                       >  2:715797/1‑150 (MQ=255)
                                                                   aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGcaca                                                                                          <  3:229020/145‑1 (MQ=255)
                                                                     cAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGAcc                                                                                  <  2:1329754/151‑1 (MQ=255)
                                                                         agaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGttt                                                                              <  2:160240/151‑1 (MQ=255)
                                                                           aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCAC‑AATTAAAAAA‑GGACCCGCCAAGCGGGATCTCT‑TCCTTACGGCTTATTCTCCTCATTACTACCggt                                                                               >  4:1089160/1‑146 (MQ=255)
                                                                             gAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCAC‑AATTAACAAC‑GGACCGGCAACGCGGGTTCTTT‑CCCTTTCGGCTTTTTCTGCACATTACGACCGgtt                                                                              >  4:823651/1‑147 (MQ=255)
                                                                              aaaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCAC‑CTTGAAAAGC‑CGACCCGCCACCCGGGTTTTCT‑CCCTTACCGCTTATCCTTCCCATTTCCACCCTTTTCg                                                                           >  4:932976/1‑149 (MQ=255)
                                                                              aaaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTa                                                                         <  4:430232/151‑1 (MQ=255)
                                                                               aaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTaa                                                                        <  4:1234880/151‑1 (MQ=255)
                                                                                aaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGGAGCACAATTGAA‑AGCGGGACGGGCACCCGGGGTTTTTC‑CCGTT‑CCGCCTATCCTGCCCACTACCACCGCTTggg                                                                           >  2:1172297/1‑145 (MQ=255)
                                                                                aaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAAc                                                                       <  4:19760/151‑1 (MQ=255)
                                                                                 aCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACg                                                                      <  3:309145/151‑1 (MQ=255)
                                                                                  cGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTGAAAAAC‑GGACCGGCCACGCGGGTTCTCT‑CCCTTACGGCTTATTCTTCCCATTACCAACGTTTGCGTaa                                                                        >  4:781778/1‑148 (MQ=255)
                                                                                   gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTTAACAAC‑GGGACGGCAACGCGGGGTCTCT‑TCGTTACGGCCTATCCTTCCCAATAAGACCGGTTGTGTAACg                                                                      >  2:1129953/1‑149 (MQ=255)
                                                                                   gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTGAACAGC‑GGACCGGCCACGCGGGTTCTCT‑CCCTTTCCGCTTATCCTTCACATTACGACCGCTTTCGTaa                                                                        >  4:137825/1‑147 (MQ=255)
                                                                                    cTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCAC‑ATTTAAAAAC‑GGGCCGGCAAAGCGGGTTTTCT‑CCCTTTCGGCTTATTCTGCCCATTACGACCGGTTGTGTGAAGAgg                                                                   >  4:1432865/1‑151 (MQ=255)
                                                                                     tGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑AATGAAAAAC‑GGGCCGGCCACGCGGGGTCTCT‑CCCTTTCGGCTTATCCTGCCCAATACGACCCGTTTCGTAACGAggg                                                                  >  3:118864/1‑151 (MQ=255)
                                                                                        cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAt                                                               <  2:226880/151‑1 (MQ=255)
                                                                                          gcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTa                                                             <  2:560189/151‑1 (MQ=255)
                                                                                               actgCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTg                                                        <  4:928314/148‑1 (MQ=255)
                                                                                                cAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTgg                                                       <  2:888355/151‑1 (MQ=255)
                                                                                                 aGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGGAC‑AATAAAAAAC‑GGGACGGCAACCCGGGGTCTCT‑CCCTTTCCGCTTTTTCTTCACCATACGACCGGTTGCGTGACGAggggggttcgtt                                                         >  2:982298/1‑140 (MQ=255)
                                                                                                  gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑AATGAAAAAC‑GGGCCGGCAAAGCGGGTTTTTT‑CCGCTACCGCTTATTCTGCACCTTACGACCGGTTTCGGAAAGAggggggttatttgggg                                                     >  4:1447428/1‑139 (MQ=255)
                                                                                                  gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAAGAC‑AATGAAAAGC‑GGACCGGCAACGCGGGGTCTCT‑CCCTTACGGCTTTTTCTGCACATTACGACCGTTTGCGTAAAGAGGGGGGTACGtttg                                                       >  3:236603/1‑147 (MQ=255)
                                                                                                  gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGGAAAAC‑ATTGTAAAGC‑GGAGCGGGCACGCGGGTTCTCTCCCTTTAGGCCTTTTCCT‑CTCAATAAGACCCTTTGGGCTAAGAGGGGGGTACGtttg                                                       >  4:1556084/1‑147 (MQ=255)
                                                                                                    aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCAC‑AATTAACAAC‑GGGCCGGCAACGCGGGTTTTCT‑CCCTTTCGGCCTTTCCTCCACCTTACGACCCGTTGTGGAAAGAGGGGGTTACGTTTgcg                                                     >  4:1248586/1‑149 (MQ=255)
                                                                                                       aTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAg                                                <  4:582288/151‑1 (MQ=255)
                                                                                                            gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCg                                            <  3:186264/150‑1 (MQ=255)
                                                                                                                 gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGc                                        <  4:876842/149‑1 (MQ=255)
                                                                                                                       cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgggg                                <  2:1661199/151‑1 (MQ=255)
                                                                                                                        gTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGCAC‑ATTTAAAACC‑GGGCCGGGAACCCCGGGTGTCT‑CCCTTACGGCGTTTCCTGTCCCTTTCCACCCGTTTTGTTACAAGGGGGTTATGTTGGGGCTAAACCCGGGGCGGGTgggg                                >  4:461985/1‑150 (MQ=255)
                                                                                                                         tGATCAGTCGCGGTAGTGGCGGGGGAGGGGGCAGCAC‑ATTGGACAAC‑GGGCCGGCAACACGGGGTCTCT‑CCCGTTCGGGTTGTTCTGGACATTAAGACCGGTTGCGCAACGAGGGGGTTACGTTTGCGCGAAGCCCggggggggcgggggg                              >  3:229882/1‑138 (MQ=255)
                                                                                                                          gATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgcgc                             <  2:822558/151‑1 (MQ=255)
                                                                                                                           aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCa                            <  2:911531/151‑1 (MQ=255)
                                                                                                                                   cgGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCAt                    <  2:1619799/151‑1 (MQ=255)
                                                                                                                                   cgGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCAt                    <  4:1499657/151‑1 (MQ=255)
                                                                                                                                       aGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTg                <  4:226487/151‑1 (MQ=255)
                                                                                                                                           gcgggcgAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAgg            <  2:1492382/151‑1 (MQ=255)
                                                                                                                                           gcgggcgAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAgg            <  4:725431/151‑1 (MQ=255)
                                                                                                                                            cgggcgAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGt           >  4:1470223/1‑151 (MQ=255)
                                                                                                                                                 gAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa      <  4:1131393/151‑1 (MQ=255)
                                                                                                                                                   gggTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAc    >  4:482849/1‑151 (MQ=255)
                                                                                                                                                     gTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAt  <  1:278308/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGC‑GGACCGGCAACGCGGATTCTCT‑CCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  NC_000913/1153423‑1153722

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