Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,896,819 +C coding (622/657 nt) ygcU ← putative oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,896,8181.C99.4% 422.6 / ‑3.7 153coding (623/657 nt)ygcUputative oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (85/67);  total (85/68)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 163 total aligned reads displayed.

GAAGAGGTTGTTTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGTGTGACGGTGCC  >  NC_000913/2896743‑2896888
                                                                            |                                                                       
gaagaGGTTGTTTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTg                                                                          >  1:5740754/1‑76 (MQ=255)
 aagaGGTTGTTTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGc                                                                        >  1:3447899/1‑76 (MQ=255)
 aagaGGTTGTTTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGc                                                                        <  2:2756563/76‑1 (MQ=255)
 aagaCGCTGTTTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGc                                                                        >  2:3506598/1‑76 (MQ=255)
  agagGTTGTTTTCCGGGGTGAATTTAAAGATTTTCACTGTTGCTTCAGTTATATAGCTCAATGCCCCTTCATTGcc                                                                       <  1:1783260/76‑2 (MQ=255)
  agagGTTGTTTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGcc                                                                       >  2:1651761/1‑75 (MQ=255)
   gagGTTGTTTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGccg                                                                      <  2:4820341/76‑3 (MQ=255)
        tgttTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCgatgat                                                                 <  2:4456037/76‑1 (MQ=255)
        tgttTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCgatgat                                                                 >  1:5254726/1‑76 (MQ=255)
        tgttTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCgatgat                                                                 <  2:3891180/76‑1 (MQ=255)
        tgttTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCgatgat                                                                 >  2:381148/1‑76 (MQ=255)
         gttTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCgatgatg                                                                <  1:4802256/76‑1 (MQ=255)
            ttCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGatgatgatg                                                             >  1:648732/1‑76 (MQ=255)
            ttCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGatgatgatg                                                             <  2:5503207/76‑1 (MQ=255)
             tCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGAtgt                                                            >  1:2469301/1‑76 (MQ=255)
             tCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGAtgt                                                            >  1:5603283/1‑76 (MQ=255)
             tCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGAtgt                                                            >  2:2741925/1‑76 (MQ=255)
             tCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGAtgt                                                            <  2:107467/76‑1 (MQ=255)
              ccGGGGTGAATTTAAAGATTCTCAATGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGAtgtg                                                           >  2:3164784/1‑76 (MQ=255)
               cGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGa                                                          <  1:2595390/76‑1 (MQ=255)
                ggggTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGAc                                                         <  1:4931730/76‑1 (MQ=255)
                ggggTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGAc                                                         >  2:727017/1‑76 (MQ=255)
                ggggTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGAc                                                         <  2:1881959/76‑1 (MQ=255)
                ggggTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGAc                                                         >  2:3101589/1‑76 (MQ=255)
                  ggTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGa                                                       >  2:4050417/1‑76 (MQ=255)
                   gTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGaa                                                      <  2:2850848/76‑1 (MQ=255)
                       aTTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTc                                                  >  1:3859158/1‑76 (MQ=255)
                       aTTTAAAGATTCTCACTGTTACGTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTc                                                  <  1:429884/76‑1 (MQ=255)
                        tttAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACg                                                        >  2:1708664/1‑69 (MQ=255)
                         ttAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCg                                                >  2:5249346/1‑76 (MQ=255)
                         ttAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCg                                                <  1:1180287/76‑1 (MQ=255)
                         ttAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCg                                                >  1:4095126/1‑76 (MQ=255)
                           aaaGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTGCgcc                                              >  1:1688417/1‑74 (MQ=255)
                            aaGATTTTCACTGTTACTTCAGTGATATGGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGc                                             <  1:5740725/76‑1 (MQ=255)
                              gATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTc                                                  >  1:2553653/1‑69 (MQ=255)
                              gATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTc                                                  <  2:2553653/69‑1 (MQ=255)
                                ttttCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGc                                         >  1:1644628/1‑76 (MQ=255)
                                ttttCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGc                                         <  1:3605953/76‑1 (MQ=255)
                                  ttCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCcg                                       <  2:565844/76‑1 (MQ=255)
                                  ttCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTACGGGCccgccg                                       >  2:4529313/1‑76 (MQ=255)
                                  ttCACTGTTACTTCAGGGATATAGCACAATGCCCCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGcca                                       <  2:5771205/76‑2 (MQ=255)
                                   tCACTGTTACTTCTGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCcgc                                      >  2:5721989/1‑76 (MQ=255)
                                   tCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCcgc                                      <  2:769874/76‑1 (MQ=255)
                                     aCTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCcgcgc                                    <  2:2081345/76‑1 (MQ=255)
                                      cTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCcgcgcg                                   >  2:3395095/1‑76 (MQ=255)
                                      cTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCcgcgcg                                   >  1:4113779/1‑76 (MQ=255)
                                      cTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCcgcgcg                                   <  2:1248475/76‑1 (MQ=255)
                                        gTTACTTCAGTGATATAGCACAATGCGCCTTCGTTGCCGCTGCTGATGTGACGAATGTCCGGGTCCGCCGCgcggc                                 <  1:723611/76‑1 (MQ=255)
                                        gTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCgcggc                                 >  1:4643777/1‑76 (MQ=255)
                                         ttACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCgcggcg                                >  1:5421789/1‑76 (MQ=255)
                                          tACTTCATTGATATAGCACAGTGCGCTTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGt                               <  1:321091/76‑1 (MQ=255)
                                          tACTTCAGTGATATAGCACAATGCCCCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGt                               <  2:3454385/76‑1 (MQ=255)
                                          tACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCgcggcgg                               >  2:1699834/1‑75 (MQ=255)
                                          tACTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGt                               <  1:3489016/76‑1 (MQ=255)
                                           aCTTCAGTGATATAGCGCAATGCGCCTTCGGTGCCGATGATGATGTGCCGACCGTCCGGTCCATCCTCTCGGCGTg                              <  1:4755959/76‑1 (MQ=255)
                                           aCTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCTGGGCCAGCCGCGCGGCGTg                              >  2:5031618/1‑76 (MQ=255)
                                           aCTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTg                              >  2:2153945/1‑76 (MQ=255)
                                           aCTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTg                              >  1:5530938/1‑76 (MQ=255)
                                           aCTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTg                              >  1:5189282/1‑76 (MQ=255)
                                           aCTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTg                              >  2:3394905/1‑76 (MQ=255)
                                           aCTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTg                              <  1:4751946/76‑1 (MQ=255)
                                            cTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTgg                             >  1:5502341/1‑76 (MQ=255)
                                            cTTCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGCCCCGCCGCGCGCCgggg                             >  2:44368/1‑76 (MQ=255)
                                             ttCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGc                            <  2:4637215/76‑1 (MQ=255)
                                              tCAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCa                           <  2:5210093/76‑1 (MQ=255)
                                               cAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCAc                          <  2:903325/76‑1 (MQ=255)
                                               cAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCAc                          <  2:3646282/76‑1 (MQ=255)
                                               cAGTGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCAc                          >  1:1109204/1‑76 (MQ=255)
                                                  tGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGtt                       <  1:5690941/76‑1 (MQ=255)
                                                  tGATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGtt                       >  1:4211125/1‑76 (MQ=255)
                                                   gATATAGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGGGGCACGttt                      >  1:2317050/1‑76 (MQ=255)
                                                     tataGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGttttt                    >  2:3679564/1‑76 (MQ=255)
                                                     tataGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGttttt                    >  1:4326593/1‑76 (MQ=255)
                                                     tataGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGttttt                    <  1:3037587/76‑1 (MQ=255)
                                                      ataGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTa                   <  2:2869129/76‑1 (MQ=255)
                                                      ataGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCACAAGCGCGGCGTGGCTCGtgttct                   >  2:676453/1‑74 (MQ=255)
                                                       taGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTaa                  >  2:2626157/1‑76 (MQ=255)
                                                       taGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTaa                  <  1:807794/76‑1 (MQ=255)
                                                        aGCACAATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAAt                 <  2:2256006/76‑1 (MQ=255)
                                                           acaATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCg              >  1:4428607/1‑76 (MQ=255)
                                                            caATGCACCTTCGTTGCCGATGATGATGTGACGACTGTCCGGGCCAGCCGCGCGGCGTGGCACGCTTTTAATGCgt             >  1:2370354/1‑76 (MQ=255)
                                                            caATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCcgcgcg                                   >  1:1357313/1‑54 (MQ=255)
                                                            caATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCcgcgcg                                   <  2:1357313/54‑1 (MQ=255)
                                                            caATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCgt             >  1:4283814/1‑76 (MQ=255)
                                                            caATGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCgt             >  1:5672643/1‑76 (MQ=255)
                                                             aaTGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCgtg            >  2:2952854/1‑76 (MQ=255)
                                                             aaTGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCgtg            >  1:2207509/1‑76 (MQ=255)
                                                             aaTGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGGTTTTAATGCgtg            >  1:197649/1‑76 (MQ=255)
                                                              aTGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCgtgt           >  2:4846302/1‑76 (MQ=255)
                                                               tGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGtgtg          >  1:1466095/1‑76 (MQ=255)
                                                               tGCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGtgtg          <  2:2469301/76‑1 (MQ=255)
                                                                gCACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGtgtgg         >  2:4440020/1‑75 (MQ=255)
                                                                 cgccTTCGTTGCCGATGATGATGTGACGAATGTCCGGCCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGTGTGAc        <  1:2909482/74‑1 (MQ=255)
                                                                 cACCTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCCCGTTTTTAATGCGTGTGAc        >  2:2573528/1‑76 (MQ=255)
                                                                    cTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGTGTGACGGt     <  2:1575745/76‑1 (MQ=255)
                                                                    cTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGTGTGACGGt     <  2:5740753/76‑1 (MQ=255)
                                                                    cTTCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGTGTGACGGt     <  1:5289234/76‑1 (MQ=255)
                                                                      tCGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCgtg            >  1:1699884/1‑67 (MQ=255)
                                                                       cGTTGCCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGTGTGACGGTGcc  >  1:3020051/1‑76 (MQ=255)
                                                                       cGTTGCTGATGATGATGTGACGAATGTCCGGTCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGTGTGACGGTGcc  <  1:3150286/76‑1 (MQ=255)
                                                                            |                                                                       
GAAGAGGTTGTTTTCCGGGGTGAATTTAAAGATTTTCACTGTTACTTCAGTGATATAGCACAATGCACCTTCGTTGCGATGATGATGTGACGAATGTCCGGGCCAGCCGCGCGGCGTGGCACGTTTTTAATGCGTGTGACGGTGCC  >  NC_000913/2896743‑2896888

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