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

Only 100 of 104 total aligned reads displayed.

CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  NC_000913/1153423‑1153722
                                                                                                                                                      |                                                                                                                                                        
cATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                                                                                                                                                          <  1:711994/151‑1 (MQ=255)
cATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                                                                                                                                                          <  1:1292957/151‑1 (MQ=255)
           aCTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAG‑ca                                                                                                                                                   >  3:508176/1‑146 (MQ=255)
            cTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAG‑cac‑aat                                                                                                                                              >  2:26123/1‑147 (MQ=255)
             tACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAG‑c                                                                                                                                                    >  2:222620/1‑143 (MQ=255)
              aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGaa                                                                                                                                           <  4:553818/150‑1 (MQ=255)
              aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGaa                                                                                                                                           <  3:608577/150‑1 (MQ=255)
              aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGaa                                                                                                                                           <  3:42600/150‑1 (MQ=255)
               cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACa                                                                                                                                         <  1:1324618/151‑1 (MQ=255)
               cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAGAG‑GAC‑ATTTaa                                                                                                                                           >  3:886936/1‑149 (MQ=255)
                aGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAg                                                                                                                                        <  1:1342200/151‑1 (MQ=255)
                aGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGGAG‑CAC‑ATTGaa                                                                                                                                           >  4:464613/1‑148 (MQ=255)
                aGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCAG‑CAC‑ATTGAACa                                                                                                                                         >  4:261730/1‑150 (MQ=255)
                 ggATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGc                                                                                                                                       <  3:56524/151‑1 (MQ=255)
                     gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGa                                                                                                                                    <  3:352257/150‑1 (MQ=255)
                       ttGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGCGGGTAAAA‑CCC‑ACTGAACAGCCGa                                                                                                                                    >  3:385801/1‑148 (MQ=255)
                        tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCg                                                                                                                                 <  3:620034/150‑1 (MQ=255)
                         gTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGc                                                                                                                               <  1:76826/151‑1 (MQ=255)
                           cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCa                                                                                                                              <  3:412618/150‑1 (MQ=255)
                            aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAG‑CAC‑Atttgaaaacgggccgg                                                                                                                                >  3:695733/1‑133 (MQ=255)
                               ggTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAG‑CAC‑ATTTAAAAGCGGaaccgcca                                                                                                                             >  3:1065945/1‑140 (MQ=255)
                                 ttGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAG‑CAC‑ATTTAAAAGCGGGCCgg                                                                                                                                >  3:1125730/1‑142 (MQ=255)
                                   gatgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGAt                                                                                                                      <  2:212483/150‑1 (MQ=255)
                                      gatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATtct                                                                                                                   <  3:495546/150‑1 (MQ=255)
                                              gtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAG‑CAC‑AATTGACAACGGGCCGGCAAAACGGGGtctcttccttt                                                                                                           >  3:267606/1‑144 (MQ=255)
                                                gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAA‑CAC‑ATTTAAAAGCGGGCCGGCAAACCGGGATCTCTCCGCTTCgg                                                                                                        >  3:323055/1‑151 (MQ=255)
                                                gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACgg                                                                                                        <  2:457076/151‑1 (MQ=255)
                                                 aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGc                                                                                                       <  3:958489/151‑1 (MQ=255)
                                                    tCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGGCG‑CAC‑ATTTAAAACCGGAGCCGCCACCCGGGTTTTCTCCCTTTCgg                                                                                                        >  3:895770/1‑147 (MQ=255)
                                                      tGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAG‑CACAATTGAAAAGGGGCCCGG‑AACGCGGGGTCTCTccccttcgggcttt                                                                                                   >  1:810869/1‑138 (MQ=255)
                                                      tGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAG‑CAA‑AATTAACAACGGGCCGGCAACGCGGGTtctcttcctttc                                                                                                          >  3:87297/1‑136 (MQ=255)
                                                       tgCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGAGAG‑ACC‑AATTTAAAGCCGACCGGGCACGCGGGGTCTCTcccctttggccttt                                                                                                   >  1:1115419/2‑137 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAG‑CAC‑AATTAACAACGGGACGGCAACGCGGGTTCTCTCCGCTACGGcctct                                                                                                   >  2:450185/1‑145 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTAt                                                                                                   <  1:809284/149‑1 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCCGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAG‑CAA‑ATTTAAAAACGGGACGGCGACGCGGGGTCTCTCCGCTTCGGccttt                                                                                                   >  1:114159/1‑145 (MQ=255)
                                                         ccTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTg                                                                                               <  3:759358/151‑1 (MQ=255)
                                                            cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGCG‑GAC‑AGTTAAAAAGGGAACGGGAACGCGGGTTTTTTCCCTTTCGGCtttt                                                                                                   >  4:180689/1‑142 (MQ=255)
                                                            cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAG‑GAC‑ATTGAAAAACGGACCGGCCACGCGGGGTCTCTCCCTTTCGGGtttt                                                                                                   >  2:31090/1‑142 (MQ=255)
                                                            cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGca                                                                                             <  3:320048/150‑1 (MQ=255)
                                                             aGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAG‑CAC‑ACTTAAAAAAGGACCGGGAACGCAGGTTCTCTCCGTTTAGGCTTATCCTTca                                                                                             >  3:1390202/1‑149 (MQ=255)
                                                             aGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGcaca                                                                                           <  2:464619/151‑1 (MQ=255)
                                                                 tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTa                                                                                        <  1:654328/150‑1 (MQ=255)
                                                                  gAACAAGAGATGAGAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGTCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTAc                                                                                       <  3:1317502/150‑1 (MQ=255)
                                                                  gAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTAc                                                                                       <  1:456783/150‑1 (MQ=255)
                                                                   aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGa                                                                                     <  2:436370/151‑1 (MQ=255)
                                                                   aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCAG‑AAC‑ATTGAACAGCGGACCGGCCACGCGGGGTCTCTCCCTTTCGGGtttttcctgaaact                                                                                         >  1:562841/1‑135 (MQ=255)
                                                                     cAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGAcc                                                                                   <  2:174862/151‑1 (MQ=255)
                                                                          gaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAG‑CAC‑AAAGAACAAAGGGCCGGCAACGCGGGGTTTCTCCCGTTCGGCCTTTCCTGcac                                                                                            >  4:449324/1‑137 (MQ=255)
                                                                            tGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAG‑CCC‑ACTTAAAAACGGGCCGGCAAAGCGGGTTTTCTCCGCTTCCGCTTTTCCTGCCCATTTCCAcccctt                                                                               >  3:1245014/1‑144 (MQ=255)
                                                                               aaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTaa                                                                         <  1:197909/151‑1 (MQ=255)
                                                                                aaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAAc                                                                        <  2:45175/151‑1 (MQ=255)
                                                                                 aCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACg                                                                       <  1:109562/151‑1 (MQ=255)
                                                                                 aCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACg                                                                       <  4:532714/151‑1 (MQ=255)
                                                                                  cGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAG‑CCA‑ATTTAAAAACGGGCCCGCAACGCGGGTTTTCTCCCTTTCGGCCTATTCTGCACATTACCACCGGTTTCGGaa                                                                         >  2:510688/1‑148 (MQ=255)
                                                                                    cTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAA‑CAC‑AATGAACAACGGACCGGCAACGCGGGTTTTCTCCCTTACGGCCTTTTCTGCACATTACGAACGTTTGTGTAACGAgg                                                                    >  3:1334366/1‑151 (MQ=255)
                                                                                     tGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAG‑CAC‑CATTAAAAACGGGACGGCAACGCGGGGTCTTTCCCTTTCGGCTTATTCTTCCCAATACCACCGgttttggtacggggg                                                                   >  1:1413212/1‑139 (MQ=255)
                                                                                     tGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAG‑CAC‑AATGTAAAGCGGGACGGGCACGCGGGGTCTTTCCCTTTCGGCCTTTTCCGCCCCTTACGACCCTTTgggggaagaagg                                                                   >  1:173417/1‑140 (MQ=255)
                                                                                     tGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAG‑CAC‑ACTTAAAAACGGACCGGCAAAGCGGGATCTCTCCCTTTCGGCGTTTCCTGCCCCCTACGACCCCTTTCGTGACGAggg                                                                   >  3:263413/1‑151 (MQ=255)
                                                                                      ggCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCCAA‑ATAAAAAAACGG‑CCGGCCACGCGGGTGCTTTCCGTTTCGGCCTATTCTTCCCCCTTAGACCGTTTTTGTGACGAgggg                                                                  >  4:413918/1‑151 (MQ=255)
                                                                                          gcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAG‑AAC‑ATTGAAAAACGGACCGGCAAAGCGGGTTTTCTCCCTTTAAGCTTTTTCTGCCCATTACCACCCGTTTTGTAACGAggggg                                                                 >  3:1013940/1‑147 (MQ=255)
                                                                                          gcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAtt                                                               <  2:346749/150‑1 (MQ=255)
                                                                                            agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACg                                                            <  4:140138/151‑1 (MQ=255)
                                                                                            agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACg                                                            <  1:553853/151‑1 (MQ=255)
                                                                                                cAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTgg                                                        <  1:959491/151‑1 (MQ=255)
                                                                                                cAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTgg                                                        <  2:436963/151‑1 (MQ=255)
                                                                                                   cAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGgcg                                                      <  2:306404/150‑1 (MQ=255)
                                                                                                     aaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTa                                                  <  1:210312/151‑1 (MQ=255)
                                                                                                      aaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTaa                                                 <  3:418871/151‑1 (MQ=255)
                                                                                                      aaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTaa                                                 <  1:1351773/151‑1 (MQ=255)
                                                                                                       aTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAg                                                <  1:410643/151‑1 (MQ=255)
                                                                                                       aTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAg                                                <  1:1415121/151‑1 (MQ=255)
                                                                                                         ggTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAG‑CAC‑ATTTAACAGCGGACCGGCAACGCGGATTCTCTCCGCTAAGGCGTATTCTGCACATTACGACCGGTTGCGTAAAGAGGGGGGTACGGTGGGG‑Ctca                                                 >  3:1367507/1‑146 (MQ=255)
                                                                                                            gtgCTGAAAGTCGTGATCAGTCGCGGTATTGGCGCGCGAGTGTACAG‑CAC‑ATTGACCAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCg                                            <  1:745108/150‑1 (MQ=255)
                                                                                                             tgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAG‑CAC‑ATTGTACAGCGGGCCCGCCACCCGGGGTTTCTCCGCTTCGGCTTATCCTGCACATTACGAACGGTTGCGTAAAGAGGGGGTTACGTTTGCG‑CTAAGccccg                                           >  2:249099/1‑148 (MQ=255)
                                                                                                              gCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGt                                          <  1:930392/150‑1 (MQ=255)
                                                                                                              gCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTg                                         <  2:258360/151‑1 (MQ=255)
                                                                                                                tGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAG‑CAC‑AATGAAAAACGGGCCGGCAAACCGGGGTCTCTTCCTTACGGGTTTTTCTTCCCATTACGACCGTTTTCGGAACGGGGGGGTTTCGCTTGGG‑CTAAAccccggg                                         >  3:1246275/1‑145 (MQ=255)
                                                                                                                    aGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAG‑CAC‑ATTTAAAAGCGGGCCGGGAACGCGGGTTCTCTCCCTTTCGGCGTTTTCTGCACATTTAGACCGCTTGCGGAAAGAGGGGGGTATGTTGGCG‑CTAAACCCGGGGCgg                                      >  4:557797/1‑148 (MQ=255)
                                                                                                                       cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTggg                                 <  3:996267/150‑1 (MQ=255)
                                                                                                                        gTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAG‑CAC‑ATTGAACAACGGACCGGCAACGCGGGTTCTCTCCGTTACGGCTTATTCTGCACATTACGACCGGTTGCGGAAAGAGGGGGGTACGTTGGCG‑CTTAGCCCGGGGCGGGTggggg                               >  3:950613/1‑150 (MQ=255)
                                                                                                                        gTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAG‑AAC‑ACTGAACAACGGGCCCGGCACGCGGGGTTTCTCCCCTTCGGCCTCTTCTCCACACTACGACCCGTTTCGGAAAGAGGGGGGTAAGTTTGCGCCTAACCCCGG‑GCGGCCggggg                               >  4:320483/1‑150 (MQ=255)
                                                                                                                           aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCa                            <  1:307919/151‑1 (MQ=255)
                                                                                                                           aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACAGTACGACCGTGTGCGTAACGAGGTGATTACGTTGGCG‑CTAAGACCGGTGCCGCTGGGgcgc                             <  1:1053597/150‑1 (MQ=255)
                                                                                                                           aTCAGTCGCGGTAGTGGCGGGCGAGGGGGAGG‑CAC‑ACTAAAAAGCGGGCCGGCGACGCGGGGTGTTTCCCTTACGGGCTAACCTTCACCTTAATACCGGTTTCGGAAAGAGGGGGTGAAGTTGGCGCCTAGACCCGGGGGGGCGgggg                                >  4:456174/1‑148 (MQ=255)
                                                                                                                           aTCAGTCGCGGTAGTGGCGGGCGAGGGGACAG‑CAC‑ATTGAAAAGAGGGCCGGCCACGCGGGGTTTTTCCGTTACGGCTTATCCTGCCCATTACGACCGTTTGCGCAACGAGGGGGGTAAGTTTGCG‑CTAAGCCCCGGGCGGCTGGGgcg                              >  1:385589/1‑149 (MQ=255)
                                                                                                                                 cgcgGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGGTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTc                      <  1:260391/151‑1 (MQ=255)
                                                                                                                                       aGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCtt                 <  3:71801/150‑1 (MQ=255)
                                                                                                                                        gTGGCGGGCGAGGGGACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCCTATCCTGCACATTACGACCGGTTTCGTAACGAGGGGATTAAGTTGGCG‑CTAAGCCCGGGGCGGCGGGGGGGCAATCCTCCTCtt                 >  1:1167805/1‑149 (MQ=255)
                                                                                                                                            cgggcgAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGt           <  1:112119/151‑1 (MQ=255)
                                                                                                                                             gggcgAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTa          <  1:886818/151‑1 (MQ=255)
                                                                                                                                             gggcgAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTa          <  1:450235/151‑1 (MQ=255)
                                                                                                                                             gggcgAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTa          <  3:636744/151‑1 (MQ=255)
                                                                                                                                              ggcgAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAAGCCTCATCTTGCAGGTAt         <  1:1029561/151‑1 (MQ=255)
                                                                                                                                               gcgAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTAtt        >  3:523287/1‑151 (MQ=255)
                                                                                                                                                 gAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATCCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCATAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa      <  4:348132/151‑1 (MQ=255)
                                                                                                                                                  aGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa      >  4:257997/1‑150 (MQ=255)
                                                                                                                                                  aGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa      >  2:103504/1‑150 (MQ=255)
                                                                                                                                                  aGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa      >  3:555940/1‑150 (MQ=255)
                                                                                                                                                    ggTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAc    >  3:824424/1‑150 (MQ=255)
                                                                                                                                                     gTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAt  >  4:108277/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAG‑CAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCG‑CTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  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: ‑