Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,828,334 A→T 12.2% I111N (ATT→AAT)  gltS ← glutamate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,828,3340AT12.2% 120.1 / 9.0 41I111N (ATT→AAT) gltSglutamate transporter
Reads supporting (aligned to +/- strand):  ref base A (17/19);  new base T (2/3);  total (19/22)
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.41e-01

ACTTCCGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAA  >  NC_000913/3828187‑3828482
                                                                                                                                                   |                                                                                                                                                    
aCTTCCGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATAccaatgccaat                                                                                                                                                     <  1:36172/149‑1 (MQ=255)
aCTTCCGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATAccaatgccaat                                                                                                                                                     <  4:43038/149‑1 (MQ=255)
  ttCCGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATgg                                                                                                                                                   >  4:18681/1‑149 (MQ=255)
  ttCCGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATgg                                                                                                                                                   >  2:31117/1‑149 (MQ=255)
      gTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCCGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCATTGGCAtt                                                                                                                                               <  4:80285/149‑1 (MQ=255)
            ttGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCAt                                                                                                                                         <  3:87857/149‑1 (MQ=255)
               gTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCAc                                                                                                                                      >  1:60656/1‑149 (MQ=255)
                      cATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGa                                                                                                                               >  4:3015/1‑149 (MQ=255)
                             ttCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCAATTTGCATCACCAACAGACcaacaa                                                                                                                        <  4:47495/149‑1 (MQ=255)
                              tCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCaacaac                                                                                                                       <  4:82310/149‑1 (MQ=255)
                                 aTGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACcac                                                                                                                    <  1:69959/149‑1 (MQ=255)
                                   gAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCATTGGCAttggtat                                                                                                                                          >  3:31275/1‑120 (MQ=255)
                                   gAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCATTGGCAttggtat                                                                                                                                          <  4:31275/125‑6 (MQ=255)
                                    aaCAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCa                                                                                                                                                 <  4:29398/117‑1 (MQ=255)
                                    aaCAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCa                                                                                                                                                 >  3:29398/1‑117 (MQ=255)
                                        aTTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCgag                                                                                                             >  1:44147/1‑149 (MQ=255)
                                         tttACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCaa                                                                                                                                   <  2:27337/126‑1 (MQ=255)
                                         tttACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCaa                                                                                                                                   >  1:27337/1‑126 (MQ=255)
                                         tttACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCaaga                                                                                                            >  4:15955/1‑149 (MQ=255)
                                          ttACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCaagaa                                                                                                           <  2:56725/149‑1 (MQ=255)
                                            aCTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCGagaaga                                                                                                         <  2:44147/149‑1 (MQ=255)
                                              tCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATg                                                                                                       >  2:28197/1‑149 (MQ=255)
                                              tCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATg                                                                                                       <  3:75758/149‑1 (MQ=255)
                                                 aCGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCATTGGCAttgg                                                                                                                                             <  4:79734/108‑3 (MQ=255)
                                                 aCGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCATTGGCAttgg                                                                                                                                             >  3:79734/1‑106 (MQ=255)
                                                        gcCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGAt                                                                                                        <  2:47633/138‑1 (MQ=255)
                                                        gcCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGAt                                                                                                        >  1:47633/1‑138 (MQ=255)
                                                              aCCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCAc                                                                                       <  3:97790/149‑1 (MQ=255)
                                                                             aaGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTgg                                                                        <  2:36953/149‑1 (MQ=255)
                                                                                agTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCaa                                                                     >  4:41974/1‑149 (MQ=255)
                                                                                    aTAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGtt                                                                 <  1:31117/149‑1 (MQ=255)
                                                                                     tAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTg                                                                <  3:5935/149‑1 (MQ=255)
                                                                                       gAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGc                                                              <  2:80679/149‑1 (MQ=255)
                                                                                                 aCAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCCATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAATCTGGCAATGTTGGCGTTCAGGCCa                                                    <  3:85404/149‑1 (MQ=255)
                                                                                                        cATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAAtggtgg                                             <  2:68332/149‑1 (MQ=255)
                                                                                                            aGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCgaa                                         <  2:50709/149‑1 (MQ=255)
                                                                                                               ggATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCgaagaa                                      >  3:51860/1‑149 (MQ=255)
                                                                                                                    aaGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCa                                 >  1:29232/1‑149 (MQ=255)
                                                                                                                      gCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGc                               >  4:48665/1‑149 (MQ=255)
                                                                                                                          tAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTa                           >  3:96279/1‑149 (MQ=255)
                                                                                                                               aGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGa                      <  4:36011/149‑1 (MQ=255)
                                                                                                                                       aTACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAg              >  2:22898/1‑149 (MQ=255)
                                                                                                                                        tACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAgg             >  3:78572/1‑149 (MQ=255)
                                                                                                                                        gcaaaatgccattgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAgg             >  4:20859/5‑149 (MQ=255)
                                                                                                                                          aaaatgccattgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGa           >  2:11851/3‑149 (MQ=255)
                                                                                                                                          aaaatgccattgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGa           >  2:60931/3‑149 (MQ=255)
                                                                                                                                           caatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGAc          <  3:15955/149‑1 (MQ=255)
                                                                                                                                           aaatgccattgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAgg                                                       <  1:92914/103‑1 (MQ=255)
                                                                                                                                           aaatgccattgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAgg                                                       >  2:92914/2‑104 (MQ=255)
                                                                                                                                            aatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACa         >  3:85982/1‑149 (MQ=255)
                                                                                                                                                 caatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCa    >  3:94732/1‑149 (MQ=255)
                                                                                                                                                   atgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCaaa  <  3:19067/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACTTCCGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAAA  >  NC_000913/3828187‑3828482

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