Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,787,427 G→A 7.1% pseudogene (1031/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,787,4270GA 7.1% 171.3 / 2.9 56pseudogene (1031/2253 nt)ygaQuncharacterized protein; putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base A (2/2);  total (26/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

AAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGG  >  NC_000913/2787279‑2787567
                                                                                                                                                    |                                                                                                                                            
aaaGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTgg                                                                                                                                              >  5:58013/1‑149 (MQ=255)
       cTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTa                                                                                                                                        <  3:79486/148‑1 (MQ=255)
        tAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAg                                                                                                                                       >  8:207048/1‑148 (MQ=255)
        tAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAg                                                                                                                                       >  8:173675/1‑148 (MQ=255)
           aaaaTTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAg                                                                                                                                   <  6:58013/149‑1 (MQ=255)
            aaaTTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAgg                                                                                                                                  >  7:367008/1‑149 (MQ=255)
                    tATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGttt                                                                                                                          >  1:216723/1‑149 (MQ=255)
                        gATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGa                                                                                                                       <  5:348059/148‑1 (MQ=255)
                         atatCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGAtt                                                                                                                     <  3:142601/149‑1 (MQ=255)
                          tatCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTc                                                                                                                    <  7:173675/149‑1 (MQ=255)
                            tCATTCCTGATTATACCTTTAACTTTATTGGAATTGTGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGAGAAGATCAGTATAGATATAGAAGGTGGTATCCATGGCGATTAGCAAGGTCAGCGTTTGCTTCca                                                                                                                  >  2:353301/1‑149 (MQ=255)
                            tCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCca                                                                                                                  >  7:442094/1‑149 (MQ=255)
                            tCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCca                                                                                                                  <  8:235710/149‑1 (MQ=255)
                               ttCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTATAGATATAGAAGGTGGTATCCCTGACTATTAGCAAGGTCAGGTTTTGAGTCcatt                                                                                                                >  8:308614/1‑148 (MQ=255)
                               gtCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGCCCTGGATTAACCATCTGATGATATACGAGCGAATATCAGTAAAGATATAGAATGTGTTATCCCTGAATATTAGCAAGTTCATGTTTTGATTAcata                                                                                                                >  8:308612/2‑147 (MQ=255)
                                  cTGATTATTCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcattcatt                                                                                                            <  3:143937/149‑1 (MQ=255)
                                  cTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcattcat                                                                                                             >  7:103256/1‑148 (MQ=255)
                                   tGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcattcatt                                                                                                            >  4:310136/1‑148 (MQ=255)
                                     aTTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTaa                                                                                                          <  7:207048/148‑1 (MQ=255)
                                          aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATc                                                                                                     >  2:249229/1‑148 (MQ=255)
                                           ccTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCaa                                                                                                                                    <  7:377277/116‑1 (MQ=255)
                                           ccTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCaa                                                                                                                                    >  8:377277/1‑116 (MQ=255)
                                           ccTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcattc                                                                                                               >  7:98711/1‑137 (MQ=255)
                                           ccTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCcattc                                                                                                               <  8:98711/137‑1 (MQ=255)
                                           ccTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCa                                                                                                   <  8:225934/149‑1 (MQ=255)
                                           ccTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCa                                                                                                   <  3:310136/149‑1 (MQ=255)
                                            cTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                                  >  7:157932/1‑149 (MQ=255)
                                             tttAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                                  >  4:206164/1‑148 (MQ=255)
                                             tttAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                                  >  3:339510/1‑148 (MQ=255)
                                             tttAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCCGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                                  >  1:37576/1‑148 (MQ=255)
                                                             ttGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTaaa                                                                                 >  3:63401/1‑149 (MQ=255)
                                                                  gttAGTGGTAAAAATGAACTGGATTATCCCTCTGCTGATATACTAGCGAAGATCAGTAAAGATATAGAATGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                            <  5:411502/147‑1 (MQ=255)
                                                                      gTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccat                                                                        <  3:393983/149‑1 (MQ=255)
                                                                       tGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAATATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                       <  1:249229/149‑1 (MQ=255)
                                                                        ggTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                       <  8:442094/148‑1 (MQ=255)
                                                                         gTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTTGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                     <  7:170344/149‑1 (MQ=255)
                                                                          tAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                     >  8:204177/1‑148 (MQ=255)
                                                                              aaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaaga                                                                <  6:354436/149‑1 (MQ=255)
                                                                              aaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaaga                                                                <  4:339510/149‑1 (MQ=255)
                                                                                   ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATAc                                                           <  2:90209/149‑1 (MQ=255)
                                                                                     cggATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAt                                                         <  2:37576/148‑1 (MQ=255)
                                                                                             cccTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTg                                                  >  4:273916/1‑148 (MQ=255)
                                                                                               ctctGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGata                                                              <  7:132598/134‑1 (MQ=255)
                                                                                               ctctGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGata                                                              >  8:132598/1‑134 (MQ=255)
                                                                                                    ctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaaa                                             <  4:63401/146‑1 (MQ=255)
                                                                                                       atatACGAGCGAATATCAGTAAAGATATATAAGGTGGTATCCCTGGCTTTTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaag                                             <  2:138233/143‑2 (MQ=255)
                                                                                                       atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAATGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACATAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACtt                                        <  8:367008/148‑1 (MQ=255)
                                                                                                       atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATATCACTGCaaag                                             >  1:138233/1‑142 (MQ=255)
                                                                                                              aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCaa                                <  8:103256/149‑1 (MQ=255)
                                                                                                                cGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                       <  2:248413/108‑1 (MQ=255)
                                                                                                                cGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                       >  1:248413/1‑108 (MQ=255)
                                                                                                                          tAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACa                    >  3:120468/1‑149 (MQ=255)
                                                                                                                                     aaGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGAt         <  2:216723/149‑1 (MQ=255)
                                                                                                                                     aaGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGAt         >  6:372712/1‑149 (MQ=255)
                                                                                                                                            tATCCCTTACTATTAGCAAGTTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATggg  <  7:308614/149‑1 (MQ=255)
                                                                                                                                            tATCCCTGACTATTAGCAAGCTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACATAAGAACGTAAATAAATCAATCCAGAAGATTTACATCTCACTGCAAAAGAGTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATggg  <  7:308612/149‑1 (MQ=255)
                                                                                                                                             aTCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATggg  >  8:201962/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AAAGATGCTAAAAAATTTGGTATCGATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGGG  >  NC_000913/2787279‑2787567

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑