Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,467,061 C→T pseudogene (1670/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,467,0610CT97.2% 117.8 / NA 36pseudogene (1670/2513 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/23);  minor base G (1/0);  total (13/23)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAA  >  NC_000913/1466929‑1467193
                                                                                                                                    |                                                                                                                                    
tCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATgg                                                                                                                                <  2:388069/139‑1 (MQ=255)
tCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATgg                                                                                                                                <  1:404750/139‑1 (MQ=255)
 cAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGgt                                                                                                                               <  2:162026/139‑1 (MQ=255)
         acgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGaa                                                                                                                       >  2:314716/1‑139 (MQ=255)
                 aCCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCAt                                                                                                                                  <  1:203333/120‑1 (MQ=255)
                 aCCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCAt                                                                                                                                  >  2:203333/1‑120 (MQ=255)
                 aCCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTg                                                                                                               >  2:218627/1‑139 (MQ=255)
                            cAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACg                                                                                                    >  2:45074/1‑139 (MQ=255)
                               tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGtt                                                                                                 <  2:28971/139‑1 (MQ=255)
                                     aaCCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaacaac                                                                                           <  1:94243/139‑1 (MQ=255)
                                         aCTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAg                                                                                       >  2:231731/1‑139 (MQ=255)
                                         aCTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAg                                                                                       >  2:470716/1‑139 (MQ=255)
                                          cTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAgg                                                                                      <  2:218970/139‑1 (MQ=255)
                                           tGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCCTTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGGGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGga                                                                                     >  1:325622/1‑139 (MQ=255)
                                                 tGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAtt                                                                               <  2:156793/139‑1 (MQ=255)
                                                         aTTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc                                                                       <  1:314716/139‑1 (MQ=255)
                                                         aTTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc                                                                       <  2:325622/139‑1 (MQ=255)
                                                                               gggCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGATAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTg                                                 <  1:25653/139‑1 (MQ=255)
                                                                                               gTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAgg                                 >  2:187490/1‑139 (MQ=255)
                                                                                                  aaTCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAggggg                              <  1:470716/139‑1 (MQ=255)
                                                                                                  aaTCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAggggg                              <  1:45074/139‑1 (MQ=255)
                                                                                                  aaTCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAggggg                              <  1:231731/139‑1 (MQ=255)
                                                                                                  aaTCTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAggggg                              <  1:165634/139‑1 (MQ=255)
                                                                                                     cTGGATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCaa                           <  1:218627/139‑1 (MQ=255)
                                                                                                     cTGGATGGTAGCCTGACTGTTACAGGCGGTGTACATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCaa                           <  2:158807/139‑1 (MQ=255)
                                                                                                       ggATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACa                         <  1:186794/139‑1 (MQ=255)
                                                                                                        gATGGTAGCGTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAg                        >  1:512126/1‑139 (MQ=255)
                                                                                                        gATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGata                                                                                 <  2:15679/82‑1 (MQ=255)
                                                                                                        gATGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGata                                                                                 >  1:15679/1‑82 (MQ=255)
                                                                                                          tGGTAGCCTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTa                      <  2:150277/139‑1 (MQ=255)
                                                                                                                ccTGACTGTTACAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAAt                >  1:53364/1‑139 (MQ=255)
                                                                                                                          aCAGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATg      <  2:316111/139‑1 (MQ=255)
                                                                                                                            aGGCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGtt    <  2:11389/139‑1 (MQ=255)
                                                                                                                              gCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaa                                                                                               >  1:73213/1‑46 (MQ=255)
                                                                                                                              gCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaa                                                                                               <  2:73213/46‑1 (MQ=255)
                                                                                                                              gCGGTGTGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTaa  >  2:231438/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                    
TCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAA  >  NC_000913/1466929‑1467193

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