Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,786,746 Δ1 bp pseudogene (350/2253 nt) ygaQ → uncharacterized protein, putative enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,786,7460C.98.2% 230.3 / ‑2.5 57pseudogene (350/2253 nt)ygaQuncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (27/29);  total (27/30)
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

TCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGACTAT  >  NC_000913/2786612‑2786876
                                                                                                                                      |                                                                                                                                  
tCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGa                                                                                                                               >  1:794472/1‑139 (MQ=255)
 cTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGAt                                                                                                                              >  2:383521/1‑139 (MQ=255)
           aCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGt                                                                                                                    <  2:794472/139‑1 (MQ=255)
              aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGAt                                                                                                                 >  2:437982/1‑139 (MQ=255)
                   aaaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                            <  2:92023/139‑1 (MQ=255)
                    aaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                                          <  2:1002542/139‑1 (MQ=255)
                     aaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                                          <  2:946368/139‑1 (MQ=255)
                     aaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCg                                                                                                          <  2:803480/139‑1 (MQ=255)
                             tttCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                                  <  1:362629/139‑1 (MQ=255)
                                      aTACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGt                                                                                              <  2:211255/134‑1 (MQ=255)
                                      aTACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGt                                                                                              >  1:211255/1‑134 (MQ=255)
                                           aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGAt                                                                                                                              <  1:925797/97‑1 (MQ=255)
                                           aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGAt                                                                                                                              >  2:925797/1‑97 (MQ=255)
                                                     agTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTc                                                                          >  2:875589/1‑139 (MQ=255)
                                                      gTTGAATAAAGCAGCACAGCGTAGCATTACCGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCgg                                                                                                         <  1:873603/107‑1 (MQ=255)
                                                      gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCgg                                                                                                         >  2:873603/1‑107 (MQ=255)
                                                      gTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCa                                                                         >  2:139442/1‑139 (MQ=255)
                                                       ttGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                                        <  1:383521/139‑1 (MQ=255)
                                                        tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       <  1:139442/139‑1 (MQ=255)
                                                        tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       <  2:780923/139‑1 (MQ=255)
                                                        tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       <  1:437982/139‑1 (MQ=255)
                                                         gAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAccc                                                                      >  2:425381/1‑139 (MQ=255)
                                                                 cagcaCCGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata                                                              <  1:946012/139‑1 (MQ=255)
                                                                 cagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtat                                                               <  2:237304/138‑1 (MQ=255)
                                                                 cagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtat                                                               >  1:237304/1‑138 (MQ=255)
                                                                      caGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGa                                                         >  1:894992/1‑139 (MQ=255)
                                                                         cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTTATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       <  1:863752/122‑1 (MQ=255)
                                                                         cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       >  2:863752/1‑122 (MQ=255)
                                                                            aGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       >  1:305816/1‑119 (MQ=255)
                                                                            aGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       <  2:305816/119‑1 (MQ=255)
                                                                            aGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAccc                                                   <  1:624944/139‑1 (MQ=255)
                                                                                   cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtat                                                               <  2:754696/120‑1 (MQ=255)
                                                                                   cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtat                                                               >  1:754696/1‑120 (MQ=255)
                                                                                   cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGAc                                                        >  2:256915/1‑127 (MQ=255)
                                                                                   cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                                            >  2:105765/1‑139 (MQ=255)
                                                                                   cAGCGACAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGAc                                                        <  1:256915/127‑1 (MQ=255)
                                                                                     gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaaa                                          <  2:199231/139‑1 (MQ=255)
                                                                                        ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGa                                       >  2:1057730/1‑139 (MQ=255)
                                                                                                     gagaAAGAAAAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCa                          >  2:138217/1‑139 (MQ=255)
                                                                                                        aaagaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCatat                       <  1:838587/139‑1 (MQ=255)
                                                                                                          agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                            >  1:971849/1‑52 (MQ=255)
                                                                                                          agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAg                                                                                                            <  2:971849/52‑1 (MQ=255)
                                                                                                          agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       >  2:990964/1‑89 (MQ=255)
                                                                                                          agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       <  2:13661/89‑1 (MQ=255)
                                                                                                          agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       <  1:990964/89‑1 (MQ=255)
                                                                                                          agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                                       >  1:13661/1‑89 (MQ=255)
                                                                                                          agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc                     >  2:1178828/1‑139 (MQ=255)
                                                                                                          agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc                     >  1:367112/1‑139 (MQ=255)
                                                                                                          agaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc                     >  1:251703/1‑139 (MQ=255)
                                                                                                           gaaaAACATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCt                    >  1:1089551/1‑139 (MQ=255)
                                                                                                             aaaaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGACAAATATGTTGCATATGCTAt                  <  1:425381/139‑1 (MQ=255)
                                                                                                               aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTa                                                                                             >  2:1115815/1‑62 (MQ=255)
                                                                                                               aaCATAAGTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTa                                                                                             <  1:1115815/62‑1 (MQ=255)
                                                                                                                      gTTTATCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTcc         <  1:1178828/139‑1 (MQ=255)
                                                                                                                           tCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata                        >  2:144597/1‑119 (MQ=255)
                                                                                                                           tCAACGAGATAAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCata                        <  1:144597/119‑1 (MQ=255)
                                                                                                                             aaCGAGATACATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCCGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGACTAt  <  1:720913/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
TCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTAACATTCCTGACTAT  >  NC_000913/2786612‑2786876

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: