Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,245,600 G→A *147* (TAG→TAA ycgY → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,245,6000GA96.8% 198.6 / ‑1.2 63coding (441/441 nt)ycgYuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/30);  minor base . (1/1);  total (32/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 = 9.99e-01

AGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAGAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACTTA  >  NC_000913/1245466‑1245730
                                                                                                                                      |                                                                                                                                  
aaccAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTaaaaaa                                                                                                                                >  2:536881/3‑139 (MQ=255)
  ccAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTaaaaaaaa                                                                                                                              >  1:961536/1‑139 (MQ=255)
     aCGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                           >  2:1264334/1‑139 (MQ=255)
     aCGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                           <  1:1506769/139‑1 (MQ=255)
      cGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCAt                                                                                                                          >  2:1458540/1‑139 (MQ=255)
      cGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCAt                                                                                                                          <  2:173894/139‑1 (MQ=255)
       gTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATa                                                                                                                         >  2:544770/1‑139 (MQ=255)
       gTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATa                                                                                                                         >  2:196131/1‑139 (MQ=255)
           gCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTa                                                                                                                     <  1:1264334/139‑1 (MQ=255)
            cGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTAt                                                                                                                    >  1:1334505/1‑139 (MQ=255)
            cGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTAt                                                                                                                    >  2:1090839/1‑139 (MQ=255)
                           gatgaAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGa                                                                                                     <  1:1090839/139‑1 (MQ=255)
                               aaaGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAAc                                                                                                 <  1:544770/139‑1 (MQ=255)
                               aaaGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAAc                                                                                                 <  1:536881/139‑1 (MQ=255)
                                   tGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTa                                                                                                                     <  1:712965/115‑1 (MQ=255)
                                   tGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTa                                                                                                                     >  2:712965/1‑115 (MQ=255)
                                    gAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCAt                                                                                            <  1:196131/139‑1 (MQ=255)
                                      aaCACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATcc                                                                                          >  1:1459059/1‑139 (MQ=255)
                                             aCCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACg                                                                                                >  2:175472/1‑126 (MQ=255)
                                             aCCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACg                                                                                                <  1:175472/126‑1 (MQ=255)
                                             aCCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAg                                                                                   >  2:700175/1‑139 (MQ=255)
                                                  tAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                           >  2:1521435/1‑94 (MQ=255)
                                                  tAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                           <  1:1521435/94‑1 (MQ=255)
                                                  tAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                      >  2:1546570/1‑99 (MQ=255)
                                                  tAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                      <  1:1546570/99‑1 (MQ=255)
                                                   aGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCAt                                                                                                                          >  1:728163/1‑94 (MQ=255)
                                                   aGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCAt                                                                                                                          <  2:728163/94‑1 (MQ=255)
                                                       aTGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACg                                                                         <  1:1458540/139‑1 (MQ=255)
                                                          aCGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGaa                                                                              <  2:1670284/131‑1 (MQ=255)
                                                          aCGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGAACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGaa                                                                              >  1:1670284/1‑131 (MQ=255)
                                                                     ggCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAaaaaaatc                                                                                                                             <  2:1015131/73‑3 (MQ=255)
                                                                     ggCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAAc                                                                                                                             >  1:1015131/1‑73 (MQ=255)
                                                                     ggCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAgc                                                           <  2:961536/139‑1 (MQ=255)
                                                                          ccACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGAAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTg                                                                      >  1:1549959/1‑123 (MQ=255)
                                                                          ccACGTATGGCAAACATTTGCGCAATACTGCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTg                                                                      <  2:1549959/123‑1 (MQ=255)
                                                                              gTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                           <  2:1032547/66‑1 (MQ=255)
                                                                              gTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                           >  1:1032547/1‑66 (MQ=255)
                                                                                   gCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGtt                                                    >  1:1092305/1‑132 (MQ=255)
                                                                                   gCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGtt                                                    <  2:1092305/132‑1 (MQ=255)
                                                                                    cAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACa                                                                                                   <  1:1484933/84‑1 (MQ=255)
                                                                                    cAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACa                                                                                                   >  2:1484933/1‑84 (MQ=255)
                                                                                      aaCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGtt                                                    <  1:287062/129‑1 (MQ=255)
                                                                                      aaCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGtt                                                    >  2:287062/1‑129 (MQ=255)
                                                                                                cAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCctct                                <  2:1459059/139‑1 (MQ=255)
                                                                                                  aTACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAACCATAGCTATGGTGTTAt                                                                                                            <  2:1571149/60‑1 (MQ=255)
                                                                                                  aTACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAACCATAGCTATGGTGTTAt                                                                                                            >  1:1571149/1‑60 (MQ=255)
                                                                                                        cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACg                                                                                                >  1:187860/1‑67 (MQ=255)
                                                                                                        cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACg                                                                                                <  2:187860/67‑1 (MQ=255)
                                                                                                        cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGt                                                                     >  1:27697/1‑94 (MQ=255)
                                                                                                        cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGt                                                                     <  2:27697/94‑1 (MQ=255)
                                                                                                        cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGa                        >  2:1594442/1‑139 (MQ=255)
                                                                                                        cccGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGa                        >  2:1684588/1‑139 (MQ=255)
                                                                                                             cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGtt                   >  2:595076/1‑139 (MQ=255)
                                                                                                             cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGtt                   >  1:950445/1‑139 (MQ=255)
                                                                                                              aGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTg                  <  2:950445/139‑1 (MQ=255)
                                                                                                                gACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCctct                                <  2:417233/123‑1 (MQ=255)
                                                                                                                gACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCctct                                >  1:417233/1‑123 (MQ=255)
                                                                                                                        aaaCCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCt        <  1:1594442/139‑1 (MQ=255)
                                                                                                                        aaaCCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTACt        <  1:369047/139‑1 (MQ=255)
                                                                                                                          aCCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGga                                                                                       >  1:883462/1‑58 (MQ=255)
                                                                                                                          aCCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGga                                                                                       <  2:883462/58‑1 (MQ=255)
                                                                                                                           ccGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGtt                                                    >  1:1201022/1‑92 (MQ=255)
                                                                                                                           ccGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGtt                                                    <  2:1201022/92‑1 (MQ=255)
                                                                                                                             gTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGgtg                                              >  2:1414524/1‑96 (MQ=255)
                                                                                                                             gTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGgtg                                              <  1:1414524/96‑1 (MQ=255)
                                                                                                                              tCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACTTa  >  1:1556736/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                  
AGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAGAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACTTA  >  NC_000913/1245466‑1245730

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