Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,432,994 (C)6→7 coding (567/1125 nt) smf ← DNA recombination‑mediator A family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,432,9881.C95.7% 151.0 / 0.1 46coding (573/1125 nt)smfDNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (22/22);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

GCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATG  >  NC_000913/3432851‑3433104
                                                                                                                                          |                                                                                                                    
gCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcgc                                                                                                                     <  1:424112/139‑1 (MQ=255)
     cACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGcccccc                                                                                                                <  2:267495/139‑1 (MQ=255)
      aCCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccccccc                                                                                                               <  1:673564/139‑2 (MQ=255)
      aCCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATAGAGGGGAAATTCCGAGACGAGAACGccacccc                                                                                                               >  1:573532/1‑138 (MQ=255)
       gcAGCGAAACACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccccccct                                                                                                              >  2:580874/2‑137 (MQ=255)
                  cTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccccccctgt                                                                                                            <  2:893459/130‑5 (MQ=255)
                  cTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccccccctgt                                                                                                            >  1:893459/1‑126 (MQ=255)
                         aaGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAg                                                                                             <  1:333024/139‑1 (MQ=255)
                         aaGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAg                                                                                             <  2:478187/139‑1 (MQ=255)
                          aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGc                                                                                            >  2:363013/1‑139 (MQ=255)
                                            tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcgg                                                                         >  1:767976/1‑139 (MQ=255)
                                                   ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                   >  2:342083/1‑139 (MQ=255)
                                                   ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                   >  2:531969/1‑139 (MQ=255)
                                                     aCTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 <  1:342083/139‑1 (MQ=255)
                                                      cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGATTGGCAGCCAGTCGGGCATGa                                                                               >  2:608679/1‑124 (MQ=255)
                                                      cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGATTGGCAGCCAGTCGGGCATGa                                                                               <  1:608679/124‑1 (MQ=255)
                                                       ttAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATgg                                                               >  2:362595/1‑139 (MQ=255)
                                                          gACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAt                                                            >  2:48662/1‑139 (MQ=255)
                                                           aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                           >  2:260300/1‑139 (MQ=255)
                                                           aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                           >  2:562220/1‑139 (MQ=255)
                                                           aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                           >  1:364429/1‑139 (MQ=255)
                                                           aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                           >  1:347224/1‑139 (MQ=255)
                                                           aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                           >  2:704297/1‑139 (MQ=255)
                                                           aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                           >  1:290921/1‑139 (MQ=255)
                                                                tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                      <  1:421725/139‑1 (MQ=255)
                                                                 gATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    <  2:767976/139‑1 (MQ=255)
                                                                     aTGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                      <  1:46547/134‑1 (MQ=255)
                                                                     aTGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                      >  2:46547/1‑134 (MQ=255)
                                                                     aTGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTc                                                 <  1:260300/139‑1 (MQ=255)
                                                                       gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGc                                                                                               <  2:926733/91‑1 (MQ=255)
                                                                       gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGc                                                                                               >  1:926733/1‑91 (MQ=255)
                                                                          aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                            <  2:347224/139‑1 (MQ=255)
                                                                          aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                            <  2:383309/139‑1 (MQ=255)
                                                                            ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGAAGAGAGCGCCCCCCCTGTTCAAGCAGACTGTCAGCCAGTCGGGCATGACGGCGGGGATGAACGGTATTAAGTCCATTCCCCAATa                                          >  2:612533/1‑139 (MQ=255)
                                                                               ttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAg                                       <  1:363013/139‑1 (MQ=255)
                                                                                 cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCa                                     <  1:531969/139‑1 (MQ=255)
                                                                                        aaTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                              <  1:562220/139‑1 (MQ=255)
                                                                                          ttGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                            <  1:431909/139‑1 (MQ=255)
                                                                                               aGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTTCATTCCCCAATACAGCAATGCTGACGCCAtt                       >  1:322762/1‑139 (MQ=255)
                                                                                                     gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 <  2:296495/91‑1 (MQ=255)
                                                                                                     gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                 >  1:296495/1‑91 (MQ=255)
                                                                                                     gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    >  2:329961/1‑104 (MQ=255)
                                                                                                     gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCa                                                    <  1:329961/104‑1 (MQ=255)
                                                                                                                cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                      >  1:817229/1‑123 (MQ=255)
                                                                                                                cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                      <  2:817229/123‑1 (MQ=255)
                                                                                                                cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCtt      >  2:289859/1‑139 (MQ=255)
                                                                                                                cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTGGGCATGACGGCGGGGATGAAt                                                                 <  2:498461/79‑1 (MQ=255)
                                                                                                                cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTGGGCATGACGGCGGGGATGAAt                                                                 >  1:498461/1‑79 (MQ=255)
                                                                                                                 gAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCttt     <  1:704297/139‑1 (MQ=255)
                                                                                                                  aGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTACGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTa    <  2:573532/139‑1 (MQ=255)
                                                                                                                    gggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTTCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATg  <  2:322762/139‑1 (MQ=255)
                                                                                                                    gggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                            <  2:124420/97‑1 (MQ=255)
                                                                                                                    gggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                            >  1:124420/1‑97 (MQ=255)
                                                                                                                                          |                                                                                                                    
GCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATG  >  NC_000913/3432851‑3433104

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