Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,608,413 T→C T204A (ACG→GCG)  yneE ← bestrophin family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,608,4130TC98.0% 184.9 / NA 51T204A (ACG→GCG) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/25);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATCCAGCGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGG  >  NC_000913/1608265‑1608552
                                                                                                                                                    |                                                                                                                                           
cAATCCAGCGAAATAAAAGTGTAAGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAACACAGATAAACCGTACGATGCAAACTCAGGGTCTAGGCCAAGGGAATTGGCGtt                                                                                                                                            <  1:600957/150‑2 (MQ=255)
          aaaTAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATgcg                                                                                                                                 >  1:721243/1‑151 (MQ=255)
            aTAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATgagc                                                                                                                                >  1:349386/1‑150 (MQ=255)
            aTAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCt                                                                                                                               >  2:500185/1‑151 (MQ=255)
                   tgtAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATcc                                                                                                                        >  2:579430/1‑151 (MQ=255)
                    gtAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCt                                                                                                                       <  1:154233/151‑1 (MQ=255)
                     tAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTg                                                                                                                      >  1:299110/1‑151 (MQ=255)
                        gAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCa                                                                                                                   <  1:279583/151‑1 (MQ=255)
                        gAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGCCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCa                                                                                                                   <  2:349673/151‑1 (MQ=255)
                          aaTCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAgg                                                                                                                 <  2:151479/151‑1 (MQ=255)
                             cAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGAcc                                                                                                              <  1:500133/151‑1 (MQ=255)
                             cAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGAc                                                                                                               <  2:216999/150‑1 (MQ=255)
                              aGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCg                                                                                                             >  1:397806/1‑151 (MQ=255)
                                cacaGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCt                                                                                                           <  1:63995/151‑1 (MQ=255)
                                   agagATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGa                                                                                                         >  2:414642/1‑150 (MQ=255)
                                           aaGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCAtt                                                                                                <  2:387908/151‑1 (MQ=255)
                                                gTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGcc                                                                                           <  2:390052/151‑1 (MQ=255)
                                                  cATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCg                                                                                          <  2:296392/150‑1 (MQ=255)
                                                       aaTGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGtt                                                                                    >  1:348016/1‑151 (MQ=255)
                                                       aaTGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGtt                                                                                    <  1:378294/151‑1 (MQ=255)
                                                          gCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTg                                                                                   >  1:200376/1‑149 (MQ=255)
                                                            aGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGc                                                                               <  1:647070/151‑1 (MQ=255)
                                                             ggTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCt                                                                              <  1:252512/151‑1 (MQ=255)
                                                                   cGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATaaa                                                                        <  2:758434/151‑1 (MQ=255)
                                                                    gACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAAc                                                                       >  1:719805/1‑151 (MQ=255)
                                                                       cAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGTCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAgg                                                                    <  1:581635/151‑1 (MQ=255)
                                                                          cgcgAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGCGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGat                                                                  >  1:64942/1‑150 (MQ=255)
                                                                                  gcagcaTGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAg                                                         <  2:23727/151‑1 (MQ=255)
                                                                                      caTGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGt                                                     >  1:687077/1‑151 (MQ=255)
                                                                                         gATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTcc                                                   >  1:270676/1‑150 (MQ=255)
                                                                                         gATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCa                                                  <  1:498644/151‑1 (MQ=255)
                                                                                          aTACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCAt                                                 <  1:309114/151‑1 (MQ=255)
                                                                                                   aGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCt                                        <  1:273976/151‑1 (MQ=255)
                                                                                                    gATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTg                                       <  1:296099/151‑1 (MQ=255)
                                                                                                       aaaCGGTACGATGCAAAATCAGGGTGTAGGCAAAAGGAATTGGCGCATAGGCAATGCACTCGCATCCTGCCAGGACCGCTAAAATATAATTAAGCCGCTCGTTGAGGCTAATAAACAGGATATCTGAAACCTGTCCATTGCGGCACTGa                                      >  1:440501/1‑149 (MQ=255)
                                                                                                        aaCGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAAcc                                   >  2:514061/1‑151 (MQ=255)
                                                                                                         aCGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCg                                  >  1:220140/1‑151 (MQ=255)
                                                                                                         aCGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCg                                  <  1:497215/151‑1 (MQ=255)
                                                                                                           ggTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGcc                                >  2:404875/1‑151 (MQ=255)
                                                                                                              aCGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAAc                             >  1:98357/1‑151 (MQ=255)
                                                                                                                 aTGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCAt                          <  2:391426/151‑1 (MQ=255)
                                                                                                                    cAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATtct                       >  1:648723/1‑151 (MQ=255)
                                                                                                                     aaaaTCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATtctc                      >  1:704986/1‑151 (MQ=255)
                                                                                                                       aaTCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTccc                    >  2:421132/1‑151 (MQ=255)
                                                                                                                       aaTCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTccc                    <  2:576343/151‑1 (MQ=255)
                                                                                                                         tCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACTGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCat                  >  2:607407/1‑151 (MQ=255)
                                                                                                                           aGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCatta                <  2:242862/151‑1 (MQ=255)
                                                                                                                           aGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCatt                 >  1:306415/1‑150 (MQ=255)
                                                                                                                            gggTGTAGGAAAAGGGAATTGGCGTGTAGGCTATGCGCCCGCATCCTGACAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGCCCAATAAACACGATATCTCAAAGCTGGCCATTGCGCCATTGACTCACCACCCATTCGcccgatat               <  1:52890/151‑6 (MQ=255)
                                                                                                                                  aGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACa          >  1:327814/1‑150 (MQ=255)
                                                                                                                                      aaaGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGAt      >  2:619563/1‑150 (MQ=255)
                                                                                                                                         gggAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACgg  <  2:538670/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CAATCCAGCGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGG  >  NC_000913/1608265‑1608552

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