Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 937,983 G→A *204* (TAG→TAA lolA → lipoprotein chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913937,9830GA98.2% 175.7 / NA 56*204* (TAG→TAAlolAlipoprotein chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/28);  minor base C (1/0);  total (28/28)
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

AGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACATTTGCTGGCTGCGGGGAAGCCGT  >  NC_000913/937850‑938114
                                                                                                                                     |                                                                                                                                   
aGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCa                                                                                                                                <  2:530208/139‑1 (MQ=255)
    tAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCtg                                                                                                                            >  1:714685/1‑139 (MQ=255)
          ggTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGAGCAACGTAAGTAAAGGCACCtgagtgag                                                                                                                      <  1:1071036/139‑1 (MQ=255)
           gTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGc                                                                                                                     >  2:414552/1‑139 (MQ=255)
                      cGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCt                                                                                                          <  1:344432/139‑1 (MQ=255)
                      cGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCt                                                                                                          <  2:520117/139‑1 (MQ=255)
                         tCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGa                                                                                                       >  2:1576980/1‑139 (MQ=255)
                             cgcAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCACGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGAtttt                                                                                                   <  1:414552/139‑1 (MQ=255)
                                  cagTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGa                                                                                              >  2:101586/1‑139 (MQ=255)
                                     ttATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGataa                                                                                           <  2:714685/139‑1 (MQ=255)
                                       aTCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCCGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATTAACGTAAGTAAAGGCACCtg                                                                                                                            <  1:1056798/104‑1 (MQ=255)
                                       aTCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTACAGGCACCtg                                                                                                                            >  2:1056798/1‑104 (MQ=255)
                                       aTCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGataata                                                                                         >  1:1767202/1‑139 (MQ=255)
                                       aTCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGataata                                                                                         <  1:101586/139‑1 (MQ=255)
                                            cTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACtttt                                                                                    >  2:20673/1‑139 (MQ=255)
                                            cTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACtttt                                                                                    <  1:20673/139‑1 (MQ=255)
                                             tGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCt                                                                                                          <  2:1001014/116‑1 (MQ=255)
                                             tGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCt                                                                                                          >  1:1001014/1‑116 (MQ=255)
                                              gAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCa                                                                                  <  1:1576980/139‑1 (MQ=255)
                                               aaaTCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCAcc                                                                                                                              <  2:460330/94‑1 (MQ=255)
                                               aaaTCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCAcc                                                                                                                              >  1:460330/1‑94 (MQ=255)
                                               aaaTCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCaa                                                                                 >  2:907862/1‑139 (MQ=255)
                                                       gCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGc                                                                         <  2:1368240/139‑1 (MQ=255)
                                                       gCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGc                                                                         >  1:1368240/1‑139 (MQ=255)
                                                       gCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGc                                                                         <  1:744376/139‑1 (MQ=255)
                                                              ggggCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGt                                                                  >  1:550366/1‑139 (MQ=255)
                                                                  cTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGc                                                              >  2:197957/1‑139 (MQ=255)
                                                                   tgtgGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCg                                                             >  2:322122/1‑139 (MQ=255)
                                                                     tgGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGa                                                                                              >  1:251965/1‑104 (MQ=255)
                                                                     tgGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGa                                                                                              <  2:251965/104‑1 (MQ=255)
                                                                                aaaTTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTa                                                <  2:550366/139‑1 (MQ=255)
                                                                                 aaTTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAg                                               <  1:322122/139‑1 (MQ=255)
                                                                                   tttACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGca                                             <  2:1721579/139‑1 (MQ=255)
                                                                                     tACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGcaca                                           <  2:1767202/139‑1 (MQ=255)
                                                                                           cACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTa                                                <  1:1659820/128‑1 (MQ=255)
                                                                                           cACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTa                                                >  2:1659820/1‑128 (MQ=255)
                                                                                           cACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTatat                                     >  1:458051/1‑139 (MQ=255)
                                                                                             ccccGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGataa                                                                                           <  2:316943/83‑1 (MQ=255)
                                                                                             ccccGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGataa                                                                                           >  1:316943/1‑83 (MQ=255)
                                                                                                cgccgcAAGGCGTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAGATTTAGCACAGTATATCGGcc                                >  2:1058680/1‑139 (MQ=255)
                                                                                                           gTCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACATTTg                     >  2:813500/1‑139 (MQ=255)
                                                                                                            tCACGGTAGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACATTTgc                    <  2:32467/139‑1 (MQ=255)
                                                                                                                   aGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCt                                                                                                          <  2:991758/46‑1 (MQ=255)
                                                                                                                   aGATGATCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCt                                                                                                          >  1:991758/1‑46 (MQ=255)
                                                                                                                    gatgatCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACATTTGCTGGCTGCg            <  1:197957/139‑1 (MQ=255)
                                                                                                                     atgatCAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACATTTGCTGGCTGCgg           >  1:719388/1‑139 (MQ=255)
                                                                                                                          cAACGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGc                                 <  1:1205432/112‑1 (MQ=255)
                                                                                                                          cAACGTAAGAAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGc                                 >  2:1205432/1‑112 (MQ=255)
                                                                                                                           aaCGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACtttt                                                                                    >  2:587462/1‑60 (MQ=255)
                                                                                                                           aaCGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACtttt                                                                                    <  1:587462/60‑1 (MQ=255)
                                                                                                                            aCGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGataat                                                                                          <  1:433466/53‑1 (MQ=255)
                                                                                                                            aCGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGataat                                                                                          >  2:433466/1‑53 (MQ=255)
                                                                                                                            aCGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACtttt                                                                                    >  2:331934/1‑59 (MQ=255)
                                                                                                                            aCGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACtttt                                                                                    <  1:331934/59‑1 (MQ=255)
                                                                                                                            aCGTAAGTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACATTTGCTGGCTGCGGGGAAGcc    >  2:1681814/1‑139 (MQ=255)
                                                                                                                              gtaagtaaAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCcg                                                                       <  2:1696012/70‑1 (MQ=255)
                                                                                                                              gtaagtaaAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCcg                                                                       >  1:1696012/1‑70 (MQ=255)
                                                                                                                              gtaagtaaAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACATTTGCTGGCTGCGGGGAAGCCGt  >  1:592942/1‑139 (MQ=255)
                                                                                                                                  gTAAAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACATTTGCTGGCTGCGGGGa        >  1:1139301/1‑129 (MQ=255)
                                                                                                                                     |                                                                                                                                   
AGTTTAGCGCGGTGGAGCAGGACGATCAGCGCAGCAGTTATCAACTGAAATCCCAGCAAAATGGGGCTGTGGATGCAGCGAAATTTACCTTCACCCCGCCGCAAGGCGTCACGGTAGATGATCAACGTAAGTAGAGGCACCTGAGTGAGCAATCTGTCGCTCGATTTTTCGGATAATACTTTTCAACCTCTGGCCGCGCGTATGCGGCCAGAAAATTTAGCACAGTATATCGGCCAGCAACATTTGCTGGCTGCGGGGAAGCCGT  >  NC_000913/937850‑938114

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