Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,154,109 G→A *270* (TAG→TAA pabC → 4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,154,1090GA100.0% 156.2 / NA 48*270* (TAG→TAApabC4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/25);  total (23/25)

TCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCTTATCAAAGAAGAGACGATATTTA  >  NC_000913/1153975‑1154238
                                                                                                                                        |                                                                                                                                 
tCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCAt                                                                                                                               >  1:638037/1‑139 (MQ=255)
  gTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGAGCGCCCGAATTAATCAt                                                                                                                               >  1:1186428/1‑139 (MQ=255)
    tGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaa                                                                                                                           >  1:135169/1‑139 (MQ=255)
       aggcagATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaaaa                                                                                                                        <  2:135169/139‑1 (MQ=255)
       aggcagATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaaaa                                                                                                                        <  2:638037/139‑1 (MQ=255)
                    ggTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGTGAGAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGAt                                                                                                             <  2:1186428/139‑1 (MQ=255)
                    ggTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATaa                                                                                                           <  1:509448/139‑1 (MQ=255)
                    ggTTATTTGTAATGCGTTAATGCCAGTGACGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATaa                                                                                                           <  1:509447/139‑1 (MQ=255)
                                     tAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaaaa                                                                                                                        >  1:1635729/1‑109 (MQ=255)
                                     tAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaaaa                                                                                                                        <  2:1635729/109‑1 (MQ=255)
                                              tGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAgtg                                                                                                                     <  2:769665/103‑1 (MQ=255)
                                              tGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAgtg                                                                                                                     >  1:769665/1‑103 (MQ=255)
                                              tGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGa                                                                                                              >  1:1587835/1‑110 (MQ=255)
                                              tGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGa                                                                                                              <  2:1587835/110‑1 (MQ=255)
                                                 tGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaaaa                                                                                                                        >  1:254021/1‑97 (MQ=255)
                                                 tGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaaaa                                                                                                                        <  2:254021/97‑1 (MQ=255)
                                                             ccTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaaaa                                                                                                                        >  2:1356377/1‑85 (MQ=255)
                                                             ccTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaaaa                                                                                                                        <  1:1356377/85‑1 (MQ=255)
                                                             ccTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGtgtg                                                                  >  1:1645637/1‑139 (MQ=255)
                                                             ccTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGtgtg                                                                  <  2:1277613/139‑1 (MQ=255)
                                                               tgtgGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTggg                                                                <  2:1217533/139‑1 (MQ=255)
                                                                 tgGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCAt                                                                                                                               >  2:783156/1‑74 (MQ=255)
                                                                 tgGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCAt                                                                                                                               <  1:783156/74‑1 (MQ=255)
                                                                 tgGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaa                                                                                                                          >  2:877871/1‑79 (MQ=255)
                                                                 tgGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGaaaa                                                                                                                          <  1:877871/79‑1 (MQ=255)
                                                                                  cGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGc                                                                                               >  2:1502365/1‑89 (MQ=255)
                                                                                  cGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGc                                                                                               <  1:1502365/89‑1 (MQ=255)
                                                                                       gCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATc                                                                              >  1:1430916/1‑101 (MQ=255)
                                                                                       gCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGa                                                        <  1:738864/123‑1 (MQ=255)
                                                                                       gCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGa                                                        >  2:738864/1‑123 (MQ=255)
                                                                                       gCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAAGCTTGTTATTGCTGGTGGTACTGGGTATc                                                                              <  2:1430916/101‑1 (MQ=255)
                                                                                                     atatTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTggg                                                                <  1:1360072/101‑1 (MQ=255)
                                                                                                     atatTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTggg                                                                >  2:1360072/1‑101 (MQ=255)
                                                                                                     atatTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGc                          >  2:1469534/1‑139 (MQ=255)
                                                                                                       atTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGAt                                                                                                             >  1:599456/1‑54 (MQ=255)
                                                                                                       atTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGAt                                                                                                             <  2:599456/54‑1 (MQ=255)
                                                                                                          tAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTAt                                                                                                                 <  1:1714302/47‑1 (MQ=255)
                                                                                                          tAGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTAt                                                                                                                 >  2:1714302/1‑47 (MQ=255)
                                                                                                           aGCCCCACTTTTT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGtt                                                                                                                   <  2:1028917/44‑1 (MQ=255)
                                                                                                           aGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGtt                                                                                                                   >  1:1028917/1‑44 (MQ=255)
                                                                                                           aGCCCCACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCTTATCa                    <  2:1645637/139‑1 (MQ=255)
                                                                                                              cccACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                            <  2:50831/64‑1 (MQ=255)
                                                                                                              cccACTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCtgg                                                                                            >  1:50831/1‑64 (MQ=255)
                                                                                                                  cTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAgtgt                                                                                                                    <  1:1641179/36‑1 (MQ=255)
                                                                                                                  cTTTGT‑‑GAGCGCCCGAATTAATCATGAAAAAAgtgt                                                                                                                    >  2:1641179/1‑36 (MQ=255)
                                                                                                                               cccGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTggg                                                                <  2:267507/77‑1 (MQ=255)
                                                                                                                               cccGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTggg                                                                >  1:267507/1‑77 (MQ=255)
                                                                                                                               cccGAATTAATCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCTTATCAAAGAAGAGACGATATTTa  >  1:588355/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                 
TCGTTGCAGGCAGATGAGATGGTTATTTGTAATGCGTTAATGCCAGTGATGCCCGTATGTGCCTGTGGCGATGTCTCCTTTTCGTCAGCAACGTTATATGAATATTTAGCCCCACTTTGT‑‑GAGCGCCCGAATTAGTCATGAAAAAAGTGTTATTGATAATCTTGTTATTGCTGGTGGTACTGGGTATCGCCGCTGGTGTGGGCGTCTGGAAGGTTCGCCATCTTGCCGACAGCAAATTGCTTATCAAAGAAGAGACGATATTTA  >  NC_000913/1153975‑1154238

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