Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 1,712,449 A→G intergenic (+291/‑320) nth → / → dtpA DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,712,4490AG98.3% 192.2 / NA 58intergenic (+291/‑320)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (28/29);  minor base C (0/1);  total (28/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

AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGA  >  NC_000913_3_pae_pgi/1712307‑1712587
                                                                                                                                              |                                                                                                                                          
aaTAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACt                                                                                                                                      <  2:715643/149‑1 (MQ=255)
aaTAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAAGTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACt                                                                                                                                      <  1:444746/149‑1 (MQ=255)
 aTAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTg                                                                                                                                     <  1:474548/149‑1 (MQ=255)
   aGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGa                                                                                                                                   <  1:1078107/149‑1 (MQ=255)
   aGCAAAATATTCTATTCATGCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCATGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGa                                                                                                                                   <  1:1167801/149‑1 (MQ=255)
      gAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTa                                                                                                                                >  1:619478/1‑149 (MQ=255)
          ttttCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTAtctc                                                                                                                            >  1:526216/1‑149 (MQ=255)
              cTATTCATTCGTGATCAAGATCACGATGCCTTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTa                                                                                                                        <  1:361688/149‑1 (MQ=255)
                 ttcattcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAGGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTa                                                                                                                     >  1:926138/1‑149 (MQ=255)
                     ttcGTTATCAAGATCACGATGCTGTAGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTaaaaa                                                                                                                 <  1:1190772/149‑1 (MQ=255)
                     ttcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTaaaaa                                                                                                                 <  2:134367/149‑1 (MQ=255)
                        gTGATCAAGCTCACGATCCTGTCGCAATGTAAAATCATTGCAAGAATCAGGTTAACATCATTTGAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGc                                                                                                              <  1:1171730/149‑1 (MQ=255)
                        gTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTaaaaa                                                                                                                 >  1:327508/1‑146 (MQ=255)
                        gTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTaaaaa                                                                                                                 <  2:327508/146‑1 (MQ=255)
                        gTGATCAAGAACACTATGCTGTAGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGc                                                                                                              <  2:926138/149‑1 (MQ=255)
                         tGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCt                                                                                                             <  1:119062/149‑1 (MQ=255)
                          gATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTa                                                                                                            >  1:612107/1‑149 (MQ=255)
                            tCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTaaa                                                                                                          <  2:370165/149‑1 (MQ=255)
                                 aTCAGGATGCTGTCGCAAGGCAAAAACATTGCAAGAATCAGGTTCAGGTGAGTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTATACTCGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCa                                                                                                     <  1:506974/149‑1 (MQ=255)
                                      gATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACa                                                                                                >  2:1220891/1‑149 (MQ=255)
                                      gATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACa                                                                                                >  1:814220/1‑149 (MQ=255)
                                               gCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGt                                                                                       >  1:871665/1‑149 (MQ=255)
                                                  aTGTAAAAACATTGCAAGAAGCATGCTAACATCATTTAAATTTACACTCGCAACAAATATTGACCTACAAAACAGTACACTGGCGATTTTTCGGAAACTGGACTAGCTGACTAGTAAAAACGCTAAATAGCAGAACATATCGCCgaact                                                                                    <  2:162519/149‑5 (MQ=255)
                                                     tAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACACTAAATAGCAGAACCTATCGCCGTTAAGcc                                                                                 >  2:636011/1‑149 (MQ=255)
                                                       aaaaCATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAt                                                                               <  1:1220891/149‑1 (MQ=255)
                                                       aaaaCAATGCAAGAATCAGGTTAAAATCAATTACATTCAACCAACAAACCCACCCAAGCGAAAAAAACATCACACCACCTATTTTCCCTTAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAt                                                                               <  2:937999/149‑1 (MQ=255)
                                                          aCATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAc                                                                            <  2:612107/149‑1 (MQ=255)
                                                           cATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcc                                                                           <  1:400693/149‑1 (MQ=255)
                                                                 aaGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCCTCAcccccccc                                                                     >  2:830173/1‑147 (MQ=255)
                                                                 aaGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAAAGCTAAATAGCAGAACATATGCCAGTTATCCCATCACCCCCctc                                                                     >  2:700023/1‑149 (MQ=255)
                                                                  aGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGGACATATTGCCGTTAAGCCATCATCCCCcttt                                                                    >  1:763683/1‑147 (MQ=255)
                                                                    aaTCAGTTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTAGTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTca                                                                  <  2:1095043/149‑1 (MQ=255)
                                                                    aaTCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTca                                                                  <  2:619478/149‑1 (MQ=255)
                                                                          gTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTcagtgcag                                                            <  1:845779/149‑1 (MQ=255)
                                                                          gTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTcagtgcag                                                            <  2:1069190/149‑1 (MQ=255)
                                                                           ttAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTcagtgcagt                                                           <  2:1181330/149‑1 (MQ=255)
                                                                             aaaaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATAACCCCCCTCTCAGTGCAGTGa                                                         >  2:683993/1‑149 (MQ=255)
                                                                                    tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGAAGTGAAAAAATc                                                  >  2:163660/1‑149 (MQ=255)
                                                                                        aaTTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAGATAGCAGAACATATCACCCTTAAGCCATCACCCCCCACTCAGTGCAGTGAACAAATATGcc                                              >  2:221580/1‑149 (MQ=255)
                                                                                        aaTTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGcc                                              <  2:599187/149‑1 (MQ=255)
                                                                                          tttACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                        <  2:1217219/121‑1 (MQ=255)
                                                                                          tttACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                        >  1:1217219/1‑121 (MQ=255)
                                                                                           ttACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATAGCCGTTAAGCCATCACCCCCCTACCAGTGCACTGAAACAATCTAccctc                                           >  2:1174468/1‑146 (MQ=255)
                                                                                                gacccaacaaAAATTGACCTACAAAACATTACAATGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACACGCATCGCAGTGTAGTGAAAAAATCTGCCGTTACGtt                                      <  2:285358/145‑1 (MQ=255)
                                                                                                  cGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAAGCTGCCGTTACGtgtt                                    >  1:274233/1‑149 (MQ=255)
                                                                                                   gCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATATCGCCAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                   >  1:309509/1‑149 (MQ=255)
                                                                                                     aacaaATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                 >  2:396811/1‑149 (MQ=255)
                                                                                                      acaaATATTGACCTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGa                                >  2:1174381/1‑149 (MQ=255)
                                                                                                                  cTACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                    >  2:1192822/1‑149 (MQ=255)
                                                                                                                   tACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                   >  1:626650/1‑149 (MQ=255)
                                                                                                                   tACAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                   <  1:1267894/149‑1 (MQ=255)
                                                                                                                     cAAAACATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTgg                 >  2:36142/1‑149 (MQ=255)
                                                                                                                       aaaCATTACACTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAGCATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGCAAAATTTAACGCTGGAt               >  1:702632/1‑149 (MQ=255)
                                                                                                                               cacTGGCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTc       >  1:46688/1‑149 (MQ=255)
                                                                                                                                   ggCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGCACTGGTCTCCTGTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg   <  1:636011/149‑1 (MQ=255)
                                                                                                                                   ggCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAAc            <  2:591259/140‑1 (MQ=255)
                                                                                                                                   ggCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAAc            >  1:591259/1‑140 (MQ=255)
                                                                                                                                    gCTATTTTTCGGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGa  >  1:734163/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                          
AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGA  >  NC_000913_3_pae_pgi/1712307‑1712587

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