Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,187,710 A→G D177D (GAT→GAC)  roxA ← ribosomal protein‑arginine oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,187,7100AG98.2% 188.3 / NA 54D177D (GAT→GAC) roxAribosomal protein‑arginine oxygenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/30);  minor base T (0/1);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACGTTGCAGCACATAATCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATATCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTCCGGTACC  >  NC_000913/1187575‑1187835
                                                                                                                                       |                                                                                                                             
ttCACGTTGCAGCACATAATCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCAc                                                                                                                           <  1:625827/140‑1 (MQ=255)
  cACGTTGCAGCACATAATCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCa                                                                                                                         <  1:332901/140‑1 (MQ=255)
     gTTGCAGCACATAATCGGCAAATCCACTAATCAGTTCCCGTGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGc                                                                                                                      <  1:330138/140‑1 (MQ=255)
       tGCAGCACATAATCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTc                                                                                                                    <  1:533235/140‑1 (MQ=255)
          agcaCATAATCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAg                                                                                                                 >  1:1213691/1‑140 (MQ=255)
               aTAATCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCtt                                                                                                            >  1:1168984/1‑140 (MQ=255)
               aTAATCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCtt                                                                                                            >  2:685674/1‑140 (MQ=255)
                 aaTCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCa                                                                                                          >  1:917752/1‑140 (MQ=255)
                 aaTCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCa                                                                                                          >  1:927986/1‑140 (MQ=255)
                 aaTCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCa                                                                                                          <  1:256228/140‑1 (MQ=255)
                 aaTCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCa                                                                                                          >  1:108145/1‑140 (MQ=255)
                   tCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATc                                                                                                        >  1:685121/1‑140 (MQ=255)
                       cAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCgatg                                                                                                    >  1:1393601/1‑140 (MQ=255)
                       cAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCgatg                                                                                                    <  1:1162381/140‑1 (MQ=255)
                         aaTCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCa                                                                                                          >  2:448691/1‑132 (MQ=255)
                         aaTCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCa                                                                                                          <  1:448691/132‑1 (MQ=255)
                            ccACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGc                                                                                               >  1:173557/1‑140 (MQ=255)
                                 aaTCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCgatga                                                                                                   >  2:572115/1‑131 (MQ=255)
                                 aaTCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCgatga                                                                                                   >  2:1053251/1‑131 (MQ=255)
                                 aaTCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCgatga                                                                                                   <  1:1053251/131‑1 (MQ=255)
                                 aaTCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCgatga                                                                                                   <  1:572115/131‑1 (MQ=255)
                                 aaTCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGa                                                                                          <  1:282118/140‑1 (MQ=255)
                                     aGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACgg                                                                                      >  2:920630/1‑140 (MQ=255)
                                       ttCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATTTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGAt                                                                                    <  2:1271728/140‑1 (MQ=255)
                                       ttCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGAt                                                                                    <  2:1168984/140‑1 (MQ=255)
                                       ttCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGAt                                                                                    <  2:108145/140‑1 (MQ=255)
                                       ttCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGAt                                                                                    >  2:691871/1‑140 (MQ=255)
                                       ttCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGAt                                                                                    <  2:927986/140‑1 (MQ=255)
                                       ttCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGAt                                                                                    <  1:602684/140‑1 (MQ=255)
                                         cccGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCg                                                                                  >  1:691864/1‑140 (MQ=255)
                                         cccGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCg                                                                                  <  2:917752/140‑1 (MQ=255)
                                            gcgTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGAcc                                                                               <  2:691864/140‑1 (MQ=255)
                                                 tttGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGAt                                                                                    <  1:690305/130‑1 (MQ=255)
                                                 tttGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGAt                                                                                    >  2:690305/1‑130 (MQ=255)
                                                  ttGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGTCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAAc                                                                         <  1:677624/140‑1 (MQ=255)
                                                  ttGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAAc                                                                         <  1:920630/140‑1 (MQ=255)
                                                  ttGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAAc                                                                         <  1:1140733/140‑1 (MQ=255)
                                                  ttGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAAc                                                                         <  2:1072915/140‑1 (MQ=255)
                                                          cgAAAGCCCACGGAATAGCTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCgg                                                                 >  2:158203/1‑140 (MQ=255)
                                                                                 cGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCAttt                                          >  2:1324109/1‑140 (MQ=255)
                                                                                 cGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCAttt                                          >  2:912502/1‑140 (MQ=255)
                                                                                 cGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCAttt                                          >  1:1298050/1‑140 (MQ=255)
                                                                                    aTTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCa                                       <  1:1324109/140‑1 (MQ=255)
                                                                                           aGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTc                                <  2:1213691/140‑1 (MQ=255)
                                                                                                         aTGCGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGc                  >  2:512072/1‑140 (MQ=255)
                                                                                                            cGGGAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGcgac               >  1:1202193/1‑140 (MQ=255)
                                                                                                             gggAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTTGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGcgacg              <  2:685121/140‑1 (MQ=255)
                                                                                                             gggAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGcgacg              <  1:389106/140‑1 (MQ=255)
                                                                                                               gAATCCTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCgacgac            <  1:838464/140‑1 (MQ=255)
                                                                                                                    cTGGCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTCCg       <  2:1393601/140‑1 (MQ=255)
                                                                                                                       gCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTCCgg      <  1:237024/138‑1 (MQ=255)
                                                                                                                       gCGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTCCgg      >  2:237024/1‑138 (MQ=255)
                                                                                                                        cGGAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTCCGGTAc   <  1:357441/140‑1 (MQ=255)
                                                                                                                         ggAATATAAAGAATGTCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTCCGGTAcc  <  2:21126/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                             
TTCACGTTGCAGCACATAATCGGCAAATCCACTAATCAGTTCCCGCGTATTTGGCGCGCGAAAGCCCACGGAATAGTTCATCGCATTTTCCAGCGCGTAGCCTTCATGCGGGAATCCTGGCGGAATATAAAGAATATCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTCCGGTACC  >  NC_000913/1187575‑1187835

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑