Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,614,614 Δ2 bp coding (681‑682/849 nt) focB → putative formate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,614,6130C.94.1% 60.9 / ‑1.8 17coding (680/849 nt)focBputative formate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (7/9);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01
*NC_0009132,614,6140T.100.0% 69.6 / NA 17coding (681/849 nt)focBputative formate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (7/10);  total (7/10)

TTATTGCGGGCGTTCGTTATGCGACAAAATCGTCGCCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCCCTCCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTG  >  NC_000913/2614476‑2614737
                                                                                                                                         ||                                                                                                                           
ttATTGCGGGCGTTCGTTATGCGACAAAATCGTCGCCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGccccc                                                                                                                             >  2:798624/1‑138 (MQ=255)
        ggCGTTCGTTATGCGACAAAATCGTCACCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAACTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCccccccccttct                                                                                                                     >  1:230065/1‑135 (MQ=255)
             tcgttATGCGACAAAATCGTCGCCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAg                                                                                                              >  2:577061/1‑139 (MQ=255)
                       aCAAAATCGTCGCCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGcaca                                                                                                    >  1:454486/1‑139 (MQ=255)
                                                  tCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGa                                                                         >  1:835810/1‑139 (MQ=255)
                                                      ccTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGt                                                                     <  2:454486/139‑1 (MQ=255)
                                                       cTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCCCCCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGt                                                                     <  1:798624/139‑1 (MQ=255)
                                                       cTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCACCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTATCGCAGACAATTTTCCGGCACTGACGGTc                                                                    <  1:577061/139‑1 (MQ=255)
                                                          tttGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGc                                                                 <  1:167089/139‑1 (MQ=255)
                                                          tttGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGc                                                                 <  2:911786/139‑1 (MQ=255)
                                                            tGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCa                                                               >  2:297856/1‑139 (MQ=255)
                                                                            gAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAAt                                               <  1:297856/139‑1 (MQ=255)
                                                                                 cTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCt                                          <  1:783273/139‑1 (MQ=255)
                                                                                  tGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTc                                         >  2:807143/1‑139 (MQ=255)
                                                                                      tCGCCATTTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCgg                                     <  2:464129/139‑1 (MQ=255)
                                                                                            aTTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGcccccc                                                                                                                            >  1:456491/1‑46 (MQ=255)
                                                                                            aTTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGcccccc                                                                                                                            <  2:456491/48‑3 (MQ=255)
                                                                                              ttgtttgtGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGCCAATTTTCCGGCCCTGACGGTCAGCCATTTTATTACCGCCAATCCGCTCCCGGTGATGcag                             >  2:141476/1‑137 (MQ=255)
                                                                                                       aTTCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTACTCGCAGCTGGCGCACAGTAGCGCAGACCATTTTCCGCCACTGACGGTCGGCCATTTTATTACCGCCAATCTGCTCCAGGTGATGCTGGGTAatagt                    <  2:230065/139‑3 (MQ=255)
                                                                                                        ttCCGTTCGCCATTGCCATTCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAAtatta                   >  2:448213/1‑139 (MQ=255)
                                                                                                                 ccattgccattCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGgcggtg          <  1:448213/139‑1 (MQ=255)
                                                                                                                         attCGCCATTTCGCCC‑‑CCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTg  <  1:141476/139‑1 (MQ=255)
                                                                                                                                         ||                                                                                                                           
TTATTGCGGGCGTTCGTTATGCGACAAAATCGTCGCCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCCCTCCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTG  >  NC_000913/2614476‑2614737

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