Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,462,645 G→C 6.3% *447S (TGA→TCA)  fadL → long‑chain fatty acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,6450GC 6.3% 222.6 / 6.1 63*447S (TGA→TCA) fadLlong‑chain fatty acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base C (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

TTACGCATTTAATAAAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACT  >  NC_000913/2462502‑2462689
                                                                                                                                               |                                             
ttACGCATTTAATAAAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAAc                                          <  2:81798/149‑1 (MQ=255)
ttACGCATTTAATAAAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAAc                                          >  3:370371/1‑149 (MQ=255)
  aCGCATTTAATAAAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAAcgc                                        <  3:791009/149‑1 (MQ=255)
  aCGCATTTAATAAAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAAcgc                                        <  4:862115/149‑1 (MQ=255)
      aTTTAATATAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTc                                    >  3:781085/1‑149 (MQ=255)
      aTTTAATAAAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTc                                    >  3:449118/1‑149 (MQ=255)
         taataaAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGcc                                 <  1:87855/149‑1 (MQ=255)
          aataaAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCt                                <  4:370371/149‑1 (MQ=255)
           ataaAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTc                                    <  2:245355/144‑1 (MQ=255)
           ataaAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTc                                    >  1:245355/1‑144 (MQ=255)
                     ttCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGa                                              >  3:319852/1‑124 (MQ=255)
                     ttCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGa                                              <  4:319852/124‑1 (MQ=255)
                     ttCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAAc                                          >  2:762153/1‑128 (MQ=255)
                     ttCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAAc                                          <  1:762153/128‑1 (MQ=255)
                     ttCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCa                     <  2:189714/149‑1 (MQ=255)
                      tCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATAAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCAc                    <  4:377757/149‑1 (MQ=255)
                           cGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCa               <  4:449118/149‑1 (MQ=255)
                            gACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTc                                    <  1:1073378/127‑1 (MQ=255)
                            gACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTc                                    >  2:1073378/1‑127 (MQ=255)
                               gTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGc          <  4:781085/149‑1 (MQ=255)
                                                gCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATaaa                         >  3:956952/1‑118 (MQ=255)
                                                gCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATaaa                         <  4:956952/118‑1 (MQ=255)
                                                gCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTATCAGAACGCGta                                     >  3:947634/1‑105 (MQ=255)
                                                gCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTATCAGAACGCGta                                     <  4:947634/106‑2 (MQ=255)
                                                  aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGAt                                             <  1:741968/96‑1 (MQ=255)
                                                  aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGAt                                             >  2:741968/1‑96 (MQ=255)
                                                  aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCt                                >  2:1248444/1‑109 (MQ=255)
                                                  aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCt                                <  1:1248444/109‑1 (MQ=255)
                                                  aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCa         <  3:351379/131‑1 (MQ=255)
                                                  aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCa         >  4:351379/1‑131 (MQ=255)
                                                               gAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTATCAGAACGCGta                                     <  2:1347917/91‑2 (MQ=255)
                                                               gAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTATCAGAACGCGta                                     >  1:1347917/1‑90 (MQ=255)
                                                                       aCGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCt                                <  2:1312738/88‑1 (MQ=255)
                                                                       aCGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCt                                >  1:1312738/1‑88 (MQ=255)
                                                                               ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGa    >  3:802560/1‑107 (MQ=255)
                                                                               ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGa    <  4:802560/107‑1 (MQ=255)
                                                                               ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa    <  2:664741/107‑1 (MQ=255)
                                                                               ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa    >  1:664741/1‑107 (MQ=255)
                                                                               ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACt  >  1:1160221/1‑109 (MQ=255)
                                                                               ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACt  >  1:884977/1‑109 (MQ=255)
                                                                               ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACt  <  2:1160221/109‑1 (MQ=255)
                                                                               ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACt  <  2:884977/109‑1 (MQ=255)
                                                                                 aTACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGAc   <  1:1415537/106‑1 (MQ=255)
                                                                                 aTACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGAc   >  2:1415537/1‑106 (MQ=255)
                                                                                  tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCt                                <  3:1190242/77‑1 (MQ=255)
                                                                                  tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCt                                >  4:1190242/1‑77 (MQ=255)
                                                                                    ccAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGt      <  1:1232990/100‑1 (MQ=255)
                                                                                              tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGt                       >  4:227943/1‑74 (MQ=255)
                                                                                              tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGt                       <  3:227943/74‑1 (MQ=255)
                                                                                              tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGa    <  2:170030/92‑1 (MQ=255)
                                                                                              tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGa    >  1:170030/1‑92 (MQ=39)
                                                                                                tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt      <  3:275882/88‑1 (MQ=255)
                                                                                                tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt      >  4:275882/1‑88 (MQ=255)
                                                                                                 gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAg        >  3:209563/1‑85 (MQ=255)
                                                                                                 gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAg        <  4:209563/85‑1 (MQ=255)
                                                                                                     gTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCt                                >  1:1353699/1‑58 (MQ=255)
                                                                                                     gTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCt                                <  2:1353699/58‑1 (MQ=255)
                                                                                                     gTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATaaa                         >  4:237264/1‑65 (MQ=255)
                                                                                                     gTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATaaa                         <  3:237264/65‑1 (MQ=255)
                                                                                                      tAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa    >  4:190841/1‑84 (MQ=255)
                                                                                                      tAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa    <  3:190841/84‑1 (MQ=255)
                                                                                                          gCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCa               >  3:1238030/1‑70 (MQ=255)
                                                                                                          gCCTGGCTGTTCGGTACTAACTTTAACTACGCGGTCTGATAACGCGTTCGCCTGGATCAAGTCACCTGCa               <  4:1238030/70‑1 (MQ=255)
                                                                                                                                               |                                             
TTACGCATTTAATAAAGATGCTTCAGTCGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACT  >  NC_000913/2462502‑2462689

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