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

CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACATC  >  NC_000913/1153423‑1153723
                                                                                                                                                      |                                                                                                                                                             
cATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                                                                                                                                                               <  1:777692/151‑1 (MQ=255)
 atatCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTa                                                                                                                                                              <  2:922944/151‑1 (MQ=255)
      cAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGca                                                                                                                                                         <  1:1181212/151‑1 (MQ=255)
      cAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCCGGGGtg                                                                                                                                                              >  1:779245/1‑145 (MQ=255)
              aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGGACACTTAAc                                                                                                                                                 >  1:80155/1‑151 (MQ=255)
               cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGcacct                                                                                                                                                      >  2:737270/1‑143 (MQ=255)
               cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACACTGaa                                                                                                                                                  >  2:138862/1‑149 (MQ=255)
               cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACa                                                                                                                                                <  2:135804/151‑1 (MQ=255)
                aGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGAAGCACAATGGACa                                                                                                                                                >  2:1167713/1‑150 (MQ=255)
                      cTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAACACAATGAACAGCgg                                                                                                                                          >  2:230032/1‑148 (MQ=255)
                       ttGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGGAGGAAAATGAACAGCgg                                                                                                                                          >  1:1334826/1‑147 (MQ=255)
                        tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGGAGCACAtttgtaataggagcgg                                                                                                                                       >  1:630660/1‑137 (MQ=255)
                        tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGAcc                                                                                                                                       <  1:1244032/149‑1 (MQ=255)
                        tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCGCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCg                                                                                                                                      <  1:544535/150‑1 (MQ=255)
                           cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAATGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCa                                                                                                                                   <  1:307765/150‑1 (MQ=255)
                            aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCaa                                                                                                                                  <  1:776266/150‑1 (MQ=255)
                            aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAAc                                                                                                                                 <  2:53415/151‑1 (MQ=255)
                               ggTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGTGCGGGGGAAGCACATTACACAGGCGAcccg                                                                                                                                     >  1:32866/1‑142 (MQ=255)
                                  tgatgaTTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGa                                                                                                                           <  1:568900/150‑1 (MQ=255)
                                   gatgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGGAACCCAATGAAAAACGGGCCGGGAAAGCgg                                                                                                                            >  2:1039769/1‑148 (MQ=255)
                                    atgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCACAATGAAAAACGGGCCGGCCAAGCgg                                                                                                                            >  2:931886/1‑147 (MQ=255)
                                    atgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATtc                                                                                                                       <  1:102929/151‑1 (MQ=255)
                                     tgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGCACAATGAAAAACGGAGCGGGaaagggg                                                                                                                             >  2:811264/1‑141 (MQ=255)
                                       atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCAAAATGGACAACGGGACCGGCacggggg                                                                                                                             >  1:417041/1‑139 (MQ=255)
                                          tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACAATTAAAAACGGGCCGGCAACGCGGGTTTTcccc                                                                                                                  >  2:253290/1‑150 (MQ=255)
                                          tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTcc                                                                                                                  <  1:666942/150‑1 (MQ=255)
                                              gtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGtt                                                                                                               <  1:368086/149‑1 (MQ=255)
                                              gtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTAc                                                                                                            <  2:935872/151‑1 (MQ=255)
                                                    tCTGGCCTCAGCTTGAACAAGAGATGAAAATGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTa                                                                                                     <  2:126484/150‑1 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTa                                                                                                     <  1:706745/148‑1 (MQ=255)
                                                       ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTAt                                                                                                    <  2:417160/149‑1 (MQ=255)
                                                        gCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCt                                                                                                 <  2:268266/151‑1 (MQ=255)
                                                         ccTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTg                                                                                                <  1:408515/151‑1 (MQ=255)
                                                           tCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACAATTAAAAGCGGACCGGCAACGCGGGTTTTCTCCCCTTCGGGCTATCCt                                                                                                 >  2:539668/1‑148 (MQ=255)
                                                              gCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGCACATTGAACAACGGGCCGGCAAAGCGGATTCTCTCCCTTTCCGCTTATTCTGccc                                                                                             >  2:27104/1‑147 (MQ=255)
                                                                ttGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAACACCCCTTGAAAAACGGGACGGCAAAGCGGGATTTTTCCCGTTCGGCCTTTCCTGCACAtt                                                                                          >  2:1072767/1‑150 (MQ=255)
                                                                ttGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCAAATTTAACACCGGACCCGCAACGCGGGTTTTTTCCCGTTCGGCCTATCCTTCACAtt                                                                                          >  2:1300985/1‑150 (MQ=255)
                                                                 tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACAtt                                                                                          <  2:504325/149‑1 (MQ=255)
                                                                         agaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGttt                                                                               <  1:1135102/151‑1 (MQ=255)
                                                                           aTGAAAGCGCTTGCAGCAGAACAGCAAAATGGTGTGCTTAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCCCATTGAACAGCGGACCGGCAACTCGGATTCTCTCCCTTACGACTTCTCCTGCCCATTCCGACCGttta                                                                              <  1:53226/150‑2 (MQ=255)
                                                                           aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGc                                                                             <  1:207240/151‑1 (MQ=255)
                                                                             gAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGCGGGGGGAGAAAAATTACAAACGGACCGGCGAAGCGGAGTTTTTTCCGTTCGGGCTTTCCTCTGCAATACGACCCttt                                                                               >  1:749649/1‑147 (MQ=255)
                                                                                  cGCCGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACTGGACATTTAAAACCGGGCCCGGAACGCGGGATCTCCCCCGTTCTGGTTTTCCTCCCCATTACGACCCCTTTGCGTa                                                                          >  1:1377414/1‑148 (MQ=255)
                                                                                   gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAg                                                                     <  1:881147/151‑1 (MQ=255)
                                                                                   gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAACACAATGAACAACGGGCCCGCAAAGCGGGTTCTCTCCCTTACGGCTTATTCTGCACATTAAGACCCTTTTCGTaa                                                                         >  2:409485/1‑147 (MQ=255)
                                                                                      ggCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAgggg                                                                  <  1:466054/151‑1 (MQ=255)
                                                                                        cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGATAGAACATTTGAAAACGGAGCGGGGAGGCGGGATCTCTCCGTTTCGGCCTATCCTGCACATTACGACCGTTTTCGTGAAGAggggg                                                                 >  2:905273/1‑149 (MQ=255)
                                                                                         agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACACTGTAAAGAGGGCCCGGAACGCGGAGGCTCCCCCTTACCGCTTTTTCTCCACAATACGACCGCTTGTGTTACCAgggg                                                                  >  1:310625/1‑148 (MQ=255)
                                                                                              aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGCACACTGAACAAAGGACCGGCAACCCCGGGTCTCCCCCGTTCGGCTTATTCTGCACATTTAGACCGTTTTTGGAAAGAggggggtatgtt                                                          >  2:757543/1‑143 (MQ=255)
                                                                                              aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGtt                                                          <  1:351784/151‑1 (MQ=255)
                                                                                               aCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACAATTAACAGCGGGCCGGCCACGCGGGGTTTCTCCCTTTCGGGTTATTCTGCACCTTACGACCGTTTTCGGAAAGAggggggttcgt                                                           >  2:640545/1‑142 (MQ=255)
                                                                                                cAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGCACAATTAACAGCGGGCCGGGCAAGCGGGTTCTCTCCGTTACGGCTTATCCTGCACAATACGACCGTTTTCGTAACGAGGGGGGTACGtt                                                          >  2:1002394/1‑149 (MQ=255)
                                                                                                 aGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGgc                                                       <  2:558529/151‑1 (MQ=255)
                                                                                                             tgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGTGGAGGGGGTAGCGCAATCAACAACGGGCCGGCAACACGGAGTCTCTCCGGTTCGGCCTTATCTTCCCAATACCACCGCCTGCGTGACTAGGGGAGTAAGTTTGGGGTAAGccc                                              >  1:1097719/1‑148 (MQ=255)
                                                                                                               cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACCGCACATTTAAAAGCGGACCGGCCACCCGGATTCTCTCCGTTTCGGCTTTTCCTGCACATTACGACCGCTTGCGTAACGAGGGGGTTACGTTTGCGCTAAGccc                                              >  2:441148/1‑146 (MQ=255)
                                                                                                                 gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCg                                        <  1:672607/150‑1 (MQ=255)
                                                                                                                   aaGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCACATTTAACAACGGGCCGGCCACGCGGGTTCTCTCCCTTACGGCTTATTCTGCACATTAAGACCGCTTTCGTAACGAGGGGGGTACGTTGGGGCTAAGCCCCGTGCgg                                       >  2:1195874/1‑149 (MQ=255)
                                                                                                                   aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGGACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCt                                     <  1:1311795/151‑1 (MQ=255)
                                                                                                                     gTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgg                                   <  1:380648/151‑1 (MQ=255)
                                                                                                                       cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTggg                                  <  2:773239/150‑1 (MQ=255)
                                                                                                                        gTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCACAATTAACAGCGGGCCGGCAAAGCGGGGTTTCCCCCTTTCGGCTTTTTCTGCACATTAAGACCGGTTTGGTAAAGAGGGGGGTACGTTGGCGCTCAGCCCGGGGCGGGCggggg                                >  2:552310/1‑150 (MQ=255)
                                                                                                                        gTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAACACATTTAAAAACGGGCCCGGCACGCGGGTTCTCTCCCTTACGGGTTATCCTGCACATTTCGACCGTTTTCGGAAAGAGGGGGTTACGTTGGGGCTAAGGCCCGGGGGGGTgggg                                 >  2:380573/1‑150 (MQ=255)
                                                                                                                         tGATCAGTCGCGGTAGTGGCGGGCGGGGGGACCGAACAATCCACAGCCGGCCCGGCACGCGGGGTTCCTCCCGTACGGGTTATCCTGCCCCGTCCGACCCGTTTGCGGAAGCGGGGGGTAAGTTTGGGCTAAGCCCCGGGGGGCTgggggg                               >  1:1275092/1‑149 (MQ=255)
                                                                                                                           aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCa                             <  1:736645/151‑1 (MQ=255)
                                                                                                                            tCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCCAAATGAAAACCCGGCCGGCCACGCGGGGTTCCTCCCTTCCCGGTTTTTCCTCCCACTACGACCGTTTGTGCAACAGGGGGGGTACGTTTGGGCTCAGCCCGCGGCGGCTgggggggc                             >  1:306640/1‑146 (MQ=255)
                                                                                                                             cAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAAt                           <  2:552897/151‑1 (MQ=255)
                                                                                                                             cAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAAt                           <  1:1303454/151‑1 (MQ=255)
                                                                                                                                 cgcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCTGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTc                       <  1:818732/151‑1 (MQ=255)
                                                                                                                                   cgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCAt                     <  2:21824/151‑1 (MQ=255)
                                                                                                                                     gTAGTGGCGGGGGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGTCCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCt                   <  1:1313629/151‑1 (MQ=255)
                                                                                                                                          ggcgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCa               <  1:40772/150‑1 (MQ=255)
                                                                                                                                            cgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGt            <  2:121512/151‑1 (MQ=255)
                                                                                                                                            cgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGt            <  1:1039003/151‑1 (MQ=255)
                                                                                                                                                 gAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa        >  2:1231856/1‑150 (MQ=255)
                                                                                                                                                 gAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa        >  1:49235/1‑150 (MQ=255)
                                                                                                                                                   gggTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAc     <  1:712999/151‑1 (MQ=255)
                                                                                                                                                    ggTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACa    <  2:1262748/151‑1 (MQ=255)
                                                                                                                                                     gTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAt   >  2:1162109/1‑151 (MQ=255)
                                                                                                                                                      tACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAtc  >  2:1306457/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                             
CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACATC  >  NC_000913/1153423‑1153723

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