Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,191,505 T→A A164A (GCA→GCT)  U712_10980 ← N‑(5'‑phosphoribosyl)anthranilate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,191,5050TA94.1% 87.6 / ‑1.7 34A164A (GCA→GCT) U712_10980N‑(5'‑phosphoribosyl)anthranilate isomerase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (16/16);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

ATCGATTCATCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCGATTCCGCTGGCAAGGTCAATTCCTGCCGGCTGATATTTCATTAAGTCTGTGATAGTGTCCGGATTCACGCCGCCTGCGATAAAACAACGTTTTCCTATAGCAGCCTGCTGGTATTCCGGCACACAGTCCCAAGAAAAAGCAACACCAGTGCCCCCTCTTGCCCCTTTGACAGATGAATCAATCACAAAGCCGTCCACATCATCTTTAAA  >  CP006881/2191369‑2191639
                                                                                                                                        |                                                                                                                                      
atCGGTTCATCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGc                                                                                                                                      >  8:48344/1‑139 (MQ=255)
atCGGTTCATCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGc                                                                                                                                      >  7:20314/1‑139 (MQ=255)
   ggttCATCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGa                                                                                                                                    <  1:45447/136‑1 (MQ=255)
   ggttCATCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGGGATGCTATCCGGATTCACGCCGCCAGCGa                                                                                                                                    <  4:275538/136‑1 (MQ=255)
    gttCATCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGAt                                                                                                                                   >  6:87644/2‑138 (MQ=255)
        aTCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATaaagc                                                                                                                              >  1:101943/1‑137 (MQ=255)
        aTCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATaaag                                                                                                                               <  1:103196/138‑2 (MQ=255)
        aTCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATaaag                                                                                                                               >  6:409972/1‑137 (MQ=255)
          ccTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAagcag                                                                                                                            >  1:14608/1‑138 (MQ=255)
              tcttctAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCg                                                                                                                          <  8:20314/137‑1 (MQ=255)
                   tAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTAcc                                                                                                                    >  2:372856/1‑138 (MQ=255)
                            catcaGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAgccgcc                                                                                                          >  1:163267/1‑139 (MQ=255)
                             atcaGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAgccgcc                                                                                                          <  1:372856/138‑1 (MQ=255)
                             atcaGATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGACGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAgccgcc                                                                                                          <  8:84025/138‑1 (MQ=255)
                                 gATTCTGATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCctgctg                                                                                                     <  8:375163/139‑1 (MQ=255)
                                    tCTGAGCTTTTTGTCCGGTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGCTGCTTTCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATa                                                                                                  >  3:45375/1‑139 (MQ=255)
                                       gATCTTTTTGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTc                                                                                               >  2:9665/1‑139 (MQ=255)
                                             tttGTCCGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGca                                                                                          >  7:97280/1‑138 (MQ=255)
                                                   cGTTTTTTTCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGcacacatt                                                                                    <  7:297269/138‑3 (MQ=255)
                                                          ttCAATTCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGa                                                                             <  1:9665/138‑1 (MQ=255)
                                                               ttCCGCTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGc                                                                       <  5:87644/139‑1 (MQ=255)
                                                                    cTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCctgctga                                                                                                    >  3:154683/1‑104 (MQ=255)
                                                                    cTGGCAAGGTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCctgctga                                                                                                    <  4:154683/105‑2 (MQ=255)
                                                                           ggTCAATTCCTTCTGGCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGCAACACCAGTTcc                                                           <  8:97280/139‑1 (MQ=255)
                                                                                         ggCTGCCATTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGCAACACCAGTTCCGCCTCTAGACCCtt                                             >  1:394073/1‑139 (MQ=255)
                                                                                               catTTCAATAGGCGTGTGATGCTATCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGCAACACCAGTTCCGCCTCTAGACCCTTTTACa                                        >  8:338961/2‑138 (MQ=255)
                                                                                                                       tCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGCAACACCAGTTCCGCCTCTAGACCCTTTTACAGATGAATCAATCACAAAGCCGt                  >  4:212251/1‑136 (MQ=255)
                                                                                                                       tCCGGATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGCAACACCAGTTCCGCCTCTAGACCCTTTTACAGATGAATCAATCACAAAGCCGt                  <  3:212251/136‑1 (MQ=255)
                                                                                                                          ggATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGCAACACCAGTTCCGCCTCTAGACCCTTTTACAGATGAATCAATCACAAAGCCGTCAAcat            <  5:179982/139‑1 (MQ=255)
                                                                                                                           gATTCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGCAACACCAGTTCCGCCTCTAGACCCTTTTACAGATGAATCAATCACAAAGCCGTCAAcatt           <  2:14608/139‑2 (MQ=255)
                                                                                                                            aTTCACGCCGCCTGCGATAAAACAACGTTTTCCTATAGCAGCCTGCTGGTATTCCGGCACACAGTCCCAAGAAAAAGCAACACCAGTGCCCCCTCTTGCCCCTTTGACAGATGAATCAATCACAAAGCCGTCcac              <  5:214856/135‑1 (MQ=255)
                                                                                                                            aTTCACGCCGCCTGCGATAAAACAACGTTTTCCTATAGCAGCCTGCTGGTATTCCGGAACACAGTCCCAAGAAAAAGCAACACCAGTGCCCCATCTTGCCCCTTTGACAGATGAATCAATCACAAAGCCGTCcac              >  6:214856/1‑135 (MQ=255)
                                                                                                                             ttCACGCCGCCAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGCAACACCAGTTCCGCCTCTAGCCCCTTTTACAGATGAATCAATCACAAAGCCGTCAAcatta          >  6:89007/1‑136 (MQ=255)
                                                                                                                                    cgccAGCGATAAAGCAGCGTTTACCAATAGCCGCCTGCTGATATTCCGGCACACATTCCCAAGAAAATGCAACACCAGTTCCGCCTCTAGACCCTTTTACAGATGAATCAATCACAAAGCCGTCAACATTATCTTTaaa  <  2:394073/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ATCGATTCATCCTTTCTTCTAAAAGCCTCATCAGATTCTGATCTTTTTGTCCGTTTTTTTCGATTCCGCTGGCAAGGTCAATTCCTGCCGGCTGATATTTCATTAAGTCTGTGATAGTGTCCGGATTCACGCCGCCTGCGATAAAACAACGTTTTCCTATAGCAGCCTGCTGGTATTCCGGCACACAGTCCCAAGAAAAAGCAACACCAGTGCCCCCTCTTGCCCCTTTGACAGATGAATCAATCACAAAGCCGTCCACATCATCTTTAAA  >  CP006881/2191369‑2191639

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