Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,865,855 G→A Q140* (CAG→TAG)  yidE ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,865,8550GA97.9% 135.2 / ‑4.2 47Q140* (CAG→TAG) yidEputative transporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (26/20);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

TCGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCG  >  NC_000913/3865713‑3865999
                                                                                                                                              |                                                                                                                                                
tcGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctaccctgc                                                                                                                                            >  2:952926/1‑149 (MQ=255)
tcGACATAGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGCCCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGctaccctgc                                                                                                                                            >  1:450428/1‑149 (MQ=255)
  gACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCtaccctgccc                                                                                                                                          >  2:431585/1‑149 (MQ=255)
      ttGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAg                                                                                                                                      <  2:233981/149‑1 (MQ=255)
           gcgGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTg                                                                                                                                 <  1:1057285/149‑1 (MQ=255)
           gcgGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTg                                                                                                                                 >  2:884612/1‑149 (MQ=255)
            cgGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTgg                                                                                                                                <  2:160881/149‑1 (MQ=255)
                    caCCCGCAACATCCACGTGGTGAGCAACATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTgg                                                                                                                        <  2:1030550/148‑1 (MQ=255)
                       ccGCAACATCCACATGGTGAACAGAATCACGCAAATGCCGAATGGATACGCCATCGCGTAACTAATCCCCATCTGATCGACCATTTCCCTTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTCCCCCCAGCGCTGGCGTATTGGTaa                                                                                                                     >  1:604037/1‑149 (MQ=255)
                             cATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGAGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCAc                                                                                                               >  1:148686/1‑149 (MQ=255)
                              aTCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCACAATCTGCTACCATGCCCCCAGCGCTGGCGTATTGGTAACCGCAcc                                                                                                              >  2:980198/1‑149 (MQ=255)
                                 cacaTGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGga                                                                                                           <  1:179405/149‑1 (MQ=255)
                                     tGGTGAACAAAATCCCGCAAATGCCGAATGGATACGGAATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGaaa                                                                                                       >  2:1000298/1‑149 (MQ=255)
                                     gggTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCAGCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGaaa                                                                                                       <  2:418651/148‑1 (MQ=255)
                                      ggTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGaaaa                                                                                                      >  1:101328/1‑149 (MQ=255)
                                          aaCAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATccc                                                                                                  <  2:1089728/149‑1 (MQ=255)
                                           aCAAAATCCCGCAAATGCCGAATGGATACGCCATCTCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATcccc                                                                                                 >  2:261802/1‑149 (MQ=255)
                                               aaTCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCTACCGTTTCCATTGGTGGACCCAGGTCGCGCAAAAGATGGTACCCTGCCCCCAGCGCTGGGGTATATGTACCCGTCTCGGACAAATTCCCCGGCa                                                                                             >  2:1065856/1‑149 (MQ=255)
                                                   ccGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCCTCTGATCGACAATTCCCACAGGAGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACtcc                                                                                         >  1:970978/1‑147 (MQ=255)
                                                             cGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCCAAATCTGCTACCCTGCCCCCAGCGCTGGCGGATTGGTAACCGCACCGGCGAAAATCCCCAGCACTACCGGCAGTGGa                                                                               >  2:45275/1‑149 (MQ=255)
                                                             cGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGa                                                                               >  1:1085110/1‑149 (MQ=255)
                                                             cGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGGATTGGTAACCGCACCGGAGGAAATCCCCAGCACTACCGGCAGTGGa                                                                               >  1:306513/1‑149 (MQ=255)
                                                                aTGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCCGCTACCCTGCCCCCCGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAata                                                                            <  1:398736/149‑1 (MQ=255)
                                                                     tatgCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCaaa                                                                       <  1:169811/146‑1 (MQ=255)
                                                                     tACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCCTTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCaaa                                                                       <  2:450428/149‑1 (MQ=255)
                                                                     tACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCaaa                                                                       >  1:983722/1‑149 (MQ=255)
                                                                            tCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTa                                                                <  1:884612/149‑1 (MQ=255)
                                                                             cgcgTAACTCATCCCCCTCTGATCGACCATTGCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCTGCGCGGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTAt                                                               <  2:968793/149‑1 (MQ=255)
                                                                              gcgTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATg                                                              <  2:688227/149‑1 (MQ=255)
                                                                                 tAACTCATCCCCATCTGATCTACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAg                                                           <  2:417263/149‑1 (MQ=255)
                                                                                     tCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCCAACAGTTTGTGCAGGagg                                                       >  2:83772/1‑147 (MQ=255)
                                                                                     tCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATg                                                       >  1:1088255/1‑149 (MQ=255)
                                                                                       aTCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCCACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGACAAAACCCCCAGCACTACCGACAGTGGAATATCACACAGTTTATCCAGACTGGc                                                     >  2:200136/1‑149 (MQ=255)
                                                                                         ccccATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAAGCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCgg                                                   >  1:126181/1‑149 (MQ=255)
                                                                                           ccATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTa                                                 >  1:260078/1‑149 (MQ=255)
                                                                                                 gATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGATTGTGCAGGATGGCGGTAACCaca                                           >  2:580624/1‑147 (MQ=255)
                                                                                                     gACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGaccac                                       >  2:864572/1‑149 (MQ=255)
                                                                                                         aTTTCCATTGGTGTACCACGGTCGCGCAATATCGGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCgat                                   <  2:175063/149‑1 (MQ=255)
                                                                                                              cATTGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGAtgatga                              >  1:969648/1‑149 (MQ=255)
                                                                                                                ttGGTGTACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGCCGGTAACCAGACCACCGATGATGACg                            >  2:947198/1‑149 (MQ=255)
                                                                                                                   gtgtACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGACGATGACGATc                         >  1:596613/1‑149 (MQ=255)
                                                                                                                    tgtACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCa                        <  1:769874/149‑1 (MQ=255)
                                                                                                                     gtACCCAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGGGGAATATCAAACAGTTTATGCCGGATGGCGGTAACCAGCCCACCGATGATGACggacgg                       >  2:6006/1‑144 (MQ=255)
                                                                                                                        cccAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACCGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAAc                    >  2:1082281/1‑149 (MQ=255)
                                                                                                                         ccAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACa                   <  2:228717/149‑1 (MQ=255)
                                                                                                                         ccAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACa                   <  2:107354/149‑1 (MQ=255)
                                                                                                                          cAGGTCGCGCAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGAGGATGACGATCAGAACAg                  <  1:431585/149‑1 (MQ=255)
                                                                                                                                   cAAAATCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGt         <  1:952926/149‑1 (MQ=255)
                                                                                                                                    aaaaTCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGAGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGtt        >  2:649038/1‑149 (MQ=255)
                                                                                                                                        tCTGCTACCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAgg    >  1:47425/1‑149 (MQ=255)
                                                                                                                                          tgctacCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGACGAGGCg  >  2:149433/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TCGACATTGACGCGGAAAATAACCCGCAACATCCACATGGTGAACAAAATCCCGCAAATGCCGAATGGATACGCCATCGCGTAACTCATCCCCATCTGATCGACCATTTCCATTGGTGTACCCAGGTCGCGCAAAATCTGCTGCCCTGCCCCCAGCGCTGGCGTATTGGTAACCGCACCGGAGAAAATCCCCAGCACTACCGGCAGTGGAATATCAAACAGTTTATGCAGGATGGCGGTAACCAGACCACCGATGATGACGATCAGAACAGCAAACAGGTTGAGGCG  >  NC_000913/3865713‑3865999

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑