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

CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  NC_000913/1153423‑1153722
                                                                                                                                                      |                                                                                                                                                        
cATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                                                                                                                                                          <  3:380187/151‑1 (MQ=255)
   atCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTAc                                                                                                                                                        <  4:1514535/150‑1 (MQ=255)
         cGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACAt                                                                                                                                                <  4:409444/151‑1 (MQ=255)
          gACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGcacaat                                                                                                                                               >  2:85970/1‑149 (MQ=255)
              aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGa                                                                                                                                             <  4:22918/149‑1 (MQ=255)
              aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAAc                                                                                                                                           <  2:1123127/151‑1 (MQ=255)
                aGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGCAGCACAATGaaa                                                                                                                                           >  1:308517/1‑148 (MQ=255)
                   aTGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGGACATTTaaaacggg                                                                                                                                      >  4:679745/1‑145 (MQ=255)
                           cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCACACTTAAccaccgaaccg                                                                                                                                 >  4:1473275/1‑138 (MQ=255)
                           cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCa                                                                                                                              <  1:436172/150‑1 (MQ=255)
                           cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCa                                                                                                                              <  4:175185/150‑1 (MQ=255)
                              cGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGAACATTTAACAACGGAcc                                                                                                                                  >  4:1649382/1‑143 (MQ=255)
                                  tgatgaTTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGa                                                                                                                       <  2:1070197/150‑1 (MQ=255)
                                    atgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGAtt                                                                                                                     <  1:189518/150‑1 (MQ=255)
                                       atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGAGGCCCATTTTACAGCGGGCCGGCAACGCgggat                                                                                                                     >  3:207733/1‑144 (MQ=255)
                                       atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAAAGCACATTTGACAGAGGACCGGCacag                                                                                                                           >  4:1475295/1‑138 (MQ=255)
                                           ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTTAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCg                                                                                                             <  1:97757/150‑1 (MQ=255)
                                           ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCg                                                                                                             <  4:71243/150‑1 (MQ=255)
                                           ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCg                                                                                                             <  3:125063/150‑1 (MQ=255)
                                             tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAAAACAGTGAACAGCGGACCGGCAACGAGGGGtctc                                                                                                                  >  3:7628/1‑144 (MQ=255)
                                             tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACATTGAACAACGGACCGGCAACGCGGGTTCTCTCCGt                                                                                                            >  4:386257/1‑149 (MQ=255)
                                              gtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTa                                                                                                          <  2:192831/150‑1 (MQ=255)
                                               tgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGGAAATTTAACAAGGGACCGGCAACGCGGCTTTTTTCCCTTAc                                                                                                         >  4:1050746/1‑150 (MQ=255)
                                                gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACgg                                                                                                       <  3:48750/151‑1 (MQ=255)
                                                 aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGCGGGGAGAACGCACTGAAAAACGGCCCGGCAACGGGGGTTCTCTCCGTTAggg                                                                                                       >  2:222377/1‑150 (MQ=255)
                                                    tCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCtt                                                                                                    <  4:597450/150‑1 (MQ=255)
                                                    tCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCt                                                                                                     <  4:142923/149‑1 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCGAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCACATTTAAAAGCCGACCGGCAACGCGGAGTTTCTCCGTTACGGCtttttc                                                                                                >  4:784692/1‑147 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGTTCAGCACACTTAAAAACGGGGCGGGCACGCGGGGTTTCTCCCTTTCGGCtttt                                                                                                  >  4:97682/1‑147 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCTGACCGGCAACGCGGATTCTCTCCGTTACGGCTTAt                                                                                                  <  2:1071064/149‑1 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTa                                                                                                   <  4:989248/148‑1 (MQ=255)
                                                                ttGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCACAATTAAAAACGGAACGGCAACGCGGGGTCTCTCCGTTTCGGCCTTTCCttccccat                                                                                        >  3:216465/1‑143 (MQ=255)
                                                                 tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCACACTGAAAAGCGGAGCGGCAACACGGGGTCTTTCCGCTTCGGCTTTTCCTTcaccatac                                                                                      >  2:135345/1‑146 (MQ=255)
                                                                 tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAGAGCAAAATGAACAGCGGACCGGCAACGCGGGGTCTCTCCCTTACCGCCTTTTCTGCACATTa                                                                                       >  2:715797/1‑150 (MQ=255)
                                                                   aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGcaca                                                                                          <  3:229020/145‑1 (MQ=255)
                                                                     cAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGAcc                                                                                  <  2:1329754/151‑1 (MQ=255)
                                                                         agaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGttt                                                                              <  2:160240/151‑1 (MQ=255)
                                                                           aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCACAATTAAAAAAGGACCCGCCAAGCGGGATCTCTTCCTTACGGCTTATTCTCCTCATTACTACCggt                                                                               >  4:1089160/1‑146 (MQ=255)
                                                                             gAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCACAATTAACAACGGACCGGCAACGCGGGTTCTTTCCCTTTCGGCTTTTTCTGCACATTACGACCGgtt                                                                              >  4:823651/1‑147 (MQ=255)
                                                                              aaaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCACCTTGAAAAGCCGACCCGCCACCCGGGTTTTCTCCCTTACCGCTTATCCTTCCCATTTCCACCCTTTTCg                                                                           >  4:932976/1‑149 (MQ=255)
                                                                              aaaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTa                                                                         <  4:430232/151‑1 (MQ=255)
                                                                               aaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTaa                                                                        <  4:1234880/151‑1 (MQ=255)
                                                                                aaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGGAGCACAATTGAAAGCGGGACGGGCACCCGGGGTTTTTCCCGTTCCGCCTATCCTGCCCACTACCACCGCTTggg                                                                           >  2:1172297/1‑145 (MQ=255)
                                                                                aaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAAc                                                                       <  4:19760/151‑1 (MQ=255)
                                                                                 aCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACg                                                                      <  3:309145/151‑1 (MQ=255)
                                                                                  cGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACATTGAAAAACGGACCGGCCACGCGGGTTCTCTCCCTTACGGCTTATTCTTCCCATTACCAACGTTTGCGTaa                                                                        >  4:781778/1‑148 (MQ=255)
                                                                                   gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACATTTAACAACGGGACGGCAACGCGGGGTCTCTTCGTTACGGCCTATCCTTCCCAATAAGACCGGTTGTGTAACg                                                                      >  2:1129953/1‑149 (MQ=255)
                                                                                   gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACATTGAACAGCGGACCGGCCACGCGGGTTCTCTCCCTTTCCGCTTATCCTTCACATTACGACCGCTTTCGTaa                                                                        >  4:137825/1‑147 (MQ=255)
                                                                                    cTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCACATTTAAAAACGGGCCGGCAAAGCGGGTTTTCTCCCTTTCGGCTTATTCTGCCCATTACGACCGGTTGTGTGAAGAgg                                                                   >  4:1432865/1‑151 (MQ=255)
                                                                                     tGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACAATGAAAAACGGGCCGGCCACGCGGGGTCTCTCCCTTTCGGCTTATCCTGCCCAATACGACCCGTTTCGTAACGAggg                                                                  >  3:118864/1‑151 (MQ=255)
                                                                                        cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAt                                                               <  2:226880/151‑1 (MQ=255)
                                                                                          gcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTa                                                             <  2:560189/151‑1 (MQ=255)
                                                                                               actgCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTg                                                        <  4:928314/148‑1 (MQ=255)
                                                                                                cAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTgg                                                       <  2:888355/151‑1 (MQ=255)
                                                                                                 aGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGGACAATAAAAAACGGGACGGCAACCCGGGGTCTCTCCCTTTCCGCTTTTTCTTCACCATACGACCGGTTGCGTGACGAggggggttcgtt                                                         >  2:982298/1‑140 (MQ=255)
                                                                                                  gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACAATGAAAAACGGGCCGGCAAAGCGGGTTTTTTCCGCTACCGCTTATTCTGCACCTTACGACCGGTTTCGGAAAGAggggggttatttgggg                                                     >  4:1447428/1‑139 (MQ=255)
                                                                                                  gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAAGACAATGAAAAGCGGACCGGCAACGCGGGGTCTCTCCCTTACGGCTTTTTCTGCACATTACGACCGTTTGCGTAAAGAGGGGGGTACGtttg                                                       >  3:236603/1‑147 (MQ=255)
                                                                                                  gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGGAAAACATTGTAAAGCGGAGCGGGCACGCGGGTTCTCTCCCTTTAGGCCTTTTCCTCTCAATAAGACCCTTTGGGCTAAGAGGGGGGTACGtttg                                                       >  4:1556084/1‑147 (MQ=255)
                                                                                                    aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACAATTAACAACGGGCCGGCAACGCGGGTTTTCTCCCTTTCGGCCTTTCCTCCACCTTACGACCCGTTGTGGAAAGAGGGGGTTACGTTTgcg                                                     >  4:1248586/1‑149 (MQ=255)
                                                                                                       aTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAg                                                <  4:582288/151‑1 (MQ=255)
                                                                                                            gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCg                                            <  3:186264/150‑1 (MQ=255)
                                                                                                                 gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGc                                        <  4:876842/149‑1 (MQ=255)
                                                                                                                       cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgggg                                <  2:1661199/151‑1 (MQ=255)
                                                                                                                        gTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGCACATTTAAAACCGGGCCGGGAACCCCGGGTGTCTCCCTTACGGCGTTTCCTGTCCCTTTCCACCCGTTTTGTTACAAGGGGGTTATGTTGGGGCTAAACCCGGGGCGGGTgggg                                >  4:461985/1‑150 (MQ=255)
                                                                                                                         tGATCAGTCGCGGTAGTGGCGGGGGAGGGGGCAGCACATTGGACAACGGGCCGGCAACACGGGGTCTCTCCCGTTCGGGTTGTTCTGGACATTAAGACCGGTTGCGCAACGAGGGGGTTACGTTTGCGCGAAGCCCggggggggcgggggg                              >  3:229882/1‑138 (MQ=255)
                                                                                                                          gATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgcgc                             <  2:822558/151‑1 (MQ=255)
                                                                                                                           aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCa                            <  2:911531/151‑1 (MQ=255)
                                                                                                                                   cgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCAt                    <  2:1619799/151‑1 (MQ=255)
                                                                                                                                   cgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCAt                    <  4:1499657/151‑1 (MQ=255)
                                                                                                                                       aGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTg                <  4:226487/151‑1 (MQ=255)
                                                                                                                                           gcgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAgg            <  2:1492382/151‑1 (MQ=255)
                                                                                                                                           gcgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAgg            <  4:725431/151‑1 (MQ=255)
                                                                                                                                            cgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGt           >  4:1470223/1‑151 (MQ=255)
                                                                                                                                                 gAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa      <  4:1131393/151‑1 (MQ=255)
                                                                                                                                                   gggTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAc    >  4:482849/1‑151 (MQ=255)
                                                                                                                                                     gTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAt  <  1:278308/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  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: