Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,001,421 2 bp→GC coding (381‑382/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,4210CG98.2% 308.2 / ‑4.5 113G128R (GGA→CGA) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base G (48/63);  total (50/63)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,001,4220GC99.1% 309.3 / ‑4.6 113D127E (GAC→GAGrarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (49/63);  total (50/63)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 114 total aligned reads displayed.

AGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCG  >  NC_000913/4001353‑4001487
                                                                     ||                                                                 
aGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAAtcgctcgcag                                                             >  2:1125437/1‑74 (MQ=255)
aGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCg                                                             >  1:5521801/1‑76 (MQ=255)
 ccccACAGCGGGACTACCCCGCCACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGa                                                            <  1:2472397/74‑1 (MQ=255)
  tCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAg                                                           >  2:4641265/1‑76 (MQ=255)
   ccACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAgg                                                          <  1:1234362/76‑1 (MQ=255)
   ccACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAgg                                                          >  1:4161841/1‑76 (MQ=255)
   ccACAGCTGAACAAACACCCCACATACCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAgg                                                          <  2:793644/76‑1 (MQ=255)
    cacaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGa                                                         >  1:2657242/1‑76 (MQ=255)
    cacaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGa                                                         >  2:4267189/1‑76 (MQ=255)
    cacaGCGGGACTACCACGCCACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGa                                                         <  1:2467994/76‑1 (MQ=255)
     ccaGCTGGACAAACCCCCCACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaa                                                        <  1:5317266/75‑1 (MQ=255)
     acaGCTGGACTAACACGCCACATATCGCCAGAACCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaa                                                        >  2:3519567/1‑76 (MQ=255)
      cggCTGGACTAACACGCACATACCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                      <  2:222186/74‑1 (MQ=255)
       aGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                      >  2:4538035/1‑76 (MQ=255)
        gCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAat                                                     >  1:5513778/1‑76 (MQ=255)
         cTGGACTAACACGCCACATATCGCCAGAAGCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatc                                                    >  2:3423180/1‑76 (MQ=255)
           gcaCTACCCCGCCACATATCCCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcat                                                  <  1:2736534/74‑1 (MQ=255)
            taCTAACCCGCCACAGAACGCGAGAACCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcatc                                                 <  2:4129919/75‑1 (MQ=255)
             aCTAACCCCCCACATATCGCCGCATTCTCCGCCAGCCATTGCATCCGGCGGAATCGCTCTCCGAGGACAATCATcc                                                <  1:3832350/76‑1 (MQ=255)
              cTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATccc                                               >  1:3646485/1‑76 (MQ=255)
               tAACACGCCATATATCGCCAGAATCACCGCCAGCCATTGCATACGGCGGAATCGCTCGCCGACGAAAATCATcccc                                              >  1:1153113/1‑76 (MQ=255)
               tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcccc                                              >  1:810870/1‑76 (MQ=255)
               aaaCACCCCACATATCCCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcccc                                              <  2:347529/75‑1 (MQ=255)
                   acGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGca                                          >  2:4432094/1‑76 (MQ=255)
                    ccccACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcac                                         <  2:3694183/74‑1 (MQ=255)
                     gCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcaca                                        >  1:3848543/1‑76 (MQ=255)
                      ccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                       >  1:1659406/1‑76 (MQ=255)
                       cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAAt                                      >  1:2723937/1‑76 (MQ=255)
                       cacaGATCGCAAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAAt                                      <  2:3942079/76‑1 (MQ=255)
                        atatatCGCTAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATg                                     <  1:2728149/74‑1 (MQ=255)
                        acaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATg                                     >  2:3852857/1‑76 (MQ=255)
                        acaTACCGCAAGATTCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATg                                     <  1:4290613/76‑1 (MQ=255)
                         caTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                    <  2:4161840/76‑1 (MQ=255)
                         caTATCGCCAGAATCACCACCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                    <  2:4054257/76‑1 (MQ=255)
                         caTATCGCCAAGACCCCCGCCAGCCATTGCATCCTGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                    <  1:3921415/76‑1 (MQ=255)
                         caTAGCTCCAGAACCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                    <  2:3513654/76‑1 (MQ=255)
                           tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTc                                  <  1:1562227/76‑1 (MQ=255)
                           tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTc                                  >  1:5547754/1‑76 (MQ=255)
                           tatCGCAAGAATAACCGACAGCCAGTGCATCCGGCGGAATCGCGCGCAGAGGAAAATCATCCCCAGCACAATggac                                  >  1:2351581/1‑73 (MQ=255)
                            atCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                 >  2:2144066/1‑76 (MQ=255)
                            atCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                 >  2:2916615/1‑76 (MQ=255)
                              cGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAcc                               <  1:4531187/76‑1 (MQ=255)
                                ccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                             >  2:4290477/1‑76 (MQ=255)
                                ccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                             >  1:758935/1‑76 (MQ=255)
                                 cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGtt                                   >  2:491497/1‑69 (MQ=255)
                                 cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGtt                                   <  1:491497/69‑1 (MQ=255)
                                 cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGc                            >  2:2626226/1‑76 (MQ=255)
                                 cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGc                            <  1:1195027/76‑1 (MQ=255)
                                 cAGAACCCCCGCCAGCCATTGCATCCGGCGGGATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGc                            <  2:2445111/76‑1 (MQ=255)
                                 cAGAACCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGc                            <  2:2948720/76‑1 (MQ=255)
                                 cAGAACCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGc                            <  1:1963660/76‑1 (MQ=255)
                                  aGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCg                           >  2:1939493/1‑76 (MQ=255)
                                  aGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCg                           >  1:2989938/1‑76 (MQ=255)
                                  aGAACCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCg                           <  1:4459784/76‑1 (MQ=255)
                                    aaccACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCggg                         <  1:3476187/73‑1 (MQ=255)
                                    aaTCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCggg                         <  2:3718935/76‑1 (MQ=255)
                                    aaTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCggg                         >  1:1084431/1‑76 (MQ=255)
                                    aaTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCggg                         >  1:4591997/1‑76 (MQ=255)
                                     aTCACCGGCAGCCATTGCAGCCGGCGGAAGCCCTCGCCGAGGAAAGTCAGACCCATCCCAATGTACACCAGCcggc                        >  2:4581146/1‑75 (MQ=255)
                                      ccACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGtt                       <  1:5751588/75‑1 (MQ=255)
                                       cgccGCCAGCGATTGCGTCCTGGTGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGTCAGCGGGTta                      <  2:3118687/74‑1 (MQ=255)
                                       cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                      <  2:2024763/76‑1 (MQ=255)
                                       cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                      <  1:2704920/76‑1 (MQ=255)
                                       cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                      >  2:1622259/1‑76 (MQ=255)
                                       cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                      <  1:3293485/76‑1 (MQ=255)
                                        cccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaa                     <  1:331392/75‑1 (MQ=255)
                                        aCCGCCAGACATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaa                     <  2:1832977/76‑1 (MQ=255)
                                         ccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaat                    >  1:1357452/1‑76 (MQ=255)
                                          cGCCAGCCATTGCATCCGGCGGAATCGCTCGCCCAGGAAAATCATCCCCAACACACTGACCACCAGCGGCCtaata                   >  2:5850548/1‑76 (MQ=255)
                                           gccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaataa                  <  2:2657242/76‑1 (MQ=255)
                                             cagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                >  2:3916838/1‑76 (MQ=255)
                                             cagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                >  2:4342367/1‑76 (MQ=255)
                                              agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAga               >  1:1090264/1‑75 (MQ=255)
                                                ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTaa             <  1:1875390/76‑1 (MQ=255)
                                                ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTaa             <  1:1900358/76‑1 (MQ=255)
                                                 caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAc            <  2:5417732/76‑1 (MQ=255)
                                                 caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAc            >  1:3422427/1‑76 (MQ=255)
                                                 caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAc            >  1:4806352/1‑76 (MQ=255)
                                                 caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAc            <  1:911035/76‑1 (MQ=255)
                                                  aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc           >  1:4728723/1‑76 (MQ=255)
                                                  aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc           >  1:259018/1‑76 (MQ=255)
                                                  aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc           <  2:3646485/76‑1 (MQ=255)
                                                   tttcATCCGGCGGAAGCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCa          <  2:2501979/73‑1 (MQ=255)
                                                    tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa         <  1:5604919/76‑1 (MQ=255)
                                                     gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg        <  1:5671972/76‑1 (MQ=255)
                                                     gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg        <  1:4865043/76‑1 (MQ=255)
                                                     gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg        <  2:2812660/76‑1 (MQ=255)
                                                     gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg        >  1:4185720/1‑76 (MQ=255)
                                                     gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg        <  1:2076264/76‑1 (MQ=255)
                                                      cATCTGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAgg       <  1:4061072/76‑1 (MQ=255)
                                                      cATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAgg       <  2:3320602/76‑1 (MQ=255)
                                                        tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt     >  2:2813677/1‑76 (MQ=255)
                                                         ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc    <  2:2448828/76‑1 (MQ=255)
                                                         ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc    <  2:5421021/76‑1 (MQ=255)
                                                         ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc    <  2:1352742/76‑1 (MQ=255)
                                                         ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc    <  1:682893/76‑1 (MQ=255)
                                                          cggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcg   >  1:3880702/1‑75 (MQ=255)
                                                          cggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcg   <  2:5547752/76‑2 (MQ=255)
                                                           ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgc  <  1:2559980/76‑3 (MQ=255)
                                                           ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgc  <  2:4438428/76‑3 (MQ=255)
                                                                     ||                                                                 
AGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCG  >  NC_000913/4001353‑4001487

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