Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,545,698 C→A 8.9% I25I (ATC→ATA)  nfuA → Fe/S biogenesis protein, putative scaffold/chaperone protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,545,6980CA 8.9% 148.1 / 5.4 45I25I (ATC→ATA) nfuAFe/S biogenesis protein, putative scaffold/chaperone protein
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base A (2/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

TCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCACCATACCTGGAAG  >  NC_000913/3545553‑3545846
                                                                                                                                                 |                                                                                                                                                    
tCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCcgc                                                                                                                                                   >  1:170581/1‑149 (MQ=255)
  gAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGt                                                                                                                                                 >  1:250053/1‑149 (MQ=255)
   aaCCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGt                                                                                                                                                 >  3:115591/1‑148 (MQ=255)
     ccTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGAAACAAATAcgcggat                                                                                                                                               >  4:322351/1‑145 (MQ=255)
        tGTAGAGCCTCGATGATGTGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTtg                                                                                                                                            >  3:319439/1‑148 (MQ=255)
        tGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTtg                                                                                                                                            <  4:319439/148‑1 (MQ=255)
          tAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGAAACAAATACGCGGATTtgtg                                                                                                                                          <  3:322351/148‑1 (MQ=255)
          tAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCCGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGa                                                                                                                                         <  1:359113/149‑1 (MQ=255)
            gagCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGAt                                                                                                                                        >  3:288450/1‑148 (MQ=255)
              gCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTaa                                                                                                                                     >  1:409297/1‑149 (MQ=255)
               ccTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAAc                                                                                                                                    >  1:184881/1‑149 (MQ=255)
                 tCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAAccc                                                                                                                                  >  1:122687/1‑149 (MQ=255)
                    atgatgGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTgg                                                                                                                               <  2:184881/149‑1 (MQ=255)
                            cGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCt                                                                                                                        >  3:379116/1‑148 (MQ=255)
                                ttatAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGc                                                                                                                   >  2:304115/1‑149 (MQ=255)
                                     aCCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTTCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTg                                                                                                                                <  3:382159/131‑1 (MQ=255)
                                     aCCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAAGCCTg                                                                                                                                >  4:382159/1‑131 (MQ=255)
                                         aCTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGa                                                                                                                >  3:376752/1‑143 (MQ=255)
                                         aCTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGa                                                                                                                <  4:376752/143‑1 (MQ=255)
                                          cTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGc                                                                                                         <  1:125432/149‑1 (MQ=255)
                                               aTAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGttt                                                                                                     >  3:14432/1‑148 (MQ=255)
                                                  gTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCtt                                                                                                  <  4:379116/148‑1 (MQ=255)
                                                  gTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTa                                                                                                 <  2:170581/149‑1 (MQ=255)
                                                     aaCTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGaa                                                                                                               <  1:113717/132‑1 (MQ=255)
                                                     aaCTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGaa                                                                                                               >  2:113717/1‑132 (MQ=255)
                                                     aaCTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTg                                                                                              >  1:251623/1‑149 (MQ=255)
                                                                  ttACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCt                                                                                 <  1:304115/149‑1 (MQ=255)
                                                                    aCTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCtg                                                                                <  4:284509/148‑1 (MQ=255)
                                                                           tCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGc                                                                         >  4:392684/1‑148 (MQ=255)
                                                                              gTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCAc                                                                      >  3:121330/1‑148 (MQ=255)
                                                                                       aTGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGcc                                                            <  2:146330/149‑1 (MQ=255)
                                                                                        tgctgcACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGcc                                                            <  3:392684/148‑1 (MQ=255)
                                                                                                 aGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGACCTGAAAttt                                                   <  3:261506/148‑1 (MQ=255)
                                                                                                  gcgcACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCTCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCGCTGAAATTTg                                                  <  4:123469/148‑1 (MQ=255)
                                                                                                  gcgcACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGa                                                 <  1:226271/149‑1 (MQ=255)
                                                                                                              aaaCTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACcg                                      >  3:215450/1‑148 (MQ=255)
                                                                                                                          aaTCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGATCGCGTATGTtgatga                         <  2:239960/149‑1 (MQ=255)
                                                                                                                          aaTCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTtgatga                         >  1:275376/1‑149 (MQ=255)
                                                                                                                                         tcacaAATACGCGGATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAGGCGCACCAt           >  4:81313/2‑148 (MQ=255)
                                                                                                                                          cacaAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGacac                                                                >  4:96622/1‑94 (MQ=255)
                                                                                                                                          cacaAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGacac                                                                <  3:96622/94‑1 (MQ=255)
                                                                                                                                            caAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAg                   <  4:107023/137‑1 (MQ=255)
                                                                                                                                            caAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAg                   >  3:107023/1‑137 (MQ=255)
                                                                                                                                            caAATACGCGGATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAGGCGCACCATAcc        <  3:81313/148‑1 (MQ=255)
                                                                                                                                                tCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCACCATACCTGGaa   <  2:250053/149‑1 (MQ=255)
                                                                                                                                                 ccGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCACCATACCTGGAAg  >  1:17630/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCCGCGTATTTGTGATTAACCCTGGCACGCCTAACGCTGAATGTGGCGTTTCTTATTGTCCGCCGGACGCTGTGGAAGCCACCGACACAGCCCTGAAATTTGACCTGCTGACCGCGTATGTTGATGAGTTAAGCGCACCATACCTGGAAG  >  NC_000913/3545553‑3545846

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