Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 685,941 (C)6→7 coding (729/741 nt) gltJ ← glutamate/aspartate ABC transporter membrane subunit GltJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913685,9351.C86.4% 60.8 / 3.8 22coding (735/741 nt)gltJglutamate/aspartate ABC transporter membrane subunit GltJ
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (5/14);  total (7/15)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CACAAAATACCAATCACTACCGCCGTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAA  >  NC_000913/685820‑686061
                                                                                                                    |                                                                                                                                 
caaaaaaTTCAAAGCACAACCGTCGTGACGGCGATTTGCCGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCg                                                                                                         <  2:223769/137‑1 (MQ=255)
     aaTACCAATCACTACCGCCGTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACtt                                                                                                    <  2:376012/140‑1 (MQ=255)
         ccAATCACTACCGCCGTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACg                                                                                                <  2:445740/140‑1 (MQ=255)
          cAATCACTACCGCCGTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTAGTTGCCCCGCCAGGTGGCCAGGAAGGCGAACTTTACGt                                                                                               >  1:200177/1‑140 (MQ=255)
             tCACTACCGCCGTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCa                                                                                            >  1:360063/1‑140 (MQ=255)
             tCACTACCGCCGTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATATGCTCCCTCATGTTGCCAGGCAGGCGAACTTTACGTTCa                                                                                            >  2:196291/1‑140 (MQ=255)
               cctcccgccgTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAAc                                                                                          <  2:444822/136‑1 (MQ=255)
                 tACCGCCGTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCa                                                                                        <  1:349292/140‑1 (MQ=255)
                 tACCGCCGTGACGGTCATTTTCCGGATGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCa                                                                                        <  1:196291/140‑1 (MQ=255)
                   ccgccgTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCACATGTTGGCAGGCAGGCGAACGTTACGTTCAACCAGc                                                                                      >  1:172630/1‑140 (MQ=255)
                          gACGGTGATTTTCAGAGTGATCGCCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGtcatca                                                                               <  1:280994/140‑1 (MQ=255)
                           aCGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCAc                                                                              <  1:534585/140‑1 (MQ=255)
                              gTGATTTTCAGGGTGATCACCATCCCGTCGACCAGCTATCGCAGGGAATGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAg                                                                           <  2:172630/140‑1 (MQ=255)
                                                         tCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATa                                                <  2:546380/140‑1 (MQ=255)
                                                            aGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCg                                                              <  2:462544/123‑1 (MQ=255)
                                                            aGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCg                                                              >  1:462544/1‑123 (MQ=255)
                                                                   atGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTAATTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATACGCCAGCGTa                                      <  2:525806/140‑1 (MQ=255)
                                                                                     tGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTcc                    >  2:354511/1‑140 (MQ=255)
                                                                                               gTCAAACTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACCCCCAGGCATggg          >  2:421357/1‑140 (MQ=255)
                                                                                               gTCAAACTCGTACATTAATTGCGCCCCCATGGTGCGAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATggg          <  2:200177/140‑1 (MQ=255)
                                                                                                  aaaCTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAAAGACTCCCCGGCATGGGGTg       >  1:492477/1‑140 (MQ=255)
                                                                                                    aCTCGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAg     <  2:395634/140‑1 (MQ=255)
                                                                                                       cGTACATTATTTGCCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTaa  <  2:453290/140‑1 (MQ=255)
                                                                                                                    |                                                                                                                                 
CACAAAATACCAATCACTACCGCCGTGACGGTGATTTTCAGGGTGATCACCAGCCCGTCGAGCAGATATGGCAGGGAAGGGACAATGGAACTCCAGTCAAACTCGTACATTATTTGCCCCCCATGTTGCCAGGCAGGCGAACTTTACGTTCAACCAGCGTCATCACCAGCATGATAAAAGCGTTAATCAACACATACGCCAGCGTAATGGCGGTAAACGACTCCCAGGCATGGGCTGAGTAA  >  NC_000913/685820‑686061

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