Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,312,212 +TC coding (538/1104 nt) ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,312,2101.T96.5% 82.2 / ‑3.1 29coding (540/1104 nt)ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (13/15);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01
*NC_0009132,312,2102.C96.5% 94.2 / ‑2.6 29coding (540/1104 nt)ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (13/15);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCT‑‑TCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAA  >  NC_000913/2312143‑2312268
                                                                    ||                                                          
cgccgACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCt                                                              <  1:570868/68‑1 (MQ=255)
cgccgACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCt                                                              <  1:744462/68‑1 (MQ=255)
     aCGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCT‑‑Tctca                                                       <  1:1325387/68‑4 (MQ=255)
                  tttGACGCTGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGc                                            <  2:397269/68‑1 (MQ=255)
                     gACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGt                                         <  1:1279298/68‑1 (MQ=255)
                      aCGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGtt                                        <  2:1286452/68‑1 (MQ=255)
                      aCGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGtt                                        <  1:153797/68‑1 (MQ=255)
                        gCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGtttt                                      >  1:1342100/1‑68 (MQ=255)
                            tGCGTCATGACCGTTGTTAGTTTCGCCACTAGTAAAGCCTTCTCACCAGTTGGGTTGCCGTTTTTAcg                                  <  2:1458819/68‑2 (MQ=255)
                            tGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTAcc                                  >  1:270540/1‑68 (MQ=255)
                                 cACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGt                             <  1:1274537/68‑1 (MQ=255)
                                 cACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGt                             <  1:1355890/68‑1 (MQ=255)
                                 cACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGt                             <  1:424707/68‑1 (MQ=255)
                                  aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTggca                            >  1:604901/1‑66 (MQ=255)
                                  aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTa                            >  2:1158549/1‑68 (MQ=255)
                                  aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTa                            >  2:1313588/1‑68 (MQ=255)
                                  aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTAACCCTGGTa                            >  2:546275/1‑68 (MQ=255)
                                    gACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTACt                          >  1:1176735/1‑68 (MQ=255)
                                      ccGTTGTTAGTTACGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGCACTGa                        <  2:604897/68‑1 (MQ=255)
                                           gttagttaCGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGc                   <  1:81585/68‑1 (MQ=38)
                                            ttagttaCGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACa                     >  2:1084040/1‑65 (MQ=38)
                                            ttagttaCGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACa                     <  1:1084044/65‑1 (MQ=38)
                                               gttaCGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAg               <  2:1243268/68‑1 (MQ=38)
                                                ttaCGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGt              <  2:1176731/68‑1 (MQ=37)
                                                ttaCGCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGt              >  1:910136/1‑68 (MQ=37)
                                                  aCGCCACTAGTAAAGCCTTCTCACAAGATGGGTTGCAGTTTTTACACTGGTACTGAACAGCAAAGTTc            >  1:877295/1‑68 (MQ=25)
                                                    gCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCg                                          >  2:517279/1‑36 (MQ=25)
                                                    gCCACTAGTAAAGCCTTCTCACCAGATGGGTTGCCg                                          <  1:517282/36‑1 (MQ=25)
                                                      cACTAGTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGc        >  1:794693/1‑68 (MQ=25)
                                                           gTAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTTCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCa   <  2:388804/68‑1 (MQ=25)
                                                            tAAAGCCTTCTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa  >  2:507078/1‑68 (MQ=25)
                                                                    ||                                                          
CGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCT‑‑TCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAA  >  NC_000913/2312143‑2312268

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