Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,070,440 3 bp→GAA intergenic (+336/‑483) unknown → / → unknown Type II secretory pathway, component PulD/Bacterial non‑heme ferritin (EC 1.16.3.2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,070,4400TG100.0% 51.3 / NA 20intergenic (+336/‑485)unknown/unknownType II secretory pathway, component PulD/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (9/11);  total (9/11)
*11,070,4410TA100.0% 52.0 / NA 20intergenic (+337/‑484)unknown/unknownType II secretory pathway, component PulD/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (9/11);  total (9/11)
*11,070,4420TA100.0% 51.3 / NA 20intergenic (+338/‑483)unknown/unknownType II secretory pathway, component PulD/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (9/11);  total (9/11)

AAATTTCCCGCCCCTGCAAATACCCCTAAAGGCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATAACGTTTAACTAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACTTATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTC  >  1/1070219‑1070627
                                                                                                                                                                                                                             |||                                                                                                                                                                                         
aaaTTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGACTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatta                                                                                                                                                                <  1:358946/251‑1 (MQ=255)
   tttCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGt                                                                                                                                                             <  1:641107/251‑1 (MQ=255)
                   aTACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATgg                                                                                                                                             <  1:334688/251‑1 (MQ=255)
                      ccccTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatt                                                                                                                                                                 <  1:332040/228‑1 (MQ=255)
                      ccccTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTaa                                                                                                                                          >  1:814500/1‑251 (MQ=255)
                                            gctgATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatt                                                                                                                                                                 >  1:591573/1‑206 (MQ=255)
                                               gATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACatatat                                                                                                                                                                  <  1:873395/202‑1 (MQ=255)
                                               gATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatta                                                                                                                                                                <  1:882711/204‑1 (MQ=255)
                                                   aGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACatatat                                                                                                                                                                  >  1:413914/1‑198 (MQ=255)
                                                             aTACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                    <  1:732398/250‑1 (MQ=255)
                                                                cGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGa                                                                                                 >  1:684168/1‑250 (MQ=255)
                                                                                                          aaTCAATATTTATTGATATTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGt                                                                                                                                                             <  1:281235/148‑1 (MQ=255)
                                                                                                               atatTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGt                                                                                                                                                             >  1:796235/1‑143 (MQ=255)
                                                                                                                            ttACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACatatat                                                                                                                                                                  >  1:294688/1‑125 (MQ=255)
                                                                                                                            ttACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatt                                                                                                                                                                 <  1:659191/126‑1 (MQ=255)
                                                                                                                            ttACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                    >  1:661888/1‑187 (MQ=255)
                                                                                                                                   gCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAAcat                                                                                                                                                                      >  1:676356/1‑112 (MQ=255)
                                                                                                                                           tcttcAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTag                     <  1:352855/251‑1 (MQ=255)
                                                                                                                                            cttcAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTaga                    >  1:313791/1‑251 (MQ=255)
                                                                                                                                                              ttAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTtc  <  1:807161/251‑1 (MQ=255)
                                                                                                                                                                                                                             |||                                                                                                                                                                                         
AAATTTCCCGCCCCTGCAAATACCCCTAAAGGCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATAACGTTTAACTAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACTTATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTC  >  1/1070219‑1070627

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