Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,743,603 Δ1 bp coding (1171/1275 nt) yiaN → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,743,6030C.96.8% 264.4 / ‑0.4 62coding (1171/1275 nt)yiaNputative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (20/40);  total (21/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATCGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTC  >  NC_000913/3743532‑3743673
                                                                       |                                                                      
tCGGTTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGc                                                                   <  1:511442/76‑1 (MQ=255)
 cGGCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCa                                                                  >  1:484951/1‑76 (MQ=255)
  ggCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAg                                                                 <  1:2311264/76‑1 (MQ=255)
     ttAATTACCCCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCa                                                              <  2:956160/76‑1 (MQ=255)
     ttAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATCACGCAGTCa                                                              >  2:737101/1‑76 (MQ=255)
     ttAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGACGATGCAGTCa                                                              <  2:947758/76‑1 (MQ=255)
      tAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCag                                                             <  2:3180498/76‑1 (MQ=255)
        aTTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCagag                                                           <  2:3859234/76‑1 (MQ=255)
         ttACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAgg                                                          <  2:635474/76‑1 (MQ=255)
          tACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGc                                                         <  2:4502110/76‑1 (MQ=255)
             acCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGtt                                                      >  1:2043463/1‑76 (MQ=255)
             acCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGtt                                                      >  1:2958715/1‑76 (MQ=255)
              ccGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGttt                                                     <  2:2836310/76‑1 (MQ=255)
               cGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGtttt                                                    <  2:4199262/76‑1 (MQ=255)
                gCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTc                                                   >  1:1029500/1‑76 (MQ=255)
                  cTATCGGTAATGTACTTAACGTTATTTCATGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTccc                                                 <  1:4682130/76‑1 (MQ=255)
                  cTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTccc                                                 >  1:2599617/1‑76 (MQ=255)
                  cTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTcca                                                 <  2:2727472/76‑2 (MQ=255)
                   tATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCt                                                <  1:3219482/76‑1 (MQ=255)
                    aTCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCtt                                               >  1:2145654/1‑76 (MQ=255)
                    aTCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCtt                                               <  1:4003344/76‑1 (MQ=255)
                     tCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTa                                              <  1:2496319/76‑1 (MQ=255)
                      cGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTAc                                             <  1:4723406/76‑1 (MQ=255)
                      cGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTAc                                             <  1:4682848/76‑1 (MQ=255)
                       ggTAATGTACTTAACGTTATTTCAGTGTTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACg                                            <  1:4519132/76‑1 (MQ=255)
                        gTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGt                                           >  2:4239848/1‑76 (MQ=255)
                        gTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGt                                           <  2:2312874/76‑1 (MQ=255)
                          aaTGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTcc                                         >  2:1279857/1‑76 (MQ=255)
                          aaTGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTcc                                         <  1:1613729/76‑1 (MQ=255)
                           aTGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGttt                                                     <  1:3348212/63‑1 (MQ=255)
                           aTGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGttt                                                     >  2:3348212/1‑63 (MQ=255)
                           aTGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCt                                        <  1:1675032/76‑1 (MQ=255)
                           aTGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCt                                        <  2:2155603/76‑1 (MQ=255)
                            tGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTg                                       <  1:3293175/76‑1 (MQ=255)
                               aCTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGtt                                    <  2:4022318/76‑1 (MQ=255)
                                  tAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTta                                 >  2:1852226/1‑76 (MQ=255)
                                    aCGTTATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTtata                               <  1:3664641/76‑1 (MQ=255)
                                       ttATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTc                            <  1:3788664/76‑1 (MQ=255)
                                        tATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCa                           <  2:1563077/76‑1 (MQ=255)
                                        tATTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCa                           <  2:879982/76‑1 (MQ=255)
                                         aTTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCat                          <  2:4728104/76‑1 (MQ=255)
                                         aTTTCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCat                          <  1:905891/76‑1 (MQ=255)
                                           ttCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCatta                        >  1:3620911/1‑76 (MQ=255)
                                            tCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGTGTTTTCCCTTACGTCCTGGTTTTATATTCattat                       <  1:3408287/76‑1 (MQ=255)
                                            tCAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCattat                       <  1:4239848/76‑1 (MQ=255)
                                             cAGGGGTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCattatt                      <  2:484951/76‑1 (MQ=255)
                                             cAGGGGCGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCattatt                      >  2:809354/1‑76 (MQ=255)
                                                gggTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGt                   >  1:490345/1‑76 (MQ=255)
                                                 ggTGGCAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTc                  >  2:66588/1‑76 (MQ=255)
                                                    ggCAAAATTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTa               <  1:2879152/76‑1 (MQ=255)
                                                      cAAAACTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTAtt             <  2:2595460/76‑1 (MQ=255)
                                                       gaaaCTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGGAttt            >  2:2579931/2‑76 (MQ=255)
                                                       aaaaCTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTAttt            >  2:1891578/1‑76 (MQ=255)
                                                       aaaaCTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTAttt            >  2:4706647/1‑76 (MQ=255)
                                                           cTCAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGttt        <  1:1040952/76‑1 (MQ=255)
                                                             cAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGttttt      <  1:2240489/76‑1 (MQ=255)
                                                             cAAATTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGttttt      >  2:4763445/1‑76 (MQ=255)
                                                              aaaTTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGtttttt     >  2:50920/1‑75 (MQ=255)
                                                               aaTTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTAt    <  2:4774963/76‑1 (MQ=255)
                                                                aTTCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTAtt   <  1:1131552/76‑1 (MQ=255)
                                                                 ttCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTc  <  1:2714979/76‑1 (MQ=255)
                                                                 ttCGATGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTc  >  1:1344448/1‑76 (MQ=255)
                                                                       |                                                                      
TCGGCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATCGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATTC  >  NC_000913/3743532‑3743673

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