Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,862,634 A→G intergenic (‑205/‑137) msrB ← / → gapA methionine sulfoxide reductase B/glyceraldehyde‑3‑phosphate dehydrogenase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,862,6340AG97.9% 156.9 / NA 47intergenic (‑205/‑137)msrB/gapAmethionine sulfoxide reductase B/glyceraldehyde‑3‑phosphate dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (21/25);  minor base T (0/1);  total (21/26)
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

CTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGA  >  NC_000913/1862512‑1862766
                                                                                                                          |                                                                                                                                    
cTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGt                                                                                                                      <  1:1336509/139‑1 (MQ=255)
               acgaAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAg                                                                                                       >  1:700262/1‑139 (MQ=255)
                        cTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGAttt                                                                                              <  2:700262/139‑1 (MQ=255)
                             ccTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGc                                                                                         >  2:259925/1‑139 (MQ=255)
                                      tgatgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                                <  2:311371/107‑1 (MQ=255)
                                      tgatgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCt                                                                                                                >  1:311371/1‑107 (MQ=255)
                                      tgatgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGa                                                                                <  2:901579/139‑1 (MQ=255)
                                       gatgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGac                                                                               >  1:1168139/1‑139 (MQ=255)
                                        atgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCACGCCCAAACGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                          >  2:952855/1‑111 (MQ=255)
                                        atgaAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGca                                                                                                          <  1:952855/111‑1 (MQ=255)
                                             aaTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGagtcag                                                                                                                            >  1:775588/1‑88 (MQ=255)
                                             aaTCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGagtcag                                                                                                                            <  2:775588/88‑1 (MQ=255)
                                                       ttATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGATTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                               <  2:1168139/139‑1 (MQ=255)
                                                               aTTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTaa                                                       >  2:284333/1‑139 (MQ=255)
                                                                        tAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGt                                              >  1:308382/1‑139 (MQ=255)
                                                                             ttCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTcgcg                                                                                                                       <  2:1302177/61‑1 (MQ=255)
                                                                             ttCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTcgcg                                                                                                                       >  1:1302177/1‑61 (MQ=255)
                                                                                  ggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGagt                                                                                                                               <  2:1240674/48‑1 (MQ=255)
                                                                                  ggCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGagt                                                                                                                               >  1:1240674/1‑48 (MQ=255)
                                                                                    cgcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGc                                  <  1:284333/139‑1 (MQ=255)
                                                                                     gcCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                                 <  1:1208897/139‑1 (MQ=255)
                                                                                         aCCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCt                             >  2:279620/1‑139 (MQ=255)
                                                                                         aCCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCt                             >  2:980215/1‑139 (MQ=255)
                                                                                               tgtgGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                       <  1:980215/139‑1 (MQ=255)
                                                                                               tgtgGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                       >  1:911140/1‑139 (MQ=255)
                                                                                               tgtgGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                       <  2:911140/139‑1 (MQ=255)
                                                                                               tgtgGTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTAtt                       <  1:279620/139‑1 (MQ=255)
                                                                                                   gTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGc                                  <  2:1144765/124‑1 (MQ=255)
                                                                                                   gTCTCAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGc                                  >  1:1144765/1‑124 (MQ=255)
                                                                                                    tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTg                                               <  2:213191/110‑1 (MQ=255)
                                                                                                    tctcAAGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTg                                               >  1:213191/1‑110 (MQ=255)
                                                                                                        aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACaa              >  1:984623/1‑139 (MQ=255)
                                                                                                        aaGCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACATTTTATTCACTAACaa              >  1:1154192/1‑139 (MQ=255)
                                                                                                          gCCCAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAAt            <  2:143714/139‑1 (MQ=255)
                                                                                                             cAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                                 >  2:710481/1‑115 (MQ=255)
                                                                                                             cAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCa                                 <  1:710481/115‑1 (MQ=255)
                                                                                                             cAAAGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGc         <  2:764279/139‑1 (MQ=255)
                                                                                                               aaGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtg       <  1:259925/139‑1 (MQ=255)
                                                                                                                aGGAAGAGTGTGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGAGTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtgg      <  2:1154192/139‑1 (MQ=255)
                                                                                                                aGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtgg      <  2:287653/139‑1 (MQ=255)
                                                                                                                aGGAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCAAAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGGTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtgg      >  1:287653/1‑139 (MQ=255)
                                                                                                                  gAAGAGTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCtggtg    <  2:308382/139‑1 (MQ=255)
                                                                                                                    agagTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGa  <  2:984623/139‑1 (MQ=255)
                                                                                                                      agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                               >  2:560942/1‑76 (MQ=255)
                                                                                                                      agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACg                                                               <  1:560942/76‑1 (MQ=255)
                                                                                                                      agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                   >  2:26538/1‑104 (MQ=255)
                                                                                                                      agTGGGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAgg                                   <  1:26538/104‑1 (MQ=255)
                                                                                                                          |                                                                                                                                    
CTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGCGAGTCAGTCGCGTAATGCTTAGGCACAGGATTGATTTGTCGCAATGATTGACACGATTCCGCTTGACGCTGCGTAAGGTTTTTGTAATTTTACAGGCAACCTTTTATTCACTAACAAATAGCTGGTGGA  >  NC_000913/1862512‑1862766

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